$D_1$ and $D_2$ resonances in coupled-channel scattering amplitudes from lattice QCD

Using lattice QCD at a pion mass of approximately 391 MeV, this study computes isospin-1/2 charmed axial-vector scattering amplitudes to identify an $S$-wave $D_1$ bound state below the $D^*\pi$ threshold, a $D$-wave $D_1'$ resonance, and a narrow tensor state, while also observing significant $S$-wave interactions and a potential additional state in the coupled $D^*\pi$-$D^*\eta$-$D_s^*\bar{K}$ channels.

The Gist

Imagine the universe is a giant, cosmic dance floor. On this floor, particles like to pair up, spin, and interact. Sometimes, they form stable couples (like a ball and a stick glued together). Other times, they bump into each other, swirl around for a split second, and then fly apart. These fleeting, swirling moments are called resonances.

This paper is like a high-tech detective story where scientists are trying to map out the dance moves of specific particles called charmed mesons (particles containing a heavy "charm" quark). Specifically, they are looking at how these particles interact with lighter particles (pions, eta mesons, etc.) to see if they form new, hidden dance partners or if they just bounce off each other.

Here is the breakdown of their adventure, explained simply:

1. The Setting: A Tiny, Trapped Universe

To study these particles, you can't just look at them in the real world; they move too fast and decay too quickly. Instead, the scientists used Lattice QCD.

• The Analogy: Imagine trying to study how a fish swims, but you can't put it in the ocean. So, you put it in a tiny, transparent fishbowl with a grid painted on the walls.
• The Catch: The fish (the particle) behaves differently in a tiny bowl than in the open ocean. The walls force the fish to swim in specific patterns.
• The Trick: By measuring exactly how the fish hits the walls and bounces back (the "energy levels"), the scientists can use complex math to figure out how the fish would swim in the open ocean (the "scattering amplitude").

In this study, they used a "fishbowl" made of supercomputers, simulating a universe where the lightest particles (pions) are heavier than in our real world. This made the math easier to solve while keeping the essential physics intact.

2. The Characters: The D-Team

The main characters in this story are the D-mesons. Think of them as heavy dancers.

• D* (D-star): A heavy dancer with a bit of spin.
• D1 and D'1: These are the "excited" versions of the heavy dancer, spinning faster or in different ways.
• The Partners: The light dancers they interact with are pions ($\pi$), eta mesons ($\eta$), and kaons ($K$).

The scientists wanted to know: When a heavy D-star dancer meets a light pion, do they just bounce off? Do they hold hands for a moment? Or do they form a new, temporary dance troupe?

3. The Discovery: Finding the Hidden Dancers

By crunching the numbers from their "fishbowl" simulation, the scientists found evidence of four distinct dance patterns (states):

A. The "Ghost" Partner (Bound State)

• What they found: A state that exists just below the energy needed for the dancers to separate.
• The Analogy: Imagine two dancers who are so attracted to each other that they can't even break apart. They are stuck in a permanent embrace, but only just barely. In physics, this is a bound state.
• The Result: They found an axial-vector state (a specific type of spin) that is tightly bound to the $D^*\pi$ channel. It's like a dancer who refuses to let go of their partner.

B. The "Flash Mob" (Narrow Resonance)

