The structure of the lightest positive-parity charmed mesons from LQCD

Through $N_f=3+1$ lattice QCD simulations, this study resolves the debate on the nature of low-lying scalar and axial-vector charmed mesons by demonstrating that the observed attractive $[6]$ and repulsive $[\overline{15}]$ $SU(3)$ multiplet interactions support a hadronic molecule structure over a compact tetraquark configuration.

The Gist

The Big Mystery: What are these "Exotic" Particles?

Imagine the universe is built out of tiny Lego bricks. For a long time, physicists thought all matter was made of simple combinations: either two bricks stuck together (a quark and an antiquark) or three bricks (like a proton). These are the "standard" particles.

However, scientists discovered some strange, heavy particles (called charmed mesons) that didn't fit the standard Lego instructions. They were lighter than expected, behaving like they were made of four bricks instead of two.

The big question was: How are these four bricks stuck together?

There were two main theories:

1. The "Compact Tetraquark" Theory: Imagine four bricks glued tightly together into a single, dense, hard ball. The glue is strong, and the bricks are all hugging each other in the center.
2. The "Hadronic Molecule" Theory: Imagine two separate Lego structures (like a small car and a small truck) driving close to each other and holding hands. They aren't glued into one solid block; they are just two separate things orbiting each other, like the Earth and the Moon.

The Experiment: A Digital Simulation

Since we can't build these particles in a lab to see what's inside, the researchers used a supercomputer to simulate the universe. This is called Lattice QCD (Quantum Chromodynamics).

Think of this simulation as a giant, digital video game where they set the rules of physics. They created a virtual world where they could build these four-brick particles and see how they behaved. They specifically looked at two different "families" of these particles (called the [6] and [15] groups) to see which theory held up.

The "Taste Test": How the Theories Predicted Different Outcomes

The researchers needed a way to tell the two theories apart. They looked at how these particles interacted with each other.

• The "Molecule" Prediction: If these are just two particles holding hands (molecules), the physics says they should feel a gentle attraction in some groups (like magnets pulling together) but a repulsion in others (like trying to push two north poles of a magnet together).
• The "Compact Tetraquark" Prediction: If these are tight, glued balls, the theory predicted something very different. It said that in one specific group (the [15] group), the particles should be deeply bound (super sticky and heavy) in the "axial-vector" category, but not in the other.

The Analogy:
Imagine you have two types of dance partners.

• Theory A (Molecule): Says that if you try to dance with Partner X, you will be pulled together. If you try to dance with Partner Y, you will push each other away.
• Theory B (Tetraquark): Says that if you dance with Partner Y, you will be locked in a tight, unbreakable embrace.

The Results: The Digital Verdict

The researchers ran the simulation and measured the energy levels. Here is what they found:

1. The [6] Group: The particles felt a gentle pull toward each other (attractive). This matched both theories, so it wasn't the deciding factor.
2. The [15] Group (The Decider):
   • In the Scalar sector, the particles pushed each other away (repulsive).
   • In the Axial-Vector sector, the particles also pushed each other away (repulsive).

The Verdict:
The "Compact Tetraquark" theory predicted that the [15] group in the axial-vector sector should be a super-tight, glued-together ball. The simulation showed the opposite: they were pushing each other apart!

The "Hadronic Molecule" theory, however, predicted exactly this: that the [15] group should be repulsive in both cases. The simulation matched this perfectly.

The Conclusion

The study concludes that these exotic particles are not tight, compact balls of four quarks glued together. Instead, they are hadronic molecules—loose associations of two separate particles orbiting each other.

In simple terms:
The universe didn't glue four bricks into a hard rock. Instead, it let two separate Lego structures hold hands and dance around each other. The supercomputer simulation proved that the "glued ball" theory was wrong and the "holding hands" theory was right.

This is a huge step forward because it helps us understand the fundamental "glue" (the strong nuclear force) that holds the universe together, showing us that sometimes, particles prefer to stay close but separate, rather than merging into a single tight knot.

Technical Summary

1. Problem Statement and Motivation

The nature of low-lying positive-parity open-charm hadronic states, specifically the scalar $D^*_{s0}(2317)$ and axial-vector $D^*_{s1}(2460)$, has been a subject of long-standing debate. These states are significantly lighter than predictions from standard quark models ($q\bar{q}$ mesons), suggesting they are exotic four-quark configurations. Two primary competing hypotheses exist:

1. Compact Tetraquarks: Bound states of a diquark and an anti-diquark ($cq\bar{q}\bar{q}$).
2. Hadronic Molecules: Loosely bound states of two mesons (e.g., $D\pi$ or $D^*\pi$).

These two models make distinct, testable predictions regarding the energy levels of specific $SU(3)$ flavor multiplets (specifically the sextet $[6]$ and the 15-plet $[15]$) in a flavor-symmetric limit. Resolving which model correctly describes the spectrum is crucial for understanding the non-perturbative dynamics of QCD.

