REVIEW 3 major objections 5 minor 2 cited by
Recovering Exact Vibrational Energies Within a Phase Space Electronic Structure Framework
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Phase-space electronic states can be corrected to give exact molecular energies.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Sec. 3.1–3.1.2, Eq. (91) and Eq. (107)] 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.
- [Sec. 3.3, Eq. (111)] 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.
- [Sec. 4.1, Eqs. (120)–(123), Table 3] 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 σ.
minor comments (5)
- [Figs. 2 and 11] 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.
- [Eq. (103)] There is a typo in Eq. (103): "muts" should be "must".
- [References] Refs. 43 and 45 appear to be the same paper (Marinica, Gaigeot, and Borgis) and should be merged.
- [Sec. 4.1] 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.
- [Notation, Sec. 3] 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.
Circularity Check
No circular derivation: the first-order PS correction is computed from the Hamiltonian itself, and the formal result is independent of the cited Gamma form; self-citations are context, not load-bearing.
full rationale
The paper's central derivation is self-contained and non-circular. Starting from a phase space electronic Hamiltonian H^PS_W (Eq. 85), the authors construct a unitary J_W by enforcing unitarity (Eq. 94) and diagonality (Eqs. 98-104) order by order. The first-order energy correction, Eq. 110, is evaluated directly from the zeroth-order eigenvectors L_W and the operators defining H^PS_W and Gamma; no eigenvalue of the full Hamiltonian is used as input. The formal calculation is explicitly independent of the specific form of Gamma ('all formal work below will not depend on the particular form of Gamma'), so the citations of Refs. 15 and 16 for Gamma are not load-bearing. The numerical model is an external benchmark (Borgis/Gross Hamiltonian), and the PS results are compared against exact diagonalization rather than fitted to it. The main weaknesses are rigor gaps, not circularity: the series in Eq. 91 is only enforced order by order, with denominators in Eq. 107 that vanish for near-degenerate electronic states, and the authors concede in Sec. 6 that Eq. 91 is 'only expected to be accurate in the limit that the ground electronic state is reasonably well separated from the higher electronic states.' Likewise, the missing sigma value in Eq. 123 is a reproducibility omission, not evidence of a fitting loop. These considerations place the paper in the 0-2 'no significant circularity' band; the score of 2 reflects the presence of numerous self-citations in the motivation and model setup, but these are not load-bearing, so the derivation itself is not circular.
Assumptions & free parameters
free parameters (1)
- sigma (Gaussian width in Gamma operator) =
not reported
assumptions (4)
- standard math Wigner-Weyl transform maps operators to phase space symbols with Moyal star product; the star product expansion in hbar is valid.
- domain assumption The electronic ground state is non-degenerate and well separated from excited states so denominators Lambda_kk - Lambda_jj in Eq. 107 do not vanish and perturbation theory is valid.
- domain assumption The particular phase space electronic Hamiltonian H_PS(R,P) (Eqs. 85-89) is a physically reasonable zeroth order; its eigenvectors are smooth in R and P.
- ad hoc to paper The choice of Gamma (Eqs. 120-123) from prior work is adequate; no uniqueness principle selects it.
Cite this review
Pith. "Pith review of Recovering Exact Vibrational Energies Within a Phase Space Electronic Structure Framework." pith.science (2026). https://pith.science/paper/CRL65QPO
@misc{pith2026250608230,
author = {Pith},
title = {Pith review of: Recovering Exact Vibrational Energies Within a Phase Space Electronic Structure Framework},
year = {2026},
howpublished = {\url{https://pith.science/paper/CRL65QPO}},
note = {Machine review of arXiv:2506.08230}
}
read the original abstract
In recent years, there has been a push to go beyond Born-Oppenheimer theory and build electronic states from a phase space perspective, i.e. parameterize electronic states by both nuclear position(R) and nuclear momentum(P). Previous empirical studies have demonstrated that such approaches can yield improved single-surface observables, including vibrational energies, electronic momenta, and vibrational circular dichroism spectra. That being said, unlike the case of BO theory, there is no unique phase space electronic Hamiltonian, nor any theory for using phase space eigenvectors (as opposed to BO eigenvectors) so as to recover exact quantum vibrational eigenvalues. As such, one might consider such phase space approaches ad hoc. To that end, here we show how to formally extract exact quantum energies from a coupled nuclear-electronic Hamiltonian using perturbation theory on top of a phase space electronic framework. Thus, while we cannot isolate an "optimal" phase space electronic Hamiltonian, this work does justify a phase space electronic structure approach by offering a rigorous framework for correcting the zeroth order phase space electronic states.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 2 Pith papers
-
Conical Intersections and Electronic Momentum As Viewed From Phase Space Electronic Structure Theory
Phase-space electronic structure theory predicts a three-dimensional branching plane and stable electronic momentum states at the BeH2 conical intersection, explaining complex Hartree-Fock instabilities.
-
The Phase-Space Way To Electronic Structure Theory and Subsequently Chemical Dynamics
The paper proposes phase-space electronic structure theory, where electronic states depend on nuclear momentum as well as position, as a general successor to the Born-Oppenheimer picture.
Reference graph
Works this paper leans on
-
[1]
Introduction to infrared and Raman spectroscopy; Elsevier, 2012
Colthup, N. Introduction to infrared and Raman spectroscopy; Elsevier, 2012
work page 2012
-
[2]
Brandt, N. C.; Keller, E. L.; Frontiera, R. R. Ultrafast Surface-Enhanced Raman Probing of the Role of Hot Electrons in Plasmon-Driven Chemistry. Journal of Physical Chemistry Letters 2016, 7, 3179--3185
work page 2016
-
[3]
o opolare Bindung nach der Quantenmechanik. Zeitschrift f \
Heitler, W.; London, F. Wechselwirkung neutraler Atome und hom \"o opolare Bindung nach der Quantenmechanik. Zeitschrift f \"u r Physik 1927, 44, 455--472
work page 1927
-
[4]
Dynamical theory of crystal lattices; Oxford university press, 1996
Born, M.; Huang, K. Dynamical theory of crystal lattices; Oxford university press, 1996
work page 1996
-
[5]
Cederbaum, L. S. In Conical Intersections: Electronic Structure, Dynamics and Spectroscopy; Domcke, W., Yarkony, D. R., Koppel, H., Eds.; World Scientific Publishing Co.: New Jersey, 2004; pp 3--40
work page 2004
-
[6]
Bernath, P. F. Spectra of atoms and molecules; Oxford university press, 2020
work page 2020
-
[7]
Lefebvre-Brion, H.; Field, R. W. The spectra and dynamics of diatomic molecules: revised and enlarged edition; Elsevier, 2004
work page 2004
-
[8]
Bunker, P.; Moss, R. The effect of the breakdown of the Born-Oppenheimer approximation on the rotation-vibration Hamiltonian of a triatomic molecule. Journal of Molecular Spectroscopy 1980, 80, 217--228
work page 1980
Show all 45 references
-
[9]
R.; Jensen, P
Bunker, P. R.; Jensen, P. Molecular symmetry and spectroscopy; NRC research press, 2006; Vol. 46853
2006
-
[10]
Schwenke, D. W. Beyond the Potential Energy Surface: Ab initio Corrections to the Born- Oppenheimer Approximation for H2O. The Journal of Physical Chemistry A 2001, 105, 2352--2360
2001
-
[11]
On the mass of atoms in molecules: Beyond the Born-Oppenheimer approximation
Scherrer, A.; Agostini, F.; Sebastiani, D.; Gross, E.; Vuilleumier, R. On the mass of atoms in molecules: Beyond the Born-Oppenheimer approximation. Physical Review X 2017, 7, 031035
2017
-
[12]
Adiabatic and diabatic representations for atom-molecule collisions: Treatment of the collinear arrangement
Baer, M. Adiabatic and diabatic representations for atom-molecule collisions: Treatment of the collinear arrangement. Chemical Physics Letters 1975, 35, 112
1975
-
[13]
A.; Truhlar, D
Mead, C. A.; Truhlar, D. G. Conditions for the definition of a strictly diabatic electronic basis for molecular systems. Journal of Chemical Physics 1982, 77, 6090--6098
1982
-
[14]
L.; Adamowicz, L
Bubin, S.; Pavanello, M.; Tung, W.-C.; Sharkey, K. L.; Adamowicz, L. Born--Oppenheimer and non-Born--Oppenheimer, atomic and molecular calculations with explicitly correlated Gaussians. Chemical reviews 2013, 113, 36--79
2013
-
[15]
Tao, Z.; Qiu, T.; Bian, X.; Duston, T.; Bradbury, N.; Subotnik, J. E. A Basis-free phase space electronic Hamiltonian that recovers beyond Born--Oppenheimer electronic momentum and current density. The Journal of Chemical Physics 2025, 162
2025
-
[16]
G.; Subotnik, J
Bian, X.; Khan, C.; Duston, T.; Rawlinson, J.; Littlejohn, R. G.; Subotnik, J. E. A Phase-Space View of Vibrational Energies without the Born–Oppenheimer Framework. Journal of Chemical Theory and Computation 0, 0, null, PMID: 40072941
-
[17]
G.; Subotnik, J
Tao, Z.; Qiu, T.; Bhati, M.; Bian, X.; Duston, T.; Rawlinson, J.; Littlejohn, R. G.; Subotnik, J. E. Practical phase-space electronic Hamiltonians for ab initio dynamics. The Journal of Chemical Physics 2024, 160
2024
-
[18]
Nafie, L. A. Adiabatic molecular properties beyond the Born--Oppenheimer approximation. Complete adiabatic wave functions and vibrationally induced electronic current density. The Journal of chemical physics 1983, 79, 4950--4957
1983
-
[19]
Electronic currents and Born-Oppenheimer molecular dynamics
Patchkovskii, S. Electronic currents and Born-Oppenheimer molecular dynamics. Journal of Chemical Physics 2012, 137, 084109
2012
-
[20]
On the molecular electronic flux: Role of nonadiabaticity and violation of conservation
Hanasaki, K.; Takatsuka, K. On the molecular electronic flux: Role of nonadiabaticity and violation of conservation. Journal of Chemical Physics 2021, 154
2021
-
[21]
G.; Subotnik, J
Qiu, T.; Bhati, M.; Tao, Z.; Bian, X.; Rawlinson, J.; Littlejohn, R. G.; Subotnik, J. E. A Simple One-Electron Expression for Electron Rotational Factors. 2024; J. Chem. Phys., in press
2024
-
[22]
Representation and conservation of angular momentum in the Born--Oppenheimer theory of polyatomic molecules
Littlejohn, R.; Rawlinson, J.; Subotnik, J. Representation and conservation of angular momentum in the Born--Oppenheimer theory of polyatomic molecules. Journal of Chemical Physics 2023, 158, 104302
2023
-
[23]
G.; Subotnik, J
Bian, X.; Tao, Z.; Wu, Y.; Rawlinson, J.; Littlejohn, R. G.; Subotnik, J. E. Total angular momentum conservation in ab initio Born-Oppenheimer molecular dynamics. Physical Review B 2023, 108, L220304
2023
-
[24]
G.; Pei, Z.; Shao, Y.; Subotnik, J
Duston, T.; Tao, Z.; Bian, X.; Bhati, M.; Rawlinson, J.; Littlejohn, R. G.; Pei, Z.; Shao, Y.; Subotnik, J. E. A phase-space electronic Hamiltonian for vibrational circular dichroism. Journal of Chemical Theory and Computation 2024, 20, 7904--7921
2024
-
[25]
Tao, Z.; Duston, T.; Pei, Z.; Shao, Y.; Rawlinson, J.; Littlejohn, R.; Subotnik, J. E. An electronic phase-space Hamiltonian approach for electronic current density and vibrational circular dichroism. The Journal of Chemical Physics 2024, 161
2024
-
[26]
Stephens, P. J. Theory of vibrational circular dichroism. Journal of Physical Chemistry 1985, 89, 748--752
1985
-
[27]
Emergent Newtonian dynamics and the geometric origin of mass
D’Alessio, L.; Polkovnikov, A. Emergent Newtonian dynamics and the geometric origin of mass. Annals of Physics 2014, 345, 141--165
2014
-
[28]
The Moving Born-Oppenheimer Approximation
Barrera, B.; Arovas, D.; Chandran, A.; Polkovnikov, A. The Moving Born-Oppenheimer Approximation. arXiv:2502.17557
-
[29]
Bloch electrons in a magnetic field
Blount, E. Bloch electrons in a magnetic field. Physical Review 1962, 126, 1636
1962
-
[30]
G.; Flynn, W
Littlejohn, R. G.; Flynn, W. G. Geometric phases and the Bohr-Sommerfeld quantization of multicomponent wave fields. Physical review letters 1991, 66, 2839
1991
-
[31]
G.; Flynn, W
Littlejohn, R. G.; Flynn, W. G. Geometric phases in the asymptotic theory of coupled wave equations. Phys. Rev. A 1991, 44, 5239--5256
1991
-
[32]
Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom
Littlejohn, R.; Rawlinson, J.; Subotnik, J. Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom . Journal of Chemical Physics 2024, 160, 114103
2024
-
[33]
Adiabatic perturbation theory in quantum dynamics; Springer Science & Business Media, 2003
Teufel, S. Adiabatic perturbation theory in quantum dynamics; Springer Science & Business Media, 2003
2003
-
[34]
Effective non-adiabatic Hamiltonians for the quantum nuclear motion over coupled electronic states
M \'a tyus, E.; Teufel, S. Effective non-adiabatic Hamiltonians for the quantum nuclear motion over coupled electronic states. The Journal of Chemical Physics 2019, 151
2019
-
[35]
C.; Fang, J
Martens, C. C.; Fang, J. Y. Semiclassical-limit molecular dynamics on multiple electronic surfaces. Journal of Chemical Physics 1997, 106, 4918--4930
1997
-
[36]
Donoso, A.; Martens, C. C. Simulation of coherent nonadiabatic dynamics using classical trajectories. Journal of Physical Chemistry A 1998, 102, 4291--4300
1998
-
[37]
Mixed quantum-classical dynamics
Kapral, R.; Ciccotti, G. Mixed quantum-classical dynamics. Journal of Chemical Physics 1999, 110, 8919--8929
1999
-
[38]
In Energy Transfer Dynamics in Biomaterial Systems; Burghardt, I., Ed.; Springer-Verlag: Berlin, 2009; p 383
Grunwald, R.; Kelly, A.; Kapral, R. In Energy Transfer Dynamics in Biomaterial Systems; Burghardt, I., Ed.; Springer-Verlag: Berlin, 2009; p 383
2009
-
[39]
Quantum-classical description of environmental effects on electronic dynamics at conical intersections
Kelly, A.; Kapral, R. Quantum-classical description of environmental effects on electronic dynamics at conical intersections. Journal of Chemical Physics 2010, 133, 084502
2010
-
[40]
Mixed quantum-classical surface hopping dynamics
Nielsen, S.; Kapral, R.; Ciccotti, G. Mixed quantum-classical surface hopping dynamics. Journal of Chemical Physics 2000, 112, 6543--6553
2000
-
[41]
Consistent schemes for non-adiabatic dynamics derived from partial linearized density matrix propagation
Huo, P.; Coker, D. Consistent schemes for non-adiabatic dynamics derived from partial linearized density matrix propagation. Journal of Chemical Physics 2012, 137, 22A535
2012
-
[42]
Semiclassical Theory of Vibrational Energy Relaxation in the Condensed Phase
Shi, Q.; Geva, E. Semiclassical Theory of Vibrational Energy Relaxation in the Condensed Phase. Journal of Physical Chemistry A 2003, 107, 9059--9069
2003
-
[43]
C.; Gaigeot, M.-P.; Borgis, D
Marinica, D. C.; Gaigeot, M.-P.; Borgis, D. Generating approximate Wigner distributions using Gaussian phase packets propagation in imaginary time. Chemical physics letters 2006, 423, 390--394
2006
-
[44]
F.; Chung, K
Davis, B. F.; Chung, K. T. Mass-polarization effect and oscillator strengths for S, P, D states of helium. Physical Review A 1982, 25, 1328
1982
-
[45]
o opolare Bindung nach der Quantenmechanik , author= Heitler, Walter and London, Fritz , journal= Zeitschrift f \
Marinica, D. C.; Gaigeot, M.-P.; Borgis, D. Generating approximate Wigner distributions using Gaussian phase packets propagation in imaginary time. Chemical physics letters 2006, 423, 390--394, Publisher: Elsevier mcitethebibliography main.tex0000664000000000000000000024473415...
2021
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.