REVIEW 2 major objections 5 minor 32 references
Born-Oppenheimer, Born-Huang, and exact factorization: quantum geometry and error in analytically transparent benchmark models
T0 review · 2 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper proves that in an exactly solvable two-oscillator model the true ground-state energy always lies between the Born–Oppenheimer and Born–Huang values, and traces both gaps to quantum geometry.
desk verdict Worth refereeing: clean closed-form benchmarks connecting DBOC, quantum metric, and nonadiabatic error, but Proposition 1's typeset inequality must be fixed before publication. 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 load-bearing object is the derivative coupling $d_{jk}(R)=\langle\phi_j|\partial_R\phi_k\rangle$, the rate at which nuclear motion rotates one clamped-electronic state toward another. The argument is carried by two spectral sums over these couplings: the quantum metric $g_k=\sum_{j\neq k}|d_{jk}|^2$, whose mass-weighted value $g_k/2M$ is exactly the diagonal Born–Huang correction; and the gap-weighted moment $G_k=\sum_{j\neq k}|d_{jk}|^2/(E_j-E_k)$, which enters the leading off-diagonal error through the integral $\int G_k(R)|\partial_R\chi_v(R)|^2\,dR$. In the exactly solvable oscillator model both sums become constants, turning the error budget into a closed identity whose two terms have opposite signs and different dependence on the vibrational quantum number. In the vibronic model the same coupling is a half-angle rotation $\vartheta'(R)/2$, giving a Lorentzian metric and a cumulative Fubini–Study length that saturates at $\pi/2$.
What would settle it
Evaluate the closed-form ground-state energies of the bilinearly coupled two-oscillator model at $M=1$ and $\beta=0.9$: if the exact energy is not strictly between the Born–Oppenheimer and Born–Huang values, Proposition 1 is false. For the error budget, compute the exact two-channel residual in the vibronic model at small $\lambda$ and compare it with the local-gap expression, the integral of $G_k$ weighted by $|\partial_R\chi_v|^2$; the comparison should degrade as vibrational spacings approach the electronic gap, revealing whether the local-resolvent assumption is the load-bearing one.
Extended reading notes
Core claim
On the paper's own terms: the three constructions should be viewed as three different reductions of the same molecular Schrödinger equation, not as rival definitions. The central quantitative claims are: (i) in the bilinearly coupled oscillator model, $E_{\mathrm{BO}} < E_{\mathrm{exact}} < E_{\mathrm{BH}}$ for the ground state for every $M>0$ and $0<|\beta|<1$, with the inequalities not extending uniformly to excited states; (ii) the diagonal Born–Huang correction is $W_n = g_n/(2M)$ with $g_n$ the quantum metric, and the leading exact-minus-Born–Oppenheimer error has the closed form $E_{0v}^{\mathrm{exact}} - E_{0v}^{\mathrm{BO}} = \frac{\beta^2}{2M}\left[\frac12 - (v+\frac12)\sqrt{\frac{1-\beta^2}{M}}\right] + O(\eta^8)$; (iii) in the two-level vibronic model the metric concentrates at the avoided crossing while the total Fubini–Study length of the lower adiabatic state is exactly $\pi/2$ for all parameters; and (iv) the exact-factorization potential is smooth for the nodeless ground state, and finite but increasingly ill-conditioned when an excited-state nuclear marginal becomes small. The paper reads these results as separating approximation error, geometric correction, and numerical conditioning.
Load-bearing premise
The general error budget assumes that the energy lost to other electronic channels can be approximated by replacing the actual nuclear vibrational energy levels with the local electronic gap; that replacement is accurate only when vibrational spacings are small compared with electronic gaps.
Editorial extensions
If this is right
- In the bilinearly coupled oscillator model, a single-surface calculation always brackets the exact ground state: $E_{\mathrm{BO}} < E_{\mathrm{exact}} < E_{\mathrm{BH}}$, for every admissible mass and coupling.
- Adding the diagonal Born–Huang correction improves the formal order of the ground-state energy, but it can worsen specific excited states, because the positive diagonal term is constant while the negative off-diagonal term grows linearly in the vibrational quantum number.
- Since the diagonal correction is $W_k = g_k/(2M)$, single-surface error grows where the electronic state changes rapidly with nuclear position; near an avoided crossing the metric sharpens as the gap closes.
- The total Fubini–Study rotation across the avoided crossing is $\pi/2$ for all parameters, so closing the gap compresses a fixed electronic-state change into a smaller interval of nuclear coordinate; the growing correction is a localization effect, not an increase in total state change.
- The exact-factorization surface can become large but remains finite and grid-converged when the nuclear marginal is small; its height is controlled by the second derivative of the marginal density, not by the small marginal alone.
Reading between the lines
- A general lesson this reading draws is that on the lowest electronic surface the two leading corrections always compete — a positive geometric term and a negative gap-weighted term — so the net single-surface error is not a property of the surface alone but depends on the nuclear state through its kinetic energy.
- The paper's own outlook suggests a direct test: replace the harmonic nuclear well in the first model with a double well while keeping the electronic geometry unchanged; any change in the ordering would isolate the role of nuclear dynamics in the error budget.
- Because the total Fubini–Study length is parameter-independent in the two-level model, cumulative state-space length is a more transferable diagnostic of nonadiabatic change than the peak metric: the peak says where change happens, the length says how much total change there is.
- The near-node analysis implies that numerical reports of exact-factorization potentials should quote the marginal amplitude alongside the potential; a large reconstructed value in a region of tiny marginal is a conditioning signal, not evidence that exact factorization has broken down.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper compares three ways of defining an effective nuclear potential from the electron-nuclear Schrödinger equation: the clamped-nucleus Born–Oppenheimer surface, the single-surface Born–Huang surface with the diagonal correction, and the state-dependent exact-factorization potential. The comparison is carried out in two solvable benchmark models. In Fernández's bilinearly coupled oscillator model the exact, Born–Oppenheimer, and Born–Huang spectra are obtained in closed form, a ground-state energy ordering is claimed (Proposition 1), excited-state counterexamples are exhibited, the diagonal Born–Huang correction is identified with the quantum metric, and a closed-form leading error budget is given (Proposition 4). In a linear vibronic-coupling model the quantum metric localizes at an avoided crossing while the total Fubini–Study length is exactly π/2 (Proposition 2), and exact-factorization reconstruction is shown to be finite but increasingly ill-conditioned near a small nuclear marginal (Proposition 3). The results are supported by exact algebra, closed forms, and grid-converged numerical diagonalization.
Significance. If the presentation issues are repaired, this is a useful benchmark contribution. Its strengths are the completely transparent closed-form treatment of the Fernández model, the direct algebraic proof of the ground-state ordering, the explicit geometric identification of the diagonal Born–Huang correction with the mass-weighted quantum metric, the clean separation in Proposition 2 between localization of the metric and total electronic-state change, and the numerical evidence for grid convergence in the vibronic model. The paper is self-contained and carefully distinguishes the three constructions, which is genuinely helpful given how often these objects are conflated. The geometric identification itself is not new—it is present in the literature and the paper says so—and the models are deliberately simple, so the contribution is primarily as an unambiguous, checkable test bed rather than as a new general theorem. With the central notation and one derivation issue fixed, it would be a solid and citable reference for nonadiabatic error analysis and exact factorization.
major comments (2)
- [Sec. 3, Eq. (22); Sec. 7; Sec. 8.1, Eq. (44)] The headline Proposition 1 is typeset as E00 < E00 < E_A00, with the exact and Born–Oppenheimer ground-state energies rendered by the same symbol. Read literally, this asserts exact < BO < BH, which contradicts the proof's own Eq. (24): (2E_exact)^2 - (2E_BO)^2 = β^2/M > 0 and (2E_BH)^2 - (2E_exact)^2 = β^2 ω/M + β^4/(4M^2) > 0, with all energies positive, give BO < exact < BH. The numerical table in Sec. 3 (0.657321, 0.657511, 0.657571 at M=10, β=0.1) confirms the intended ascending order as BO, exact, BH. The same symbol collision affects the vibronic ordering statement in Sec. 7 and Proposition 4's Eq. (44), where E_{0v} appears on both sides of the equation. The algebra is clearly correct, but a reader cannot verify the paper's central claim without guessing which E is which; the manuscript needs distinct symbols for the exact, Born–Oppenheimer, and Born–Huang energies throughout, including in tables and equation displays.
- [Sec. 8.1, Eq. (42); Appendix B, Eq. (45)] Equation (42) is presented as the leading off-diagonal correction obtained by keeping the first-derivative coupling and closing over the upper vibrational manifold in the local electronic-resolvent approximation. As the derivation in Appendix B shows, the relevant off-diagonal operator is Λ_{kj} = -(1/M)d_{kj}∂_R - (1/(2M))d'_{kj} in the two-state case. After integration by parts and closure over the vibrational basis, the summed squared matrix element is proportional to ∫(d χ'_v + (1/2)d' χ_v)^2, not to ∫|d|^2 |χ'_v|^2. The Fernández model has d' = 0, so Proposition 4 and the closed-form budget (44) are unaffected. But as a general geometric error budget, Eq. (42) drops a d'-dependent term without stating or justifying the approximation. The text should either include the additional term, or explicitly restrict Eq. (42) to locally constant derivative couplings and show that the omitted term is higher order in the intended expansion.
minor comments (5)
- [Sec. 3, Table] The table header uses identical E symbols for the Born–Oppenheimer and exact columns; even with the surrounding text, the column order is ambiguous. Please label the columns explicitly, e.g., E_{nv}^{BO}, E_{nv}^{exact}, E_{nv}^{BH}, and use the same labels in the 'comparison' column.
- [Sec. 4, Eq. (31) vs Eq. (34)] Equation (31) writes d_{+-} = (1/2)ϑ'(R), while Eq. (34) gives d_{+-} = -κλ/(2(κ^2R^2+λ^2)) for the vibronic model. The sign convention for ϑ should be stated explicitly so that the two formulas are compatible, or the sign in the explicit expression should be reconciled with the definition of ϑ.
- [Sec. 5] The statement that a gauge with vanishing vector potential is chosen for real stationary states should include a brief justification: for a real wavefunction, the conditional electronic factor Φ_R can be chosen real, making the vector potential zero. As written, the reader may wonder whether this gauge choice is always available.
- [Abstract and Sec. 8.1] The phrase 'sixth-order result recast as a leading geometric error budget' would be more precise as 'the leading terms of the sixth-order result,' because Eq. (44) carries an O(η^8) remainder, and the quoted ratio of about 0.77 at v=0 shows that the remainder is not numerically negligible at M=10.
- [Sec. 7, Fig. 3(c)] The statement that all results are stable under grid refinement is supported by the figure, but the text would benefit from a brief description of which quantity is plotted in Fig. 3(c) and the numerical tolerance used to declare convergence.
Circularity Check
No significant circularity: all load-bearing results are derived by direct algebra from the stated model Hamiltonians, with no fitted parameters and no self-citation chain.
full rationale
The paper's central claims—Proposition 1 (ground-state ordering), Proposition 2 (Fubini–Study length π/2), Proposition 3 (near-node conditioning), and Proposition 4 (closed-form error budget)—are all obtained by explicit calculation from the two model Hamiltonians. The exact spectrum, clamped-nucleus surfaces, diagonal Born–Huang correction, derivative couplings, and exact-factorization potential are each computed in closed form from the defining operators; none is assumed from or fitted to the target result. The identification of the diagonal Born–Huang correction with the mass-weighted quantum metric is derived in Section 4 from the overlap expansion and completeness (Eqs. 27-28), then independently evaluated in the models, rather than being imported as an input. The recasting of Fernández's sixth-order result as a geometric budget in Proposition 4 is an ex-post reinterpretation of an independently re-derived perturbation expansion, and the paper explicitly acknowledges the provenance of the result. There are no self-citations by the author, and the cited external results (Hunter, Cederbaum, Abedi-Maitra-Gross, Provost-Vallee, etc.) are standard and not used to force any model-specific conclusion. The main presentation blemish is the symbol collision in Proposition 1 and surrounding text, where the exact and Born–Oppenheimer energies both appear as E00, making the printed inequality ambiguous; however, the proof's explicit algebra (Eqs. 23-24) and the numerical table resolve the intended ordering, so this is a typesetting/verifiability issue rather than a circular argument. No fitted parameter is renamed as a prediction, no uniqueness theorem is invoked to forbid alternatives, and no ansatz is smuggled in via citation. The paper is therefore self-contained and non-circular.
Assumptions & free parameters
assumptions (4)
- domain assumption The finite-dimensional model Hamiltonians (two coupled oscillators, two-level vibronic) have complete adiabatic bases, so resolutions of identity and perturbation theory are exact in these spaces.
- domain assumption For real stationary states, a gauge exists in which the exact-factorization vector potential vanishes, so the scalar potential alone defines the nuclear equation.
- domain assumption The local electronic-resolvent approximation, replacing exact vibrational denominators by local electronic gaps in the second-order off-diagonal correction, yields the leading error budget.
- standard math Standard quantum-mechanical results: Rayleigh-Ritz variational principle, virial theorem, Hellmann-Feynman relation, second-order perturbation theory, and the Fubini-Study metric definition.
Cite this review
Pith. "Pith review of Born-Oppenheimer, Born-Huang, and exact factorization: quantum geometry and error in analytically transparent benchmark models." pith.science (2026). https://pith.science/paper/UUNZ7UOG
@misc{pith2026260808668,
author = {Pith},
title = {Pith review of: Born-Oppenheimer, Born-Huang, and exact factorization: quantum geometry and error in analytically transparent benchmark models},
year = {2026},
howpublished = {\url{https://pith.science/paper/UUNZ7UOG}},
note = {Machine review of arXiv:2608.08668}
}
abstract
The phrase "potential energy surface" refers to several distinct objects. The Born-Oppenheimer construction gives a clamped-nucleus electronic eigenvalue, the single-surface Born-Huang construction adds the diagonal correction, and exact factorization gives a state-dependent exact scalar potential. These constructions answer different questions and should not be regarded as competing definitions of one universal surface. We compare them in two analytically transparent benchmark models. For Fern\'andez's bilinearly coupled oscillators, the exact molecular spectrum, the Born-Oppenheimer and Born-Huang spectra, and the ground-state exact-factorization surface are obtained in closed form. We prove the ground-state ordering of these energies for every admissible mass ratio and coupling and show why it does not extend uniformly to excited states. The diagonal Born-Huang correction is identified with the mass-weighted quantum metric, and Fern\'andez's sixth-order result is recast as a leading geometric error budget involving the metric and a gap-weighted spectral moment of the same derivative couplings. In a linear vibronic-coupling model the metric localizes at an avoided crossing while the total Fubini-Study length remains $\pi/2$, separating the localization of electronic-state change from its total magnitude. Exact factorization is smooth for the nodeless ground state but becomes increasingly ill-conditioned, while remaining finite, when an excited-state nuclear marginal becomes small. These models separate approximation error, geometric correction, and conditioning in a form that can be checked directly.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
M. Born and R. Oppenheimer,Zur Quantentheorie der Molekeln, Ann. Phys. (Leipzig)389, 457 (1927)
work page 1927
-
[2]
M. Born and K. Huang,Dynamical Theory of Crystal Lattices(Oxford University Press, Oxford, 1954), Appendix VIII
work page 1954
-
[3]
Kutzelnigg,The adiabatic approximation I
W. Kutzelnigg,The adiabatic approximation I. The physical background of the Born–Handy ansatz, Mol. Phys.90, 909 (1997)
work page 1997
-
[4]
N. C. Handy, Y. Yamaguchi, and H. F. Schaefer III,The diagonal correction to the Born– Oppenheimer approximation: Its effect on the singlet–triplet splitting of CH 2 and other molecular effects, J. Chem. Phys.84, 4481 (1986)
work page 1986
-
[5]
F. M. Fern´ andez,Corrections to the Born–Oppenheimer approximation by means of perturbation theory, Phys. Rev. A50, 2953 (1994)
work page 1994
-
[6]
G. I. Kerley,On corrections to the Born–Oppenheimer approximation, Kerley Technical Services Report (2011); arXiv:1306.6574 (2013)
work page Pith review arXiv 2011
-
[7]
C. W. Duncan and A. del Campo,Shortcuts to adiabaticity assisted by counterdiabatic Born– Oppenheimer dynamics, New J. Phys.20, 085003 (2018)
work page 2018
-
[8]
C. A. Arango,The use of the Born–Oppenheimer factorization in the phase-space representa- tion of the time-independent Schr¨ odinger equation for bilinearly coupled harmonic oscillators, arXiv:2108.11500 (2021)
work page Pith review arXiv 2021
Show all 32 references
-
[9]
V. F. Brattsev,The ground-state energy of a molecule in the adiabatic approximation, Dokl. Akad. Nauk SSSR160, 570 (1965) [Sov. Phys. Dokl.10, 44 (1965)]
1965
-
[10]
S. T. Epstein,Ground-state energy of a molecule in the adiabatic approximation, J. Chem. Phys.44, 836–837 (1966); erratum,44, 4062 (1966). 30
1966
-
[11]
C. W. Patterson,Corrections to the Born–Oppenheimer approximation for a harmonic oscillator, Phys. Rev. A47, 879 (1993)
1993
-
[12]
J. H. Van Vleck,On the isotope corrections in molecular spectra, J. Chem. Phys.4, 327 (1936)
1936
-
[13]
Fock,Bemerkung zum Virialsatz, Z
V. Fock,Bemerkung zum Virialsatz, Z. Phys.63, 855 (1930)
1930
-
[14]
Hunter,Conditional probability amplitudes in wave mechanics, Int
G. Hunter,Conditional probability amplitudes in wave mechanics, Int. J. Quantum Chem.9, 237 (1975)
1975
-
[15]
L. S. Cederbaum,The exact molecular wavefunction as a product of an electronic and a nuclear wavefunction, J. Chem. Phys.138, 224110 (2013)
2013
-
[16]
Abedi, N
A. Abedi, N. T. Maitra, and E. K. U. Gross,Exact factorization of the time-dependent electron-nuclear wave function, Phys. Rev. Lett.105, 123002 (2010)
2010
-
[17]
Abedi, N
A. Abedi, N. T. Maitra, and E. K. U. Gross,Correlated electron-nuclear dynamics: Exact factorization of the molecular wavefunction, J. Chem. Phys.137, 22A530 (2012)
2012
-
[18]
Abedi, F
A. Abedi, F. Agostini, Y. Suzuki, and E. K. U. Gross,Dynamical steps that bridge piecewise adiabatic shapes in the exact time-dependent potential energy surface, Phys. Rev. Lett.110, 263001 (2013)
2013
-
[19]
Requist, F
R. Requist, F. Tandetzky, and E. K. U. Gross,Molecular geometric phase from the exact electron-nuclear factorization, Phys. Rev. A93, 042108 (2016)
2016
-
[20]
Jecko, B
T. Jecko, B. T. Sutcliffe, and R. G. Woolley,On factorization of molecular wavefunctions, J. Phys. A: Math. Theor.48, 445201 (2015)
2015
-
[21]
Jecko, B
T. Jecko, B. T. Sutcliffe, and R. G. Woolley,Corrigendum: On factorization of molecular wavefunctions, J. Phys. A: Math. Theor.51, 149501 (2018)
2018
-
[22]
B. T. Sutcliffe and R. G. Woolley,On the quantum theory of molecules, J. Chem. Phys.137, 22A544 (2012)
2012
-
[23]
Shin and H
S. Shin and H. Metiu,Nonadiabatic effects on the charge transfer rate constant: A numerical study of a simple model system, J. Chem. Phys.102, 9285 (1995)
1995
-
[24]
K¨ oppel, W
H. K¨ oppel, W. Domcke, and L. S. Cederbaum,Multimode molecular dynamics beyond the Born–Oppenheimer approximation, Adv. Chem. Phys.57, 59 (1984)
1984
-
[25]
A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Dynamics of the dissipative two-state system, Rev. Mod. Phys.59, 1 (1987)
1987
-
[26]
Braak,Integrability of the Rabi model, Phys
D. Braak,Integrability of the Rabi model, Phys. Rev. Lett.107, 100401 (2011)
2011
-
[27]
Hellmann,Einf¨ uhrung in die Quantenchemie(Franz Deuticke, Leipzig, 1937)
H. Hellmann,Einf¨ uhrung in die Quantenchemie(Franz Deuticke, Leipzig, 1937)
1937
-
[28]
R. P. Feynman,Forces in molecules, Phys. Rev.56, 340 (1939)
1939
-
[29]
J. P. Provost and G. Vallee,Riemannian structure on manifolds of quantum states, Commun. Math. Phys.76, 289 (1980)
1980
-
[30]
M. V. Berry,Quantal phase factors accompanying adiabatic changes, Proc. R. Soc. Lond. A 392, 45 (1984). 31
1984
-
[31]
Herzberg and H
G. Herzberg and H. C. Longuet-Higgins,Intersection of potential energy surfaces in polyatomic molecules, Discuss. Faraday Soc.35, 77 (1963)
1963
-
[32]
C. A. Mead and D. G. Truhlar,On the determination of Born–Oppenheimer nuclear motion wave functions including complications due to conical intersections and identical nuclei, J. Chem. Phys.70, 2284 (1979). 32
1979
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.