⚛️ lattice

Properties of the positive and negative parity charm-strange and bottom-strange mesons $D_s$ , $D_s^$ , $D_{s0}^$ , $D_{s1}$ , $B_s$ , $B_s^$ , $B_{s0}^$ , $B_{s1}$ from lattice QCD: masses, decay constants, and compositeness

This paper presents a comprehensive lattice QCD study using RBC/UKQCD ensembles to determine the masses, decay constants, and compositeness parameters of positive- and negative-parity charm-strange and bottom-strange mesons, providing the first lattice results for $f_{B^*_{s0}}$ and $f_{B_{s1}}$ and confirming that the positive-parity states are predominantly molecular in nature.

Original authors: Forrest Guyton, Stefan Meinel

Published 2026-09-11

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Forrest Guyton, Stefan Meinel

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

In the subatomic world, matter is built from a small family of fundamental particles called quarks. Most of the matter we see around us, including the protons and neutrons in our own bodies, is made of quarks bound together in groups of three. However, nature also allows quarks to pair up, forming short-lived particles known as mesons. For decades, physicists have used a framework called the quark model to predict how these pairs should behave, assuming they are simple, compact objects held together by a single force. But in 2003, an experiment discovered a meson that did not fit this simple picture. This particle, named $D^*_{s0}(2317)$ , was found to be lighter than theory predicted and sitting just below the energy threshold where it should have broken apart into two other particles. This anomaly suggested that the particle might not be a simple pair of quarks at all, but rather a loose, fragile molecule made of two distinct pairs of quarks orbiting each other. Understanding whether these strange particles are compact dots or loose molecules is crucial, as it tests the limits of our understanding of the strong force that binds the universe together.

A team of researchers at the University of Arizona has now taken a deep dive into this mystery, using a powerful computational technique called lattice quantum chromodynamics to simulate the behavior of these particles from first principles. Instead of observing particles in a physical accelerator, the scientists built a virtual grid of space and time on supercomputers, filling it with the mathematical rules that govern quarks and gluons. They focused on a specific family of particles containing a heavy quark (either charm or bottom) and a strange quark. Their goal was to measure the masses and internal structures of the lightest versions of these particles, including both the standard, stable ones and the more exotic, positive-parity states that had caused so much confusion. By running these simulations across seven different virtual environments with varying sizes and levels of detail, the team was able to extrapolate their findings to the real world, effectively removing the artificial constraints of their computer grid to see what nature actually looks like.

The results provide a detailed map of these heavy-strange mesons, confirming the existence of the elusive states and offering a new perspective on their composition. For the standard, negative-parity particles, the team calculated their decay constants, which are essentially measures of how tightly the quarks are bound and how likely the particle is to decay. These values matched well with previous world averages, validating the reliability of their method. However, the real breakthrough came in the positive-parity sector. The researchers found that the binding energies of the exotic states—how much energy is required to pull them apart—were consistent with the particles being just barely held together. Specifically, the charm-strange particle $D^*_{s0}$ was found to be bound by about 48 MeV, and the bottom-strange partner $B^*_{s0}$ by about 69 MeV. These numbers place the particles just below the point where they would naturally fall apart, a hallmark of a molecular structure rather than a tight, compact one.

To determine the true nature of these particles, the team applied a mathematical test known as Weinberg's compositeness criterion. This test analyzes the relationship between the particle's binding energy and its size to estimate how much of the particle is made of the two pairs of quarks versus a single, compact core. The simulations yielded a very small value for the compact core component, suggesting that these particles are indeed predominantly molecular. In other words, the data points to a structure where two separate meson pairs are orbiting each other in a loose embrace, rather than a single, dense blob of quarks. This finding supports the idea that the $D^*_{s0}(2317)$ and its relatives are exotic molecules, a conclusion that aligns with their proximity to the energy thresholds where they could split apart.

