Quark-Mass Dependence of Light-Nuclei Masses from Lattice QCD and Trace-Anomaly Contributions to Nuclear Bindings

This paper presents lattice QCD calculations of light-nuclei masses across various quark masses to establish first-principles constraints on nuclear interactions and reveals that, while quark-mass contributions to binding energy are small and additive, the dominant contribution arises from the gluonic component via the QCD trace anomaly.

The Gist

Imagine the universe as a giant, intricate LEGO castle. For a long time, scientists knew that the bricks (protons and neutrons) were held together by a mysterious "glue" to form the walls of the castle (atomic nuclei). But they didn't fully understand what that glue was made of, or how changing the recipe for the bricks themselves would change the strength of the castle.

This paper is like a team of master builders using a super-powerful, digital simulation (called Lattice QCD) to rebuild the universe's smallest castles from scratch, right down to the fundamental ingredients.

Here is the story of what they found, explained simply:

1. The Ingredients: Quarks and Glue

Inside every proton or neutron, there are tiny particles called quarks. Think of quarks as the "flour" in a cake. But here's the twist: the weight of the cake (the proton) doesn't come mostly from the flour. It comes from the energy of the mixing and baking process (the gluons, or "glue").

The scientists wanted to know: If we change the "flour" (make the quarks heavier or lighter), how does that change the weight of the cake, and how does it change the strength of the glue holding the nuclear "castles" together?

2. The Experiment: Building Castles in a Digital Box

The researchers built a virtual box (a computer simulation) and created four specific nuclear structures:

• The Deuteron: A pair of protons/neutrons stuck together (like a small tent).
• The Dineutron: Two neutrons trying to stick together (like two magnets that usually repel).
• Helium-3 & Helium-4: Slightly larger, more complex structures (like a small house).

They ran these simulations with different "recipes." Sometimes they used heavy quarks (heavy flour), and sometimes they used light quarks (light flour), eventually reaching the physical point—the exact recipe nature uses in our real universe.

3. The Big Discovery: The "Glue" Does the Heavy Lifting

The most surprising finding is about why these nuclei hold together.

• The Old Guess: Scientists thought that changing the quark mass (the flour) would be the main reason the binding energy (the glue strength) changed.
• The Reality: The paper found that the quark mass contribution is actually tiny. It's like adding a pinch of salt to a soup; it changes the flavor slightly, but it doesn't change the volume.
• The Real Hero: The gluonic component (the energy of the glue itself) is the dominant force. When the nuclei get bigger (from 2 particles to 4), the "glue" gets stronger and more complex. The binding energy grows because of the collective dance of the gluons, not just because of the quarks.

The Analogy: Imagine a group of people holding hands to form a circle.

• The quarks are the people's hands.
• The gluons are the tension in their arms and the energy they spend holding on.
• The paper shows that if you change the size of the people's hands (quark mass), the circle doesn't change much. But if you change how hard they squeeze or how many people join in (gluonic dynamics), the circle becomes incredibly strong.

4. The "Unbound" Neighbor

One of the most fascinating results concerns the dineutron (two neutrons).

• In the simulations with heavy quarks (heavy flour), the two neutrons stuck together tightly. They were a happy couple.
• But as they tuned the simulation to the real world (physical quark mass), the couple broke up. The dineutron became unbound. It couldn't hold itself together.
• Meanwhile, the deuteron (proton + neutron) stayed together, but just barely (a "shallow" bond).

This tells us that the universe is in a very delicate balance. If the quarks were just a little bit heavier, the dineutron would exist, and the chemistry of the universe (and life itself) would be completely different.

5. Why This Matters

This paper is a "first-principles" proof. It doesn't rely on guessing or approximations; it calculates the laws of physics directly from the bottom up.

• It solves a mystery: It confirms that the mass of the universe (and the stability of matter) comes mostly from the trace anomaly—a fancy term for the energy generated by the gluon field, not the mass of the particles themselves.
• It sets the rules: By showing how nuclear binding changes with quark mass, they have provided a strict rulebook for other theories (like Effective Field Theories) to follow.
• It explains the "Why": It tells us that the reason atoms hold together is primarily a gluonic phenomenon. The quarks are just the passengers; the gluons are the engine.

In a Nutshell

The universe is held together not by the weight of its smallest parts, but by the energy of the invisible glue that binds them. This paper proves that if you tweak the ingredients slightly, the glue behaves differently, but the glue itself is the true hero of nuclear stability. Without this specific "glue" behavior, the stars, planets, and us would simply fall apart.

Technical Summary

1. Problem Statement

The origin of visible mass in the universe is rooted in Quantum Chromodynamics (QCD). While hadron masses are predominantly generated by gluonic interactions (via the trace anomaly) rather than bare quark masses, the specific mechanism by which nuclear binding energy emerges from QCD dynamics remains a central open question.

• The Gap: Previous lattice QCD studies of multi-nucleon systems were largely performed at unphysically heavy pion masses ($m_\pi > 400$ MeV), where binding patterns often differ qualitatively from nature. There was a lack of first-principles calculations at physical quark masses to determine how nuclear binding depends on the light-quark mass.
• The Goal: To compute the masses of light nuclei ($^2$H, $nn$, $^3$He, $^4$He) using physical sea quarks and varying valence quark masses, determine the nuclear sigma terms (quantifying the quark-mass contribution), and decompose the nuclear binding energy into quark-mass and gluonic components using the QCD trace anomaly.

2. Methodology

A. Lattice Setup

