{"id":"6ba3a626-6cf6-4fce-9079-c36976437cee","arxiv_id":"2506.08230","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A Wigner-Weyl perturbative expansion around momentum-dependent electronic states recovers exact coupled nuclear-electronic eigenvalues and improves model vibrational energies.","lead":"This paper derives a perturbation theory that starts from electronic states depending on both nuclear position and momentum, instead of the standard Born-Oppenheimer states, and corrects them step by step toward the exact quantum energy. The framework gives phase-space electronic structure calculations a formal justification and improves computed vibrational energies on a model hydrogen bond.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central exactness claim rests on the convergence of the perturbative series in Eq. 91; the paper enforces unitarity and diagonality only order by order, and near degeneracies the denominators in Eq. 107 (acknowledged in Sec. 6) prevent convergence, so 'recover exact energies' is not established.","rationale":"I checked the algebra of the central first-order derivation. The expansion of L^dagger_W * H_W * L_W around H^PS_W is consistent: the difference term i hbar Gamma P / M contributes only its direct term at O(hbar), the unitarity condition fixes L^(1)_W via Eq. 97, and the diagonality condition yields A_ij = -B_ij / (Lambda_ii - Lambda_jj) and the first-order diagonal correction hbar B00. This part is sound. The numerical results are also credible evidence that the first-order correction helps in the well-separated model. The load-bearing weakness is the step from a formal order-by-order series to the claim of exactness. The series in Eq. 91 is never shown to converge or to define a genuine unitary operator after resummation; the denominators in Eq. 107 make the construction singular at degeneracy, and the authors themselves limit the validity to well-separated ground states in Sec. 6. Because the abstract and conclusions use 'exact' and 'rigorous framework,' the paper should either prove convergence or summability for the excluded cases or explicitly frame the result as an asymptotic expansion. The reader's weakest assumption identified the same issue, and the conditional verdict is appropriate.","tokens_in":19960,"tokens_out":19297,"duration_ms":235023,"concrete_test":"For a two-surface linear vibronic coupling model with an exactly solvable nuclear Hamiltonian and tunable electronic gap Delta, compute the phase-space series terms E^(0), E^(1), and the full E^(2) from Eq. 108 (not just the Gamma^2 term Z00) at fixed reduced mass, and compare with exact diagonalization as Delta is decreased. If the sequence of partial sums does not approach the exact energy for small Delta, or if including E^(2) increases the error, the exactness claim fails in the regime the paper excludes; if the series converges for Delta above a finite threshold, the conditional claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The formal claim is that the order-by-order Moyal construction of J_W in Eq. 92 sums to a unitary operator that diagonalizes the full Hamiltonian and thereby recovers exact eigenvalues. What is actually shown is weaker: Eqs. 94 and 104 enforce unitarity and block diagonality order by order in hbar, with the first-order generator L^(1)_W explicitly containing denominators (Lambda^PS_kk - Lambda^PS_jj) in Eq. 107. No convergence, summability, or error bound for the infinite series in Eq. 91 is provided. For near-degenerate or multiconfigurational electronic states these denominators are small or zero, and the expansion is not a controlled perturbation; the authors concede this limitation in Sec. 6. Since Littlejohn-Flynn-type expansions are generically asymptotic rather than convergent in hbar, the phrase 'recover exact quantum energies' overstates what the derivation supports: the paper demonstrates an asymptotic correction scheme that improves energies when electronic gaps are large, not a route to exact eigenvalues in general. The numerical validation does not resolve this, because the model has a well-separated ground surface and only first-order (plus a truncated second-order) terms are tested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Wigner–Weyl/Moyal perturbative framework for correcting phase-space electronic Hamiltonian eigenstates so as to converge toward the exact eigenvalues of a coupled nuclear-electronic Hamiltonian. The authors first review the Born–Oppenheimer and Littlejohn–Flynn expansions, then construct a unitary transformation J_W = L_W + ℏ L_W^(1) + ... that diagonalizes the full Hamiltonian order by order in ℏ, deriving a first-order energy correction ℏB00 (Eq. 110) and a truncated second-order term (Eq. 111). The formalism is tested on a one-dimensional hydrogen-bond model, with numerical results showing that the first-order correction improves vibrational energy gaps by roughly an order of magnitude relative to zeroth-order phase-space and Born–Oppenheimer results. The central algebraic derivation is explicit and does not involve fitting of final energies, but the claimed exactness rests on an order-by-order expansion whose convergence is not established, and the numerical validation is limited to a single well-separated model with a partially implemented second-order term.","tokens_in":20253,"tokens_out":4963,"duration_ms":62276,"significance":"If the formal expansion is understood as an asymptotic series that can be controlled in favorable regimes, this is a useful and conceptually important contribution: it provides a systematic way to improve phase-space electronic structure calculations and partially answers the criticism that such Hamiltonians are ad hoc. The derivation is independent of the specific form of Γ, which is a genuine strength, and the first-order energy correction has a simple closed form. The numerical results demonstrate a practically meaningful improvement over both Born-Oppenheimer and zeroth-order phase-space approaches. However, the significance is tempered by the absence of convergence or error estimates, the use of a single model with a well-separated ground state, and the fact that the reported second-order results do not implement the full second-order theory.","major_comments":[{"comment":"The central claim that the construction recovers exact eigenvalues is supported only order by order in ℏ. The unitarity condition (Eq. 94) and block-diagonality condition (Eq. 104) fix the perturbative generators, but no convergence, summability, or error bound is provided for the infinite series in Eq. 91, whose first-order generator contains the energy-gap denominators (Λ_k^PS - Λ_j^PS) in Eq. 107. The acknowledgment in Sec. 6 that the expansion is perturbative and requires a well-separated ground electronic state means that the phrase \"recover exact quantum eigenvalues\" is stronger than what is demonstrated. Please either add a remainder estimate or a rigorous statement of the asymptotic nature of the series, or temper the title/abstract wording accordingly.","section":"Sec. 3.1–3.1.2, Eq. (91) and Eq. (107)"},{"comment":"The numerical E^(2)_PS2 results are based on a \"very naive\" second-order correction that keeps only (L† Γ² L / 2M)_00 and drops all other terms of Eq. 108, including all terms involving L^(1)_W. Consequently, the observed accuracy and µ-scaling of E^(2)_PS2 in Figs. 2, 7, 10, and 11 do not validate the full second-order theory. Please either implement the complete second-order correction or present E^(2)_PS2 solely as an ad hoc benchmark and base the numerical case for the perturbative framework on the first-order results.","section":"Sec. 3.3, Eq. (111)"},{"comment":"The Gaussian width σ appearing in the Γ operator (Eq. 123) is never specified. Since all phase-space calculations use this operator, the numerical results in Figs. 1–11 cannot be reproduced without knowing σ. Please report the value used and, ideally, examine the sensitivity of the first-order correction to σ.","section":"Sec. 4.1, Eqs. (120)–(123), Table 3"}],"minor_comments":[{"comment":"The legends in Figs. 2 and 11 contain garbled labels such as \"ΔE(0) ̃S1\" and \"ΔE(2) ̃S2\"; these should be corrected to the notation used in Table 2.","section":"Figs. 2 and 11"},{"comment":"There is a typo in Eq. (103): \"muts\" should be \"must\".","section":"Eq. (103)"},{"comment":"Refs. 43 and 45 appear to be the same paper (Marinica, Gaigeot, and Borgis) and should be merged.","section":"References"},{"comment":"Please state the basis size or grid parameters used for the exact diagonalization reference, and add a data/code availability statement, so that the numerical results can be independently reproduced.","section":"Sec. 4.1"},{"comment":"The boldface notation for nuclear operators is not consistently rendered in the typeset equations; for example, Eq. (85) mixes P as an operator and as a phase-space variable. Please ensure consistent formatting.","section":"Notation, Sec. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is heavily self-referential (Refs. 15, 16, 24, 25), which is natural for a group developing its own theoretical framework but should be kept in mind when assessing novelty. The main technical derivation appears sound and the first-order correction is a useful contribution. The major concern is that the title and abstract overstate the exactness result relative to the order-by-order construction, and the missing σ parameter prevents reproduction of the numerical study. These issues are fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe key result is Eqs. 91–110: a first-order Moyal/Littlejohn-Flynn correction built on top of a phase-space electronic Hamiltonian. Unlike the BO case, where the first-order diagonal correction vanishes, here you get a nonzero hbar B00 term, and the paper shows numerically on the Borgis model that this buys roughly an order of magnitude in vibrational energy accuracy. The construction is genuinely new as far as I can tell — previous phase-space work was empirical, and Littlejohn-Flynn theory targeted a BO reference. The algebra in Sec. 3 is explicit and I didn't find a sign or factor error in the parts I checked. The self-citation to Refs. 15 and 16 is appropriate: the Gamma operator comes from there, and the main derivation is independent of its specific form. Credit also for not fitting the final energies; sigma enters through Gamma, not as a fit to target eigenvalues.\n\nSoft spots, in order of significance. First, the phrase 'recover exact quantum energies' is too strong. What is actually shown is an order-by-order block diagonalization. The denominators in Eq. 107 are electronic energy gaps, so the series is asymptotic at best, and no convergence or error bound is given. The authors do acknowledge this in Sec. 6 — they explicitly restrict to well-separated ground states and exclude near-degenerate and multiconfigurational cases. So the concern is less that they hid the limitation and more that the abstract and conclusion overstate the formal status. That should be fixed in revision.\n\nSecond, the value of sigma in Eq. 123 is never stated. That's a reproducibility gap for the numerics. Easy to fix, but without it the figures can't be independently reproduced. Third, the second-order correction is truncated to a single term, and the authors say so plainly. Fine for a benchmark, but it means the E(2) results are not a test of the full second-order theory. Fourth, only one model, one dimension. Acceptable for a first demonstration, but it limits how much weight the numerical evidence can carry.\n\nBottom line: the central construction is credible and the paper is a real step toward putting phase-space electronic structure on a systematic footing. It deserves a serious referee — the convergence question and the missing parameter need addressing, but the derivation is worth engaging with.\n\nBest.","headline":"New first-order perturbative correction makes phase-space electronic structure systematically improvable; the 'exact' wording overstates a formal asymptotic scheme, but the core construction is sound and worth refereeing.","tokens_in":20747,"tokens_out":3098,"would_cite":true,"duration_ms":33796,"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":"Phase-space electronic states can be corrected to give exact molecular energies.","keywords":["phase space electronic structure","Born-Oppenheimer","Wigner-Weyl transform","vibrational energies","perturbation theory","nonadiabatic coupling","quantum chemistry","Weyl transform"],"falsifier":"Compute the first-order corrected energy with Eq. 110 on a model where the ground and excited electronic states become nearly degenerate, and compare with exact diagonalization; if the correction diverges or fails to approach the exact eigenvalue as the nuclear basis is enlarged, the formal claim of recovering exact energies breaks down.","tokens_in":1680,"feed_emoji":"⚛️","tokens_out":4079,"duration_ms":106416,"temperature":0.7,"pith_summary":"This paper makes the case that phase-space electronic structure theory, where electronic states are parameterized by both nuclear position and nuclear momentum, is not an ad hoc procedure. The authors show that, starting from the eigenvectors of any phase-space electronic Hamiltonian, one can construct a nearby unitary transformation order by order in ℏ so that the transformed Hamiltonian has the same eigenvalues as the full coupled nuclear-electronic Hamiltonian. The first-order correction is a single matrix element, ℏB00, which can be evaluated without computing the first-order eigenvector corrections. In a one-dimensional hydrogen-bond model, adding this correction improves vibrational energies by roughly an order of magnitude over zeroth-order phase-space results and by much more over standard Born-Oppenheimer results. The formal result is perturbative and is expected to be accurate when the ground electronic state is well separated from excited states.","feed_headline":"Phase-space states yield exact vibrational energies","feed_subtitle":"A first-order correction makes phase-space electronic Hamiltonians formally exact.","key_machinery":"The machinery is the Wigner-Weyl transform and its star product, which allows a quantum operator to be represented by a function of phase-space variables (R,P) while preserving the full quantum structure through an ℏ expansion. On top of this, the paper uses a perturbative unitary transformation: it starts with the phase-space diagonalizing matrix L_W, adds an ℏ-dependent correction L_W^(1), and enforces both unitarity and block-diagonality order by order. The key simplification is that the first-order correction to the energy can be evaluated directly from the matrix element B00, without explicitly constructing L_W^(1), because the diagonal part of the commutator term vanishes. This machinery converts a potentially ad hoc phase-space electronic Hamiltonian into one with a well-defined perturbative route to exact eigenvalues.","core_discovery":"The central claim is that phase-space electronic eigenstates can serve as a rigorous zeroth-order starting point for exact quantum nuclear-electronic calculations. Concretely, the paper constructs a formal operator J_W = L_W + ℏL_W^(1) + ... such that $W^{{-1}}$(J_W^† * H_W * J_W) is diagonal in the electronic subspace and reproduces the exact eigenvalues of the full Hamiltonian. The paper shows that this is possible by enforcing unitarity of J_W order by order in ℏ and then choosing the remaining freedom to make the transformed Hamiltonian block-diagonal. The first-order diagonal matrix element is (H_FPS^(1)W)_00 = ℏB00, given by Eq. 110, which depends only on the zeroth-order phase-space eigenvectors and their derivatives. The paper demonstrates numerically that this correction substantially improves vibrational energies, and it argues that the construction removes the previous objection that phase-space approaches lack a systematic way to approach the exact answer.","pith_inferences":["A similar unitary-reconstruction scheme could be applied to extract exact electronic momentum or current density from phase-space eigenstates, not just vibrational energies.","The first-order correction may be combined with machine-learned phase-space surfaces to avoid explicit derivatives with respect to R and P, making the correction practical for larger molecules.","The requirement of well-separated electronic states points toward the need for a quasi-degenerate or resummed version of this perturbation theory for conical intersections and strongly multiconfigurational systems.","The harmonic-limit finding suggests that the corrected phase-space method is most valuable precisely when anharmonicity is strong; one could test this on a library of small molecules with known anharmonic spectra."],"forward_implications":["Phase-space electronic structure methods can now be systematically improved beyond the zeroth-order Weyl-transformed surface.","Vibrational energies computed from phase-space surfaces can be corrected toward exact coupled nuclear-electronic results by adding the first-order term ℏB00.","The formalism works for any choice of the phase-space operator Γ, so it justifies a family of phase-space Hamiltonians rather than a single preferred one.","The improvement is especially significant outside the harmonic limit, where anharmonic couplings play a larger role.","The perturbative series provides a practical route to non-Born-Oppenheimer corrections without treating the full nonadiabatic problem explicitly."],"supporting_citations":[{"why":"Introduces the phase-space electronic Hamiltonian used as the zeroth-order reference.","marker":"[15]"},{"why":"Provides previous empirical evidence that diagonalizing a Weyl-transformed phase-space surface improves vibrational energies.","marker":"[16]"},{"why":"Supplies the Littlejohn-Flynn perturbative expansion that the paper adapts to the phase-space setting.","marker":"[31]"},{"why":"Demonstrates Moyal perturbation theory for diagonalizing the Born-Oppenheimer Hamiltonian, the direct methodological precursor.","marker":"[32]"},{"why":"Provides the three-particle hydrogen-bond model used for all numerical tests.","marker":"[43]"},{"why":"Recent study of the same model that supplies the baseline comparison for exact and Born-Oppenheimer results.","marker":"[11]"}],"fun_headline_variants":["First-order fix makes phase-space electronic states exact","Exact vibrational energies from phase-space electronic states","Phase-space eigenstates made exact by first-order correction","Formal route to exact energies in phase-space electronic structure","Perturbation theory recovers exact energies in phase-space framework"],"cache_read_input_tokens":22912,"weakest_assumption_plain":"The perturbative expansion in ℏ is valid, meaning the off-diagonal couplings divided by electronic energy gaps are small enough for the series to converge or behave asymptotically.","fun_headline_variants_meta":{"raw":{"variants":["First-order fix makes phase-space electronic states exact","Exact vibrational energies from phase-space electronic states","Phase-space eigenstates made exact by first-order correction","Formal route to exact energies in phase-space electronic structure","Perturbation theory recovers exact energies in phase-space framework"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1522,"prompt_tokens":910,"completion_tokens":612,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":544}},"tokens_in":526,"tokens_out":612,"duration_ms":6269,"temperature":1.0,"reasoning_tokens":544,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:17:40.695406+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the first-order corrected energy with Eq. 110 on a model where the ground and excited electronic states become nearly degenerate, and compare with exact diagonalization; if the correction diverges or fails to approach the exact eigenvalue as the nuclear basis is enlarged, the formal claim of recovering exact energies breaks down.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the phase-space electronic Hamiltonian used as the zeroth-order reference."},{"cited_title":"G.; Subotnik, J","cited_arxiv_id":null,"evidence_quote":"Provides previous empirical evidence that diagonalizing a Weyl-transformed phase-space surface improves vibrational energies."},{"cited_title":"G.; Flynn, W","cited_arxiv_id":null,"evidence_quote":"Supplies the Littlejohn-Flynn perturbative expansion that the paper adapts to the phase-space setting."},{"cited_title":"Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom","cited_arxiv_id":null,"evidence_quote":"Demonstrates Moyal perturbation theory for diagonalizing the Born-Oppenheimer Hamiltonian, the direct methodological precursor."},{"cited_title":"C.; Gaigeot, M.-P.; Borgis, D","cited_arxiv_id":null,"evidence_quote":"Provides the three-particle hydrogen-bond model used for all numerical tests."},{"cited_title":"On the mass of atoms in molecules: Beyond the Born-Oppenheimer approximation","cited_arxiv_id":null,"evidence_quote":"Recent study of the same model that supplies the baseline comparison for exact and Born-Oppenheimer results."}],"review_version":1}