The $T_{bc}$ tetraquarks near the $B\bar{D}$ threshold

Using the dynamic diquark model with a Born-Oppenheimer potential derived from lattice QCD, this study predicts that the doubly heavy scalar tetraquark $T_{bc}^{(0)}$ lies near the $B\bar{D}$ threshold as a potential bound state or narrow resonance, while the $T_{bc}^{(1)}$ axial-vector state is a compact $S$-wave resonance situated approximately 23–28 MeV above the $B^{*}\bar{D}$ threshold.

The Gist

Imagine the universe is made of tiny, fundamental Lego bricks called quarks. Normally, these bricks assemble in two standard ways: either two bricks form a "meson" (like a cousin of the proton) or three bricks form a "baryon" (like a proton or neutron). For decades, physicists believed these were the only ways to build stable structures.

However, in the last twenty years, scientists have begun to find "exotic" structures made of four glued-together bricks. These are called tetraquarks. It is like finding a stable Lego house made of four bricks instead of the usual two or three.

This work by Halil Mutuk is a theoretical investigation of a very specific, rare type of these four-brick houses. Here is a breakdown of what they did and what they found, using simple analogies.

1. The special "heavy" tetraquark

Most exotic particles found so far consist of light bricks. This work examines a "heavy" version called $T_{bc}$.

• The ingredients: Imagine a heavy brick made of two antiquarks (a bottom and a charm quark) acting as a heavy core, plus a light pair of quarks (up and down) serving as a light shell.
• The structure: The scientists modeled this as a heavy "antidiquark" (the heavy core) and a light "diquark" (the light shell) holding hands.

2. The method: The "slow and fast" dance

To determine how heavy this particle is and how it behaves, the authors used a method called the Born-Oppenheimer approximation.

• The analogy: Think of a heavy elephant (the heavy quarks) walking slowly through a field, while a swarm of faster, buzzing bees (the light quarks and gluons) zip around it at lightning speed.
• How it works: Since the elephant moves so slowly, the bees adjust almost instantaneously to its position. The bees create an invisible "force field" (a potential) that determines how the elephant can move. The scientists calculated the energy of this dance to predict the weight of the resulting particle.

3. The two predicted particles

The study predicts two specific versions of this $T_{bc}$ particle, differing in how their internal "spins" (a quantum property like a tiny magnet) are arranged:

• The scalar state ($0^+$): This is the "quiet" version.

  • The prediction: It weighs about 7.14 to 7.16 GeV.
  • The location: It sits almost exactly at the "edge" of a cliff called the $B\bar{D}$ threshold.
  • What this means: It is so close to the edge that it is hard to say whether it is a stable, bound particle (sitting securely on the ground) or a fleeting "resonance" (a temporary wobble right at the edge). If it is stable, it would be incredibly long-lived because it cannot easily decay into lighter parts.

• The axial vector state ($1^+$): This is the "rotating" version.

  • The prediction: It weighs about 7.22 GeV.
  • The location: It sits clearly above another threshold ($B^*\bar{D}$) but below yet another ($B\bar{D}^*$).
  • What this means: It is definitely a "resonance." It is like a ball rolling in a shallow depression just above a hill. It will exist for a short time and then decay (fall apart) into other particles. The work predicts it will appear as a distinct bump in experimental data, but its shape will be distorted because it lies so close to the edge of the hill.

4. How tight is the grip?

The scientists calculated the size of these particles.

• The result: They are very small and compact, with a radius of about 0.45 femtometers (a femtometer is one quadrillionth of a meter).
• The analogy: This is much smaller than a "loose molecule," where two separate particles hold hands only from a distance. Instead, these four bricks are firmly fused into a single, dense lump. It is like a tightly packed suitcase compared to two suitcases tied together with a long rope.

5. Why the difference?

The work explains that the weight difference between the "quiet" and "rotating" versions comes from two things:

1. Mass difference: The heavy core is slightly heavier when the spins are aligned one way or the other.
2. Magnetic interaction: The quarks have tiny magnetic properties. When they interact, this adds a small amount of energy. The study found that the "rotating" version is about 60 to 80 MeV heavier than the "quiet" one.

6. The big picture

The authors compare their results with other current studies (such as those using supercomputers called lattice QCD).

• Agreement: Their predictions fit well within the range of other theories.
• Deviation: Their "rotating" particle is predicted to be somewhat heavier (by about 30–70 MeV) than some current supercomputer calculations suggest. The authors propose this might be because their model treats the particles as a single solid unit, whereas supercomputer models may capture subtle, long-range interactions between the particles that their model simplifies.

Conclusion

In short, this work predicts that nature has two new, heavy, four-quark particles waiting to be discovered.

• One is a tightly bound, compact object sitting right at the edge of stability and possibly very difficult to detect because it hardly decays.
• The other is a short-lived resonance that should show up clearly in particle accelerators like the LHCb or Belle-II experiments.

The authors essentially say: "If you look at the data around 7.15 GeV and 7.22 GeV, you should see these specific patterns. Their detection would prove that these four bricks can indeed stick together in a firm, compact knot."

Technical Summary

Problem Description
The existence and nature of doubly heavy open-flavor tetraquarks, specifically the $T_{bc}$ system composed of a light diquark ($ud$) and a heavy antidiquark ($\bar{b}\bar{c}$), remain subjects of theoretical debate. While recent Lattice QCD studies have provided evidence for poles near the threshold in the $I(J^P) = 0(1^+)$ channel, results regarding binding energy, the existence of the $0^+$ partner, and the distinction between genuine bound states, virtual states, or resonances are inconsistent across various calculations. Furthermore, theoretical predictions from quark models and QCD sum rules vary widely, ranging from unbound to deeply bound. The $T_{bc}$ system is particularly critical as it lies at the boundary between bound and unbound configurations, offering a stringent test of the dynamics responsible for binding exotic hadrons.

