Binding energy of the $T_{bb}$ tetraquark from lattice QCD with relativistic and nonrelativistic heavy-quark actions

This paper presents a new lattice QCD determination of the $T_{bb}$ tetraquark binding energy using relativistic heavy-quark actions, yielding results of approximately $-79$ MeV that are consistent with prior NRQCD findings but exhibit reduced magnitude due to the exclusive use of symmetric correlation matrices with local four-quark operators.

The Gist

Imagine the universe is built out of tiny, invisible Lego bricks called quarks. Usually, these bricks snap together in pairs (like a proton and an antiproton) or triplets (like a proton or neutron). But physicists have long wondered: What if we could snap four of them together? Specifically, what if we took two heavy "bottom" quarks and two light "up/down" quarks?

This paper is about building a super-precise digital simulation to see if such a four-brick structure, called a Tbb tetraquark, can actually exist and stay stuck together without falling apart.

Here is the story of how the authors solved this puzzle, explained in everyday terms.

1. The Challenge: The "Heavy" Problem

The main characters in this story are the bottom quarks. They are incredibly heavy—about as heavy as a gold atom! In the world of quantum physics, heavy things are tricky to simulate on a computer.

Think of the computer simulation as a video game.

• The Old Way (NRQCD): Previously, scientists played this game using a "cheat code" or a simplified physics engine. They assumed the heavy quarks moved so slowly they could be treated as stationary statues. This made the game run faster, but it meant the physics wasn't 100% real. It was like playing a racing game where the cars don't actually accelerate; they just teleport.
• The New Way (RHQ): In this new study, the authors built a brand-new physics engine. They programmed the heavy quarks to move and behave exactly like real heavy objects, even though they are still heavy. This is like playing the racing game with full, realistic physics.

2. The Experiment: Building the Digital House

The team ran their simulation on seven different "worlds" (computer grids).

• Some worlds had "rough" terrain (coarse grids), and some had "smooth" terrain (fine grids).
• Some worlds had "heavy" pions (a type of particle acting like the air in the room), and some had "light" pions (closer to our real universe).

By running the simulation in all these different worlds, they could mathematically smooth out the rough edges and figure out what would happen in our real, physical universe.

3. The Big Discovery: Is it Stuck?

The goal was to measure the binding energy. Think of this as the "glue strength."

• If the glue is negative (in physics terms), it means the four bricks are stuck together tightly. They are a stable molecule.
• If the glue is zero or positive, the structure falls apart immediately.

The Result: The simulation showed a strong negative number. The Tbb tetraquark does exist and is stable! It's like finding a new, rare Lego castle that doesn't fall apart when you shake the table.

4. The Twist: Why the Numbers Changed

Here is the most interesting part of the paper. The authors compared their new, high-fidelity simulation (the "Real Physics" engine) with a re-analysis of an old simulation (the "Cheat Code" engine) using the exact same data.

• The Old Result: The old study suggested the glue was very strong (a binding energy of about -100 MeV).
• The New Result: The new study found the glue is still strong, but not as strong as previously thought (about -79 MeV).

Why the difference?
Imagine you are trying to hear a whisper in a noisy room.

• The Old Method: They used a microphone that picked up a lot of background noise (excited states). They thought the whisper was louder than it actually was because the noise made it sound deeper.
• The New Method: They used a high-tech noise-canceling microphone. They realized the "whisper" (the true energy of the particle) was actually a bit quieter than they thought.

The authors realized that the old method had a slight bias that made the particle look more tightly bound than it really is. By using better math and cleaner data, they got a more accurate, slightly "looser" measurement.

5. The Takeaway

This paper is a victory for precision.

1. Confirmation: It confirms that this exotic four-quark particle is real and stable.
2. Correction: It corrects previous estimates, showing that while the particle is stable, it isn't quite as tightly bound as we first thought.
3. Methodology: It proves that using the "Real Physics" engine (RHQ) gives results that match the "Cheat Code" engine (NRQCD) surprisingly well, giving us confidence that our understanding of the heavy quark world is solid.

In a nutshell: The scientists built a better microscope, looked at a strange new particle made of four quarks, and confirmed it exists. They also realized their old microscope was slightly zoomed in too much, making the particle look a bit more "glued together" than it actually is. Now, we have a clearer picture of this tiny, exotic piece of the universe.

Technical Summary

1. Problem Statement

The existence of a QCD-stable doubly-bottom tetraquark ($T_{bb}$) with quantum numbers $I(J^P) = 0(1^+)$ (composition $\bar{b}\bar{b}ud$) is theoretically well-established. However, previous lattice QCD calculations have yielded significant uncertainties in the binding energy ($\Delta E = m_{T_{bb}} - m_B - m_{B^*}$). Two primary sources of systematic error have plagued these calculations:

1. Heavy-quark discretization errors: Arising from the difficulty of simulating $b$-quarks on the lattice, often requiring non-relativistic approximations (NRQCD) or specific relativistic tuning.
2. Ground-state extraction: The difficulty of isolating the true ground-state energy from correlation functions, particularly due to excited-state contamination. Previous studies often used asymmetric correlation matrices or "wall" sources with point sinks, which can introduce a negative bias in the extracted energy (overestimating the binding magnitude) if the fit time separation is insufficient.

