Heavy hadron spectrum from 2+1+1 flavor MILC lattices

This paper presents a study of the mass spectra and mass differences of heavy hadrons containing bottom quarks using MILC's 2+1+1 flavor HISQ gauge ensembles, employing a mixed-action approach with NRQCD for bottom quarks, anisotropic Clover for charm, and O(a)-improved Wilson-Clover for light and strange quarks.

The Gist

Imagine the universe is built out of tiny, invisible LEGO bricks. In the world of particle physics, these bricks are called quarks. Most of the time, these bricks snap together to form protons and neutrons (the stuff inside atoms). But sometimes, they snap together in exotic ways to form "heavy hadrons"—particles that contain super-heavy bricks called bottom quarks and charm quarks.

This paper is a report from a team of scientists who are trying to build a perfect digital model of these heavy particles to understand how they work. Here is the story of what they did, explained simply.

1. The Challenge: Building with Heavy Bricks

Think of the bottom quark as a massive, heavy bowling ball, while the up and down quarks (the normal ones) are like ping-pong balls.

• The Problem: If you try to simulate a bowling ball moving on a computer grid that is too "chunky" (coarse), the simulation breaks. The ball is too heavy for the grid to handle accurately using standard physics rules.
• The Solution: The scientists used a special set of rules called NRQCD (Non-Relativistic Quantum Chromodynamics) for the heavy bowling balls. It's like saying, "We know this ball is heavy and moves slowly, so let's use a simplified rulebook just for it." For the lighter ping-pong balls, they used standard rules.

2. The Laboratory: A Digital Grid

The scientists didn't use a physical lab; they used a supercomputer to create a virtual 3D grid (a lattice).

• They used three different sizes of grids: a coarse one, a medium one, and a fine one.
• Think of this like taking a photo of an object with three different camera lenses: one blurry, one clear, and one ultra-sharp. By comparing the results from all three, they can figure out the "true" shape of the object without the blur of the camera messing things up.

3. Tuning the Engine: The "Dial" Process

Before they could build the heavy particles, they had to "tune" their simulation so it matched reality. This is like tuning a radio to get a clear signal.

• The Charm Quark (The Middleweight): They adjusted the settings until their digital "Charm Meson" weighed exactly the same as the real one measured in experiments (about 2.076 GeV).
• The Bottom Quark (The Heavyweight): They adjusted the settings for the bottom quark until their digital "Bottom Meson" matched the real-world weight (about 5.367 GeV).
• The Strange/Light Quarks: They tuned these to match the mass of a specific theoretical particle called the $\eta_s$.

Once these "dials" were set correctly, their computer simulation was ready to build new, never-before-seen particles.

4. The Construction: Building Heavy Baryons

With the settings tuned, they started building baryons (particles made of three quarks). They mixed and matched the heavy bottom quarks with charm, strange, and light quarks in every possible combination.

• Triple-Bottom: Three bottom quarks stuck together (like three bowling balls glued together).
• Double-Bottom: Two bottom quarks and one light quark.
• Single-Bottom: One bottom quark and two light quarks.

They created "interpolating operators," which are essentially mathematical blueprints or molds used to force the computer to form these specific shapes.

5. The Results: Weighing the Invisible

Once the particles were built in the simulation, the scientists measured their "effective mass" over time.

• The Plateau: Imagine watching a scale settle. At first, the number jumps around, but eventually, it flattens out into a steady line (a plateau). When the line flattens, that's the true weight of the particle.
• The Check: They tested their system by building a particle made of both a bottom and a charm quark ($B_c$ meson). It worked perfectly, proving their mixed-rule system was accurate.
• The Comparison: They compared their new weights for these heavy particles with previous studies and experimental data. The results matched up beautifully. Their digital models predicted the masses of these exotic particles with high precision.

6. Why Does This Matter?

You might ask, "Who cares about particles made of three bottom quarks?"

• Testing the Rules of the Universe: These heavy particles are like stress tests for the Standard Model of physics. If our calculations don't match reality, it means our understanding of the fundamental forces of nature is wrong.
• Future Experiments: As new particle accelerators come online, they will likely discover these heavy particles in the real world. This paper provides a "cheat sheet" for experimentalists, telling them exactly what mass to look for.
• Precision: By using three different grid sizes and carefully tuning the quark masses, the scientists are reducing the "fuzziness" in their calculations, moving us closer to a perfect understanding of how the universe is built.

Summary

In short, this paper is about a team of scientists using a supercomputer to build a digital LEGO set of the heaviest particles in the universe. They carefully calibrated their tools to match real-world weights, built complex structures made of heavy bottom quarks, and confirmed that their digital models are accurate. This work helps us understand the deep, non-perturbative (messy and complex) dynamics of the strong force that holds matter together.

Technical Summary

Here is a detailed technical summary of the paper "Heavy hadron spectrum from 2+1+1 flavor MILC lattices" by Sabiar Shaikh et al.

1. Problem Statement

The paper addresses the need for high-precision, first-principles calculations of the mass spectra and mass differences of heavy hadrons containing one or more bottom ($b$) quarks. While Lattice QCD has been successful in $B$-physics (decay constants, mixing parameters), the spectroscopy of heavy baryons (containing $b$, $c$, $s$, and light $u/d$ quarks) remains challenging due to:

• Systematic Uncertainties: Controlling discretization errors, particularly for heavy quarks on relatively coarse lattices.
• Action Compatibility: The difficulty of simultaneously simulating bottom quarks (non-relativistic), charm quarks (relativistic but heavy), and light quarks within a single framework.
• Experimental Context: Ongoing experimental progress in heavy-flavor physics requires precise theoretical benchmarks to interpret new data.

