{"id":"19a3e7d4-5588-462c-a3b2-942883dad2dd","arxiv_id":"2411.12885","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Time-resolved Coulomb explosion imaging measures the oscillating internuclear distribution of vibrational wave packets in triplet-state K2 and Rb2 on helium nanodroplets.","lead":"Vibrational wave packets in potassium and rubidium dimers on helium nanodroplets were imaged in real time using Coulomb explosion. The experiment tracks the bond-length distribution as it oscillates, confirming a two-level wave packet and a slow decay attributed to interactions with the droplet.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"K2 amplitude decay time is likely misstated by a factor of two: the SFT fit is to spectral power (τ=260±30 ps), so the oscillation amplitude lifetime is ≈520 ps, contradicting the abstract's '~260 ps' and the observed 0.035→0.020 Å over 300 ps.","rationale":"I read the paper as making two connected claims: first, that Coulomb explosion imaging can map the measured kinetic-energy release to a time-dependent internuclear distribution P(R,t) and thereby image vibrational wave packets on helium nanodroplets; second, that the K2 wave-packet oscillation amplitude decays with a time constant of about 260 ps due to coupling to the droplet. The imaging claim is well supported: the static P(R) for the v=0 triplet state matches known distributions, the measured oscillation frequency 611.6 GHz agrees with the calculated ν1,0 for 39K2, and the early-time P(R,t) evolution agrees qualitatively with the dynamic Stark simulation. The probe-induced distortion invoked in Sec. V.C is a legitimate concern for the early transient, but the authors exclude t≤40 ps from the decay fit, and a time-independent probe shift would not create the observed oscillatory structure; I therefore do not treat the probe mapping as the most decisive weakness. The factor-of-two power/amplitude confusion is instead concrete, internal, and checkable. The SFT analysis fits spectral power, whose decay constant is half the amplitude decay constant, yet the abstract reports the result as an amplitude decay time. The observed amplitudes themselves, 0.035 Å early and 0.020 Å near 300 ps, select τ_amp≈520 ps, not 260 ps. This error changes a headline number and weakens the comparison to Grüner et al. and the droplet-relaxation interpretation. Because the imaging demonstration remains intact and the decay issue is quantitative rather than foundational, the reader's conditional verdict is appropriate and should not be altered. My agreement with the reader is partial: their weakest_assumption focuses on the Coulomb mapping, but their rationale already lists the power/amplitude conflation; my stress-test picks that conflation as the most load-bearing specific defect.","tokens_in":18200,"tokens_out":12456,"duration_ms":139062,"concrete_test":"Reanalyze the K2 ⟨R⟩(t) trace by extracting the instantaneous oscillation envelope, e.g. via a Hilbert transform after bandpass filtering around 612 GHz or by fitting the detrended trace's local maxima and minima, for t=10–300 ps while excluding t≤40 ps as in the paper. Fit A(t)=A0 exp(−t/τ_amp) to this envelope and compare τ_amp with both 260 ps and 520 ps. Also refit the SFT power decay and verify whether τ_power≈τ_amp/2. If τ_amp≈520 ps, the abstract, conclusion, and comparison to Ref. 35 must be corrected; if τ_amp≈260 ps, the current wording is consistent and no correction is needed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim about long-term dynamics is internally inconsistent. In Sec. V.D the authors fit an exponential decay ∝ exp(−t/τ) to the spectral power |SFT(t_i,ν)|² extracted from 20 ps sliding windows, obtaining τ=260±30 ps. But spectral power is proportional to the square of the oscillation amplitude, so an amplitude decay exp(−t/τ_amp) produces a power decay exp(−2t/τ_amp). The fitted value is therefore τ_power=260±30 ps, implying τ_amp≈520 ps. The abstract and conclusion nevertheless state that 'the decay time of the amplitude is ∼260 ps' and compare this with the 0.3 ns amplitude lifetime reported by Grüner et al. in Ref. 35. The directly reported amplitudes confirm the factor-of-two correction: with τ_amp=260 ps, 0.035 Å at early times would decay to about 0.011 Å at 300 ps, whereas the measured value is 0.020 Å; with τ_amp≈520 ps the expected endpoint is 0.020 Å. Thus the headline decay time is off by a factor of two, and the asserted closeness to Ref. 35 plus the attribution to droplet-induced vibrational relaxation rest on the wrong time constant. This does not invalidate the core imaging demonstration, but it changes a specific quantitative result and its physical interpretation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports time-resolved Coulomb explosion imaging of vibrational wave packets in K2 and Rb2 in their lowest triplet state on helium nanodroplets. A nonresonant pump pulse creates a coherent superposition of v=0 and v=1 (with a small v=2 contribution) via the dynamic Stark effect, and a delayed intense probe pulse doubly ionizes the dimer; the kinetic energy of the fragment ions is converted to an internuclear distance distribution P(R,t) using a Coulomb repulsion law. The authors observe oscillatory P(R,t) and ⟨R⟩(t) for 300 ps (K2) and 100 ps (Rb2), identify the main oscillation frequencies as the v=0–1 vibrational coherences, compare them with independent calculations, and reproduce the early dynamics with a one-dimensional Schrödinger-equation simulation based on literature potentials and polarizabilities. They also report a gradual decrease of the K2 oscillation amplitude and extract a decay time τ=260±30 ps from a sliding-window Fourier transform, attributing the decay to coupling with the helium droplet.","tokens_in":18510,"tokens_out":3232,"duration_ms":37012,"significance":"If the central claim holds, the work provides a direct, time-resolved structural measurement of vibrational wave packets in molecules on helium nanodroplets, going beyond previous yield-based probes and extending Coulomb explosion imaging to a new class of weakly bound, droplet-embedded systems. The strength of the paper lies in the direct comparison of measured beat frequencies with independently calculated vibrational level spacings, and in the use of a dynamic Stark simulation with no fitted parameters aimed at reproducing the measured dynamics. The observed agreement with the simulated ⟨R⟩(t) and the matching of the main spectral peak at 611.7 GHz versus the calculated 611.8 GHz (K2) and 396.8 GHz versus 398.3/396.1 GHz (Rb2) are convincing evidence that the observed oscillations are the anticipated v=0–1 vibrational wave packets. The main quantitative weakness is the treatment of the decay time of the oscillation amplitude, which is computed from the decay of spectral power rather than amplitude.","major_comments":[{"comment":"The extracted decay constant is inconsistent with its interpretation. The sliding Fourier transform yields power spectra |SFT(t_i,ν)|², and the fit in Fig. 8(b) is to the spectral power; an exponential amplitude decay exp(−t/τ_amp) produces a power decay exp(−2t/τ_amp). Therefore the fit result τ=260±30 ps is the power decay time, implying an amplitude lifetime of about 520 ps. This is directly contradicted by the reported amplitude values: with τ_amp=260 ps, an initial amplitude of 0.035 Å would fall to 0.011 Å after 300 ps, whereas the measured value is 0.020 Å; with τ_amp≈520 ps the expected endpoint is 0.020 Å. The statements in the abstract ('decay time of the amplitude is ~260 ps'), in Sec. V.D, and in the conclusion ('on a time scale of ~0.3 ns') therefore misstate the central long-term quantitative result, and the comparison with the 0.3 ns lifetime of Grüner et al. (Ref. 35) as well as the attribution to droplet-induced vibrational relaxation rest on the wrong time constant. The authors should redo the analysis and correct all statements that quote the 260 ps value as an amplitude decay time.","section":"Sec. V.D, Fig. 8(b) and Abstract/Conclusion"},{"comment":"The paper's own explanation of the early transient modulation of ⟨R⟩(t) as a probe-induced distortion of the neutral 1³Σu⁺ potential, dependent on molecular alignment via Eq. (2), means that the measured P(R,t) is not strictly the field-free wave packet probability |Ψ(R,t)|². If this probe distortion is appreciable at other delays, the central imaging claim in the abstract ('P(R,t) ... represents the modulus square of the wave packet within the accuracy of the experiment') needs qualification. The authors should either estimate the size of the probe-induced potential distortion (and hence the systematic shift in R at the moment of ionization) or explicitly restrict the claim of mapping the field-free wave packet to the delays where the distortion is negligible, and state the resulting uncertainty in P(R,t).","section":"Sec. V.C and Sec. II (Eq. 1)"}],"minor_comments":[{"comment":"The caption of Fig. 8(b) says 'Exponential fit used to determine the decay time τ of the oscillation amplitude', but the fit is performed on spectral power; the caption and text should state that the fit is to the power and clarify the relation between the power decay time and the amplitude decay time.","section":"Sec. V.D, Fig. 8 caption"},{"comment":"The phrase 'the v = 1–2 correspondence' in the comparison of the simulated and measured spectra should read 'the v = 1–2 coherence' or 'the v = 1–2 transition'.","section":"Sec. V.C, last paragraph"},{"comment":"The conversion from ion velocity to kinetic energy uses a calibration constant k chosen to match previous works; the resulting systematic uncertainty in the absolute R scale is not propagated into the reported ⟨R⟩ values and should be stated.","section":"Sec. V.A and Sec. II"},{"comment":"The figure caption and text state that the vertical red line indicates ⟨R⟩(t), but the main text refers to the black dotted line as ⟨R⟩; the consistent notation used in the text should be reflected in the figure for clarity.","section":"Sec. V.B, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The central imaging result is sound and the frequency analysis is convincing. The main outstanding issue is the factor-of-two error in the decay time, which is a quantitative result with a physical interpretation attached to it. This is fixable by re-analyzing the data and revising the affected text, so I recommend major revision rather than rejection. The probe-distortion concern should also be addressed explicitly, but it does not invalidate the main demonstration."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the paper on time-resolved Coulomb explosion imaging of vibrational wave packets in K2 and Rb2 on helium nanodroplets. The core result is solid and genuinely new: they turn P(Ekin) from Coulomb explosion into a time-dependent P(R,t), so for the first time you get direct bond-length distributions of a vibrating molecule on a droplet instead of just ionization yields. The measured vibrational frequencies (611.7 GHz for K2, 396.8 GHz for Rb2) match independent calculations from literature potentials, and the dynamic Stark simulation reproduces the early-time response and the two-state composition of the wave packet. That part deserves a serious referee.\n\nThe soft spot is quantitative and it is in the headline number. In Sec. V.D they fit an exponential decay to the spectral power |SFT(t_i, ν)|² from 20-ps sliding windows and get τ = 260 ± 30 ps. They then call this the 'decay time of the amplitude' in the abstract and conclusion. But spectral power is the square of the amplitude, so an amplitude decay exp(−t/τ_amp) gives a power decay exp(−2t/τ_amp). The fitted τ is the power lifetime; the amplitude lifetime is about 520 ps. Their own data confirm this: 0.035 Å at early times falls to 0.020 Å at 300 ps, which is a decay factor of 0.57 over 300 ps, consistent with τ_amp ≈ 520 ps, not 260 ps. So the abstract's '~260 ps' and the comparison to Grüner et al.'s 0.3 ns lifetime rest on a factor-of-two error. The physical mechanism (droplet-induced vibrational relaxation) might still be right, but the quantitative claim needs correction.\n\nA second, minor caveat: the paper acknowledges a probe-induced distortion of the neutral potential during the first ~30 ps and excludes those points from the decay fit. That is reasonable, but it leaves a small question about whether P(R,t) at other delays is truly the field-free wave packet squared. The authors address it by correlating the transient with alignment dynamics and excluding affected data, so this does not sink the imaging claim.\n\nThe paper does not ship code or data, but the key results are reproducible from the described methods, and the fit error is the main issue. I'd send it to review with a request to fix the amplitude versus power lifetime, and to state the uncertainty on the early-time amplitude decay explicitly.\n\nRecommendation: engage with it—the imaging capability is worth citing, and the factor-of-two issue is correctable.","headline":"First time-resolved CEI of vibrational wave packets on He droplets is a real advance, but the ~260 ps amplitude decay time is off by a factor of two because it is fit to spectral power, not amplitude.","tokens_in":19029,"tokens_out":3074,"would_cite":true,"duration_ms":25761,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A delayed femtosecond probe that Coulomb-explodes alkali dimers on helium nanodroplets maps a vibrational wave packet's bond-length distribution over time, for K2 across more than 180 vibrational periods.","keywords":["time-resolved Coulomb explosion imaging","vibrational wave packet","alkali dimers","helium nanodroplets","dynamic Stark effect","stimulated Raman scattering","velocity map imaging","wave packet decoherence"],"falsifier":"Compute the $K2^{2}$+ and $Rb2^{2}$+ potential curves with a correlated electronic-structure method over the R range 3.5–7 Å: if the potential deviates from the 1/R Coulomb form by more than the VMI energy resolution at any R, the one-to-one mapping in Eq. (1) fails for that portion of P(R,t). Or, in a separate experiment, image the same K2 wave packet with a probe that does not Stark-shift the neutral potential, such as femtosecond electron diffraction; agreement would confirm that P(R,t) is the field-free |Ψ(R,t)|^2.","tokens_in":18002,"feed_emoji":"⚛️","tokens_out":9393,"duration_ms":84563,"temperature":0.7,"pith_summary":"Vibrational wave packets in the lowest triplet state of potassium and rubidium dimers, sitting on helium nanodroplets, are created by a non-resonant femtosecond pump pulse and then imaged by a delayed 50-fs probe that doubly ionizes the dimer and sends the two alkali ions flying apart. The paper's central claim is that the measured kinetic-energy distribution of the fragment ions can be converted, at each delay, into the time-dependent internuclear distance distribution P(R,t), which represents the squared wave packet within experimental accuracy. The resulting P(R,t) oscillates with the vibrational period for the full observation window, and for K2 the oscillation survives more than 180 periods while its amplitude decays from 0.035 Å to 0.020 Å on a ~260 ps timescale, attributed to weak coupling to the droplet. This matters because it turns Coulomb explosion imaging into a structural movie of vibrational motion on a helium droplet, going beyond earlier work that only recorded time-dependent ionization yields.","feed_headline":"Coulomb explosion tracks vibrating dimers for 300 ps","feed_subtitle":"A 50-fs laser pulse maps bond length in K2 and Rb2 on helium droplets through 180 vibrations","key_machinery":"The central object is the Coulomb explosion imaging relation Ekin = 7.2 eV / R (R in Å), which gives a one-to-one mapping between the kinetic energy of each Ak+ fragment and the internuclear distance at the instant of double ionization. Combined with the Jacobian transformation from P(Ekin) to P(R), it turns velocity-map images recorded at many pump–probe delays into the time-dependent distribution P(R,t). The wave-packet creation is described by the dynamic Stark effect, where the pump pulse transiently deepens and shifts the 13Sigma+u potential through the polarizability interaction, launching the coherent superposition of v=0 and v=1 states.","core_discovery":"On its own terms, the paper demonstrates that timed Coulomb explosion imaging can recover the time-dependent internuclear separation distribution of a vibrational wave packet in K2 and Rb2 on helium nanodroplets. The pump pulse creates a coherent superposition dominated by the v=0 and v=1 vibrational states of the 13Sigma+u state; the probe pulse doubly ionizes the dimer, and Eq. (1), Ekin = 7.2 eV / R, converts each fragment kinetic energy into the bond length at the moment of explosion. Fourier analysis of the mean bond length <R>(t) yields the v=0–1 beat frequency, 611.7 GHz for K2 and 396.8 GHz for Rb2, consistent with calculated nu_1,0 for the free dimers. For K2 the oscillations persist for 300 ps with a gradually decreasing amplitude, and a sliding-window spectral analysis gives a decay time of 260 ± 30 ps, which the paper ascribes to vibrational relaxation caused by the nearby droplet.","pith_inferences":["An extension the paper does not make: if the probe-induced distortion is modeled quantitatively, the early transient in <R>(t) should be removable, yielding a cleaner field-free P(R,t).","Because P(R,t) contains full shape information, dephasing and revival of the wave packet should be visible as broadening and re-narrowing of the R-distribution, not just as amplitude decay of the mean.","The same Coulomb explosion protocol applied to a dimer cation such as K2+ would image vibrational relaxation directly, expected to be much faster because the charged dimer couples more strongly to the helium droplet.","The droplet-coupling interpretation predicts that the 260 ps decay time should vary with droplet size and temperature; scanning those parameters would test whether the mechanism is vibrational relaxation rather than pure dephasing."],"forward_implications":["The method provides direct structural information—bond-length distributions—rather than only time-dependent ionization yields, so it can show how the shape of the wave packet evolves, not just when it returns to the Franck–Condon region.","For K2, the ~260 ps decay of the oscillation amplitude, attributed to weak coupling with the helium droplet, implies the droplet acts as a relaxation bath for the vibrating triplet-state dimers.","The same Coulomb explosion approach is expected to work for the other homonuclear alkali dimers (Li2, Na2, Cs2) and for heteronuclear alkali dimers, including in the singlet ground state with shorter pump pulses.","The dominant v=0–1 coherence produces cosine-like <R>(t) oscillations, with the weaker v=1–2 contribution visible as a small satellite peak in the power spectrum, matching dynamic Stark effect simulations.","The early transient modulation of <R>(t) is explained as a probe-induced distortion of the neutral potential that scales with the degree of alignment, meaning the probe is not completely passive while it images."],"supporting_citations":[{"why":"Establishes the velocity-map imaging procedure and the P(Ekin)-to-P(R) transformation used here, along with the K+ kinetic-energy calibration.","marker":"[49]"},{"why":"Identifies the outermost and innermost Coulomb explosion channels in K2 on helium droplets and assigns the 13Sigma+u channel used for the wave-packet analysis.","marker":"[38]"},{"why":"Reports the ~0.3 ns lifetime of Rb2 triplet-state vibrational wave packets on helium droplets, the comparison basis for the ~260 ps decay measured here.","marker":"[35]"},{"why":"Supplies the 13Sigma+u potential of K2 used for the dynamic Stark effect simulations and for computing the nu_1,0 frequency.","marker":"[40]"},{"why":"Supplies the parallel polarizability component alpha_parallel(R) used in the dynamic Stark effect model of the pump-induced potential.","marker":"[59]"},{"why":"Provides the dynamic Stark effect description of the non-resonant Raman pump used to create the wave packets.","marker":"[46]"},{"why":"Provides the Polar Onion Peeling Abel-inversion algorithm used to convert the 2D velocity images into kinetic energy distributions.","marker":"[69]"},{"why":"Shows that the same pump pulses induce rotational wave packets in the dimers, supporting the alignment-dependent probe-effect explanation for the early transients.","marker":"[71]"}],"fun_headline_variants":["Coulomb explosion images dimer vibrations on helium droplets","180 vibrations of K2 and Rb2 captured in 300 ps","Laser pulses trace bond-length dance in alkali dimers","Time-resolved explosion reveals wave packet oscillations","Vibrating dimers on droplets imaged via Coulomb explosion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire kinetic-energy-to-bond-length mapping assumes the exploding dication feels a pure Coulomb repulsion and that the 50-fs probe does not distort the neutral molecule's potential, yet the paper itself invokes probe-induced distortion to explain the early oscillations in <R>(t).","fun_headline_variants_meta":{"raw":{"variants":["Coulomb explosion images dimer vibrations on helium droplets","180 vibrations of K2 and Rb2 captured in 300 ps","Laser pulses trace bond-length dance in alkali dimers","Time-resolved explosion reveals wave packet oscillations","Vibrating dimers on droplets imaged via Coulomb explosion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000959,"raw_usage":{"total_tokens":4155,"prompt_tokens":1083,"completion_tokens":3072,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":699,"completion_tokens_details":{"reasoning_tokens":2993}},"tokens_in":699,"tokens_out":3072,"duration_ms":23250,"temperature":1.0,"reasoning_tokens":2993,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:04:44.435428+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the $K2^{2}$+ and $Rb2^{2}$+ potential curves with a correlated electronic-structure method over the R range 3.5–7 Å: if the potential deviates from the 1/R Coulomb form by more than the VMI energy resolution at any R, the one-to-one mapping in Eq. (1) fails for that portion of P(R,t). Or, in a separate experiment, image the same K2 wave packet with a probe that does not Stark-shift the neutral potential, such as femtosecond electron diffraction; agreement would confirm that P(R,t) is the field-free |Ψ(R,t)|^2.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the velocity-map imaging procedure and the P(Ekin)-to-P(R) transformation used here, along with the K+ kinetic-energy calibration."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the outermost and innermost Coulomb explosion channels in K2 on helium droplets and assigns the 13Sigma+u channel used for the wave-packet analysis."},{"cited_title":"Grüner , author M","cited_arxiv_id":null,"evidence_quote":"Reports the ~0.3 ns lifetime of Rb2 triplet-state vibrational wave packets on helium droplets, the comparison basis for the ~260 ps decay measured here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the 13Sigma+u potential of K2 used for the dynamic Stark effect simulations and for computing the nu_1,0 frequency."},{"cited_title":"Deiglmayr , author M","cited_arxiv_id":null,"evidence_quote":"Supplies the parallel polarizability component alpha_parallel(R) used in the dynamic Stark effect model of the pump-induced potential."},{"cited_title":"\\ Shu , author E","cited_arxiv_id":null,"evidence_quote":"Provides the dynamic Stark effect description of the non-resonant Raman pump used to create the wave packets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Polar Onion Peeling Abel-inversion algorithm used to convert the 2D velocity images into kinetic energy distributions."},{"cited_title":"Kranabetter , author H","cited_arxiv_id":null,"evidence_quote":"Shows that the same pump pulses induce rotational wave packets in the dimers, supporting the alignment-dependent probe-effect explanation for the early transients."}],"review_version":1}