Preliminary study of the $H$ dibaryon in $N_{\rm f}=2+1$ lattice QCD

This paper presents preliminary lattice QCD results for the $I=0$, $S=-2$ $H$ dibaryon using $N_{\rm f}=2+1$ flavors at unphysical quark masses, employing the distillation technique and multiple momentum frames to determine the interacting spectrum and scattering amplitude as a step toward establishing the particle's existence at physical quark masses.

The Gist

The Search for the "Super-Particle": A Layman's Guide to the H-Dibaryon

Imagine the atomic nucleus as a bustling city made of tiny building blocks called quarks. Usually, these blocks stick together in small, stable groups: three quarks make a proton or a neutron (the "families" of the city), and two quarks make a pion (the "glue" that holds them together).

But what if six quarks decided to hold hands and form a single, super-tight unit? That is the H-dibaryon.

This paper is a report from a team of scientists (the "BaSc" collaboration) who are trying to find out if this six-quark "super-family" actually exists, or if it's just a theoretical ghost.

Here is the story of their hunt, explained without the heavy math.

---

1. The Mystery Guest: The H-Dibaryon

In 1977, a physicist named Jaffe predicted that six quarks (two up, two down, and two strange) could bind together so tightly that they would be lighter than two separate particles (called Lambda baryons) floating apart.

• The Analogy: Imagine you have two heavy suitcases. Usually, carrying them separately is hard. But if you could magically fuse them into one super-suitcase that weighs less than the two separate ones, that would be the H-dibaryon.
• The Problem: Despite decades of searching in particle accelerators (like the ones at KEK in Japan or ALICE at CERN), no one has ever seen this "super-suitcase" in real life. It might be too heavy to exist, or it might be so unstable it falls apart instantly.

2. The Laboratory: A Digital Universe

Since we can't easily catch this particle in the real world, the scientists built a digital universe inside a supercomputer. This is called Lattice QCD.