2. Methodology

The authors employed Lattice Quantum Chromodynamics (LQCD) to calculate the energy levels of these states directly from first principles.

• Simulation Setup:

  • Ensemble: $N_f = 3+1$ dynamical clover fermion ensembles with physical charm mass and three degenerate light quark flavors ($u, d, s$).
  • Parameters: Lattice coupling $\beta = 3.6$, lattice size $64^4$, and six iterations of stout smearing.
  • Tuning: The light quark mass was tuned to a pion mass range of $600 < M_\pi < 700$ MeV (specifically $M_\pi = 613 \pm 1$ MeV). The charm quark mass was tuned to reproduce the physical hyperfine splitting ratio of $J/\psi$ and $\eta_c$. The lattice spacing was determined to be $a \approx 0.27$ GeV$^{-1}$.
  • Statistics: A 9,960-trajectory ensemble was generated.

• Operator Construction and Contractions:

  • The study focused on the $SU(3)$ multiplets $[6]$ and $[15]$ in both scalar ($0^+$) and axial-vector ($1^+$) sectors.
  • The $[3]$ multiplet was neglected as it requires disconnected diagrams.
  • Correlators: The authors constructed four-quark interpolating operators. The contraction for the $[6]$ multiplet involves specific trace structures of light ($S$) and charm ($C$) propagators. The contraction for the $[15]$ is identical to the $[6]$ but includes a relative minus sign between terms, reflecting the different symmetry properties of the flavor wavefunctions.
  • Sources/Sinks: Smeared sources were used with both smeared and point sinks to improve signal-to-noise ratios.

• Analysis Strategy:

  • Correlators were fitted using a multi-exponential ansatz (including a term for the scattering threshold artifact where a pion propagates forward and a $D/D^*$ meson backward).
  • Fits were performed with $N=2, 3, 4$ states.
  • Model Averaging: The Akaike Information Criterion (AIC) was used to weight different fit models ($t_{min}$ and number of states) to determine the ground-state energy expectation values, minimizing model bias.
  • Error Analysis: A binned jackknife procedure was used to estimate statistical uncertainties.

3. Key Contributions

• First Direct Lattice Calculation: This work provides the first lattice calculation of the $SU(3)$ $[6]$ and $[15]$ multiplets for both scalar and axial-vector charmed mesons in a flavor-symmetric setting.
• Discrimination of Models: By calculating the mass shifts relative to the scattering thresholds ($D\pi$ for scalar, $D^*\pi$ for axial-vector), the authors directly test the competing predictions of the tetraquark vs. molecule models.
• Methodological Rigor: The use of AIC-based model averaging and the explicit handling of scattering artifacts in the correlator fits ensures robust extraction of ground-state energies.

4. Results

The authors calculated the mass shifts ($\Delta E$) relative to the relevant two-meson thresholds:

Sector	Multiplet	Mass Shift ($\Delta E$)	Interpretation
Scalar ($0^+$)	$[6]$	$-13 \pm 2$ MeV	Attractive (Bound state)
Scalar ($0^+$)	$[15]$	$+12 \pm 1$ MeV	Repulsive
Axial-Vector ($1^+$)	$[6]$	$-7 \pm 4$ MeV	Attractive (Bound state)
Axial-Vector ($1^+$)	$[15]$	$+11 \pm 3$ MeV	Repulsive

Key Finding: In both the scalar and axial-vector sectors, the $[6]$ multiplet exhibits attractive interactions (negative mass shift), while the $[15]$ multiplet exhibits repulsive interactions (positive mass shift).

5. Significance and Conclusion

The results decisively favor the hadronic molecule picture and rule out the compact tetraquark model for these states:

1. Consistency with Molecules: Unitarized Chiral Perturbation Theory (UChPT) predicts attractive interactions for the $[6]$ and repulsive interactions for the $[15]$ in both sectors. The lattice results match this pattern perfectly.
2. Inconsistency with Tetraquarks:
   • The tetraquark model predicts that the scalar $[15]$ should be repulsive (consistent with results).
   • However, the tetraquark model predicts that the axial-vector $[15]$ should be a deeply bound state (attractive), as it can be formed by combining a "good" diquark and a "bad" diquark with roughly equal energy costs.
   • The lattice result shows the axial-vector $[15]$ is repulsive ($\Delta E \approx +11$ MeV), directly contradicting the tetraquark prediction.

Conclusion: The study concludes that compact diquark-antidiquark tetraquarks are not a viable description for the lightest positive-parity open-charm states. Instead, the spectrum is wholly consistent with a hadronic molecule description, where these states arise from meson-meson interactions. This resolves a long-standing debate in hadron spectroscopy using first-principles QCD calculations.