• What they found: A state that appears right at the energy where the dancers usually separate, but it has a very specific, sharp signature.
• The Analogy: Think of a flash mob. A group of people suddenly starts dancing in perfect sync for a few seconds, then disappears. It's very distinct and short-lived.
• The Result: They found a narrow resonance (called $D'_1$) that is very picky. It only really dances with the "D-wave" (a specific way of spinning) and doesn't mix much with other moves. It's a precise, high-energy move.

C. The "Wild Card" (Broad Resonance)

• What they found: A state that is very messy and energetic, appearing at higher energies where many different dance partners are available.
• The Analogy: Imagine a chaotic mosh pit. Everyone is bumping into everyone else. It's hard to tell who is dancing with whom because the energy is so high and the interactions are so strong.
• The Result: They found a broad resonance (called $D_1$) that is very "fuzzy." It couples strongly to multiple channels ($D^*\pi$, $D^*\eta$, $D^*_s \bar{K}$). It's like a dancer who is good at everything but doesn't specialize in one move. This explains why experimentalists have had a hard time pinning down its exact mass in the real world—it's a "fuzzy" particle.

D. The "Twin" (Tensor State)

• What they found: A state with a different kind of spin (spin-2).
• The Analogy: If the previous dancers were spinning on one foot, this one is doing a complex pirouette.
• The Result: They found a narrow, sharp resonance in the "tensor" channel, acting like a well-defined, stable flash mob.

4. The Big Mystery: Why is the Heavy One Lighter?

One of the most surprising findings is about the "Wild Card" (the broad resonance).

• The Expectation: In the real world, experiments see a particle called $D_1(2430)$ that is quite heavy.
• The Discovery: The scientists found their version of this particle to be lighter than the experimental one, even though they used heavier "ingredients" (quarks) in their simulation.
• The Explanation: This suggests that the way we usually measure these particles (looking for a simple peak in the data) is misleading. The "peak" we see in experiments is actually a distorted shadow of the real particle. The real particle is actually lower in mass, but because it interacts so strongly with other particles, it gets "stretched out" and looks heavier and fuzzier. It's like looking at a reflection in a funhouse mirror; the image is distorted, but the object is real.

5. The "Six-Pack" Theory (SU(3) Symmetry)

The paper also touches on a deep theory called SU(3) Flavor Symmetry.

• The Analogy: Imagine if you could swap a "light" dancer for a "strange" dancer, and the dance moves would stay exactly the same. In reality, they aren't exactly the same, but they are very similar.
• The Result: The scientists found evidence for a "six-pack" of particles (a theoretical group called a sextet). They found a pole (a mathematical signature of a particle) that seems to belong to this group. This supports the idea that there is a hidden structure to how these particles organize themselves, beyond just simple pairs.

Summary

This paper is a triumph of "digital archaeology." By simulating a tiny, artificial universe on a supercomputer, the scientists were able to:

1. Map the dance floor: They calculated exactly how heavy charmed particles interact with light ones.
2. Find the ghosts: They identified a bound state that is stuck together.
3. Spot the flash mobs: They found sharp, narrow resonances that dance briefly.
4. Uncover the chaos: They found a broad, messy resonance that explains why experiments have been confused about the mass of certain particles.
5. Fix the mirror: They showed that the "fuzzy" nature of these particles makes them look heavier in experiments than they really are.

In short, they used math and supercomputers to decode the secret choreography of the subatomic world, revealing that the universe's dance floor is more complex, and more beautiful, than we previously thought.

Technical Summary

1. Problem Statement

The paper addresses the spectroscopy of excited charmed mesons, specifically the lightest open-charm states with valence quark content $c\bar{q}$ (where $q$ is a light or strange quark). While ground-state mesons ($D, D^*$) are well-described by the quark model, excited states like $D_1(2430)$, $D_1(2420)$, and $D^*_2(2460)$ present significant challenges:

• Broad vs. Narrow: Some states (e.g., $D_1(2430)$ and the scalar $D^*_0(2300)$) appear as broad enhancements above strong-decay thresholds, while others (e.g., $D_1(2420)$) are narrow.
• Quark Model Limitations: Simple quark-potential models struggle to explain the mass ordering and decay widths, particularly the discrepancy where the broad $D^*_0(2300)$ is measured heavier than the narrow $D^*_{s0}(2317)$, contrary to expectations based on light vs. strange quark mass differences.
• Theoretical Ambiguity: Unstable resonances are rigorously defined by pole singularities in the complex energy plane of scattering amplitudes, not just by Breit-Wigner peak positions. The interplay between quark-model configurations ($c\bar{q}$), meson-meson components, and potential tetraquark structures requires a non-perturbative treatment that accounts for coupled-channel dynamics and unitarity.

2. Methodology

The authors employ Lattice Quantum Chromodynamics (LQCD) to compute scattering amplitudes from first principles, avoiding model-dependent assumptions.

• Lattice Setup:

  • Ensembles: Three ensembles of gauge configurations with anisotropic lattices (finer temporal spacing) to improve energy resolution.
  • Quark Masses: Simulations are performed at a light-quark mass corresponding to a pion mass of $m_\pi \approx 391$ MeV. This unphysical mass is chosen to stabilize the $D^*$ meson (preventing $D^* \to D\pi$ decay) and push three-hadron thresholds ($D\pi\pi$) high enough to treat the system using only two-hadron channels.
  • Volumes: Three spatial volumes ($1.9, 2.4, 2.9$ fm) to extract finite-volume spectra.

• Operator Basis & Spectra:

  • A large basis of interpolating operators is used, including single-meson-like ($c\bar{q}$) and two-meson-like ($D^{(*)}\pi, D^{(*)}\eta, D^{(*)}_s\bar{K}$) operators.
  • Variational Analysis: The Generalized Eigenvalue Problem (GEVP) is solved on correlation matrices to extract a discrete spectrum of energy levels ($E_n$) in finite volume.

• Scattering Formalism:

  • Lüscher's Method: The finite-volume energy levels are mapped to infinite-volume scattering amplitudes using the quantization condition.
  • Coupled-Channel K-Matrix: The scattering amplitudes are parametrized using a K-matrix formalism to ensure unitarity and analyticity. The analysis covers:
    • $J^P = 1^+$ (Axial-vector): Coupled channels $D^*\pi$, $D^*\eta$, and $D^*_s\bar{K}$ in S-wave, with mixing to $D^*\pi$ D-wave.
    • $J^P = 2^+$ (Tensor): Coupled channels $D\pi$ and $D^*\pi$ in D-wave.
    • Backgrounds: Negative parity channels ($1^-, 2^-$) are treated as non-resonant backgrounds.
  • Global Fit: A global $\chi^2$ minimization is performed to fit the K-matrix parameters to 107 lattice energy levels across multiple irreducible representations (irreps) and moving frames.

• Pole Extraction:

  • The fitted amplitudes are analytically continued to the complex energy plane ($s = E_{cm}^2$).
  • Pole singularities are searched for on various Riemann sheets (defined by the signs of the imaginary parts of momenta in open channels) to identify bound states and resonances.
  • Systematic uncertainties are estimated by varying hadron masses, anisotropy, and the functional form of the K-matrix (polynomial orders).

3. Key Contributions

• First Coupled-Channel Axial-Vector Study: This is the first lattice QCD calculation of coupled $D^*\pi - D^*\eta - D^*_s\bar{K}$ scattering amplitudes in the $J^P=1^+$ channel, extending previous work that was limited to elastic $D^*\pi$.
• High-Energy Reach: The analysis extends to energies above the $D^*\eta$ and $D^*_s\bar{K}$ thresholds, probing the region where multiple open channels exist.
• SU(3) Flavor Symmetry Test: The authors test amplitudes constrained by $SU(3)_F$ flavor symmetry to investigate the existence of a "sextet" pole predicted by unitarized chiral perturbation theory.
• Rigorous Pole Analysis: Instead of relying on peak positions, the paper provides a rigorous extraction of pole locations and residues, distinguishing between bound states, narrow resonances, and broad structures.

4. Key Results

The analysis identifies four distinct pole singularities:

$J^P = 1^+$ (Axial-Vector) States:

1. State (I) - Bound State:

   • Location: A shallow bound state below the $D^*\pi$ threshold.
   • Mass: $2395.6 \pm 1.4$ MeV (in lattice units).
   • Coupling: Strongly coupled to $D^*\pi$ S-wave; decoupled from D-wave.
   • Interpretation: Likely the lattice counterpart to the broad experimental $D_1(2430)$. The fact that it is a bound state at $m_\pi=391$ MeV but expected to become a broad resonance at physical pion mass explains the experimental broadness.

2. State (II) - Narrow Resonance:

   • Location: Just above the $D^*\pi$ threshold, below the $D^*\eta$ threshold.
   • Mass: $2478.0 \pm 3.2$ MeV.
   • Width: Very narrow ($< 0.0013$ in lattice units, corresponding to a few MeV).
   • Coupling: Dominantly coupled to $D^*\pi$ D-wave ($^3D_1$).
   • Interpretation: Matches the experimental $D_1(2420)$. The mass is $\approx 55$ MeV higher than experiment, consistent with the heavier pion mass used.

3. State (III) - Broad Resonance:

   • Location: Deep in the complex plane, above the $D^*_s\bar{K}$ threshold.
   • Mass: $2737 \pm 79$ MeV.
   • Width: Large ($\approx 156$ MeV).
   • Coupling: Strongly coupled to $D^*\pi$ S-wave and heavier channels.
   • Interpretation: A new state not previously identified in this channel. It is consistent with the flavor-sextet pole predicted by $SU(3)$ chiral perturbation theory, supporting the "two-pole structure" hypothesis for the S-wave amplitude.

$J^P = 2^+$ (Tensor) State:

4. State (IV) - Narrow Resonance:
   • Location: Near the $D\pi$ and $D^*\pi$ thresholds.
   • Mass: $2524.6 \pm 2.2$ MeV.
   • Width: Very narrow ($\approx 1.6$ MeV).
   • Coupling: Strongly coupled to $D\pi$ ($^1D_2$) and $D^*\pi$ ($^3D_2$).
   • Interpretation: Corresponds to the experimental $D^*_2(2460)$. The mass difference ($\approx 64$ MeV) is consistent with the heavy-light quark mass scaling.

5. Significance and Implications

• Resolution of the $D_1(2430)$ Puzzle: The results support the hypothesis that the experimental $D_1(2430)$ is a broad resonance whose pole lies significantly lower in mass than the Breit-Wigner peak suggests. The lattice bound state at $m_\pi=391$ MeV is expected to evolve into this broad resonance as the pion mass is lowered to the physical point.
• Validation of Heavy Quark Symmetry (HQSS): The study confirms strong signatures of HQSS. The S-wave amplitudes for $D\pi$ (scalar) and $D^*\pi$ (axial-vector) show similar behaviors, and the narrow resonances in D-waves ($D_1(2420)$ and $D^*_2(2460)$) appear as degenerate doublets as predicted.
• Evidence for Exotic/Non-Quark-Model Dynamics: The existence of the broad State (III) and the two-pole structure in the S-wave supports theoretical predictions that these states are not simple $c\bar{q}$ excitations but involve significant meson-meson components (molecular-like or dynamically generated states) arising from $SU(3)$ flavor symmetry breaking.
• Methodological Benchmark: The work demonstrates the capability of modern lattice QCD to handle complex, multi-channel scattering problems with high precision, providing a rigorous framework for future studies of exotic hadrons.

In conclusion, this paper provides the first ab initio determination of coupled-channel axial-vector and tensor scattering amplitudes involving charmed mesons, identifying key resonant structures that align with experimental observations while offering a deeper theoretical understanding of their nature through pole analysis.