Analysis of the hadronic molecules $DK$, $D^*K$, $DK^*$ and their bottom analogs with QCD sum rules

This study employs QCD sum rules with color-singlet currents and dimension-12 vacuum condensates to predict the masses of charm-strange and bottom-strange hadronic molecules, finding that the predicted masses for $DK$, $D^*K$, and $DK^*$ states align well with experimental data for $D_{s0}(2317)$, $D_{s1}(2460)$, and $D_{s1}(2536)$, while the $BK^*$ state is predicted to be a bound molecular state below its threshold.

The Gist

Imagine the universe as a giant, bustling construction site. For decades, physicists have been trying to understand the blueprints of the smallest building blocks: quarks. Usually, these quarks team up in simple pairs (like a proton and an electron) or triplets (like a proton made of three quarks).

But sometimes, nature gets creative and builds "weird" structures where four quarks huddle together. These are called tetraquarks. Think of them not as a tight-knit family of four, but more like a dance couple holding hands with another dance couple. They are loosely bound together, forming a "molecule" made of particles rather than atoms.

This paper is a detective story about finding and measuring these exotic "four-quark molecules."

The Mystery: The "Too Light" Particles

In the early 2000s, scientists discovered some strange particles named $D_s$ mesons (specifically $D_{s0}(2317)$, $D_{s1}(2460)$, and $D_{s1}(2536)$).

Here's the problem: When physicists used their standard "mathematical rulers" (the Quark Model) to predict how heavy these particles should be, the rulers said they should be much heavier. But the experiments showed they were surprisingly light—about 100 units of weight lighter than expected.

It was like ordering a standard brick and receiving a feather. This confused everyone. Was the brick broken? Was the ruler wrong? Or was the brick actually a hollow shell made of something else?

The Theory: The "Molecular" Hypothesis

The authors of this paper propose a solution: Maybe these aren't solid bricks at all. Maybe they are molecules.

Imagine a heavy truck (a charm quark) and a light car (a strange quark).

• The Old View: They are glued together into a single, solid block.
• The New View: The truck and the car are driving side-by-side, holding hands loosely. They are two separate vehicles moving as a unit.

In the language of physics, the authors suggest that the mysterious light particles are actually molecules formed by a D-meson (a charm quark + a light quark) and a K-meson (a strange quark + a light quark) sticking together.

The Tool: The "QCD Sum Rule" Calculator

How do you prove a particle is a molecule without seeing it directly? You can't just look under a microscope; these things are too small and exist for a split second.

The authors use a powerful mathematical technique called QCD Sum Rules.

• The Analogy: Imagine you are trying to figure out the weight of a hidden object inside a sealed box. You can't open the box, but you can shake it, listen to the sound it makes, and feel how it vibrates.
• The Math: The authors write down complex equations that describe how these particles "vibrate" based on the fundamental forces of nature (Quantum Chromodynamics). They calculate what the mass should be if the particle is a "molecule" of specific ingredients.

They didn't just look at the obvious ingredients; they accounted for the "background noise" of the vacuum (empty space) up to a very high level of complexity (dimension 12), ensuring their calculation was incredibly precise.

The Results: A Perfect Match

The authors calculated the masses of these theoretical molecules:

1. $DK$ Molecule: Predicted mass $\approx$ 2.322 GeV.
   • Real World Match: This matches the mysterious $D_{s0}(2317)$ perfectly.
2. $D^*K$ Molecule: Predicted mass $\approx$ 2.457 GeV.
   • Real World Match: This matches $D_{s1}(2460)$ perfectly.
3. $DK^*$ Molecule: Predicted mass $\approx$ 2.538 GeV.
   • Real World Match: This matches $D_{s1}(2536)$ perfectly.

The Verdict: The "molecule" theory works! The math predicts that if these particles are loose molecules, they will weigh exactly what the experiments measured. This solves the "too light" mystery.

The Crystal Ball: Predicting the Future

The authors didn't stop at explaining the past. They used their same "molecule calculator" to predict what happens with bottom quarks (which are even heavier than charm quarks).

They predicted three new "bottom-strange" molecules:

• $BK$ and $B^*K$: They predict these exist, but they are likely resonances (unstable vibrations that fall apart quickly), not stable bound states.
• $BK^*$: This is the exciting one. The math predicts a mass of 6.158 GeV.
  • Why it matters: This mass is lower than the energy needed to break it apart. This implies it is a stable bound state.
  • Real World Connection: The LHCb experiment recently saw a particle called $B_{sJ}(6158)$. The authors' prediction matches this observation almost perfectly, suggesting this new particle is indeed a $BK^*$ molecule.

Summary

Think of this paper as a master architect who finally figured out the blueprint for a weird, floating house.

1. The Problem: Some houses (particles) were too light to be solid bricks.
2. The Solution: They realized the houses were actually two smaller houses tied together with a rope (molecules).
3. The Proof: The math of the "rope-tied" houses matched the weight of the real houses exactly.
4. The Future: They used this new blueprint to predict where to find a new, stable house made of heavier materials, and it seems we've already found it!

This work helps us understand that the universe is full of complex, molecular structures made of the tiniest building blocks, not just simple, solid ones.

Technical Summary

1. Problem Statement

The paper addresses the long-standing "low-mass issue" in heavy-light meson spectroscopy. Specifically, the experimentally observed masses of the charm-strange states $D^*_s0(2317)$, $D_{s1}(2460)$, and $D_{s1}(2536)$ are approximately 100 MeV lower than predictions from conventional quark models and early lattice QCD simulations. This discrepancy has led to various theoretical interpretations, including compact tetraquarks, conventional $q\bar{q}$ states, mixed states, and hadronic molecules.