2. Methodology

The authors performed a controlled comparison using two different heavy-quark actions on the same set of gauge configurations to isolate discretization errors and refine ground-state extraction techniques.

• Gauge Configurations: The study utilized 2+1 flavor domain-wall fermion configurations generated by the RBC/UKQCD collaborations. Seven ensembles were used, covering lattice spacings from $0.114$ fm to $0.073$ fm and pion masses from $431$ MeV down to the physical point ($139$ MeV).
• Heavy-Quark Actions:
  • Relativistic Heavy Quark (RHQ): A new calculation using a nonperturbatively tuned three-parameter anisotropic-clover action (Fermilab method) for the $b$-quarks. This allows simulations at the physical $b$-quark mass with arbitrary lattice spacing.
  • NRQCD: A reanalysis of data from a previous study (Ref. [25]) using a tadpole- and Symanzik-improved NRQCD action for the $b$-quarks.
• Interpolating Operators and Correlation Matrices:
  • The authors focused exclusively on local four-quark operators ($3 \times 3$ correlation matrices).
  • They explicitly excluded "scattering operators" (involving $B$ and $B^*$ mesons) from the source and sink, relying on recent findings (Ref. [48]) that scattering operators are crucial for excited states but do not alter the ground-state energy extracted via the Generalized Eigenvalue Problem (GEVP) when local operators are used.
  • This approach contrasts with previous works that used asymmetric matrices (local source, scattering sink), which were suspected of causing negative energy biases.
• Analysis Strategy:
  • GEVP: The ground-state energy was extracted by solving the GEVP for the $3 \times 3$ correlation matrices.
  • Model Averaging: To handle systematic uncertainties from fit ranges ($t_{min}$ to $t_{max}$), the authors employed Bayesian model averaging (Ref. [99]) rather than selecting a single "best" fit range.
  • Extrapolations: Combined chiral (pion mass) and continuum (lattice spacing) extrapolations were performed for the RHQ data, while chiral-only extrapolations were used for the NRQCD data (as NRQCD does not have a continuum limit in the same sense).

3. Key Contributions

• First Physical-Point RHQ Determination: This is the first lattice QCD determination of the $T_{bb}$ binding energy at the physical pion mass and in the continuum limit using a relativistic heavy-quark action.
• Controlled Comparison of Actions: By reanalyzing NRQCD data with the exact same fit methodology and operator basis as the new RHQ data, the authors isolated the impact of the heavy-quark action. They found remarkable agreement between the two approaches, suggesting that heavy-quark discretization errors are under control in both methods when tuned correctly.
• Resolution of Ground-State Bias: The paper demonstrates that previous discrepancies in binding energy magnitudes (specifically the larger binding found in Ref. [25]) were likely due to the use of asymmetric correlation matrices and scattering operators at the sink, which introduced a negative bias. The current symmetric $3 \times 3$ analysis yields a smaller (less deeply bound) magnitude.

4. Results

The authors obtained the binding energy relative to the $B B^*$ threshold ($m_{T_{bb}} - m_B - m_{B^*}$) for both actions:

• RHQ Result:
  $$(m_{T_{bb}} - m_B - m_{B^*})_{RHQ} = (-79 \pm 23) \text{ MeV}$$
  Derived from combined chiral and continuum extrapolations of seven ensembles.

• NRQCD Result:
  $$(m_{T_{bb}} - m_B - m_{B^*})_{NRQCD} = (-74 \pm 17 \pm 10) \text{ MeV}$$
  Derived from chiral-only extrapolations of five ensembles, with an additional $10$ MeV systematic uncertainty estimate for discretization/NRQCD truncation errors.

• Comparison with Previous Work:
  The results are significantly smaller in magnitude than the original analysis in Ref. [25] (which reported $\approx -100$ to $-120$ MeV depending on the ensemble) and the earlier prediction by Francis et al. (2016). The authors attribute this reduction to the removal of excited-state contamination biases present in previous asymmetric matrix analyses.

5. Significance

• Confirmation of Stability: The results confirm the existence of a QCD-stable $T_{bb}$ tetraquark, as the binding energy is significantly negative (below the $B B^*$ threshold) in both relativistic and non-relativistic frameworks.
• Refined Prediction: The study provides the most precise lattice determination to date, suggesting a binding energy of approximately $-79$ MeV with a total uncertainty of roughly $25$ MeV. This indicates a deeply bound state, but with a smaller magnitude than the earliest lattice predictions.
• Methodological Impact: The work highlights the critical importance of using symmetric correlation matrices and local operators for ground-state extraction in tetraquark systems. It validates the use of the RHQ action for heavy-quark physics at the physical mass and demonstrates that NRQCD, when carefully reanalyzed, yields consistent results.
• Theoretical Implications: The smaller binding energy compared to early predictions may have implications for the production rates and decay modes of this state in future experiments (e.g., at LHCb or Belle II), as the phase space for decay into $B B^*$ is larger than previously thought, though the state remains stable against strong decay.