Recognition: no theorem link
Assessing the validity of the Born-Oppenheimer approximation in potential models for doubly heavy hadrons
Pith reviewed 2026-05-16 05:24 UTC · model grok-4.3
The pith
The Born-Oppenheimer approximation deviates from benchmark solutions in potential models for doubly heavy hadrons as the heavy quark mass increases.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In potential models for doubly heavy hadrons, the Born-Oppenheimer approximation yields results close to those from the Gaussian expansion method when the heavy-quark mass is relatively small. As the heavy-quark mass increases, calculations employing Slater-type functions yield larger values than those from the Gaussian expansion method, whereas those using Gaussian-type functions lead to smaller ones. The underestimation observed in Born-Oppenheimer approximation calculations with Gaussian-type functions primarily stems from the neglect of non-adiabatic corrections, and Slater-type functions generally lead to an enhanced binding energy.
What carries the argument
Born-Oppenheimer approximation applied to potential models for doubly heavy hadrons, benchmarked against Gaussian expansion method results and using Slater-type versus Gaussian-type trial wave functions.
If this is right
- Born-Oppenheimer results remain close to the benchmark at small heavy quark masses.
- Slater-type trial functions systematically increase the predicted binding energy relative to the benchmark.
- Gaussian-type trial functions systematically decrease the predicted binding energy, mainly because non-adiabatic corrections are omitted.
- The size of the discrepancy grows with rising heavy quark mass.
- The comparison clarifies the practical limits of the Born-Oppenheimer treatment for doubly heavy hadron structure.
Where Pith is reading between the lines
- For bottom-quark systems the size of non-adiabatic corrections may exceed the accuracy needed for reliable mass predictions.
- The same mass-dependent trend could appear in other quarkonium or tetraquark models that rely on the Born-Oppenheimer separation.
- Testing the pattern across several different potential shapes would show whether the observed deviations are model-independent.
Load-bearing premise
The Gaussian expansion method itself supplies an accurate, converged reference solution to the potential-model equations.
What would settle it
Exact numerical solution of the same potential-model Schrödinger equation for a sequence of increasing heavy-quark masses, followed by direct comparison with the Born-Oppenheimer values obtained using Gaussian-type functions.
Figures
read the original abstract
The Born-Oppenheimer approximation is widely used to investigate the properties of hydrogen-like systems and doubly heavy hadrons. However, the extent to which this approximation captures the features of such systems within potential models remains an open question. In this work, we adopt the results obtained with the Gaussian expansion method as a benchmark to assess the validity of the Born-Oppenheimer approximation within potential models for hadronic systems. We also investigate the dependence of the Born-Oppenheimer approximation results on the choice of trial wave functions. A comprehensive study of the Born-Oppenheimer approximation is carried out by performing calculations using Slater-type functions and Gaussian-type functions as trial wave functions, and by comparing the resulting predictions with those obtained from the Gaussian expansion method. We find that the calculations performed within the Born-Oppenheimer approximation are close to those obtained with the Gaussian expansion method when the heavy-quark mass is relatively small. However, as the heavy-quark mass increases, calculations employing Slater-type functions yield larger values than those from the Gaussian expansion method, whereas those using Gaussian-type functions lead to smaller ones. The use of Slater-type functions generally leads to an enhanced binding energy. The underestimation observed in Born-Oppenheimer approximation calculations with Gaussian-type functions primarily stems from the neglect of non-adiabatic corrections. This comparative study provides deeper insight into the structure of doubly heavy hadrons and helps clarify the applicability and limitations of the Born-Oppenheimer treatment within potential models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript assesses the validity of the Born-Oppenheimer (BO) approximation for doubly heavy hadrons in potential models by benchmarking BO calculations (using Slater-type and Gaussian-type trial wave functions) against the Gaussian expansion method (GEM). It reports good agreement at small heavy-quark masses but increasing deviations at larger masses, with Slater-type functions producing larger binding energies than GEM and Gaussian-type functions producing smaller ones, attributing the underestimation to neglected non-adiabatic corrections.
Significance. If the GEM benchmark is shown to be converged, the work provides useful quantitative insight into the mass-dependent applicability of the BO approximation in hadronic potential models and the sensitivity to trial-function choice. It supplies concrete evidence that non-adiabatic effects become relevant for heavier systems and could inform improved modeling of doubly heavy baryons and tetraquarks.
major comments (2)
- [Methods / GEM implementation] The central claim that deviations from GEM reflect BO limitations (rather than numerical artifacts) rests on GEM being a converged reference. No systematic convergence tests with respect to basis size, number of Gaussians, or variational parameters are reported, nor are stability checks shown as heavy-quark mass increases. This is load-bearing because the reported discrepancies grow with mass, exactly the regime where basis incompleteness could produce the observed pattern (Slater over, Gaussian under).
- [Results] Quantitative uncertainties (error bars, basis truncation estimates) are absent from all reported energies and differences. Without them, the statistical significance of the mass-dependent deviations cannot be assessed, weakening the attribution to non-adiabatic effects.
minor comments (3)
- The specific potential (Cornell, linear-plus-Coulomb, etc.) and its parameters should be stated explicitly in the text or a table, together with the precise range of heavy-quark masses examined.
- [BO calculations] Reproducibility would be aided by listing the explicit functional forms and variational parameters of the Slater-type and Gaussian-type trial functions used in the BO calculations.
- Figure captions and table headings should clarify whether the plotted/ tabulated quantities are binding energies, total energies, or mass differences, and should indicate the sign convention for binding.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments, which help clarify the presentation of our results on the validity of the Born-Oppenheimer approximation. We address each major comment below and indicate the revisions planned for the updated version.
read point-by-point responses
-
Referee: [Methods / GEM implementation] The central claim that deviations from GEM reflect BO limitations (rather than numerical artifacts) rests on GEM being a converged reference. No systematic convergence tests with respect to basis size, number of Gaussians, or variational parameters are reported, nor are stability checks shown as heavy-quark mass increases. This is load-bearing because the reported discrepancies grow with mass, exactly the regime where basis incompleteness could produce the observed pattern (Slater over, Gaussian under).
Authors: We agree that explicit demonstration of GEM convergence is essential to attribute the observed deviations to limitations of the Born-Oppenheimer approximation rather than numerical artifacts. In our calculations, the GEM employs up to 25 Gaussian basis functions per relative coordinate in the three-body system, selected to achieve variational convergence based on prior benchmarks for similar potentials. However, we acknowledge that systematic tests varying the basis size and stability with increasing heavy-quark mass were not reported. In the revised manuscript, we will add convergence tables and a brief discussion showing binding energies for representative masses (e.g., 1.5 GeV and 5 GeV) as the number of Gaussians increases from 10 to 30, confirming stabilization to within ~1-2 MeV. These checks will rule out basis incompleteness as the source of the mass-dependent pattern, while preserving the interpretation that non-adiabatic effects grow with mass. revision: yes
-
Referee: [Results] Quantitative uncertainties (error bars, basis truncation estimates) are absent from all reported energies and differences. Without them, the statistical significance of the mass-dependent deviations cannot be assessed, weakening the attribution to non-adiabatic effects.
Authors: We concur that the absence of quantitative uncertainties limits the ability to assess the significance of the deviations. The reported GEM energies are obtained from well-converged variational calculations, but truncation errors were not explicitly estimated or displayed. In the revised manuscript, we will include basis-truncation uncertainty estimates (derived from the energy difference between calculations with N and N+5 Gaussians) as error bars or parenthetical values in the tables and figures. This will allow direct evaluation of whether the mass-dependent discrepancies between BO (Slater/Gaussian) and GEM exceed the numerical precision, thereby strengthening the link to neglected non-adiabatic corrections. revision: yes
Circularity Check
Independent GEM benchmark avoids circularity in BO validity assessment
full rationale
The paper performs a direct numerical comparison of Born-Oppenheimer results (using Slater-type and Gaussian-type trial functions) against Gaussian expansion method energies treated as an external benchmark. No derivation step reduces to its own inputs by construction, no parameters are fitted from the target data and then relabeled as predictions, and no load-bearing uniqueness theorem or ansatz is imported via self-citation. The central claims rest on this explicit method-to-method comparison rather than self-referential logic. Minor self-citations for context or potential models may exist but are not load-bearing for the validity assessment.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Potential models accurately describe the interactions in doubly heavy hadrons.
- domain assumption The Gaussian expansion method provides a reliable benchmark for the exact numerical solution of the potential model.
Reference graph
Works this paper leans on
-
[1]
First Observation of the Doubly Charmed Baryon Xi_cc^+
M. Mattsonet al.(SELEX), First Observation of the Dou- bly Charmed BaryonΞ + cc, Phys. Rev. Lett.89, 112001 (2002), arXiv:hep-ex/0208014
work page internal anchor Pith review Pith/arXiv arXiv 2002
-
[2]
S. P. Ratti, New Results on c-Baryons and a Search for cc- Baryons in FOCUS, Nucl. Phys. B Proc. Suppl.115, 33 (2003)
work page 2003
-
[3]
Search for Doubly Charmed Baryons Xi_cc^+ and Xi_cc^++ in BABAR
B. Aubertet al.(BaBar), Search for doubly charmed baryons Ξ+ cc andΞ ++ cc in BABAR, Phys. Rev. D74, 011103 (2006), arXiv:hep-ex/0605075
work page internal anchor Pith review Pith/arXiv arXiv 2006
-
[4]
Chistovet al.(Belle), Observation of new states decaying intoΛ + c K−π+ andΛ + c K0 S π−, Phys
R. Chistovet al.(Belle), Observation of new states decaying intoΛ + c K−π+ andΛ + c K0 S π−, Phys. Rev. Lett.97, 162001 (2006), arXiv:hep-ex/0606051
-
[5]
Aaijet al.(LHCb), Precision measurement of theΞ ++ cc mass, JHEP02, 049, arXiv:1911.08594 [hep-ex]
R. Aaijet al.(LHCb), Precision measurement of theΞ ++ cc mass, JHEP02, 049, arXiv:1911.08594 [hep-ex]
-
[6]
Measurement of the lifetime of the doubly charmed baryon $\Xi_{cc}^{++}$
R. Aaijet al.(LHCb), Measurement of the Lifetime of the Dou- bly Charmed BaryonΞ ++ cc , Phys. Rev. Lett.121, 052002 (2018), arXiv:1806.02744 [hep-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[7]
First observation of the doubly charmed baryon decay $\Xi_{cc}^{++}\rightarrow \Xi_{c}^{+}\pi^{+}$
R. Aaijet al.(LHCb), First Observation of the Doubly Charmed Baryon DecayΞ ++ cc →Ξ + c π+, Phys. Rev. Lett.121, 162002 (2018), arXiv:1807.01919 [hep-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[8]
R. Aaijet al.(LHCb), Observation of the doubly charmed baryon decayΞ ++ cc →Ξ ′+ c π+, JHEP05, 038, arXiv:2202.05648 [hep-ex]
-
[9]
R. Aaijet al.(LHCb), Observation of the doubly-charmed- baryon decayΞ ++ cc →Ξ 0 cπ+π+, JHEP10, 136, arXiv:2504.05063 [hep-ex]
-
[10]
R. Aaijet al.(LHCb), Observation of an exotic narrow doubly charmed tetraquark, Nature Phys.18, 751 (2022), arXiv:2109.01038 [hep-ex]
-
[11]
R. Aaijet al.(LHCb), Study of the doubly charmed tetraquark T + cc, Nature Commun.13, 3351 (2022), arXiv:2109.01056 [hep- ex]
-
[12]
H.-X. Chen, W. Chen, X. Liu, and S.-L. Zhu, The hidden-charm pentaquark and tetraquark states, Phys. Rept.639, 1 (2016), arXiv:1601.02092 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[13]
Pentaquark and Tetraquark states
Y .-R. Liu, H.-X. Chen, W. Chen, X. Liu, and S.-L. Zhu, Pen- taquark and Tetraquark states, Prog. Part. Nucl. Phys.107, 237 (2019), arXiv:1903.11976 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2019
- [14]
-
[15]
M.-Z. Liu, Y .-W. Pan, Z.-W. Liu, T.-W. Wu, J.-X. Lu, and L.- S. Geng, Three ways to decipher the nature of exotic hadrons: Multiplets, three-body hadronic molecules, and correlation functions, Phys. Rept.1108, 1 (2025), arXiv:2404.06399 [hep- ph]
-
[16]
X. Wang, X. Liu, and Y . Gao, Colloquium: Hadron production in open-charm meson pairs at e+e- colliders, Rev. Mod. Phys. 98, 021001 (2026), arXiv:2502.15117 [hep-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[17]
M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln, Annalen Phys.389, 457 (1927)
work page 1927
- [18]
-
[19]
Born-Oppenheimer Approximation for the XYZ Mesons
E. Braaten, C. Langmack, and D. H. Smith, Born-Oppenheimer Approximation for the XYZ Mesons, Phys. Rev. D90, 014044 (2014), arXiv:1402.0438 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[20]
P. Bicudo, K. Cichy, A. Peters, B. Wagenbach, and M. Wagner, Evidence for the existence ofud ¯b¯band the non-existence of ss ¯b¯bandcc ¯b¯btetraquarks from lattice QCD, Phys. Rev. D92, 014507 (2015), arXiv:1505.00613 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[21]
P. Bicudo, J. Scheunert, and M. Wagner, Including heavy spin effects in the prediction of a¯b¯budtetraquark with lattice QCD potentials, Phys. Rev. D95, 034502 (2017), arXiv:1612.02758 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[22]
R. Bruschini, Heavy-quark spin symmetry breaking in the Born-Oppenheimer approximation, JHEP08, 219, arXiv:2303.17533 [hep-ph]
-
[23]
M. Berwein, N. Brambilla, A. Mohapatra, and A. Vairo, Hy- brids, tetraquarks, pentaquarks, doubly heavy baryons, and quarkonia in Born-Oppenheimer effective theory, Phys. Rev. D 110, 094040 (2024), arXiv:2408.04719 [hep-ph]
-
[24]
E. Braaten and R. Bruschini, Model-independent predictions for decays of hidden-heavy hadrons into pairs of heavy hadrons, Phys. Rev. D109, 094051 (2024), arXiv:2403.12868 [hep-ph]
-
[25]
E. Braaten and R. Bruschini, Exotic hidden-heavy hadrons and where to find them, Phys. Lett. B863, 139386 (2025), arXiv:2409.08002 [hep-ph]
-
[26]
L. Maiani, A. D. Polosa, and V . Riquer, Hydrogen bond of QCD, Phys. Rev. D100, 014002 (2019), arXiv:1903.10253 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[27]
B. Grinstein, L. Maiani, and A. D. Polosa, Radiative decays of X(3872) discriminate between the molecular and compact inter- pretations, Phys. Rev. D109, 074009 (2024), arXiv:2401.11623 [hep-ph]
- [28]
-
[29]
D. Germani, B. Grinstein, and A. D. Polosa, Tetraquarks in the Born-Oppenheimer approximation, JHEP04, 004, arXiv:2501.13249 [hep-ph]
- [30]
-
[31]
L. Pauling, The application of the quantum mechanics to the structure of the hydrogen molecule and hydrogen molecule-ion and to related problems., Chemical Reviews5, 173 (1928)
work page 1928
-
[32]
A. Szabo and N. S. Ostlund,Modern quantum chemistry: intro- duction to advanced electronic structure theory(Dover Publi- cations, 1996)
work page 1996
- [33]
-
[34]
E. Hiyama, Gaussian expansion method for few-body systems and its applications to atomic and nuclear physics, PTEP2012, 01A204 (2012)
work page 2012
- [35]
- [36]
-
[37]
H. Zhou, S.-Q. Luo, and X. Liu, Triply heavy baryon spectroscopy revisited, Phys. Rev. D112, 074007 (2025), arXiv:2507.10243 [hep-ph]
- [38]
-
[39]
Z.-L. Zhang and S.-Q. Luo, Spectroscopic properties of 1F- wave singly bottom baryons, Phys. Rev. D112, 074020 (2025), arXiv:2504.17507 [hep-ph]
- [40]
- [41]
-
[42]
L. Meng, Y .-K. Chen, Y . Ma, and S.-L. Zhu, Tetraquark bound states in constituent quark models: Benchmark test calcula- tions, Phys. Rev. D108, 114016 (2023), arXiv:2310.13354 [hep-ph]
- [43]
- [44]
- [45]
- [46]
-
[47]
Y .-X. Peng, S.-Q. Luo, and X. Liu, Refining radiative decay studies in singly heavy baryons, Phys. Rev. D110, 074034 (2024), arXiv:2405.12812 [hep-ph]
- [48]
- [49]
-
[50]
Weinberg,Lectures on Quantum Mechanics(Cambridge University Press, 2015)
S. Weinberg,Lectures on Quantum Mechanics(Cambridge University Press, 2015)
work page 2015
-
[51]
V . F. Bratsev, The ground state energy of a molecule in adia- batic approximation, Dokl. Akad. Nauk SSSR160, 570 (1965), [Proc. Acad. Sci. USSR]
work page 1965
-
[52]
T. K. Das, H. T. Coelho, and V . P. Brito, Comparison of Born- Oppenheimer and hyperspherical adiabatic approximations in the trinucleon problem, Phys. Rev. C48, 2201 (1993)
work page 1993
- [53]
- [54]
- [55]
- [56]
- [57]
-
[58]
J. L. Richardson, The Heavy Quark Potential and the Upsilon, J/ψSystems, Phys. Lett. B82, 272 (1979)
work page 1979
-
[59]
Hadronic matrix elements and radiative $B\to K^{*}\gamma$ decay
J. Tang, J.-H. Liu, and K.-T. Chao, Hadronic matrix elements and radiativeB→K ∗γdecay, Phys. Rev. D51, 3501 (1995), arXiv:hep-ph/9502411
work page internal anchor Pith review Pith/arXiv arXiv 1995
-
[60]
G.-L. Wang, T.-F. Feng, and Y .-Q. Wang, Mass spectra and wave functions of toponia, Phys. Rev. D111, 096016 (2025), arXiv:2411.17955 [hep-ph]
-
[61]
Charmonium potential from full lattice QCD
T. Kawanai and S. Sasaki, Charmonium potential from full lat- tice QCD, Phys. Rev. D85, 091503 (2012), arXiv:1110.0888 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[62]
Navaset al.(Particle Data Group), Review of particle physics, Phys
S. Navaset al.(Particle Data Group), Review of particle physics, Phys. Rev. D110, 030001 (2024)
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.