• Ensemble: The authors utilized a $2+1$-flavor Highly Improved Staggered Quark (HISQ) gauge ensemble generated by the HotQCD Collaboration.
• Parameters:
  • Sea quarks: Physical strange mass and light sea quarks corresponding to $m_\pi \approx 140$ MeV.
  • Lattice spacing: $a \approx 0.076$ fm on a $64^4$ lattice (physical volume $\approx (5 \text{ fm})^3$).
  • Valence quarks: Tree-level tadpole-improved Wilson–clover action with HYP-smeared gauge links.
  • Valence pion masses: Three sets corresponding to $m_\pi \approx 140$ MeV (physical), $300$ MeV, and $700$ MeV.
• Operators: Interpolating operators were constructed for single nucleons and light nuclei ($^2$H, $nn$, $^3$He, $^4$He) with well-defined spin and isospin quantum numbers. Both symmetric smeared-smeared (SS) and asymmetric smeared-point (SP) source-sink setups were employed to control systematic errors.

B. Energy Extraction

• Techniques: The finite-volume energy levels ($E_A$) were extracted using two complementary methods to ensure robustness:
  1. Multi-state fits: Standard exponential fits to the two-point correlation functions.
  2. Transfer Matrix Generalized Eigenvalue Problem (TGEVP): A recently developed method to handle excited-state contamination and spurious eigenvalues.
• Spurious Mode Removal: Spurious eigenvalues in the TGEVP were identified and removed using adaptive Kernel Density Estimation (KDE) on bootstrapped eigenvalue distributions.
• Binding Energy: The energy splitting $\Delta E_A = E_A - A \cdot m_N$ was calculated, where $m_N$ is the nucleon mass.

C. Theoretical Frameworks

• Feynman-Hellmann Theorem: Used to extract the nuclear sigma term ratio $\sigma_A / \sigma_N = \partial E_A / \partial m_N$ from the slope of the energy vs. nucleon mass curve.
• QCD Trace Anomaly: The energy-momentum tensor trace relation was used to decompose the nuclear binding energy ($\Delta E_A$) into:
  • Quark-mass contribution ($\Delta E^\sigma_A$): Arising from explicit chiral symmetry breaking.
  • Gluonic contribution ($\Delta E^{F^2}_A$): Arising from the trace anomaly (gluon dynamics).
  • Decomposition performed at the renormalization scale $\mu = 2$ GeV.

3. Key Results

A. Nuclear Binding Patterns

• Physical Point ($m_\pi \approx 140$ MeV):
  • Deuteron ($^2$H): The lowest energy level is consistent with a shallow bound state within uncertainties.
  • Dineutron ($nn$): The energy shift changes sign as $m_\pi$ approaches the physical value, indicating an unbound system (consistent with nature).
  • Helium Isotopes: At the physical point, correlation functions for $^3$He and $^4$He were noise-dominated, preventing precise extraction, but results at $m_\pi \approx 300$ MeV showed bound states consistent with previous lattice studies.
• Heavy Quark Masses ($m_\pi \approx 700-800$ MeV): Both deuteron and dineutron channels exhibit deeply bound states, confirming the trend that nuclear binding strengthens significantly at heavier quark masses.
• Operator Sensitivity: At heavy masses, SS setups yielded smaller (less negative) energy shifts than SP setups, highlighting sensitivity to operator structure when scattering states are dense. At lighter masses, SS and SP results were consistent.

B. Nuclear Sigma Terms

• The ratio of nuclear to nucleon sigma terms ($\sigma_A / \sigma_N$) was determined via linear fits of $E_A(m_N)$.
• Findings: The ratios are approximately equal to the mass number $A$ (e.g., $\approx 2$ for $^2$H, $\approx 3$ for $^3$He, $\approx 4$ for $^4$He).
• Implication: This suggests that the quark-mass contribution to nuclear binding is approximately additive with respect to the number of nucleons, with only small non-additive (many-body) corrections within current precision.

C. Trace-Anomaly Decomposition

• Using the calculated sigma terms and physical nuclear masses, the binding energy was decomposed at $\mu = 2$ GeV.
• Quark-Mass Component: The contribution from the quark-mass term ($\Delta E^\sigma_A$) to the binding energy per nucleon is small and roughly constant across different nuclei.
• Gluonic Component: The gluonic contribution ($\Delta E^{F^2}_A$) provides the dominant share of the binding energy.
• Mass Number Dependence: The growth of binding energy with increasing mass number ($A$) is driven primarily by the increase in the gluonic component, suggesting that nuclear binding is governed by gluonic many-body dynamics associated with the QCD trace anomaly.

4. Significance and Impact

1. First-Principles Constraints: This work provides the first lattice QCD constraints on the quark-mass dependence of two- and three-nucleon interactions using physical sea quarks. The results align with expectations from low-energy effective field theories (EFT) and phenomenological potentials (like AV18).
2. Origin of Nuclear Mass: The study quantitatively demonstrates that while the mass of individual nucleons is largely gluonic, the binding energy holding nuclei together is also predominantly a gluonic phenomenon. The quark-mass contribution is a subleading correction.
3. Validation of EFT: The observed mass dependence validates the use of pionless EFT and chiral EFT frameworks in describing the quark-mass dependence of nuclear forces, providing a bridge between lattice QCD and nuclear phenomenology.
4. Methodological Advancement: The successful application of TGEVP with adaptive KDE for spurious mode removal in multi-nucleon systems at physical pion masses sets a new standard for extracting energy levels in noisy, multi-hadron channels.

Conclusion

The paper establishes a direct link between nuclear binding and the QCD trace anomaly. It concludes that the stability of light nuclei and the increase in binding energy with mass number are driven primarily by gluonic dynamics, while the explicit chiral symmetry breaking (quark-mass effects) plays a minor, approximately additive role. This resolves a long-standing question regarding the QCD origin of nuclear binding and provides a rigorous benchmark for future nuclear physics calculations.