$T_{cc}$ pole trajectory

Using MILC's $N_f=2+1+1$ HISQ gauge ensembles with anisotropic Clover and $O(a)$-improved Wilson--Clover actions, this study investigates the $T_{cc}$ tetraquark spectrum across varying quark masses and operator bases, employing a modified Lüscher method to address non-analyticity near the Left Hand Cut.

The Gist

The Big Picture: Hunting for a "Ghost" Particle

Imagine the universe is a giant construction site filled with tiny, invisible LEGO bricks called quarks. Usually, these bricks snap together in pairs (like protons and neutrons) or triplets. But physicists have been wondering: Can four bricks snap together to form a stable, new shape?

Specifically, they are hunting for a "doubly charmed tetraquark" (called $T_{cc}$). Think of this as a rare, four-brick structure made of two heavy "Charm" bricks and two light "Up/Down" bricks.

For a long time, scientists weren't sure if this structure could actually hold together or if it would immediately fall apart. The paper you shared is a progress report from a team of scientists who are using a supercomputer to build this structure virtually and see if it sticks.

The Tools: The Digital Sandbox

To do this, they can't just use real bricks; they need a Lattice QCD simulation.

• The Lattice: Imagine a giant, 3D grid (like a giant wireframe cube) representing space and time.
• The Ensembles: They are using pre-made "worlds" (called MILC ensembles) generated by a massive collaboration. Think of these as different "weather conditions" or "terrain types" in their simulation.
• The Action: They use specific rules (called "actions") to tell the bricks how to move and interact.
  • For the light bricks, they use a standard rulebook (Wilson-Clover).
  • For the heavy "Charm" bricks, they use a special, faster rulebook (Anisotropic Clover) because heavy things move differently and need more precise timing to track.

The Strategy: Trying Different Recipes

The scientists are trying to figure out exactly how heavy the "Charm" bricks need to be for the structure to hold.

1. Varying the Weight: They didn't just test the real Charm mass. They tested a range of weights, from "light Charm" to "heavy Bottom." It's like testing a cake recipe by changing the amount of sugar to see when the batter finally sets.
2. The "Left Hand Cut" Problem: As they get closer to the real physical world (using lighter pions), the math gets weird and "bumpy" (this is the Left Hand Cut). To handle this, they used a special mathematical trick (a modified Lüscher method) to smooth out the bumps and get a clear reading.

The Ingredients: The Operator Basis

To find the particle, they need to look for it in the right way. In the past, some scientists only looked for the particle in one specific shape (like a molecule). This team decided to look at three different shapes at once to be sure they didn't miss it:

1. The Diquark-Antidiquark: Two heavy bricks glued together, and two light bricks glued together, then those two pairs stuck to each other.
2. The Molecular: Two separate pairs (like a $D$ meson and a $D^*$ meson) floating near each other.
3. The Scattering: The bricks just bouncing off each other.

They built a giant "correlation matrix" (a 5x5 grid of data) to see how these different shapes interact over time.

The Results: It's Holding Together!

Here is what they found in their preliminary data:

• The Setup: They simulated a box of space and watched how the energy levels of the bricks changed.
• The Discovery: When they used a specific setting (light quark mass $\kappa = 0.12566$), they found a "ground state" (the most stable version of the particle) that sits below the energy level where the bricks would naturally fly apart.
• The Analogy: Imagine two magnets. If you pull them apart, they snap back together. If the energy level is lower than the point where they separate, it means there is an attractive force holding them together.
• The Conclusion: The $T_{cc}$ particle seems to be a real, bound state, but it's a very "loose" hug. It's only bound by a tiny amount of energy (about 0.36 MeV), which is why it's so hard to find in real life. It's like a house of cards that is barely holding its shape.

What's Next?

This is a "status report." The team is currently:

1. Running more simulations on larger grids (243 x 64) to get clearer pictures.
2. Testing even lighter pion masses (getting closer to the real world).
3. Trying to map out the exact "trajectory" of the particle's pole (its stability) as they change the mass of the heavy quarks.

In a nutshell: The scientists built a digital universe, mixed different types of quark "dough," and found strong evidence that a rare, four-quark particle exists. It's a fragile, barely-holding-together structure, but it is there. They are now fine-tuning their recipe to see exactly how it behaves in the real world.

Technical Summary

1. Problem Statement

The paper addresses the theoretical uncertainty regarding the existence and binding energy of the doubly charmed tetraquark state, $T_{cc}$ ($cc\bar{u}\bar{d}$), with quantum numbers $I(J^P) = 0(1^+)$.

• Context: While lattice QCD strongly predicts a deeply bound doubly bottom tetraquark ($T_{bb}$), the status of $T_{cc}$ is less certain. Experimental and theoretical estimates suggest $T_{cc}$ lies very close to the $D^0 D^{*+}$ threshold (approx. 0.36 MeV below), making it a shallow bound state or a virtual state.
• Challenges:
  1. Operator Basis: Previous lattice studies often excluded diquark-antidiquark operators, relying solely on molecular/scattering operators. However, heavy quark symmetry suggests diquark operators are crucial for double-heavy systems.
  2. Left-Hand Cut (LHC): As the light quark mass approaches the physical pion mass, the scattering amplitude develops non-analyticities due to the Left-Hand Cut, complicating the extraction of pole positions using standard Lüscher formalism.
  3. Mass Variation: To understand the binding mechanism, it is necessary to vary the heavy quark mass (between charm and bottom) and the light quark mass (pion mass).