• The Analogy: Think of the computer as a giant, 3D grid (like a massive Rubik's cube). The scientists simulate the laws of physics on this grid. They place "quarks" on the grid and watch how they interact.
• The Catch: To make the simulation run fast enough, they can't use the exact weight of real quarks. It's like trying to simulate a hurricane, but you have to use "heavy" air molecules to make the math work. In this study, the "pions" (the glue particles) are about 280 MeV heavy, whereas in our real universe, they are lighter (about 140 MeV). They are simulating a universe where the rules are slightly "off," hoping to see if the H-dibaryon appears there first.

3. The Detective Work: Listening for Echoes

How do you find a particle that might not even exist? You don't look for it directly; you listen for the echoes of its presence.

• The Method: The scientists create a "correlation matrix." Imagine shouting into a cave. If the cave is empty, the echo is simple. If there are hidden chambers (particles) inside, the echo gets complex and changes pitch.
• The Channels: They are looking at three specific "rooms" where particles might hang out:
  1. $\Lambda\Lambda$ (Two Lambda particles)
  2. $N\Xi$ (A Nucleon and a Xi particle)
  3. $\Sigma\Sigma$ (Two Sigma particles)
• The Distillation: To hear these echoes clearly, they use a technique called "distillation." It's like using a high-quality noise-canceling microphone to filter out the static and hear the faintest whisper of the six-quark state.

4. The Findings: A Glimmer of Hope (But Not Proof Yet)

The team ran their simulation and found the energy levels of these particle pairs.

• The Result: They found that the particles are interacting. The data shows a "spectrum" (a list of energy levels) that is slightly different from what you would expect if the particles were just floating past each other.
• The Interpretation: They used a mathematical tool called the Lüscher quantization condition.
  • Analogy: Imagine a guitar string. If you pluck it, it vibrates at a specific note. If you put a small weight on the string, the note changes. By measuring exactly how the "note" (energy level) changes in their digital cave, they can calculate how strongly the particles are pulling on each other.
• The Current Status: This is a preliminary study. They found that the particles are interacting, but they haven't yet confirmed if the H-dibaryon is a "bound state" (stuck together forever) or just a "resonance" (a temporary hug that breaks apart).

5. Why This Matters

The paper highlights a major headache in physics: Discretization Errors.

• The Analogy: Imagine trying to draw a perfect circle using only square pixels. If the pixels are huge, the circle looks blocky and wrong. If you make the pixels tiny, the circle looks smooth.
• The Issue: In previous studies, when scientists changed the "pixel size" (the lattice spacing) of their computer grid, the answer changed drastically. Sometimes the H-dibaryon looked heavy; sometimes light. This paper is a step toward solving that puzzle. They are working to ensure that when they finally lower the "pixel size" to match the real world, the answer remains consistent.

Summary: What's Next?

The scientists are essentially saying:

"We built a digital model of a universe with slightly heavy particles. We listened to the echoes of six quarks trying to hold hands. We heard something interesting! The particles are definitely interacting. But we need to refine our 'microphone' (better math), check our 'pixels' (smaller grid sizes), and lower the 'weight' of our particles to the real world to see if the H-dibaryon is truly there."

If they find it, it would rewrite our understanding of how matter holds together, potentially explaining how the cores of neutron stars (which are basically giant balls of neutrons) behave under extreme pressure. It's a hunt for the "Holy Grail" of nuclear physics, one digital simulation at a time.

Technical Summary

1. Problem Statement and Motivation

The existence of the $H$ dibaryon, a six-quark state ($uuddss$) with quantum numbers $I=0, S=-2, J^P=0^+$, remains an open question in nuclear physics.

• Theoretical Context: First proposed by Jaffe (1977) as a deeply bound state, experimental constraints (e.g., the Nagara event) suggest a much shallower binding energy ($B_H \lesssim 7$ MeV) or non-existence.
• Lattice QCD Discrepancies: Previous lattice calculations have yielded conflicting results regarding the binding energy, often dependent on the pion mass ($m_\pi$) and lattice spacing ($a$).
  • $SU(3)$-symmetric calculations often find significant binding.
  • Calculations near the physical point show weak interactions or vanishing binding.
  • Critical Issue: Significant discretization errors have been observed, where binding energies vary drastically with lattice spacing (e.g., a factor of 7 difference at $a \approx 0.1$ fm vs. the continuum limit).
• Goal: This study aims to provide preliminary results on the $H$ dibaryon using $N_f = 2+1$ dynamical quarks at a pion mass of $m_\pi \approx 280$ MeV, investigating the interacting spectrum across multiple momentum frames to determine if a bound state exists before extrapolating to physical quark masses.

2. Methodology

2.1 Lattice Setup

• Ensemble: The calculation uses a single CLS (Coordinated Lattice Simulations) ensemble (D251) with $N_f = 2+1$ flavors.
• Parameters:
  • Pion mass: $m_\pi \approx 280$ MeV (heavier than physical).
  • Lattice spacing: $a \approx 0.064$ fm.
  • Volume: $64^3 \times 128$ sites ($m_\pi L \approx 5.91$, suppressing finite-volume exponential corrections).
  • Action: $O(a^2)$ improved Lüscher-Weisz gauge action and $O(a)$ improved Wilson fermions with non-perturbatively tuned $c_{sw}$.
• Statistics: 209 configurations with all-to-all smeared propagators computed using sources on every time slice.

2.2 Operator Construction

• Basis: A large basis of bilocal baryon-baryon interpolating operators was constructed using the distillation technique.
• Channels: The study investigates the three relevant coupled channels: $\Lambda\Lambda$, $N\Xi$, and $\Sigma\Sigma$.
• Momentum Frames: Calculations were performed in multiple momentum frames ($\vec{P}^2 = (2\pi/L)^2 n$ with $n=0,1,2,3,4$) to access different energy levels and partial waves.
• Hexaquark Operators: Local hexaquark operators were excluded from the basis. Previous studies suggested they have small overlap with the flavor-singlet ground state, whereas bilocal baryon operators are more effective for constraining the finite-volume spectrum.

2.3 Spectrum Extraction

• Correlation Matrices: $N \times N$ correlation matrices $C_{ij}(t)$ were built and solved using the Generalized Eigenvalue Problem (GEVP).
• Fitting Strategy: To avoid "false plateaus," the authors employed multi-exponential fits directly to the diagonal elements of the rotated correlation matrix.
• Model Averaging: An information-theory-based criterion (Akaike Information Criterion, AIC) was used to select the optimal model (number of exponentials, priors, $t_{min}$). This approach averages over different models to reduce systematic bias from excited state contamination.
• Single Hadron Masses: Single baryon masses ($N, \Lambda, \Sigma, \Xi$) and meson masses were extracted to define thresholds.

2.4 Scattering Analysis (Lüscher Method)

• Formalism: The infinite-volume scattering amplitude is extracted via the Lüscher quantization condition, relating finite-volume energy levels to the $K$-matrix.
• Coupled Channels: The analysis considers the mixing of $S$ and $P$ waves in moving frames and the coupling between $\Lambda\Lambda$, $N\Xi$, and $\Sigma\Sigma$.
• Parameterization: A minimal 5-parameter model was used for the inverse $K$-matrix, separating $SU(3)$-symmetric terms and $SU(3)$-breaking terms. The flavor-singlet channel ($c_1$) was allowed to have linear energy dependence, while mixing terms were simplified based on chiral effective theory expectations.

3. Key Results

• Energy Spectrum: The authors extracted the finite-volume energy spectrum for the lowest-lying states in the rest frame and moving frames.
  • Figure 3 (in the paper) displays the extracted energy levels (blue points) compared to non-interacting levels (gray dashed lines) and physical thresholds.
• Preliminary Binding Analysis:
  • A preliminary analysis was performed on the six lowest energy levels in the rest frame.
  • The analysis neglected partial wave mixing and left-hand cuts (due to one-pion exchange) to test a simplified parameterization.
  • The model successfully described the lowest-lying flavor-symmetric spectrum, with the $SU(3)$-breaking term ($c_{1 \leftrightarrow 27}$) found to be crucial for a good fit quality.
• Limitations Identified:
  • The simple parameterization failed to describe higher-lying energy levels, likely due to the large $SU(3)$ breaking ($M_K^2 - M_\pi^2$) on this ensemble.
  • Left-Hand Cuts: The standard Lüscher formalism used does not yet incorporate the effects of left-hand cuts (specifically from one-pion exchange and two-pion exchange in the $\Sigma\Sigma$ channel), which are present in the kinematic region between thresholds. This is a known source of systematic error that requires extended formalism.

4. Significance and Future Outlook

• Contribution: This work represents a step toward resolving the $H$ dibaryon puzzle by utilizing a modern $N_f=2+1$ ensemble with improved actions and a rigorous operator basis (distillation) at a pion mass closer to the physical point than many previous $SU(3)$-symmetric studies.
• Systematic Control: The study highlights the critical importance of:
  1. Discretization Errors: Emphasizing that simulations at different lattice spacings are necessary to control $O(a)$ errors, given the large variations seen in previous works.
  2. Coupled Channels: The necessity of treating $\Lambda\Lambda$, $N\Xi$, and $\Sigma\Sigma$ as a coupled system with $S$-$P$ wave mixing.
  3. Left-Hand Cuts: The need to extend the Lüscher formalism to account for long-range interactions (pion exchange) which break down the standard formalism near thresholds.
• Next Steps: The collaboration plans to:
  • Perform simultaneous fits of single and two-baryon correlators to better capture correlations.
  • Implement model-averaging for two-baryon correlators to quantify systematic uncertainties.
  • Incorporate left-hand cut effects into the quantization condition.
  • Conduct simulations at different lattice spacings to extrapolate to the continuum limit.

In conclusion, while this is a preliminary study that does not yet definitively confirm or rule out the existence of the $H$ dibaryon at physical masses, it establishes a robust framework for future high-precision calculations and identifies the specific systematic effects (discretization, left-hand cuts) that must be addressed to reach a conclusive answer.