{"id":"bfe1d066-b9ca-4f3d-be4c-d19ea2eaaf9b","arxiv_id":"2501.16137","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"A screened-potential variant of time-dependent quantum Monte Carlo reproduces exact one-dimensional helium dynamics, but the polynomial-scaling claim is supported only by limited tests.","lead":"This paper extends a time-dependent quantum Monte Carlo scheme to fermionic atoms by replacing explicit Slater determinants with screened Coulomb potentials whose exchange-hole size is set by hand. For a one-dimensional helium model, the method matches exact dipole and ionization curves, but the headline claim of polynomial-time scaling rests on a small scaling test rather than a proof.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fitted exchange-hole parameter rs is not derived or shown transferable, and no convergence study supports the controlled-accuracy polynomial-time claim.","rationale":"The paper's central claim is that TDQMC reduces exponential many-body cost to polynomial time with controlled accuracy. For this to hold, two things are required: (i) a representation whose error can be systematically reduced by increasing computational resources, and (ii) a proof or demonstration that the required resources stay polynomial in system size and accuracy. Neither is provided. The fermionic representation of Eq. (8) discards antisymmetrization and instead inserts a fixed screened Coulomb interaction (Eq. (10)) with hole size rs. The value rs = 0.67 a.u. is selected solely to match the exact ground-state energy of the single 1D ortho-helium test case; it is not derived from the local walker density or from the antisymmetric nodal structure. Consequently, the dynamics in Figs. 4-5 could simply reflect a refitted interaction rather than a successful approximation to the full many-body wavefunction. The reported scaling tests are not a complexity analysis: they measure parallel speed-up for a fixed 1D problem, not the growth of total cost with N, M, or required error tolerance, and the statement that Monte Carlo sums can be truncated is not substantiated. The paper is clear that this is an approximate method, but 'controlled accuracy' is precisely what is not shown. A transferability test on a second few-electron system, or a walker-convergence study, would settle whether the central claim is supportable. As it stands, the rejection with moderate confidence is appropriate; no change to the reader's verdict is warranted.","tokens_in":9192,"tokens_out":6072,"duration_ms":57132,"concrete_test":"Apply TDQMC to a second two-electron system, e.g., Li+ in the same 1D model, using rs = 0.67 a.u. and also rs re-optimized to reproduce the exact ground-state energy of Li+. Compare the predicted strong-field dipole and ionization with exact TDSE results. If the re-optimized rs changes by more than ~10% from 0.67 a.u., or if the fixed-rs dynamics deviate substantially from exact results, the exchange-hole parameter is not transferable and the central claim of a controlled-accuracy polynomial-time fermionic method is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of polynomial-time, controlled-accuracy fermionic dynamics rests on two replacements: Eq. (8) substitutes the antisymmetrized many-body wavefunction with a simple product of single-particle guide waves, and Eq. (10) encodes all exchange physics through a static screened Coulomb potential with one hole-size parameter rs. In the Results section, rs = 0.67 a.u. is selected so that the TDQMC ground-state energy of 1D ortho-helium reproduces the exact diagonalization result, and the same rs is then used for the time-dependent response. No derivation of rs from first principles is offered, and no evidence shows that this parameter is transferable across systems, field strengths, or observables. Consequently, the agreement in Figs. 4-5 for ortho-helium may reflect a refitted effective interaction rather than a robust approximation to the full many-body wavefunction. The scaling discussion in the Conclusions is a heuristic strategy statement, not a complexity analysis; the reported 128-512 processor test in one dimension measures parallel speed-up for a fixed two-electron problem, not the growth of total computational cost with particle number, walker number, or desired accuracy. Thus the load-bearing premise that a finite set of single-particle guide waves and walkers with a single fixed rs yields a controlled-accuracy polynomial-time method for fermions is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a time-dependent quantum Monte Carlo (TDQMC) method in which ensembles of Monte Carlo walkers and single-particle guide waves in physical space evolve concurrently, with the electron-electron interaction replaced by an effective potential that includes a same-spin screened Coulomb term intended to account for exchange. The method is applied to one-dimensional para- and ortho-helium, computing ground-state energies, electron-pair density functions, and the dipole and ionization response to an intense few-femtosecond pulse. The authors claim that this approach reduces the exponential scaling of the many-body problem to polynomial scaling while retaining controlled accuracy for fermions.","tokens_in":9426,"tokens_out":6694,"duration_ms":61610,"significance":"If the polynomial-scaling and controlled-accuracy claims were established, the method would be a significant practical tool for strong-field multi-electron dynamics, an area where exact grid methods scale exponentially and TDDFT suffers from exchange-correlation inaccuracies. The paper ships concrete comparisons against exact two-dimensional TDSE solutions for two-electron one-dimensional atoms, which is a valuable testbed and a strength. However, the central claims rest on a tunable exchange-hole parameter rs and on a scaling discussion that is not backed by convergence or complexity measurements. As presented, the evidence supports the method only as a promising heuristic that needs further validation, and the strong claims in the abstract and conclusions are not yet justified.","major_comments":[{"comment":"The exchange-hole size rs is not derived from first principles; it is chosen as rs=0.67 a.u. so that the TDQMC ground-state energy of ortho-helium reproduces the exact diagonalization value. The same fitted rs is then used for the time-dependent dipole and ionization curves in Fig. 5. Consequently, the agreement in Fig. 5 is an in-sample consistency result, not an independent prediction of the exchange model. The statement that rs 'can be estimated from the corresponding TDHF solution' is not carried out, and no transferability tests (e.g., other intensities, wavelengths, or atoms) are given. This is load-bearing because it is the only evidence that Eq. (10) correctly captures fermionic exchange dynamics. The authors should either derive rs from TDHF or a similar first-principles procedure, or demonstrate that a single rs predicts multiple observables and conditions out of sample.","section":"Screened potential formulation, Eq. (10); Results and discussion (rs fitting)"},{"comment":"The central claim of a 'controlled-accuracy approximate solution which scales as a low-order polynomial in time' is not supported by the presented evidence. The parallel test is a fixed two-electron problem in one spatial dimension, reporting only an increase in processor count from 128 to 512 with a statement of 'lower than quadratic' scaling in processors. There is no scaling study with respect to walker number M, particle number N, or desired accuracy, and no convergence study of the observables with respect to M, kernel bandwidth, time step, or complex-time step. The phrase 'controlled accuracy' is never given an error bound or a convergence threshold. As it stands, the scaling claim is a heuristic strategy statement, not a measured result. The authors should present scaling data in terms of total computational cost versus N, M, and error tolerance.","section":"Results and discussion (parallel scaling); Conclusions"},{"comment":"The ground-state calculations use a=1 a.u. and b=0.5 a.u. in Eqs. (11) and (12), while the time-dependent simulations in Figs. 4 and 5 use a=1 a.u. and b=1.5 a.u. This changes the Hamiltonian between the preparation stage and the dynamics. The manuscript does not state whether the TDQMC ground state is re-prepared with b=1.5 before the laser pulse is applied, nor whether the exact TDSE reference also uses b=1.5. If the b values differ, the comparison in Figs. 4 and 5 may not be between identical physical systems. The authors should justify the parameter change or use one consistent parameter set for both ground state and dynamics.","section":"Results and discussion (parameters a and b)"},{"comment":"The replacement of the antisymmetrized many-body wavefunction by a simple product in Eq. (8) removes the explicit antisymmetric nodal structure, and all fermionic exchange is relegated to the same-spin screened potential in Eq. (10). The paper does not provide a derivation or a quantitative error estimate for this replacement. In particular, the statement in the Conclusions that 'all the information about the nodal regions in case of fermions is carried by the evolving guide waves' is difficult to reconcile with the product form, which has no exchange-induced sign change. Even with the Gram-Schmidt orthogonalization mentioned for ortho-helium, the guide waves are not antisymmetrized. A more careful justification of this central approximation is needed, or at least a numerical comparison between the product form with screening and a full Slater-determinant version of the method.","section":"Screened potential formulation, Eq. (8) and Conclusions"}],"minor_comments":[{"comment":"The estimate 'workload proportional to M3K' is unclear; presumably it should be K^{3N} for N particles on a grid with K points per dimension, and the symbol M is later reused for the number of walkers. Please clarify the notation.","section":"Introduction, scaling estimate"},{"comment":"The screened potential expression involving 'erf' is garbled in the typeset equation. Please define the modified e-e potential explicitly, e.g., V(r) erfc(r / rs) / r or an equivalent form, and ensure that Eq. (10) matches the description in the text.","section":"Eq. (10) and surrounding text"},{"comment":"The caption contains a typo: 'exact resut' should be 'exact result'.","section":"Figure 5 caption"},{"comment":"The adaptive kernel exponent α=0.2 is introduced without justification, and no sensitivity analysis is provided for α or the pilot bandwidth σ. Since the effective potential in Eq. (6) depends on these choices, a brief study of their influence on the reported observables would strengthen the claims.","section":"Results and discussion, adaptive kernel"},{"comment":"No statistical error bars are reported for the Monte Carlo observables (dipole moment, ionization probability, pair density). For a stochastic method, error bars are essential to assess whether the agreement with the exact results is meaningful.","section":"Results and discussion, Monte Carlo statistics"},{"comment":"Reference [13] contains five separate papers; it would be helpful to cite the specific paper that describes the complex-time ground-state preparation procedure and the random component in the walker motion.","section":"References"},{"comment":"The statement that known quantum Monte Carlo techniques 'do not allow treatment of time-dependent processes' is too strong; time-dependent diffusion Monte Carlo and path-integral approaches exist. A more nuanced statement with appropriate citations would be more accurate.","section":"Introduction, QMC statement"},{"comment":"The sentence introducing the pair density function contains a stray 'e.g. 3'; this appears to be a missing citation or equation and should be corrected.","section":"Results and discussion, pair density paragraph"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a single-author paper that builds on a series of prior TDQMC papers by the same author, but the novelty relative to those earlier works is not clearly delineated. The strongest concern, reflected in the major comments, is that the central polynomial-scaling claim and the fermionic exchange model are validated only through a circular fitting procedure and a scaling test that is not a complexity measurement. These issues could in principle be addressed by additional numerical experiments, so I do not recommend rejection at this stage, but I cannot accept the paper in its current form. I would also encourage the editor to ensure that the claimed 'controlled accuracy' is defined operationally before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, the screened-exchange-hole formulation (Eq. 10) is a real modification of the author's earlier TDQMC work, and the 1D ortho-helium dynamics with rs=0.67 a.u. do reproduce the exact dipole and ionization curves quite well. Second, the headline claim of polynomial-time controlled-accuracy fermionic dynamics is not backed by the evidence. There is no complexity analysis, no walker-convergence study, and the only scaling data is a 128-to-512-processor weak-scaling test on a fixed two-electron problem.\n\nWhat is genuinely new and good: replacing Slater determinants with a product of guide waves plus a static screened Coulomb potential restricted to same-spin walkers is a clean idea that simplifies the guiding equations and avoids the sign problem in a pragmatic way. The para-helium case, where no screening is used, gives a good match to exact strong-field ionization. The ground-state pair density functions also show the expected Coulomb and exchange features. For a one-dimensional model, the method performs well as an approximate heuristic.\n\nThe soft spots are substantial. rs is chosen so that the TDQMC ground-state energy of ortho-helium matches exact diagonalization. That same rs is then used for the time-dependent response, so the agreement in Figs. 4-5 is not an independent prediction. No evidence is given that rs transfers to other systems, field strengths, or observables. The scaling discussion in the conclusions is a strategy statement, not a quantitative argument; the parallel speed-up measurement says nothing about how the total cost grows with electron number, walker number, or desired accuracy. The paper also never reports error bars or a sensitivity study for the softening parameters a and b, the kernel exponent alpha, or the walker count. Without those, \"controlled accuracy\" is an assertion, not a demonstrated property.\n\nI think the author is honestly engaging with the literature, and the method may be useful for qualitative strong-field simulations after more work. But the paper's load-bearing claim is the polynomial scaling, and that claim is unsupported. I would recommend reject in present form, with a clear path to resubmission: a systematic walker-convergence study, a derivation or at least a scaling argument that separates errors from walker number and time step, and a 3D few-electron demonstration with rs determined without referencing the target observable. If the venue is willing to referee a method paper that needs major revision, send it; a good referee can articulate exactly what is missing and the new formulation deserves that input.","headline":"Genuine screened-exchange-hole extension of TDQMC with a nice 1D helium demo, but the polynomial-time scaling claim is not supported and the key fermionic parameter is fitted, so the paper as written does not establish its central thesis.","tokens_in":9976,"tokens_out":2300,"would_cite":false,"duration_ms":25916,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims a hybrid wave–walker Monte Carlo scheme cuts exponential many-body scaling to polynomial time.","keywords":["time-dependent quantum Monte Carlo","guide waves","exchange hole","screened Coulomb potential","electron-pair density","strong-field ionization","helium atom","polynomial-time scaling"],"falsifier":"Fix $r_s=0.67$ a.u. using only the ortho-helium ground-state energy, then compute the dipole response and ionization for a range of pulse intensities and wavelengths and compare with exact TDSE results; the central claim survives only if those predictions stay accurate without refitting $r_s$ for each observable.","tokens_in":8895,"feed_emoji":"⚛️","tokens_out":9932,"duration_ms":86744,"temperature":0.7,"pith_summary":"This paper claims that a time-dependent quantum Monte Carlo method can reduce the exponential computational cost of the quantum many-body problem to polynomial time by replacing the configuration-space wavefunction with a set of single-particle guide waves and moving sample points (walkers) in physical space. Each guide wave obeys a coupled Schrödinger equation, each walker follows a first-order guiding equation, and the electron–electron interaction is built from Gaussian-smoothed walker distributions whose widths set the nonlocal correlation length. For fermionic states, exchange is encoded by screening the Coulomb repulsion with an error-function factor of width $r_s$ restricted to same-spin walkers, so the exchange hole is explicit and no fermion sign problem arises. For one-dimensional para- and ortho-helium, the method reproduces the ground-state electron-pair density and gives dipole moments and ionization probabilities close to exact TDSE results under a strong ultrashort pulse. The author concludes that this is a controlled-accuracy approximate solution whose cost scales as a low-order polynomial in time for both bosons and fermions.","feed_headline":"Many-body quantum cost drops to polynomial time in a guided Monte Carlo scheme","feed_subtitle":"Guide waves plus walkers match exact helium dipole and ionization under strong laser pulses.","key_machinery":"The central machinery is the coupled pair of evolution equations for guide waves and walkers. Each walker moves by a first-order de Broglie–Bohm guiding equation evaluated from its own guide wave, while each guide wave evolves by a Schrödinger equation in physical space whose electron–electron potential is a Monte Carlo sum over Gaussian-smoothed walker distributions. The smoothing width is the nonlocal correlation length, set by adaptive kernel density estimation. For fermions the Coulomb interaction is multiplied by $\\mathrm{erf}(|\\mathbf r_i-\\mathbf r_j|/r_s)$ restricted to same-spin walkers, which opens an exchange hole around each electron and removes the Coulomb cusp; $r_s$ is the size of that hole. This lets a few thousand walkers per electron sample the probability distributions provided by the guide waves, which is what keeps the computational cost polynomial.","core_discovery":"The central claim is that all of the information carried by the many-body wavefunction in configuration space can be carried, for computational purposes, by finite ensembles of single-particle guide waves in physical space coupled through effective potentials. The paper replaces the Slater-determinant structure for fermions by a simple product of guide waves and modifies the electron–electron Coulomb potential by a screening factor $\\mathrm{erf}(|\\mathbf r_i-\\mathbf r_j|/r_s)$ with a Kronecker-delta restriction to same-spin walkers. With $r_s=0.67$ a.u., the ortho-helium ground-state energy comes within 5% of the exact value, and the time-dependent dipole and ionization in a strong laser pulse agree with the exact TDSE results, while setting $r_s=0$ changes the response qualitatively. The paper takes this as evidence that dynamic exchange effects are captured and that the method provides a polynomial-time approximate solution for both bosonic and fermionic systems.","pith_inferences":["A decisive stress test the paper leaves implicit is to fix $r_s$ once from a ground-state energy and then predict several independent strong-field observables for the same atom; if each observable requires its own $r_s$, the method is a fitting scheme rather than a predictive theory.","The exchange-hole size is treated as a fixed constant, but during strong ionization the electron cloud is strongly deformed; allowing $r_s$ to become a dynamical, locally density-dependent quantity is a direct extension that could be tested against the same exact TDSE results.","If the polynomial scaling survives beyond two electrons, the method would open a class of classically accessible correlated multi-electron strong-field problems where exact TDSE is unavailable; benchmarking against beryllium or small molecules would settle that question."],"forward_implications":["Correlated multi-electron dynamics in strong fields would become tractable on classical hardware, with cost growing as a low-order polynomial in time rather than exponentially.","Fermionic exchange is essential in this model: the ortho-helium dipole and ionization match exact results with hole size $r_s=0.67$ a.u., while $r_s=0$ gives a qualitatively different response.","The walker ensembles give direct access to statistical observables such as the electron-pair density without computing configuration-space integrals.","Parallel scaling tests reported for one dimension find work below quadratic in the number of walkers and close to linear in the number of electrons."],"supporting_citations":[{"why":"Defines the TDQMC algorithm of concurrent propagation of walker ensembles and guide waves in complex then real time, on which the whole paper builds.","marker":"[13]"},{"why":"Supplies the first-order de Broglie–Bohm guiding equation used to move the walkers.","marker":"[14, 15]"},{"why":"Provides the kernel density estimation used to smooth walker distributions and set the nonlocal correlation lengths in the effective potentials.","marker":"[19, 20]"},{"why":"Proposes screened Coulomb potentials to represent fermion exchange, adapted here as the error-function exchange hole of size $r_s$.","marker":"[21]"},{"why":"Provides the soft-core one-dimensional potentials used to avoid the Coulomb singularity in the helium model.","marker":"[25]"},{"why":"Supplies accurate helium pair-density results used to validate the walker-derived electron-pair density function.","marker":"[26]"},{"why":"Documents the exponential-cost fermion sign problem that the paper claims to bypass with positive-definite walker densities and guide waves.","marker":"[27]"}],"fun_headline_variants":["Quantum dynamics speedup: polynomial time with guided Monte Carlo","Polynomial-time quantum dynamics via guided Monte Carlo","Guided Monte Carlo drops quantum many-body cost to polynomial","Exchange-hole Monte Carlo achieves polynomial-time quantum dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a finite set of single-particle guide waves in physical space, coupled through effective potentials built from smoothed walker distributions, captures the full many-body dynamics, with all fermion exchange represented by one fixed hole size $r_s$ tuned to the ground-state energy.","fun_headline_variants_meta":{"raw":{"variants":["Quantum dynamics speedup: polynomial time with guided Monte Carlo","Polynomial-time quantum dynamics via guided Monte Carlo","Guided Monte Carlo drops quantum many-body cost to polynomial","Exchange-hole Monte Carlo achieves polynomial-time quantum dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00053,"raw_usage":{"total_tokens":2527,"prompt_tokens":890,"completion_tokens":1637,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":1574}},"tokens_in":506,"tokens_out":1637,"duration_ms":10291,"temperature":1.0,"reasoning_tokens":1574,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:42:14.530538+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fix $r_s=0.67$ a.u. using only the ortho-helium ground-state energy, then compute the dipole response and ionization for a range of pulse intensities and wavelengths and compare with exact TDSE results; the central claim survives only if those predictions stay accurate without refitting $r_s$ for each observable.","supporting_citations":[],"review_version":1}