The study also provided the first-ever lattice QCD calculations for the decay constants of the bottom-strange positive-parity mesons, $B^*_{s0}$ and $B_{s1}$ . These values are essential for understanding how these particles interact and decay, filling a gap in the theoretical knowledge of the bottom quark sector. While the results for the charm particles showed a clear molecular character, the bottom particles presented a slightly more complex picture. The data still favored a molecular interpretation, but the uncertainties were larger, leaving room for the possibility that these heavier states might have a more significant compact core than their lighter cousins. The researchers noted that their analysis had limitations, particularly regarding certain complex interactions that become important at very low energies, which means the molecular picture, while strongly supported, is not yet proven beyond all doubt.

Ultimately, this work represents a significant step forward in mapping the landscape of heavy mesons. By combining high-precision simulations with rigorous statistical methods, the team has moved beyond speculation to provide concrete numbers that describe the mass, binding energy, and internal makeup of these particles. The findings suggest that the subatomic world is richer and more varied than the simple quark model once implied, hosting particles that behave like fragile molecules held together by the delicate balance of the strong force. As experimentalists continue to refine their measurements of these particles, the precise numbers provided by this study will serve as a critical benchmark, helping to confirm whether nature truly builds these exotic states from loose pairs of quarks or if a more complex story is still waiting to be uncovered.

Technical Summary: Properties of Positive and Negative Parity Charm-Strange and Bottom-Strange Mesons from Lattice QCD

Problem Statement
Since the 2003 discovery of the $D^*_{s0}(2317)^\pm$ by the BaBar collaboration, there has been significant interest in the nature of positive-parity heavy-strange mesons. The $D^*_{s0}(2317)^\pm$ lies approximately 45 MeV below the $DK$ kinematic threshold, a position far lower than the $^3P_0$ $q\bar{q}$ state predicted by potential models. This discrepancy has prompted interpretations of these states as exotic, potentially molecular structures (four-quark states consisting of two shallowly bound color-singlet pairs). While experimental data supports this view through features like proximity to thresholds and strong coupling to meson-meson channels, a rigorous theoretical determination of their properties—specifically their masses, decay constants, and compositeness—remains challenging. Furthermore, the positive-parity bottom-strange sector ( $B^*_{s0}$ and $B_{s1}$ ) is less explored experimentally, with ground states not yet observed, necessitating lattice QCD predictions.

Methodology
The authors perform a comprehensive lattice QCD calculation using seven ensembles of gauge configurations generated by the RBC/UKQCD collaboration. The study employs:

Fermion Actions: Domain-wall fermions for light and strange quarks, and anisotropic clover actions for charm and bottom quarks. The heavy-quark actions are nonperturbatively tuned to remove discretization errors, allowing work at physical charm and bottom masses across all lattice spacings.
Ensembles: Seven ensembles with pion masses ranging from 431 MeV to 139 MeV (near-physical) and lattice spacings from 0.114 fm to 0.073 fm. This allows for combined chiral and continuum extrapolations.
Operator Basis: To resolve the low-lying states, the authors utilize a correlation matrix constructed from three types of interpolating operators:
1. Standard $q\bar{q}$ operators.
2. Operators with covariant derivatives.
3. Meson-meson scattering operators (e.g., $D K$ , $D^* K$ ) at both source and sink.
Analysis Techniques:
- Generalized Eigenvalue Problem (GEVP): Used to extract finite-volume energy levels and optimized eigenvectors from the correlation matrices.
- Lüscher's Method: Applied to the finite-volume spectra to extract infinite-volume scattering phase shifts and determine bound-state masses. The authors parametrize the scattering amplitudes using Effective-Range Expansions (ERE).
- Decay Constants: Extracted using optimized operators derived from the GEVP eigenvectors.
- Compositeness: The Weinberg compositeness criterion is applied to evaluate the "molecular" nature of the states, calculating the field-strength renormalization $Z$ (where $1-Z$ represents the probability of the state being a composite molecule).

Key Contributions