While the molecular picture (specifically $DK$, $D^*K$, and $DK^*$) has been proposed to explain these states, there is a lack of consensus among different theoretical approaches. Furthermore, the nature of their bottom-strange analogs ($BK$, $B^*K$, $BK^*$) remains largely unexplored experimentally, despite recent LHCb observations of structures near the $BK^*$ threshold. The authors aim to resolve the structural nature of these states and predict the properties of their bottom counterparts using a rigorous non-perturbative framework.

2. Methodology

The authors employ QCD Sum Rules (QCDSR), a powerful non-perturbative method connecting hadronic observables to QCD parameters via the Operator Product Expansion (OPE).

• Interpolating Currents: They construct color-singlet-color-singlet type currents to represent the hadronic molecular states with quantum numbers $J^P = 0^+$ and $1^+$.
  • $J_0(x)$ for scalar states ($DK/BK$).
  • $J_{1\mu}(x)$ and $J_{2\mu}(x)$ for axial-vector states ($D^*K/B^*K$ and $DK^*/BK^*$).
  • The currents are formed by products of light quark ($u, d, s$) and heavy quark ($c, b$) fields.
• Correlation Functions: Two-point correlation functions are calculated on two sides:
  1. Phenomenological Side: Saturated by the ground state hadrons and a continuum, parameterized by masses ($M$) and pole residues ($\lambda$).
  2. QCD Side: Calculated using Wick's theorem and quark propagators, expanded in terms of vacuum condensates.
• OPE Precision: A key methodological feature is the inclusion of vacuum condensates up to dimension 12. This high-order expansion is crucial for improving the convergence of the OPE and the reliability of the sum rules.
• Borel Transformation: The sum rules are derived by matching the two sides and applying a Borel transformation to suppress higher-energy continuum contributions and enhance the ground state signal.
• Stability Criteria: The analysis strictly adheres to three criteria for selecting the working window:
  1. Pole Dominance: The pole contribution must exceed 50% (typically found between 40–60% in the window).
  2. OPE Convergence: Higher-dimensional condensate contributions must be small (verified up to dimension 12).
  3. Energy Scale Formula: The optimal energy scale $\mu$ is determined using the formula $\mu = \sqrt{M^2 - M_{c/b}^2 - k M_s}$, where $M$ is the hadron mass and $k$ is the number of strange quarks.

3. Key Contributions

• High-Precision OPE: The calculation of spectral densities up to dimension 12 condensates represents a significant improvement over previous studies that often truncated at dimension 6 or 8.
• Systematic Bottom Analog Prediction: While many studies focus on charm, this work provides a systematic, parallel analysis of the bottom-strange sector ($BK, B^*K, BK^*$), offering predictions for states not yet fully characterized experimentally.
• Refined Energy Scale: The application of a specific energy scale formula tailored for tetraquark/molecular states ensures the renormalization group evolution of parameters is handled consistently.

4. Results

The authors extracted the masses ($M$) and pole residues ($\lambda$) for six molecular states. The results are summarized below:

State	$J^P$	Predicted Mass (GeV)	Experimental Candidate	Consistency
**$DK$** |$0^+$|$2.322^{+0.066}_{-0.072}$|$D^*_s0(2317)$	Excellent
$D^*K$	$1^+$	$2.457^{+0.064}_{-0.068}$	$D_{s1}(2460)$	Excellent
$DK^*$	$1^+$	$2.538^{+0.059}_{-0.062}$	$D_{s1}(2536)$	Excellent
**$BK$** |$0^+$|$5.970^{+0.061}_{-0.064}$	N/A	Above Threshold
$B^*K$	$1^+$	$6.050^{+0.062}_{-0.064}$	N/A	Above Threshold
$BK^*$	$1^+$	$6.158^{+0.061}_{-0.063}$	$B_{sJ}(6158)$	Consistent

• Charm Sector: The predicted masses for $DK$, $D^*K$, and $DK^*$ align remarkably well with the experimental masses of $D^*_s0(2317)$, $D_{s1}(2460)$, and $D_{s1}(2536)$. This strongly supports the interpretation of these particles as hadronic molecules rather than compact tetraquarks or conventional $q\bar{q}$ states.
• Bottom Sector:
  • The predicted masses for $BK$ and $B^*K$ are higher than their respective $BK$ and $B^*K$ thresholds. This implies they are likely resonance states rather than bound states.
  • The predicted mass for $BK^*$ is $6.158$ GeV, which is lower than the $BK^*$ threshold. This suggests $BK^*$ is a bound hadronic molecular state. This prediction is consistent with the LHCb observation of the $B_{sJ}(6158)$ structure, which decays via the $B^*K$ channel.

5. Significance

• Validation of Molecular Picture: The study provides robust theoretical evidence supporting the hadronic molecular interpretation for the low-mass charm-strange states, resolving a decade-old discrepancy with quark models.
• Guidance for Future Experiments: The prediction that $BK^*$ is a bound state with a mass of $\approx 6.158$ GeV offers a specific target for future experimental searches at facilities like LHCb and Belle II. Conversely, the prediction that $BK$ and $B^*K$ are resonances above threshold guides the search for broad structures rather than narrow peaks.
• Methodological Benchmark: By including dimension-12 condensates, the paper sets a new standard for precision in QCD sum rule analyses of exotic hadrons, demonstrating that higher-order terms are necessary for stable and reliable mass predictions in the heavy-flavor sector.
• Future Directions: The authors note that while masses are consistent, further work is needed to calculate decay widths and partial widths (e.g., radiative vs. strong decays) to fully confirm the molecular nature and distinguish it from other interpretations. The pole residues calculated here serve as essential inputs for such future three-point sum rule calculations.