{"id":"2b801318-3f2a-4c6d-a705-76ecdeb1e5d3","arxiv_id":"2412.06680","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"TDDFT simulations of Coulomb explosions produce lower ion kinetic energies and broader angular distributions than classical fixed-charge models, and the authors attribute the experimentally observed broadening to these quantum electronic effects.","lead":"This paper uses time-dependent density functional theory to simulate Coulomb explosions of acetylene, butane, and isoxazole, and compares the results with classical and semi-classical models. The simulations show that explicitly including electrons lowers the kinetic energies of the ejected ions and broadens their angular distributions, which the authors say matches the broader momentum patterns seen in experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Classical baseline conflates gradual ionization charging with quantum electron dynamics; semi-classical data locate the effect before t*, so a time-dependent-charge classical control is needed before attributing broadening to quantum effects.","rationale":"The paper's core claim is that the discrepancy between classical Coulomb-explosion simulations and experiment arises from quantum effects, specifically electron dynamics. The evidence is that TDDFT runs produce lower ion kinetic energies and broader angular distributions than a classical model with fixed post-ionization charges. The reader correctly identifies the confound: the classical model also starts full Coulomb repulsion at t=0, while TDDFT charges rise over the pulse. The present stress test confirms this is the most load-bearing issue. The semi-classical simulations, which use TDDFT until t*=25 fs and classical afterwards, reproduce the quantum angular distributions, so the decisive physics is in the ionization window. The text itself (Section IV) says the kinetic-energy reduction is 'primarily due to ... the gradual increase in ion charge during ionization.' Thus the only missing control is a classical simulation that includes that gradual charge increase but no other electron dynamics. The proposed test would discriminate between two explanations: (1) the broader distributions are caused by the time-dependent charging alone, in which case the attribution to quantum effects is an overstatement; or (2) the broader distributions require the explicit electron density dynamics of TDDFT, in which case the claim is supported. Because the paper currently provides no such control, the conditional verdict is appropriate.","tokens_in":19239,"tokens_out":4235,"duration_ms":45299,"concrete_test":"For C2H2, define a 'time-dependent-charge classical' model: use Eq. (12) with charges q_i(t) = Z_i * [N(0)-N(t)]/[N(0)-N(t*)] for t<t* and q_i(t)=final average charge for t>=t*, with the same 300 K Boltzmann initial conditions and the same number of trajectories (200). Compare the final velocity and angular distributions (Figs. 6-10) to the quantum and semi-classical results. If the time-dependent-charge classical distributions match the quantum ones (e.g., angular spread within ~2-3 degrees and kinetic-energy spreads within ~0.5 eV), then the broadening is due to finite ionization time, not to quantum electron dynamics, and the central claim must be revised. If they remain narrow like the fixed-charge classical results, the quantum attribution is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the experimentally observed broadening is due to 'quantum effects' rests on comparing TDDFT runs to a classical baseline in which every atom carries its final average charge and full Coulomb repulsion acts from t=0. The quantum runs instead accumulate charge during the ~12-25 fs ionization window. The paper itself (Section IV) attributes the lower kinetic energies to 'the gradual increase in ion charge during ionization,' and the semi-classical runs, which switch to classical at t*=25 fs, closely reproduce the quantum angular distributions (Figs. 8-10). This shows the relevant physics is established during the ionization phase, before t*. But within that phase, 'gradual charging' and 'quantum electron dynamics' are not separated. A classical simulation with fixed final charges but delayed onset would already lower kinetic energies; a simulation with charges ramped following the TDDFT ionization curve would test whether the broader angular distributions also follow from charging dynamics alone. Without such a control, the phrase 'quantum effects' overstates what is demonstrated: the results establish that a fixed-charge, t=0 full-repulsion classical model is inadequate, not that electron dynamics beyond a time-dependent charge state are responsible.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares three levels of simulation for laser-induced Coulomb explosion of acetylene, butane, and isoxazole: full real-time TDDFT with Ehrenfest ion dynamics, fixed-charge classical dynamics, and a semi-classical hybrid that switches from TDDFT to classical at a time t* after ionization. The central claim is that quantum electron dynamics lower the final kinetic energies of all ions and broaden the angular and momentum distributions relative to fixed-charge classical simulations, and that this broadening is the origin of the discrepancy between classical Coulomb-explosion simulations and experimental ion-momentum patterns from Ref. [58]. The TDDFT method, classical baseline construction, and semi-classical switching procedure are described in Section II; results for the three molecules are presented in Section III; Section IV summarizes the claims.","tokens_in":19528,"tokens_out":4321,"duration_ms":52114,"significance":"If established, the result is valuable for the Coulomb-explosion-imaging community because it would show that standard fixed-charge classical simulations omit a physically important ionization-phase effect, and it would provide an interpretation of the experimentally observed broad momentum distributions. The paper has clear strengths: it uses explicit real-time electron dynamics, reports consistent trends across three molecules and two field strengths, introduces a semi-classical control that localizes the differences to the ionization phase, and does not fit any target quantity to produce the quantum results. The main limitation is that the comparison does not cleanly separate 'time-dependent charging' from 'quantum electron dynamics,' and the paper does not quantitatively compare its quantum distributions with the experiment it cites.","major_comments":[{"comment":"The classical baseline assigns each atom a fixed charge equal to its average post-ionization value and applies full Coulomb repulsion from t=0, whereas the TDDFT runs accumulate ion charge gradually over the approximately 12–25 fs ionization window. Section IV itself states that the kinetic-energy reduction is 'primarily due to ... the gradual increase in ion charge during ionization.' Therefore, a classical control simulation with charges ramped according to the TDDFT ionization curve, or at least with the same delayed onset of repulsion, is required before the lower kinetic energies and broader angular distributions can be attributed to quantum electron dynamics. As written, the comparison demonstrates that the fixed-charge/t=0 classical model is inadequate, not that electron dynamics beyond time-dependent charging are responsible for the broadening.","section":"§II, Eq. (12) and §IV"},{"comment":"The paper claims that quantum TDDFT results 'align closely' with experimental observations and that the experimentally observed broader momentum distributions are due to the quantum effects identified in this work, but no quantitative comparison with the experimental data of Ref. [58] is presented. The figures compare classical, semi-classical, and quantum simulations only with each other. A direct comparison, such as measured versus simulated Newton-plot widths or kinetic-energy distributions for isoxazole, is needed to support the experimental-broadening claim; absent that, the statement overreaches what the simulations alone can establish.","section":"Introduction and §IV, Ref. [58]"}],"minor_comments":[{"comment":"The phrase 'lower kinetic energies all ions' is missing a preposition; it should read 'lower kinetic energies for all ions.'","section":"Abstract"},{"comment":"The 'classical (with pulse)' runs are compared with the no-pulse classical results, but Eq. (12) contains no laser term and the text does not define how the pulse is added in those runs; please specify the modified equation or clarify the setup.","section":"§III A, Figs. 2 and 3"},{"comment":"The angle convention in Fig. 10d is described with 'when when' and the treatment of Quadrant III as 180–270 degrees is physically equivalent to negative angles between -180 and -90 degrees; please clarify the plotting convention so the reader can interpret the x-axis unambiguously.","section":"§III B, Fig. 10d"},{"comment":"The sentence 'The choice of isoxazole as the target molecule was motivated by previous by the aforementioned recent study [58]' contains a doubled phrase and should be rewritten.","section":"§III D"},{"comment":"The transition time t* is selected using the electron count from a single quantum trajectory; if the semi-classical simulations all use the same t*, please state whether the spread in ionization completion times across trajectories was considered, since per-trajectory variations could affect the comparison.","section":"§II, t* definition"}],"recommendation":"major_revision","confidential_remarks":"The central scientific question is well posed and the simulations appear internally consistent, but the attribution of the observed broadening to 'quantum effects' hinges on a control that is currently missing: a classical simulation with a time-dependent charge ramp extracted from the TDDFT ionization dynamics. This is a fixable gap, but it is load-bearing for the paper's main claim and for the claimed connection to Ref. [58]. I also note that the title and abstract use 'quantum effects' broadly while the nuclei are treated classically; the claim is specifically about electron dynamics, which should be reflected in the wording."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper compares classical, semi-classical, and TDDFT Coulomb explosion simulations for acetylene, butane, and isoxazole. The genuinely new thing is the semi-classical decomposition: run TDDFT up to a transition time t* (when ionization finishes), then hand off to classical point-charge dynamics. The authors show that the angular broadening of the full quantum runs is already present at t*, which localizes the physics to the ionization phase. That is a useful framing. The internal comparisons are consistent across three molecules, the statistics are adequate for the qualitative claims, and the classical charges are taken from the quantum averages rather than fitted, so there is no circularity.\n\nThe soft spot is the attribution of the effect to 'quantum effects.' The classical baseline uses fixed final-average charges with full Coulomb repulsion from t=0, while the TDDFT runs charge up gradually during the pulse. The paper itself (Section IV) credits both electron density attraction and gradual charging for the lower kinetic energies. So the comparison conflates two differences: electron dynamics beyond the charge state, and the temporal charging profile. The semi-classical results show the effect is locked in before t*, but they do not separate these two factors. A classical control with charges ramped following the TDDFT ionization curve would settle whether the broadening is a genuine electron-dynamics effect or just a delayed-onset Coulomb explosion. Without that, 'quantum effects' overstates what has been demonstrated.\n\nSecond, the comparison with experiment [58] is only qualitative. The authors say the quantum distributions resemble the measured Newton plots, but they do not overlay data, quantify widths, or compute any discrepancy metric. For a paper motivated by explaining an experimental discrepancy, that is a notable gap.\n\nNone of this sinks the paper. The semi-classical decomposition is clean, the direction of the effect is robust, and the conclusion that fixed-charge classical models miss something real is well supported. I would send this to peer review, but I would require either the time-dependent-charge classical control or a more careful wording of the claim, plus a quantitative experimental comparison.","headline":"Good semi-classical decomposition, but the main claim overreaches until a time-dependent-charge classical control is tested.","tokens_in":19991,"tokens_out":3312,"would_cite":true,"duration_ms":34650,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.15.ee","33.80.Rv"],"model":"deepseek-v4-flash","headline":"Quantum electron dynamics, not classical point charges, are what makes experimental Coulomb explosion ion momenta spread as widely as they do.","keywords":["Coulomb explosion","time-dependent density functional theory","quantum effects","ion momentum distributions","Newton plots","angular distribution","kinetic energy","molecular imaging"],"falsifier":"Run a classical explosion model in which each ion's charge follows the same time-dependent ionization curve produced by the TDDFT runs, with the electron clouds still omitted. If that delayed-charging classical model already reproduces the quantum kinetic energies, then the energy lowering is a charging-schedule effect rather than a quantum electron effect, and only the angular broadening, if it persists as a quantum feature, remains a genuine electron-dynamics signature. The predicted difference is quantitative and directly checkable against the same Newton plots.","tokens_in":19087,"feed_emoji":"⚛️","tokens_out":6134,"duration_ms":61827,"temperature":0.7,"pith_summary":"This paper argues that the missing ingredient in classical Coulomb explosion simulations is the electrons themselves. Using time-dependent density functional theory, the authors find that when electron density is propagated in real time, all ions end up with lower kinetic energies and their trajectories spread over wider angles than the fixed-charge classical models predict. The broader, more diverse distributions match the patterns seen in experimental ion-momentum (Newton plot) measurements on acetylene, butane, and isoxazole. A sympathetic reader would care because Coulomb explosion imaging is used to reconstruct molecular structures, and the paper claims that classical analyses overstate how sharply that structure is defined.","feed_headline":"Quantum effects widen Coulomb explosion ion patterns","feed_subtitle":"In TDDFT, electrons and gradual ionization slow ions and broaden their trajectories, matching experiments classical models miss.","key_machinery":"The load-bearing mechanism is real-time TDDFT with Ehrenfest dynamics: Kohn–Sham orbitals are propagated in time on a real-space grid, the laser field is applied in the dipole approximation, ions are moved classically under forces from the laser, Coulomb repulsion, and the gradient of the electron–ion interaction energy, and a complex absorbing potential removes ionized density. The comparison set is what isolates quantum effects: purely classical simulations that assign each atom its average post-ionization charge from the start, and semi-classical simulations that run TDDFT through the ionization phase and then switch to classical repulsion at t* when ionization ends. The semi-classical results sit close to the quantum ones, which locates the decisive quantum influence in the close-proximity, charge-evolving ionization phase rather than in the later separated-fragment dynamics.","core_discovery":"The central discovery is that quantum effects, meaning the effects of explicitly propagated electron density, change the Coulomb explosion in a systematic way: lower final speeds for every ion and a broader, more varied set of final velocity vectors. In the quantum simulations the ions only acquire their full charge gradually as the laser pulse strips electrons away, and the remaining electron clouds exert attractive forces on the nuclei, so the repulsion is weaker and varies from run to run. The classical fixed-charge model, which switches on the full post-ionization charges at the first time step and has no electron density, produces faster, more collimated ions. Because the quantum Newton plots are visibly broader, the authors conclude that this broadening, rather than experimental noise or alignment effects, is what makes experimental ion momentum distributions look wider than the classical predictions.","pith_inferences":["A classical simulation with time-dependent charges matched to the TDDFT electron-loss curve might recover much of the kinetic-energy deficit, suggesting the 'quantum effect' label partly covers a delayed-ionization effect that a classical model could emulate.","The paper's comparisons are visual; quantifying the overlap between simulated and experimental Newton plots with a statistical measure would convert the qualitative match into a testable claim about the origin of the broadening.","If the quantum spread is as large as the paper indicates, Coulomb explosion imaging analyses should treat the classical-model angular precision as an upper bound on structural resolution, not the expected experimental uncertainty."],"forward_implications":["Classical fixed-charge Coulomb explosion models systematically overestimate final ion kinetic energies; quantum simulations give lower energies for all three molecules studied.","Newton plots from quantum simulations are broader than classical ones and closer to experiment, so structural reconstruction from Coulomb explosion imaging should use quantum or semi-classical distributions rather than classical point-charge spreads.","The semi-classical method—TDDFT through ionization, then classical—reproduces the quantum distributions to good approximation, suggesting a practical route for larger molecules where full quantum propagation is expensive.","Doubling the laser field strength increases ionization and ion speeds in both models, but the quantum-versus-classical differences persist, so the effect is not merely an artifact of low charge states."],"supporting_citations":[{"why":"provides the experimental ion-momentum patterns that appear broader than classical simulations and supplies the isoxazole molecule and laser parameters","marker":"[58]"},{"why":"earlier trajectory-surface-hopping study of quantum contributions to Coulomb explosion imaging that this work positions its TDDFT approach against","marker":"[67]"},{"why":"prior TDDFT study of strong-field ionization of hydrocarbon molecules that validates the Ehrenfest/TDDFT treatment of Coulomb explosion","marker":"[70]"},{"why":"TDDFT study of alignment-dependent ionization of acetylene and ethylene, the method used here for acetylene","marker":"[71]"},{"why":"TDDFT study of ionization and fragmentation of C2H2 by circularly polarized pulses, informing the ionization and fragmentation analysis","marker":"[72]"},{"why":"TDDFT strong-laser-pulse ionization and Coulomb explosion of hydrocarbon molecules, the methodological basis for the quantum simulations","marker":"[74]"},{"why":"supplies the norm-conserving pseudopotentials used for the ion–electron interaction","marker":"[76]"},{"why":"provides the Perdew–Zunger parameterization used for the exchange-correlation potential in the TDDFT Hamiltonian","marker":"[77]"},{"why":"provides the complex absorbing potential that removes ionized electron density and defines the electron-loss measurement","marker":"[78]"}],"fun_headline_variants":["Quantum effects slow and scatter Coulomb explosion ions","Electron dynamics broaden ion patterns in quantum Coulomb blasts","Quantum simulations yield slower, wider ion trajectories","TDDFT shows quantum electrons alter explosion ion speeds","Gradual ionization and electrons defy classical Coulomb predictions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The attribution of the broader experimental ion distributions to quantum effects assumes that the classical baseline differs from the quantum runs only by missing electron dynamics, and not also because the classical model starts with full charges at t = 0 whereas the quantum model builds up charge gradually during the pulse.","fun_headline_variants_meta":{"raw":{"variants":["Quantum effects slow and scatter Coulomb explosion ions","Electron dynamics broaden ion patterns in quantum Coulomb blasts","Quantum simulations yield slower, wider ion trajectories","TDDFT shows quantum electrons alter explosion ion speeds","Gradual ionization and electrons defy classical Coulomb predictions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1303,"prompt_tokens":822,"completion_tokens":481,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":438,"completion_tokens_details":{"reasoning_tokens":409}},"tokens_in":438,"tokens_out":481,"duration_ms":6206,"temperature":1.0,"reasoning_tokens":409,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:25:27.227570+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a classical explosion model in which each ion's charge follows the same time-dependent ionization curve produced by the TDDFT runs, with the electron clouds still omitted. If that delayed-charging classical model already reproduces the quantum kinetic energies, then the energy lowering is a charging-schedule effect rather than a quantum electron effect, and only the angular broadening, if it persists as a quantum feature, remains a genuine electron-dynamics signature. The predicted difference is quantitative and directly checkable against the same Newton plots.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the experimental ion-momentum patterns that appear broader than classical simulations and supplies the isoxazole molecule and laser parameters"},{"cited_title":"Singh, C","cited_arxiv_id":null,"evidence_quote":"earlier trajectory-surface-hopping study of quantum contributions to Coulomb explosion imaging that this work positions its TDDFT approach against"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"prior TDDFT study of strong-field ionization of hydrocarbon molecules that validates the Ehrenfest/TDDFT treatment of Coulomb explosion"},{"cited_title":"Russakoff, S","cited_arxiv_id":null,"evidence_quote":"TDDFT study of alignment-dependent ionization of acetylene and ethylene, the method used here for acetylene"},{"cited_title":"Russakoff and K","cited_arxiv_id":null,"evidence_quote":"TDDFT study of ionization and fragmentation of C2H2 by circularly polarized pulses, informing the ionization and fragmentation analysis"},{"cited_title":"Bubin, M","cited_arxiv_id":null,"evidence_quote":"TDDFT strong-laser-pulse ionization and Coulomb explosion of hydrocarbon molecules, the methodological basis for the quantum simulations"},{"cited_title":"Troullier and J","cited_arxiv_id":null,"evidence_quote":"supplies the norm-conserving pseudopotentials used for the ion–electron interaction"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the Perdew–Zunger parameterization used for the exchange-correlation potential in the TDDFT Hamiltonian"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the complex absorbing potential that removes ionized electron density and defines the electron-loss measurement"}],"review_version":1}