Assessing the validity of the Born-Oppenheimer approximation in potential models for doubly heavy hadrons

This study assesses the validity of the Born-Oppenheimer approximation in potential models for doubly heavy hadrons by benchmarking it against Gaussian expansion method results, revealing that while the approximation is accurate for lighter heavy-quark masses, its performance with different trial wave functions (Slater vs. Gaussian) diverges at higher masses due to the neglect of non-adiabatic corrections.

The Gist

The Heavy-Handed Approximation: A Story of Quarks, Molecules, and Guessing Games

Imagine you are trying to understand how a complex machine works. You have two heavy, slow-moving gears (the heavy quarks) and a swarm of tiny, hyperactive bees buzzing around them (the light quarks). To figure out how the whole machine moves, you have two main strategies.

This paper is essentially a "stress test" for one of those strategies, called the Born-Oppenheimer Approximation (BOA), to see if it works well for the strange world of "doubly heavy hadrons" (particles made of two heavy quarks and some light ones).

Here is the breakdown in simple terms:

1. The Two Strategies

• The "Slow-Motion" Strategy (Born-Oppenheimer Approximation):
  Imagine the heavy gears are so massive and slow that they barely move at all. You freeze them in place. Then, you watch how the bees buzz around the frozen gears. You calculate the energy of the bees, and then you let the gears move, using the bees' energy as a guide.

  • Why use it? It's much easier to calculate. It's like taking a snapshot of a moving car to study its wheels.
  • The Problem: It assumes the heavy parts are infinitely heavy compared to the light parts. In a hydrogen molecule (a proton and an electron), this works perfectly because the proton is 1,800 times heavier than the electron. But in a particle with two heavy quarks, the heavy quark is only about 5 times heavier than the light one. That's like comparing a bowling ball to a basketball—not quite "infinitely heavy."

• The "Full-Blow" Strategy (Gaussian Expansion Method - GEM):
  This is the "gold standard." Instead of freezing anything, you let the heavy gears and the light bees move all at once, interacting with each other in real-time. You solve the math for the whole chaotic system together.

  • Why use it? It's incredibly accurate but computationally expensive (like simulating every single atom in a car crash).
  • The Paper's Role: The authors use this method as the "Truth Benchmark." They want to see how close the "Slow-Motion" strategy gets to the "Full-Blow" truth.

2. The Experiment: Testing the Guess

The authors ran simulations on two types of systems:

1. Hydrogen Molecules (The Control Group): They tested their methods on real atoms. As expected, when the heavy part is much heavier than the light part, the "Slow-Motion" guess was almost perfect. But as the heavy part got lighter (closer to the light part's weight), the guess started to drift away from the truth.
2. Doubly Heavy Hadrons (The Real Test): They applied this to particles like the $\Xi_{cc}$ (two charm quarks) and $\Xi_{bb}$ (two bottom quarks).

3. The Twist: It Depends on Your "Lens"

Here is where it gets interesting. When using the "Slow-Motion" strategy, you have to choose a mathematical "lens" (a trial wave function) to describe how the light quarks behave. The authors tested two lenses:

• Lens A: Slater-Type Functions (STFs)

  • Analogy: Think of this lens as a flashlight with a sharp, focused beam. It's great at seeing details close up (short distances) but gets fuzzy far away.
  • The Result: When the heavy quarks got very heavy (like the bottom quark), this lens overestimated the binding energy. It made the particle look like it was holding on too tightly. It was like the flashlight missed the fact that the "glue" (confinement) gets weaker at long distances.

• Lens B: Gaussian-Type Functions (GTFs)

  • Analogy: Think of this lens as a soft, diffused fog. It spreads out nicely but blurs the sharp edges up close.
  • The Result: When the heavy quarks got very heavy, this lens underestimated the binding energy. It made the particle look like it was holding on too loosely. The authors found this happened because this lens ignores a subtle effect called "non-adiabatic corrections"—basically, it forgets that the heavy gears aren't perfectly frozen; they wiggle a little, and that wiggle matters.

4. The Big Takeaway

The paper concludes that the "Slow-Motion" strategy (Born-Oppenheimer) is a good qualitative tool (it gives you the right general idea) for particles with charm quarks.

However, as the heavy quarks get heavier (like bottom quarks), the strategy starts to break down quantitatively:

• If you use the "Sharp Beam" lens, you think the particle is too stable.
• If you use the "Foggy" lens, you think the particle is too unstable.

The Final Verdict:
If you want a quick, rough estimate for heavy particles, the Born-Oppenheimer approximation is fine. But if you want high-precision predictions (like predicting the exact mass of a new particle to match an experiment), you cannot rely on the approximation. You must use the "Full-Blow" method (GEM) that treats all the moving parts together without freezing anything.

In short: The approximation is like using a map of a city to navigate a single room. It works if the room is huge and simple, but if you need to find a specific coffee cup in a cluttered room, you need to look at the room itself, not the map.

Technical Summary

1. Problem Statement

The Born-Oppenheimer Approximation (BOA) is a standard tool in molecular physics, relying on a large mass hierarchy between heavy nuclei and light electrons to separate degrees of freedom. In Quantum Chromodynamics (QCD), BOA has been widely applied to study doubly heavy hadrons (doubly heavy baryons and tetraquarks), treating heavy quarks as static sources for light quark dynamics.

However, the validity of this approximation in hadronic systems is questionable for two main reasons:

1. Mass Hierarchy: Unlike the proton-electron mass ratio ($\sim 1800$), the heavy quark ($c$ or $b$) to light quark ($u, d, s$) mass ratio is much smaller (e.g., $m_c/m_{u/d} \approx 5$). This weak scale separation challenges the core assumption of BOA.
2. Basis Function Dependence: Previous studies using BOA in potential models have yielded varying results depending on the choice of trial wave functions (Slater-type vs. Gaussian-type), without a rigorous benchmark to determine which is more accurate or if the approximation itself introduces systematic errors.

The paper aims to quantify the validity of BOA within potential models by comparing it against a fully dynamical method and analyzing the impact of trial wave function choices.

2. Methodology

The authors employ a comparative approach using three distinct methods within a unified potential model framework:

• Benchmark Method: Gaussian Expansion Method (GEM)

  • Solves the multi-body Schrödinger equation without assuming a separation of heavy and light degrees of freedom.
  • Uses a complete set of Gaussian basis functions to capture full correlations between constituent quarks.
  • Serves as the "exact" numerical benchmark for the study.

• Approximation Method: Born-Oppenheimer Approximation (BOA)

  • Decomposes the Hamiltonian into heavy ($H_{heavy}$) and light ($H_{light}$) components.
  • Solves the light quark dynamics for fixed heavy quark positions to obtain an effective potential $\epsilon(r_{AB})$.
  • Solves the heavy quark motion within this effective potential.
  • Variational Wave Functions: The authors test two types of trial wave functions for the light quarks within the BOA framework:
    1. Slater-Type Functions (STFs): Known for accurate short-distance behavior (finite slope) and exact hydrogen atom solutions.
    2. Gaussian-Type Functions (GTFs): Preferred in hadron physics for rapid long-range suppression.

• Systematic Study:

  • Hydrogen-like Systems: Validated the methods on hydrogen molecular ions and molecules to establish the relationship between mass ratios and BOA accuracy.
  • Hadronic Systems: Applied the methods to $S$-wave doubly heavy baryons ($QQq$) and tetraquarks ($QQ\bar{q}\bar{q}$).
  • Hamiltonian: Used a general potential model including Coulombic, linear confinement, and hyperfine interactions. Parameters were fixed by reproducing known meson masses ($D, B$ mesons) using GEM.

3. Key Contributions

• Quantitative Benchmarking: The first comprehensive comparison of BOA-STFs, BOA-GTFs, and GEM results for doubly heavy hadrons within the same potential model, using GEM as a rigorous benchmark.
• Mass Dependence Analysis: Demonstrated that the accuracy of BOA is highly sensitive to the heavy-quark mass ($m_Q$).
• Wave Function Sensitivity: Revealed that the choice of trial wave function in BOA leads to divergent trends as $m_Q$ increases:
  • BOA-STFs tend to overestimate binding energy (predicting lower masses/higher binding).
  • BOA-GTFs tend to underestimate binding energy (predicting higher masses/lower binding).
• Physical Interpretation: Identified the physical origins of these deviations:
  • GTF Underestimation: Primarily due to the neglect of non-adiabatic corrections inherent in the BOA.
  • STF Overestimation: Arises because STFs are not well-suited to describe the long-range behavior induced by the confining potential in hadronic systems.

4. Key Results

A. Hydrogen-like Systems (Validation)

• When the mass ratio $m_{heavy}/m_{light}$ is large, BOA results (both STF and GTF) converge with GEM.
• As the mass ratio decreases (approaching hadronic scales), BOA predictions deviate significantly, generally yielding deeper binding energies than GEM.

B. Doubly Heavy Baryons ($QQq$) and Tetraquarks ($QQ\bar{q}\bar{q}$)

• Small Heavy Quark Mass ($m_Q \approx m_c$): All three methods (BOA-STF, BOA-GTF, GEM) yield comparable results.
• Large Heavy Quark Mass ($m_Q \approx m_b$ or higher): Divergence becomes pronounced:
  • BOA-STF: Predicts eigenvalues and masses that are larger (more bound) than GEM.
  • BOA-GTF: Predicts eigenvalues and masses that are smaller (less bound) than GEM.
  • GEM: Provides the intermediate, most reliable value.

Specific Examples (from Tables IV & V):

• For the $\Xi_{cc}$ (doubly charmed baryon):
  • GEM Mass: $\sim 3566$ MeV.
  • BOA-STF Mass: $\sim 3550$ MeV (Over-binding).
  • BOA-GTF Mass: $\sim 3525$ MeV (Under-binding).
• For the $T_{cc}$-like tetraquark ($cc\bar{n}\bar{n}$):
  • GEM Mass: $\sim 4051$ MeV.
  • BOA-STF Mass: $\sim 4088$ MeV.
  • BOA-GTF Mass: $\sim 4043$ MeV.

C. Trend Reversal

The study highlights a critical finding: as the heavy quark mass increases, the error direction flips between the two basis choices. STFs become increasingly inaccurate by over-binding, while GTFs become inaccurate by under-binding due to non-adiabatic effects.

5. Significance and Conclusion

• Limitations of BOA: While BOA offers an intuitive physical picture for doubly heavy hadrons, its quantitative accuracy is limited, especially for bottom ($b$) and heavier systems. The approximation breaks down as the mass hierarchy weakens or when specific basis functions fail to capture long-range confinement dynamics.
• Recommendation: For qualitative studies of charm systems, BOA may provide reasonable estimates. However, for quantitative predictions (especially for bottom systems or precision spectroscopy), fully dynamical treatments like the Gaussian Expansion Method (GEM) are indispensable.
• Model Dependence: The authors clarify that their conclusions are specific to the potential model framework used. They do not invalidate QCD-based BOA studies (e.g., lattice QCD or effective field theories) which may handle non-adiabatic effects differently, but they emphasize that within potential models, the choice of basis function and the neglect of non-adiabatic terms are major sources of uncertainty.

In summary, this work provides a critical "reality check" for potential model calculations of doubly heavy hadrons, establishing GEM as the necessary benchmark and warning against the uncritical application of BOA with specific trial wave functions for high-precision mass predictions.