First Lattice Determination of $B^*_{s0}$ and $B_{s1}$ Decay Constants: The paper provides the first lattice QCD calculations for the decay constants of the positive-parity bottom-strange mesons.
Comprehensive Positive-Parity Spectrum: The study simultaneously analyzes the scalar ( $0^+$ ) and axial-vector ( $1^+$ ) states for both charm-strange ( $D^*_{s0}, D_{s1}$ ) and bottom-strange ( $B^*_{s0}, B_{s1}$ ) systems, including their binding energies and decay constants.
Systematic Treatment of Finite Volume and Chiral Limits: By utilizing multiple ensembles with varying pion masses and lattice spacings, the authors perform robust chiral-continuum extrapolations. They also address finite-volume effects on decay constants using chiral perturbation theory-inspired corrections.
Evaluation of Exotic Nature: The authors rigorously test the molecular hypothesis using Weinberg's criterion, calculating the coupling constants ( $g$ ) and compositeness parameters ( $Z$ ) for all four positive-parity states.

Results

Negative-Parity Mesons: The calculated decay constants for the ground states are:
- $f_{D_s} = 251.4(2.1)(0.4)(2.5)$ MeV
- $f_{D^*_s} = 272.5(4.3)(1.0)(2.7)$ MeV
- $f_{B_s} = 228.4(5.8)(0.5)(2.3)$ MeV
- $f_{B^*_s} = 229.2(4.4)(0.8)(2.3)$ MeV
  These results are consistent with FLAG world averages. The ratios $f_{D^*_s}/f_{D_s} = 1.086(14)(11)$ and $f_{B^*_s}/f_{B_s} = 1.003(22)(10)$ are also determined.
Positive-Parity Binding Energies: After extrapolation to the physical point, the binding energies ( $m_{H^*} - m_{H^{(*)}} - m_K$ $m_{H^{*}} - m_{H^{(*)}} - m_{K}$ ) are:
- $D^*_{s0}$ : $-48(14)(4)$ MeV
- $D_{s1}$ : $-61(15)(2)$ MeV
- $B^*_{s0}$ : $-69(13)(4)$ MeV
- $B_{s1}$ : $-77(10)(5)$ MeV
  The charm-strange results are consistent with experimental values.
Positive-Parity Decay Constants:
- $f_{D^*_{s0}} = 136.6(8.0)(4.0)(1.4)$ MeV
- $f_{D_{s1}} = 200(33)(24)(2)$ MeV
- $f_{B^*_{s0}} = 207(12)(8)(2)$ MeV
- $f_{B_{s1}} = 196(16)(11)(2)$ MeV
Compositeness and Molecular Interpretation: The analysis yields small values for the field-strength renormalization $Z$ (e.g., $Z \approx 0.03$ for $D^*_{s0}$ in the zeroth-order ERE), implying large compositeness ( $1-Z \approx 0.97$ ). The scattering lengths satisfy $1/a_0 \sim -\gamma$ , and effective ranges are small and positive. These findings are consistent with an interpretation of the positive-parity states as predominantly molecular. However, the authors note that uncertainties in the first-order ERE fits are large, particularly for the bottom sector, and corrections to the weak-binding limit may be significant.

Significance and Claims
The paper claims to provide the first lattice QCD determination of the decay constants for the positive-parity bottom-strange mesons ( $B^*_{s0}$ and $B_{s1}$ ). It offers a unified, high-precision study of the heavy-strange spectrum that confirms the molecular interpretation of the $D^*_{s0}(2317)$ and $D_{s1}(2460)$ states, extending this analysis to the bottom sector. The results for binding energies in the charm sector align with experimental data, while the bottom sector predictions await experimental confirmation. The authors explicitly state that while their analysis supports a molecular interpretation, limitations exist: left-hand cuts in scattering amplitudes (relevant at low pion mass) are neglected, and corrections to the weak-binding limit in Weinberg's criterion are expected to be large in some cases. Consequently, while the results "hint at a considerable molecular content," particularly for $D^*_{s0}$ , the situation for the $B_{s1}$ remains less clear due to large uncertainties in the first-order ERE analysis.

More like this