2. Methodology

The authors employ Lattice QCD simulations to investigate the spectrum of $T_{cc}$ across varying quark masses.

A. Simulation Setup

• Ensembles: Used MILC $N_f = 2+1+1$ Highly Improved Staggered Quark (HISQ) gauge configurations at two lattice spacings ($a \approx 0.15$ fm and $0.12$ fm).
• Quark Actions:
  • Light Quarks (u/d): $O(a)$-improved Wilson–Clover fermion action.
  • Strange Quark: Wilson–Clover action, tuned to reproduce the $\eta_s$ mass (688.5 MeV).
  • Charm Quark: Anisotropic Clover-improved Wilson action (Relativistic Heavy Quark or RHQ action) to handle discretization errors and improve temporal resolution.
• Parameter Tuning:
  • RHQ parameters ($m_0, c_P, \nu$) were tuned by matching spin-averaged $D_s$ masses, hyperfine splittings ($D^*_s - D_s$), and ensuring the relativistic dispersion relation ($c^2=1$) holds.
  • Light quark masses were varied to achieve pion masses ranging from $\sim 688$ MeV down to $400$ MeV.

B. Operator Basis

To capture the full dynamics of the tetraquark, the authors constructed a 5 $\times$ 5 correlation matrix using three distinct types of interpolating operators:

1. Diquark-Antidiquark ($D$): $\epsilon_{abc} [c^a T_C \gamma_k c^b] [\bar{u}^a C \gamma_5 \bar{d}^b T]$. Based on heavy quark-diquark symmetry, these are expected to be significant.
2. Molecular ($M_1, M_2$): Local meson-meson operators resembling $D D^*$ structures.
3. Scattering ($S$): Non-local operators representing $D$ and $D^*$ mesons with back-to-back momenta ($\vec{p}_1 + \vec{p}_2 = 0$).

C. Analysis Techniques

• GEVP: The Generalized Eigenvalue Problem (GEVP) was applied to the 5 $\times$ 5 correlator matrix to extract energy levels.
• Handling LHC: The authors utilized the standard Lüscher method when far from the Left-Hand Cut and a modified Lüscher method when close to it to handle non-analyticities.
• Computational Tricks: For the scattering operator $S$, the "one-end trick" with complex $Z(2) \times Z(2)$ random noise was employed to efficiently compute the necessary propagators.

3. Key Contributions

• Inclusion of Diquark Operators: This work explicitly includes diquark-antidiquark operators in the basis for $T_{cc}$ analysis, addressing a gap in previous lattice studies that often omitted them despite heavy quark symmetry arguments.
• Systematic Mass Variation: The study varies both heavy quark masses (interpolating between charm and bottom) and light quark masses (from heavy to near-physical pion masses) to map the pole trajectory of the $T_{cc}$.
• Robust Operator Basis: The construction of a 5 $\times$ 5 matrix combining diquark, molecular, and scattering operators provides a more comprehensive view of the state's nature (compact vs. molecular).
• LHC Handling: The application of modified Lüscher formalism near the Left-Hand Cut ensures more reliable extraction of scattering parameters in the physically relevant region.

4. Results

• Tuning Validation: The RHQ action parameters were successfully tuned, reproducing experimental $D_s$ and $D^*_s$ masses and hyperfine splittings within errors. The dispersion relations for $D$ and $D^*$ mesons confirmed $c^2 \approx 1$.
• Energy Spectrum (Preliminary):
  • At the heaviest pion mass ($m_\pi \approx 688$ MeV, $\kappa = 0.12566$), the GEVP analysis on the $16^3 \times 48$ lattice revealed three distinct energy levels below the inelastic $D^* D^*$ threshold.
  • Ground State: The lowest energy state lies below the $D^* D$ threshold, indicating an attractive interaction and a bound state.
  • Excited States: The first and second excited states appear between the $D^* D$ and $D^* D^*$ thresholds.
• Threshold Behavior: The analysis confirms that the simulation volume and momentum setup are sufficient to resolve the relevant energy levels for Lüscher analysis. The energy shifts from non-interacting thresholds clearly demonstrate the presence of interactions between $D$ and $D^*$ mesons.

5. Significance and Outlook

• Confirmation of Binding: The preliminary results at heavy pion masses confirm the existence of an attractive potential leading to a bound $T_{cc}$ state, consistent with the expectation that binding energy increases as the heavy quark mass increases.
• Future Directions: The authors plan to extend this analysis to lighter pion masses ($m_\pi \approx 400$ MeV) and heavier heavy-quark mass points (approaching the bottom sector).
• Scientific Impact: By mapping the pole trajectory, this work aims to clarify the nature of the $T_{cc}$ (whether it is a compact tetraquark or a loosely bound molecule) and resolve the theoretical uncertainties regarding its binding energy near the physical point. The inclusion of diquark operators is a critical step in distinguishing between these structural interpretations.

In summary, this paper presents a rigorous, multi-operator lattice QCD study of the $T_{cc}$ tetraquark, establishing a robust framework for extracting its pole trajectory by systematically varying quark masses and employing advanced analysis techniques to handle complex scattering dynamics.