⚛️ phenomenology

Singly heavy tetraquarks

This paper employs a semi-relativistic hybrid quark potential model to systematically investigate the spectra and fall-apart decays of singly-heavy tetraquarks, concluding that while most predicted states are compact and lie above meson-meson thresholds, specific experimental resonances like $D_{s0}(2317)$ and $T_{b\bar{s}}(5568)$ cannot be explained as compact tetraquarks, whereas others such as $T_{\bar{c}\bar{s}0}(2870)$ and $T_{c\bar{s}0}(2900)$ are consistent with the tetraquark interpretation.

Original authors: Jun-Jie Liu, Zhi-Biao Liang, Feng-Xiao Liu, Mu-Yang Chen, Xian-Hui Zhong, Qiang Zhao

Published 2026-07-29

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Jun-Jie Liu, Zhi-Biao Liang, Feng-Xiao Liu, Mu-Yang Chen, Xian-Hui Zhong, Qiang Zhao

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is built from tiny, invisible LEGO bricks called quarks. Usually, these bricks snap together in very predictable ways: two bricks make a "meson" (like a proton's cousin), and three bricks make a "baryon" (like a proton or neutron). But sometimes, nature gets creative and snaps four bricks together to form a "tetraquark." For decades, physicists have been hunting for these exotic four-brick structures, wondering if they are just loose clusters of two mesons bumping into each other, or if they are tight, compact little balls of four quarks glued together. The big question is: what holds these four bricks together? Is it just the standard "glue" of the strong force, or do other invisible forces, like exchanging tiny messenger particles, play a bigger role? Understanding these strange particles helps scientists decode the secret rules of how matter holds itself together, a puzzle that has stumped physicists for years.

In this new study, a team of researchers decided to build a massive digital simulation to map out the entire family of "singly-heavy" tetraquarks. Think of these as exotic four-brick structures where one brick is a heavy, stubborn "heavy quark" (either a charm or a bottom quark), and the other three are lighter, more energetic "light quarks." The team didn't just look at the mass (how heavy the particle is); they also calculated how these particles might fall apart. They used a sophisticated "hybrid" model, which is like a recipe that mixes two different types of glue: the standard "one-gluon exchange" (the usual strong force) and "one-boson exchange" (forces carried by particles like pions and rho mesons). By running these calculations, they predicted the exact masses and decay patterns for every possible arrangement of these four-quark systems.

The results of their simulation paint a very clear picture. First, they found that these tetraquarks are indeed compact, tight little balls, not loose clouds. They sit at high energy levels, far above the point where they would naturally break apart into two separate mesons. Because they are so compact, the researchers predict they are surprisingly stable, with "fall-apart" widths (a measure of how quickly they decay) ranging from a tiny 1 MeV up to about 120 MeV. This suggests that if we build these particles in a lab, they might stick around long enough to be spotted by detectors.

However, the simulation also acted as a detective, ruling out some popular suspects. The researchers found that several famous particles reported by experiments—specifically the $D_{s0}(2317)$ , $D_{s1}(2460)$ , $T_{b\bar{s}}(5568)$ , and $T_{c\bar{s}}(2327)$ —cannot be explained as these compact tetraquarks. Their predicted masses are simply too high to match the experimental data for these specific particles. In other words, if these particles exist, they are likely something else entirely, not the tight four-quark balls the team modeled.

On the flip side, the study offers some exciting new leads. The researchers suggest that two recently observed resonances, $T_{\bar{c}\bar{s}0}(2870)$ and $T_{c\bar{s}0}(2900)$ , are very likely to be these compact tetraquarks. Their simulations predict masses and decay patterns that line up well with what experiments have seen. Specifically, they propose that $T_{\bar{c}\bar{s}0}(2870)$ is a particle with specific quantum numbers (spin and isospin) and that $T_{c\bar{s}0}(2900)$ is its partner with slightly different properties.

The team also mapped out a whole zoo of other potential tetraquarks that haven't been seen yet. They predict that many of these new states should be observable in specific decay channels, such as particles breaking apart into a $D$ meson and a pion, or a $B$ meson and a kaon. For instance, they suggest looking for a narrow state around 2.2 GeV in the $D\pi$ channel, or a bottom-quark partner around 5.6 GeV in the $B\pi$ channel. While the paper doesn't claim these have been found yet, it provides a detailed "treasure map" for experimentalists at facilities like the LHC, telling them exactly where to look and what signals to expect. The study emphasizes that including the "boson exchange" forces was crucial; without them, the predictions would have been completely off, showing that these extra forces are essential for understanding how these exotic particles behave.

Technical Summary: Systematic Study of Singly-Heavy Tetraquark Spectra and Decays

Problem Statement
The nature of exotic hadrons, particularly singly-heavy tetraquarks ( $Qq\bar{q}\bar{q}$ ), remains a subject of intense theoretical and experimental investigation. While numerous candidates have been observed (e.g., $D_{s0}(2317)$ , $D_{s1}(2460)$ , $X(5568)$ , $T_{\bar{c}\bar{s}0}(2870)$ , $T_{c\bar{s}0}(2900)$ ), their internal structure—whether compact tetraquarks, hadronic molecules, or kinematic effects—is often ambiguous due to significant model dependencies in theoretical predictions. Existing literature frequently focuses on specific systems, neglects the full spectrum of singly-heavy tetraquarks, or omits decay properties alongside mass spectra. Furthermore, many studies rely solely on One-Gluon Exchange (OGE) potentials, potentially missing crucial dynamics arising from light meson exchanges in systems containing light quarks.

Methodology
The authors perform a systematic study of the $1S$ -wave states for the entire family of singly-heavy tetraquark systems: $Qn\bar{n}\bar{n}$ , $Qs\bar{n}\bar{n}$ , $Qn\bar{s}\bar{n}$ , $Qs\bar{s}\bar{n}$ , $Qn\bar{s}\bar{s}$ , and $Qs\bar{s}\bar{s}$ (where $Q=c, b$ and $n=u, d$ ). The study employs a semi-relativistic hybrid quark potential model and a quark-exchange model for decays.

Hybrid Potential Model: The Hamiltonian includes relativistic kinetic energy terms for light quarks. The interaction potential $V_{ij}$ combines:
- One-Gluon Exchange (OGE): Includes linear confinement, color-Coulomb, and color-magnetic spin-spin interactions.
- One-Boson Exchange (OBE): Incorporates exchanges of pseudoscalar mesons ( $\pi, K, \eta, \eta'$ ), scalar mesons ( $\sigma$ ), and vector mesons ( $\rho, \omega, K^*, \phi$ ). This is a key improvement over previous works that often neglected OBE or treated it simplistically.
- Configuration Space: Wave functions are constructed based on flavor, spin, and color symmetries ( $SU(3)_F \otimes SU(2)_S \otimes SU(3)_C$ ), ensuring total antisymmetry for identical quarks/antiquarks. Two color configurations, $|6\bar{6}\rangle$ and $|\bar{3}3\rangle$ , are considered, allowing for configuration mixing.
Numerical Solution: The four-body Schrödinger equation is solved using the Explicitly Correlated Gaussian (ECG) method with a variational approach. The spatial wave functions are expanded using non-diagonal Gaussian basis functions in Jacobi coordinates.
Decay Analysis: Fall-apart decays ( $A \to B + C$ ) are evaluated using the quark-exchange model. The decay amplitudes are calculated using the genuine tetraquark wave functions obtained from the mass spectrum calculations, rather than simplified Gaussian approximations. The partial widths are derived from the OGE-induced quark rearrangement mechanism.

Key Contributions

Unified Framework: Provides the first systematic calculation of both mass spectra and fall-apart decay widths for the entire set of singly-heavy tetraquark systems within a single, consistent hybrid potential model.
Inclusion of OBE: Demonstrates that One-Boson Exchange potentials (specifically $\pi, \rho, \omega, K^*$ ) play a crucial role in determining the spectrum, particularly for states with specific isospin and color configurations, leading to non-degenerate spectra for different isospin states within the same flavor representation.
Relativistic and Wave Function Improvements: Utilizes a semi-relativistic kinetic term and genuine tetraquark wave functions for decay calculations, enhancing the reliability of width predictions compared to previous non-relativistic or single-Gaussian approaches.

Results

General Characteristics: All calculated $1S$ -wave states are compact (RMS radii $\sim 0.4-0.7$ fm) and lie significantly above the lowest meson-meson dissociation thresholds.
- $J^P = 2^+$ states exhibit extremely narrow fall-apart widths ( $\sim 1$ MeV).
- $J^P = 0^+, 1^+$ states generally have narrow widths ( $\sim 1-120$ MeV), suggesting they are potentially observable.
- Significant configuration mixing occurs between $|6\bar{6}\rangle$ and $|\bar{3}3\rangle$ states, driven by color-magnetic interactions, resulting in large mass splittings ( $\sim 200-600$ MeV) for mixed doublets.
Specific System Findings:
- $Qn\bar{n}\bar{n}$ : Two low-mass states, $T^{1/2}_{(cn[\bar{n}\bar{n}])0^+}(2212)$ and $T^{1/2}_{(bn[\bar{n}\bar{n}])0^+}(5606)$ , are predicted. Their masses overlap with conventional orbitally excited $D$ and $B$ mesons. The authors note that while $T^{1/2}_{(cn[\bar{n}\bar{n}])0^+}(2212)$ shares mass properties with $D_0(2300)$ , its predicted narrow width ( $\sim 41$ MeV) contradicts the broad experimental width of $D_0(2300)$ , suggesting $D_0(2300)$ is not a compact tetraquark.
- $Qn\bar{s}\bar{n}$ : The lowest states are predicted at $\sim 2527$ MeV ( $cn\bar{s}\bar{n}$ ) and $\sim 5850$ MeV ( $bn\bar{s}\bar{n}$ ). These masses are too high to explain $D_{s0}(2317)$ , $D_{s1}(2460)$ , or the $Tb\bar{s}(5568)$ resonance. The $T_{c\bar{s}0}(2900)$ resonance is favored as the $I=1$ scalar state $T^1_{(cn\{\bar{s}\bar{n}\})0^+}(2922)$ (mass $\sim 2922$ MeV), consistent with its observed quantum numbers and decay modes, though the predicted width is slightly narrower than data.
- $Qs\bar{n}\bar{n}$ : The $I^J P = 00^+$ state $T^0_{(\bar{c}\bar{s}[ud])0^+}(2919)$ is identified as a strong candidate for the $T_{\bar{c}\bar{s}0}(2870)$ observed in $D^-K^+$ . Its predicted mass ( $\sim 2919$ MeV) and width ( $\sim 71$ MeV) align well with experimental data. A lower-mass partner $T^0_{(\bar{c}\bar{s}[ud])0^+}(2510)$ is also predicted.
- Other Systems: Predictions are provided for $Qs\bar{s}\bar{n}$ , $Qn\bar{s}\bar{s}$ , and $Qs\bar{s}\bar{s}$ systems, identifying various narrow states with specific dominant decay channels (e.g., $D_s K$ , $B_s \pi$ , $D^*_s \phi$ ).

Significance and Claims
The paper claims that the inclusion of OBE potentials is essential for accurately describing the spectra of singly-heavy tetraquarks, resolving degeneracies found in pure OGE models and shifting mass predictions significantly. The authors assert that their results rule out compact tetraquark interpretations for $D_{s0}(2317)$ , $D_{s1}(2460)$ , $Tb\bar{s}(5568)$ , and $Tc\bar{s}(2327)$ due to mass mismatches. Conversely, they provide strong theoretical support for the tetraquark nature of $T_{\bar{c}\bar{s}0}(2870)$ and $T_{c\bar{s}0}(2900)$ , offering specific mass, width, and decay channel predictions to guide future experimental searches at facilities like LHCb, Belle II, and BESIII. The study emphasizes that many predicted states possess narrow widths, making them viable targets for observation in dominant decay channels.

More like this