Methodology
The authors investigate the states $T_{bc}^{(0)}$ ($J^P=0^+$) and $T_{bc}^{(1)}$ ($J^P=1^+$) using the dynamic diquark model within the Born-Oppenheimer (BO) approximation.

• System Description: The tetraquark is modeled as a bound state of a colored light diquark ($\delta = [ud]$) and a colored heavy antidiquark ($\bar{\delta}' = [\bar{b}\bar{c}]$).
• Potential: The interaction between the clusters is governed by the static $\Sigma_g^+(1S)$ BO potential, extracted from Lattice QCD calculations of the static potential in the fundamental representation.
• Hamiltonian: The effective Hamiltonian includes a kinetic term, the BO potential, and a spin-dependent term of the chromomagnetic interaction ($H_\kappa$). The model accounts for the specific spin configurations of the light diquark (bound to isospin $I=0$ for the "good" diquark) and the heavy antidiquark (allowing both scalar $s_{\bar{b}\bar{c}}=0$ and axial-vector $s_{\bar{b}\bar{c}}=1$ configurations).
• Numerical Solution: The radial Schrödinger equation is solved numerically on a discretized grid to obtain eigenvalues ($E_n$). Physical masses are calculated by adding the constituent masses and the expectation values of the chromomagnetic interaction.
• Parameters: The study uses inputs from QCD sum rules and meson phenomenology for diquark masses and the coupling of the light diquark ($\kappa_{ud}$). The chromomagnetic coupling of the heavy antidiquark ($\kappa_{\bar{b}\bar{c}}$) is varied over three sets (10, 15, and 20 MeV) to examine sensitivity.

Main Contributions and Results

1. Mass Predictions:

   • The scalar state $T_{bc}^{(0)}$ is predicted in the narrow range 7.143–7.158 GeV.
   • The axial-vector state $T_{bc}^{(1)}$ is predicted in the range 7.217–7.222 GeV.
   • The total spread caused by the variation of $\kappa_{\bar{b}\bar{c}}$ is small ($\sim 15$ MeV for $0^+$ and $\sim 5$ MeV for $1^+$), indicating that the spectrum is dominated by constituent masses and the BO potential.

2. Hyperfine Splitting:

   • The hyperfine splitting $\Delta_{HF} = M(T_{bc}^{(1)}) - M(T_{bc}^{(0)})$ is calculated across parameter sets as 59–79 MeV.
   • The splitting is driven primarily by the kinematic mass difference between the symmetric and antisymmetric configurations of the heavy antidiquark ($\sim 40$ MeV), with a linear chromomagnetic contribution ($\partial \Delta_{HF} / \partial \kappa_{\bar{b}\bar{c}} = 2$).

3. Spatial Structure:

   • The mean cluster separation is $\langle r \rangle \simeq 0.45\text{--}0.46$ fm, and the reciprocal mean radius is $\langle 1/r \rangle^{-1} \simeq 0.33\text{--}0.34$ fm.
   • These values are weakly dependent on parameters and significantly smaller than typical hadronic molecular scales ($\sim 1$ fm), supporting a compact diquark-antidiquark interpretation rather than a loosely bound molecular structure.
   • The ratio $\langle 1/r \rangle^{-1} / \langle r \rangle \simeq 0.73\text{--}0.74$ suggests a wavefunction form lying between pure Coulomb and harmonic oscillator regimes, implying comparable contributions from short-range gluon exchange and long-range confinement effects.

4. Phenomenology Relative to Thresholds:

   • Scalar State ($0^+$): The predicted mass spans the $B\bar{D}$ threshold (7149 MeV). Depending on the specific parameter set, the state lies either slightly above or slightly below the threshold. Consequently, it may appear either as a weakly decaying bound state (if below threshold) or as a narrow resonance near the threshold.
   • Axial-Vector State ($1^+$): The state is robustly predicted 23–28 MeV above the $B^*\bar{D}$ threshold and ~70 MeV below the $B\bar{D}^*$ threshold. This positions it as an S-wave resonance in the $B^*\bar{D}$ channel, where the line shape is expected to be strongly influenced by the nearby threshold cusp.

Significance and Claims
The work claims that the dynamic diquark model, supplemented by Lattice QCD BO potentials, provides a coherent and falsifiable framework for the $T_{bc}$ spectrum with minimal parametric freedom.

• Predictive Power: Once the mass gap of the heavy antidiquark is fixed, the linear response of the hyperfine splitting to the chromomagnetic coupling offers a potential method to utilize future experimental data to constrain $\kappa_{\bar{b}\bar{c}}$, a parameter currently inaccessible to Lattice QCD or phenomenology.
• Experimental Signatures: The authors suggest that a resonant structure near 7.15 GeV in the $B\bar{D}$ invariant mass (within a few MeV of the threshold) would be a candidate for the scalar state, while a structure near 7.22 GeV in the $B^*\bar{D}$ channel (approx. 25 MeV above threshold) would signal the axial-vector state. The predicted mass hierarchy and specific branching asymmetry (favoring $B^*\bar{D}$ over $B\bar{D}^*$ for the $1^+$ state) serve as distinguishing features against molecular interpretations.
• Limitations: The authors acknowledge that their single-channel BO treatment neglects long-range one-pion exchange and explicit dynamical coupling channels, which could explain the $\sim 30\text{--}70$ MeV discrepancy between their predictions and some recent Lattice QCD results on binding energy. They regard their prediction as the value for the constituent $\delta\text{--}\bar{\delta}'$, noting that a comparison with Lattice results with an accuracy better than $O(\Lambda_{QCD})$ requires the inclusion of these omitted effects.