{"id":"c5f21aa1-702c-4ed1-a214-c0da7105da79","arxiv_id":"2608.08668","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"In two solvable molecular models, the Born-Oppenheimer, Born-Huang, and exact-factorization energy surfaces are compared in closed form, yielding a ground-state ordering, a quantum-geometric error budget, and a conditioning result.","lead":"This paper compares three different ways to define a molecule's energy surface using two simple, exactly solvable model molecules. It shows which energy ordering holds for the ground state, why it breaks for excited states, and how geometry and numerical conditioning shape the corrections.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"As printed, Proposition 1 conflates the exact and Born–Oppenheimer energy symbols, so the headline ordering cannot be checked; the proof and Table instead support E_BO < E_exact < E_BH.","rationale":"I checked the mathematics behind the central ordering and error-budget claims before settling on a concern. The proof of Proposition 1 is internally consistent once the two energy symbols are distinguished: Eq. (24) gives exact > BO, and the variational argument gives exact ≤ BH, so the intended chain is BO < exact < BH. The numerical table supports this chain. The local-resolvent approximation in Appendix B is explicitly identified by the paper, its O(η^8) remainder is stated, and the model-specific ratios are given, so I do not regard it as a hidden flaw. The gauge choice for the exact-factorization scalar potential is legitimate for the real stationary states treated here. The genuinely load-bearing issue is the notation collision in the statement of the headline result: as printed, Proposition 1's inequality appears to assert the opposite of what the proof and table establish. This is a presentation flaw, but it is load-bearing because the ordering is the paper's central theorem and cannot be checked without guessing which E is which. The reader identified this issue in the rationale but did not make it the weakest assumption; I agree that the local-resolvent approximation is also an assumption, but it is clearly disclosed and tested. My recommendation is therefore to keep the conditional verdict, pending an unambiguous restatement of Proposition 1 and the other affected energy comparisons.","tokens_in":24058,"tokens_out":31026,"duration_ms":329266,"concrete_test":"Re-typeset every occurrence of the distinct energy symbols in Proposition 1, Eqs. (23)–(24), the Sec. 3 table, and the Sec. 7 vibronic sentence with unambiguous glyphs (e.g., E for BO, \\mathcal{E} for exact, E^A for Born–Huang). Then re-check the inequalities: the proof of Eq. (24) plus positivity requires BO < exact, while the variational argument requires exact < BH; if the corrected chain BO < exact < BH is consistent with both the table and the vibronic values, the theorem stands. Also verify Proposition 4's exact-minus-BO sign by evaluating Eq. (44) against the exact oscillator energies for M = 10, β = 0.3, v = 0..2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Proposition 1 (Eq. 22) is the paper's headline result, but the symbol rendering makes it impossible to tell which of the two identical-looking E symbols is the exact molecular energy and which is the Born–Oppenheimer energy. Read literally, the inequality E00 < E00 < EA00 asserts exact < BO < BH. That contradicts the proof's own Eq. (24): (2E_exact)^2 − (2E_BO)^2 = β²/M > 0 with all energies positive forces BO < exact, and the Sec. 3 table (0.657321, 0.657511, 0.657571 at M = 10, β = 0.1) shows the ascending order is BO, exact, BH. The same symbol collision affects the Sec. 7 vibronic ordering statement and Proposition 4's error formula. The underlying derivations appear correct, because the proof and the variational upper bound fix the intended inequality, but the central claim as typeset is ambiguous, so a reader cannot verify the paper's main theorem without guessing which symbol is which.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":24189,"tokens_out":17926,"duration_ms":190250,"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":[{"comment":"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.","section":"Sec. 3, Eq. (22); Sec. 7; Sec. 8.1, Eq. (44)"},{"comment":"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.","section":"Sec. 8.1, Eq. (42); Appendix B, Eq. (45)"}],"minor_comments":[{"comment":"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.","section":"Sec. 3, Table"},{"comment":"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 ϑ.","section":"Sec. 4, Eq. (31) vs Eq. (34)"},{"comment":"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.","section":"Sec. 5"},{"comment":"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.","section":"Abstract and Sec. 8.1"},{"comment":"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.","section":"Sec. 7, Fig. 3(c)"}],"recommendation":"major_revision","confidential_remarks":"The central derivations appear to be sound, and the intended ordering BO < exact < BH is fully recoverable from Eq. (24) and the tables. The notation collision in Proposition 1 is conspicuous enough that I suspect a LaTeX macro error in the source; it nevertheless must be fixed because it makes the paper's main theorem unreadable as printed. The derivation issue in Eq. (42) is the other substantive point: it is a genuine gap in the general formulation, even though Proposition 4 is protected by d'=0 in the Fernández model. I do not see a deeper mathematical flaw, and the paper is well within the scope of a molecular-physics or quantum-chemistry journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWhat you need to know: this paper delivers exactly what it promises—two analytically solvable models where the Born–Oppenheimer, Born–Huang, and exact-factorization surfaces can be compared in closed form. The line that carries the paper is Proposition 1: in the Fernández model the ground-state energies are strictly ordered for every M>0 and 0<|β|<1. The proof in Eq. (24) is correct and the numerics in Sec. 3 confirm the ascending order BO, exact, BH. There is a real flaw in the typesetting: Eq. (22) prints E00 < E00 < E_A00, with two identical symbols where one is the exact energy and the other the BO energy. Read literally it says exact < BO < BH, which contradicts the proof. The fix is just renaming one symbol, but it has to happen; the headline claim is uncheckable as printed.\n\nWhat is genuinely new: the all-parameter ordering proof tied to the Brattsev/Epstein inequalities, the excited-state reversals (especially (0,2), where BO beats BH), the Fubini–Study length π/2 independence in the vibronic model, and the closed-form error budget identifying the diagonal correction with the mass-weighted quantum metric and the leading nonadiabatic term with a gap-weighted spectral moment. These are not cosmetic repackagings. The derivations are transparent and the numerical tables match the formulas.\n\nThe soft spots are minor relative to the core. The error budget in Eq. (42) uses a local electronic-resolvent approximation—replace exact vibrational denominators by the local gap and close over the upper manifold. The paper flags this and verifies the remainder explicitly in the Fernández model, so it is not a hidden assumption. The gauge choice, real stationary states with vanishing vector potential, is appropriate for the models and stated. No code ships, but the methods are specified closely enough to reimplement.\n\nWho it is for: anyone doing nonadiabatic dynamics, quantum chemistry method development, or exact-factorization theory who wants a test bed that separates approximation error from geometric correction. It deserves a serious referee; the only blocking issue is the symbol collision in Proposition 1 and the parallel statements in Sec. 7 and Proposition 4. Fix that and I would be happy to see it published.","headline":"Worth refereeing: clean closed-form benchmarks connecting DBOC, quantum metric, and nonadiabatic error, but Proposition 1's typeset inequality must be fixed before publication.","tokens_in":24761,"tokens_out":1456,"would_cite":true,"duration_ms":15335,"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":"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.","keywords":["Born–Oppenheimer approximation","Born–Huang expansion","exact factorization","quantum metric","Fubini–Study geometry","nonadiabatic coupling","potential energy surface","avoided crossing"],"falsifier":"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.","tokens_in":23776,"feed_emoji":"⚛️","tokens_out":15141,"duration_ms":141795,"temperature":0.7,"pith_summary":"This paper tries to establish that the 'potential energy surface' is not one thing: the clamped-nucleus Born–Oppenheimer surface, the Born–Huang surface with its diagonal correction, and the state-dependent exact-factorization potential answer different questions, so asking which is the correct surface is the wrong question. To make that distinction checkable, it solves two transparent benchmark models. In the bilinearly coupled two-oscillator model it proves, for every mass ratio and admissible coupling, that the true ground-state energy lies strictly between the Born–Oppenheimer value (below) and the Born–Huang value (above), and it shows by exact computation why that order does not survive for all excited states. It then identifies the diagonal Born–Huang correction with the mass-weighted quantum metric and writes the leading nonadiabatic error as a gap-weighted spectral moment of the same derivative couplings. A reader should care because nearly every molecular calculation runs on one of these surfaces, and the paper shows precisely which couplings each one keeps or discards.","feed_headline":"Ground-state energy always sits between two molecular surfaces","feed_subtitle":"In a solvable two-oscillator model, the two approximations bracket the exact value; the gap is quantum geometry.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"It supplies the Born–Huang channel expansion and the diagonal derivative correction used throughout.","marker":"[2]"},{"why":"It provides the bilinearly coupled oscillator benchmark and the sixth-order perturbation framework that Proposition 4 recasts geometrically.","marker":"[5]"},{"why":"It supplies the virial theorem used to evaluate the derivative norm that gives the diagonal correction in closed form.","marker":"[13]"},{"why":"It introduces the exact factorization of the molecular wavefunction into a nuclear amplitude and a conditional electronic state.","marker":"[14]"},{"why":"It gives the variational formulation of exact factorization that underlies the state-dependent surface.","marker":"[15]"},{"why":"It establishes the exact nuclear equation and uniqueness framework for the factorization potentials.","marker":"[16]"},{"why":"It analyzes regularity of the quotient in exact factorization on nodal sets, framing the conditioning discussion.","marker":"[20]"},{"why":"It is the corrigendum that sharpens the nodal-set regularity analysis for the exact-factorization quotient.","marker":"[21]"},{"why":"It supplies the Hellmann–Feynman force-matrix form that underlies the derivative-coupling formula.","marker":"[27]"},{"why":"It defines the quantum (Fubini–Study) metric that the paper identifies with the diagonal Born–Huang correction.","marker":"[29]"}],"fun_headline_variants":["Quantum metric sets the gap between two molecular surfaces","Exact ground state sits between Born-Oppenheimer and Born-Huang","Three molecular surfaces, one exact, in the middle","Metric localizes at avoided crossing, total length stays pi/2","Born-Huang correction equals quantum metric in solvable model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Quantum metric sets the gap between two molecular surfaces","Exact ground state sits between Born-Oppenheimer and Born-Huang","Three molecular surfaces, one exact, in the middle","Metric localizes at avoided crossing, total length stays pi/2","Born-Huang correction equals quantum metric in solvable model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001225,"raw_usage":{"total_tokens":5120,"prompt_tokens":1115,"completion_tokens":4005,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":731,"completion_tokens_details":{"reasoning_tokens":3922}},"tokens_in":731,"tokens_out":4005,"duration_ms":30708,"temperature":1.0,"reasoning_tokens":3922,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:29:31.256416+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Born and K","cited_arxiv_id":null,"evidence_quote":"It supplies the Born–Huang channel expansion and the diagonal derivative correction used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the bilinearly coupled oscillator benchmark and the sixth-order perturbation framework that Proposition 4 recasts geometrically."},{"cited_title":"Fock,Bemerkung zum Virialsatz, Z","cited_arxiv_id":null,"evidence_quote":"It supplies the virial theorem used to evaluate the derivative norm that gives the diagonal correction in closed form."},{"cited_title":"Hunter,Conditional probability amplitudes in wave mechanics, Int","cited_arxiv_id":null,"evidence_quote":"It introduces the exact factorization of the molecular wavefunction into a nuclear amplitude and a conditional electronic state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the variational formulation of exact factorization that underlies the state-dependent surface."},{"cited_title":"Abedi, N","cited_arxiv_id":null,"evidence_quote":"It establishes the exact nuclear equation and uniqueness framework for the factorization potentials."},{"cited_title":"Jecko, B","cited_arxiv_id":null,"evidence_quote":"It analyzes regularity of the quotient in exact factorization on nodal sets, framing the conditioning discussion."},{"cited_title":"Jecko, B","cited_arxiv_id":null,"evidence_quote":"It is the corrigendum that sharpens the nodal-set regularity analysis for the exact-factorization quotient."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines the quantum (Fubini–Study) metric that the paper identifies with the diagonal Born–Huang correction."}],"review_version":1}