2. Methodology

The authors employ a mixed-action approach on $N_f = 2+1+1$ flavor Highly Improved Staggered Quark (HISQ) gauge ensembles generated by the MILC Collaboration. The study utilizes three lattice spacings (though results are currently presented for two) to control discretization effects.

A. Lattice Ensembles

• Source: MILC $N_f = 2+1+1$ HISQ gauge configurations.
• Lattice Spacings ($a$): Three ensembles are used:
  • Coarse: $a \approx 0.15$ fm ($16^3 \times 48$)
  • Medium: $a \approx 0.12$ fm ($24^3 \times 64$)
  • Fine: $a \approx 0.09$ fm ($32^3 \times 96$) (Production runs ongoing; results currently limited to the first two).
• Sea Quarks: Includes up, down, strange, and charm quarks.

B. Valence Quark Actions

Different lattice actions are employed for different quark flavors to optimize accuracy and computational efficiency:

1. Bottom Quarks ($b$): NRQCD (Non-Relativistic QCD) action.
   • Justification: The bottom quark is highly non-relativistic ($v^2 \sim 0.1$). NRQCD expands the Lagrangian in powers of velocity ($v$) up to $O(v^6)$.
   • Tuning: The bare mass ($m_b$) and hyperfine coefficient ($c_4$) are tuned to match the $B_s$ meson kinetic mass and hyperfine splitting ($\Delta M_{B_s}$) to PDG values.
2. Charm Quarks ($c$): Anisotropic Clover (Relativistic Heavy Quark - RHQ) action.
   • Justification: Introduces distinct spatial ($a_s$) and temporal ($a_t$) lattice spacings to improve temporal resolution for heavy quarks.
   • Tuning: Parameters ($m_0, c_P, \nu$) are tuned using the spin-averaged $D_s$ mass and the $D_s$-$D^*_s$ hyperfine splitting.
3. Strange and Light ($s, u, d$) Quarks: $O(a)$-improved Wilson–Clover action.
   • Tuning: The hopping parameter ($\kappa$) is tuned to reproduce the $\eta_s$ meson mass (688.5 MeV) and the physical pion mass. The clover coefficient ($c_{SW}$) is determined via tadpole improvement.

C. Operator Construction

• Interpolating Operators: Constructed for bottom mesons and baryons containing all combinations of $b, c, s, u, d$ quarks.
• Baryon Focus: Positive-parity states with spins $J=1/2$ and $J=3/2$.
• Structure: Operators respect color antisymmetry and specific spin structures (e.g., using charge-conjugation matrix $C$ and gamma matrices $\gamma_k, \gamma_5$).

3. Key Contributions

• Mixed-Action Framework: Successfully implements a hybrid setup (NRQCD + RHQ + Clover) on $2+1+1$ flavor HISQ ensembles, allowing for a consistent treatment of bottom, charm, and light quarks.
• Comprehensive Spectroscopy: Constructs interpolating operators for a wide range of heavy baryons, including triple-bottom ($\Omega_{bbb}$), double-bottom ($\Xi_{bb}, \Omega_{bb}$), and single-bottom ($\Lambda_b, \Xi_b, \Omega_b$) states, as well as mixed heavy-flavor states like $B_c$ mesons.
• Unitarity Check: Validates the four-component NRQCD formulation by comparing forward and backward propagation of $\eta_b$ correlators, confirming unitarity is preserved within statistical uncertainties.
• Tuning Validation: Demonstrates robust tuning of the bottom and charm sectors by reproducing experimental $B_s$ and $D_s$ masses and splittings across different lattice spacings.

4. Results

• Ground-State Extraction: Effective mass plots for $B_c$ mesons and various bottom baryons ($\Omega^*_{bbb}, \Omega^*_{cbb}, \tilde{\Omega}_{ccb}$) show clear, stable plateaus on the $16^3 \times 48$and$24^3 \times 64$ ensembles, allowing for reliable ground-state mass extraction.
• Mass Spectra:
  • The computed mass spectra for single, double, and triple-bottom baryons exhibit the expected mass hierarchies.
  • Results are in good agreement with previous lattice calculations (e.g., Brown et al., Burch, Mathur et al.) and experimental data where available.
• Consistency: The mass differences and splittings are consistent across the two lattice spacings used, indicating controlled discretization errors.
• Status: The study is presented as a status update; the finest lattice ($0.09$ fm) is currently being processed to enable continuum extrapolation.

5. Significance

• Precision Benchmarking: Provides a critical theoretical benchmark for current and future heavy-flavor experiments (e.g., LHCb, Belle II) by predicting masses for states that are difficult to observe or have large experimental uncertainties.
• Methodological Validation: Validates the mixed-action approach (NRQCD for $b$, RHQ for $c$, Clover for light) as a viable and precise method for heavy-hadron spectroscopy, bridging the gap between non-relativistic and relativistic regimes.
• Future Directions: This work lays the foundation for calculating hyperfine splittings, mass differences, and excited states, which are essential for understanding non-perturbative QCD dynamics in the heavy-quark sector. The inclusion of the finer lattice spacing will further reduce systematic errors, moving toward percent-level precision.

In summary, this paper represents a significant step toward precision lattice QCD calculations of the heavy-hadron spectrum, successfully combining advanced lattice actions to tackle the complex dynamics of multi-heavy-quark systems.