{"id":"5a38126e-ea7b-4dac-920e-fa48db951214","arxiv_id":"2501.16132","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A quantum Monte Carlo molecular dynamics scheme represents both electrons and nuclei by guided walker ensembles and is benchmarked on laser-driven H2 against Ehrenfest and Hartree-Fock.","lead":"This paper proposes a computational scheme called time-dependent quantum Monte Carlo for molecular dynamics, where both electrons and nuclei are represented by swarms of trajectories guided by quantum waves. The method is tested on a laser-driven hydrogen molecule and compared with Ehrenfest dynamics and Hartree-Fock.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (16)-(19) do not follow from substituting smoothed densities (12)-(13) into the Hartree equations (5)-(6); the effective potentials replace a convolution with an ad hoc weighted sum, so the 'ab initio' derivation is unsupported.","rationale":"The reader's weakest_assumption identifies the kernel-smoothed effective potentials as the place where the 'ab initio' claim fails. My stress-test sharpens this: the effective potentials are not merely heuristically motivated, they are not the mathematical consequence of the stated substitution. The convolution-versus-weighted-sum mismatch means Eqs. (16)-(19) introduce an uncontrolled modification of the Hartree potential, and the adaptive bandwidths add further freedom. Since the numerical results are single runs without error bars or convergence tests, and the benchmark is itself approximate Ehrenfest dynamics, the central claim that TDQMC captures quantum correlations self-consistently is not supported. The proposed test directly checks whether the reported agreement is a consequence of the model's stated equations or of the ad hoc smoothing. If the test showed the convolution form gives the same dynamics, the concern would be weakened; but as written, the burden is on the derivation. Therefore the reader's REJECT verdict remains appropriate, and no verdict change is needed.","tokens_in":11757,"tokens_out":7096,"duration_ms":74586,"concrete_test":"Re-derive Eq. (16) from Eqs. (5)-(6) with Eq. (12): for the paper's 1D H2 ground-state walker configuration, compute the Hartree integral H(r_i)=∫dr_j V_ee(r_i-r_j)ρ_j^smooth(r_j) by high-resolution quadrature and compare it pointwise with the right-hand side of Eq. (16) at the same r_i and walker positions. If the two differ by more than a few percent for the finite σ values used, the derivation step is invalid. Then rerun the Fig. 2 ionization dynamics using the convolution form and check whether the reported ~5% agreement with Ehrenfest survives; if it does not, the agreement is an artifact of the un-derived Gaussian weighting.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2's transition from the Hartree equations (5)-(6) to the TDQMC equations (14)-(19) is the load-bearing step for the paper's central claim, but it is not a valid derivation. The text states that substituting the smoothed walker densities (12)-(13) into (5)-(6) yields the effective potentials (16)-(19). However, those effective potentials are not the integrals of the pair interaction against the smoothed density. For example, inserting (12) into the electron-electron term of (6) gives the convolution integral ∫dr_j V_ee(r_i-r_j) [Σ_l exp(-(r_j-r_j^l)^2/σ^2)]/z_j, whereas Eq. (16) evaluates the bare potential directly at the walker separation r_i-r_j^l and multiplies it by a Gaussian in that same separation: Σ_l V_ee(r_i-r_j^l) exp(-(r_i-r_j^l)^2/σ^2)/z_j. These expressions differ for any finite bandwidth; no limit, quadrature, or other controlled approximation is provided to connect them. The adaptive bandwidths (22)-(23) are standard density-estimation choices, not determined by the Schrödinger dynamics. Thus the only mechanism by which correlations beyond Hartree enter is an ad hoc modification of the mean-field potential. Agreement with the 'exact' Ehrenfest benchmark - itself an approximation that treats nuclei classically and neglects nuclear quantum correlations - cannot establish that the method is ab initio.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a time-dependent quantum Monte Carlo (TDQMC) method for molecular dynamics in which both electrons and nuclei are represented by ensembles of Bohmian walkers and guiding waves. The formalism is derived from a single Hartree product ansatz by replacing the electronic and nuclear densities in the Hartree equations with kernel-smoothed walker densities, yielding coupled Schr\\\"odinger equations for the guiding waves and effective nonlocal potentials. The method is applied to a one-dimensional model of H2 in an intense laser pulse, with comparisons to \"exact\" Ehrenfest molecular dynamics, time-dependent Hartree-Fock, and ground-state energy estimates; a two-dimensional orientation-dynamics example is also presented. The central claim is that TDQMC is an ab initio method that captures local and nonlocal quantum correlations self-consistently.","tokens_in":12149,"tokens_out":4447,"duration_ms":41437,"significance":"If the derivation were rigorous and the numerical evidence compelling, the method could be a significant step toward scalable correlated molecular dynamics beyond mean-field theory, with potential impact on attosecond and strong-field physics. The manuscript does present a concrete algorithm and reports agreement with an established approximate benchmark, but the central derivation is not mathematically valid, and the numerical validation is insufficient to support the \"ab initio\" claim. The paper contains no convergence studies, no error bars, and no independent confirmation of the effective-potential approximation, so the current contribution does not meet the standard for a methods paper in this field.","major_comments":[{"comment":"The claim that substituting the smoothed walker densities (12)-(13) into the Hartree equations (5)-(6) yields the effective potentials (16)-(19) is not correct. For example, inserting (12) into the electron-electron term of (6) gives a convolution of the bare potential against the Gaussian kernel, namely \\int dr_j V_ee(r_i-r_j) \\sum_l \\exp(-(r_j-r_j^l)^2/\\sigma^2)/z_j, whereas Eq. (16) evaluates the potential directly at the walker separation and multiplies by a Gaussian in that same separation. These expressions differ for any finite bandwidth, and no limit, quadrature, or other controlled approximation is provided to connect them. Since this substitution is the only mechanism by which correlations beyond the single Hartree product are introduced, the paper's central claim that the method is ab initio is unsupported.","section":"Section 2, Eqs. (12)-(19)"},{"comment":"The comparison benchmark, labeled \"exact\" Ehrenfest molecular dynamics, is not an exact solution of the full many-body Schr\\\"odinger equation: it treats the nuclei as classical particles and neglects nuclear quantum correlations. Agreement to within 5% for the ionization and internuclear distance therefore cannot validate the claim that TDQMC captures local and nonlocal quantum correlations for both electrons and nuclei self-consistently. The manuscript also reports no convergence study with respect to the number of walkers M or the kernel bandwidths, so it is unclear whether the observed agreement is robust or fortuitous.","section":"Section 5, benchmark comparison"},{"comment":"The adaptive bandwidths \\sigma_j^k and \\Sigma_J^k are introduced using the standard density-estimation formulas of Abramson and Silverman, but they are not derived from the Schr\\\"odinger dynamics. They therefore constitute an uncontrolled approximation inside the effective potentials, and the paper provides no evidence that the results are insensitive to the choice of pilot bandwidths \\sigma and \\Sigma or to the initial wavefunction widths \\sigma_e and \\Sigma_n. This issue is load-bearing because the effective potentials are the sole device for modeling nonlocal correlations, and the absence of a derivation or sensitivity analysis undermines the ab initio claim.","section":"Section 2, Eqs. (22)-(23)"}],"minor_comments":[{"comment":"There are several typographical errors, including \"nuc lear\" in the abstract, \"ultashort\" in the abstract, and \"Metopolis\" in Section 5 (should be \"Metropolis\").","section":"Abstract and Section 5"},{"comment":"The notation in Eq. (7) refers to an antisymmetrized product for electrons, but the later discussion of H2 in Section 5 mentions a symmetrized product for the two electronic walkers; the treatment of spin statistics in the one-dimensional model should be clarified.","section":"Section 2, Eq. (7)"},{"comment":"The captions for Figures 2-5 do not specify the color/line conventions in a way that is self-explanatory without the main text, and Figure 1's caption lacks a description of the pulse envelope that is supposedly shown.","section":"Figure captions"},{"comment":"The energy expression in Eq. (28) requires an estimate of the many-body quantum potential using adaptive kernel density estimation, but no specific algorithm or bandwidth choice is given for the multivariate case, making the reported ground-state energy difficult to reproduce.","section":"Section 4"}],"recommendation":"reject","confidential_remarks":"The manuscript relies heavily on the author's prior TDQMC publications (refs 24-26) and does not clearly delineate the novel contributions of this work relative to those papers. The central derivation defect and the use of an approximate benchmark as 'exact' raise concerns about the validity of the main claims; these issues are not minor and would require a fundamental reworking of the paper's theoretical basis and numerical validation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one with a red pen on Section 2. The genuinely new piece is the quantum-quantum extension of your earlier TDQMC work: coupled ensembles of nuclear and electronic walkers, with guiding waves for both, plus the classical-nucleus limit. The 1D H2 test is sensible, and the ground-state energy (-1.661 a.u.) sits close to the direct Ehrenfest number (-1.665), with the time-dependent ionization within 5%. That part is fine, and the limiting-case discussion around Eqs. (24)-(25) is clear and honest.\n\nThe problem is the derivation of the effective potentials, Eqs. (16)-(19). The paper says they come from substituting the smoothed walker densities (12)-(13) into the Hartree equations (5)-(6). That is not what substitution gives. You would get a convolution of the interaction potential with the Gaussian mixture; what the paper writes is the bare potential evaluated at the walker separation, multiplied by a Gaussian in that same separation. These are different expressions, and no limit, quadrature, or other controlled approximation is offered to connect them. Since those potentials are the sole place where correlations beyond Hartree enter, the 'ab initio' claim does not follow. The benchmark itself is approximate—Ehrenfest MD keeps nuclei classical—so agreement with it can't rescue the derivation. And there are no convergence tests in the number of walkers, grid spacing, or time step, no code, and the soft-core parameters and initial widths are chosen without a sensitivity study.\n\nThat is a load-bearing flaw, not a cosmetic one. But it's a fixable framing issue as much as a mathematical one. If the author stops calling the effective potentials a consequence of the Hartree equations and instead presents them as an independent approximation—a trajectory-based, mean-field-with-nonlocal-corrections scheme—then the 5% agreement becomes a plausible demonstration of the heuristic's behavior, not a proof of ab initio status. A serious referee should ask for that reframing, plus convergence data and at least one stronger benchmark.\n\nFor whom: strong-field/attosecond chemists interested in trajectory-based methods. This deserves peer review rather than desk rejection—the extension is real, and the numerical comparison is informative—but it needs major revision before publication.","headline":"A real extension of TDQMC to nuclei, but the effective-potential step is not a derivation—'ab initio' is unsupported and needs reframing.","tokens_in":12578,"tokens_out":3965,"would_cite":false,"duration_ms":37346,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.25.-v","02.70.Ss"],"model":"deepseek-v4-flash","headline":"This paper argues that molecular dynamics can be solved ab initio by time-dependent quantum Monte Carlo, with both electrons and nuclei represented by guided walker ensembles, and demonstrates about 5 percent agreement with exact…","keywords":["time-dependent quantum Monte Carlo","molecular dynamics","Bohmian trajectories","guiding waves","quantum correlations","strong-field ionization","hydrogen molecule","Ehrenfest molecular dynamics"],"falsifier":"A decisive test would be to solve the same one-dimensional two-electron $\\mathrm{H}_2$ model, with the same softened Coulomb potentials and the same 335 nm pulse at $9\\times10^{14}$ W/cm$^2$, by a numerically exact time-dependent Schrödinger calculation and compare the survival probability and internuclear distance over the full pulse; if the exact curves differ from TDQMC by much more than the claimed 5 percent, the effective-potential approximation has failed. A second check is convergence in the number of walkers: if increasing $M$ systematically changes the ionization dynamics by more than the statistical error, the result is not converged and the ab initio claim is not met.","tokens_in":11553,"feed_emoji":"⚫️","tokens_out":7547,"duration_ms":69040,"temperature":0.7,"pith_summary":"This paper proposes an ab initio molecular-dynamics method, time-dependent quantum Monte Carlo (TDQMC), that treats both electrons and nuclei quantum mechanically. Each physical particle is replaced by an ensemble of Monte Carlo walkers and guiding waves in physical space, with the walkers propagated by de Broglie-Bohm velocities and the guiding waves by coupled linear Schrödinger equations. Nonlocal quantum correlations are introduced through effective potentials obtained by adaptive kernel smoothing of the walker densities, avoiding both quantum potentials and explicit exchange-correlation integrals. Applied to one- and two-dimensional $\\mathrm{H}_2$ in intense ultrashort pulses, TDQMC gives ionization and internuclear-distance evolution within about 5 percent of exact Ehrenfest molecular dynamics, while time-dependent Hartree-Fock deviates by about 20 percent and does not dissociate the molecule.","feed_headline":"Walker-based dynamics tracks laser-driven H2 within 5 percent","feed_subtitle":"Both electrons and nuclei move as quantum walkers, capturing correlations that mean-field Hartree-Fock misses.","key_machinery":"The central object is the walker–guiding-wave pair. For each electron $i$ and nucleus $I$, an ensemble of $M$ walkers with positions $\\mathbf{r}_i^{(k)}$ and $\\mathbf{R}_I^{(k)}$ moves under the de Broglie-Bohm guiding equations (8)–(9) and a random drift, while separate guiding waves $\\phi_i^{(k)}$ and $\\Phi_I^{(k)}$ satisfy the coupled linear Schrödinger equations (14)–(15). The exact densities in the Hartree equations are replaced by kernel density estimates (12)–(13), and the resulting effective potentials (16)–(19) encode nonlocal correlations with adaptive bandwidths $\\sigma$ and $\\Sigma$ set by formulas (22)–(23). The key move is that no quantum potential is computed and correlated interactions are reduced to Monte Carlo sums over walkers, which replaces the costly integrals of multi-configuration methods.","core_discovery":"The paper's central claim is that a fully quantum treatment of molecular dynamics is possible with TDQMC, in which both electrons and nuclei are described on equal footing by ensembles of walkers and guiding waves in physical space. The guiding waves obey coupled linear Schrödinger equations, and the walkers move by the de Broglie-Bohm guidance equations plus a random diffusive drift, so that the walker density reproduces the time-dependent probability density of each particle. Quantum correlations, including nonlocal effects, are encoded through effective potentials built from adaptive kernel density estimates of the walker distributions, with no quantum potentials and no overlap, exchange, or correlation integrals. The author shows for a one-dimensional hydrogen molecule in a strong 335 nm pulse that the TDQMC ionization and internuclear-distance evolution stay within about 5 percent of the exact Ehrenfest molecular dynamics benchmark, whereas time-dependent Hartree-Fock underestimates ionization by about 20 percent and does not dissociate the molecule.","pith_inferences":["Editorial inference: the same effective-potential construction could be tested on exactly solvable two-electron models with known correlation energies, such as a Hooke's-law atom, to isolate what the adaptive kernel smoothing adds beyond mean-field theory.","Editorial inference: because the method avoids overlap and exchange integrals, its scaling in particle number is likely its main advantage; a quantitative complexity study for three or more electrons would show whether the approach is practical beyond the two-electron molecule demonstrated here.","Editorial inference: the 5 percent figure is established for a softened one-dimensional model; extension to three dimensions would need to show that the adaptive bandwidth formulas remain stable when the Coulomb potentials are not softened and the walker density varies in all directions.","Editorial inference: the quantum-quantum and quantum-classical variants suggest a natural place to add trajectory surface hopping or ring-polymer corrections to restore nuclear quantum effects beyond the de Broglie-Bohm guidance, though the paper does not pursue this."],"forward_implications":["Correlated electron–nuclear dynamics can be propagated in physical space without solving the full many-body Schrödinger equation, using explicit Coulomb potentials and no exchange-correlation integrals.","The quantum-classical limit gives a Newtonian nuclear trajectory bundle with self-consistent back-reaction, recovering previous mixed quantum-classical methods and the Ehrenfest mean-field limit as special cases.","For the one-dimensional $\\mathrm{H}_2$ benchmark, TDQMC predicts final ionization only about 5 percent above exact Ehrenfest, whereas time-dependent Hartree-Fock is about 20 percent below and fails to dissociate the molecule.","The nonlocal correlation lengths $\\sigma$ and $\\Sigma$ grow as the molecule ionizes, providing a dynamical diagnostic of how electron–electron and electron–nuclear correlations spread during dissociation.","The two-dimensional orientation simulation indicates the method can follow molecular rotation and dissociation in an external field, not just fixed-nucleus electron dynamics."],"supporting_citations":[{"why":"Supplies the exact Ehrenfest molecular dynamics benchmark against which TDQMC ionization and internuclear-distance evolution are compared.","marker":"[20]"},{"why":"Introduces the time-dependent quantum Monte Carlo method with walkers and guiding waves that this paper extends to nuclear degrees of freedom.","marker":"[24]"},{"why":"Establishes walker propagation, thermalization by importance sampling, and adaptive kernel bandwidth formulas used in the correlated ensembles.","marker":"[25]"},{"why":"Provides kernel density estimation of particle densities for constructing the nonlocal effective potentials.","marker":"[26]"},{"why":"Source of the de Broglie-Bohm guidance equations used to move the walkers and of the limiting interpretation for vanishing quantum potential.","marker":"[32]"},{"why":"Supplies the adaptive kernel bandwidth formula underlying equations (22)–(23) for the electron and nuclear correlation lengths.","marker":"[33]"},{"why":"Supports the adaptive bandwidth selection and the relation of correlation lengths to walker covariance matrices.","marker":"[34]"},{"why":"Provides the softened Coulomb potentials used in the one-dimensional $\\mathrm{H}_2$ simulations to avoid singularities at the nuclei.","marker":"[39]"}],"fun_headline_variants":["Walkers simulate laser-driven H2 within 5 percent","Quantum walkers capture molecular dynamics accurately","TDQMC: electrons and nuclei as quantum walkers","Laser pulse dissociates H2 in walker-based simulation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assumption that replacing the true electron and nuclear densities in the Hartree equations by broadened, walker-sampled densities with adaptive widths faithfully encodes the quantum correlations that the simple product wavefunction leaves out; if that smoothing distorts the correlations, the claimed agreement with exact Ehrenfest dynamics would not be guaranteed and the ab initio claim would fail.","fun_headline_variants_meta":{"raw":{"variants":["Walkers simulate laser-driven H2 within 5 percent","Quantum walkers capture molecular dynamics accurately","TDQMC: electrons and nuclei as quantum walkers","Laser pulse dissociates H2 in walker-based simulation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000608,"raw_usage":{"total_tokens":2796,"prompt_tokens":871,"completion_tokens":1925,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":487,"completion_tokens_details":{"reasoning_tokens":1860}},"tokens_in":487,"tokens_out":1925,"duration_ms":14076,"temperature":1.0,"reasoning_tokens":1860,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:42:30.232924+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to solve the same one-dimensional two-electron $\\mathrm{H}_2$ model, with the same softened Coulomb potentials and the same 335 nm pulse at $9\\times10^{14}$ W/cm$^2$, by a numerically exact time-dependent Schrödinger calculation and compare the survival probability and internuclear distance over the full pulse; if the exact curves differ from TDQMC by much more than the claimed 5 percent, the effective-potential approximation has failed. A second check is convergence in the number of walkers: if increasing $M$ systematically changes the ionization dynamics by more than the statistical error, the result is not converged and the ab initio claim is not met.","supporting_citations":[],"review_version":1}