Lattice determination of the QCD low-energy constant $\ell_{\scriptscriptstyle{7}}$

This paper presents a non-perturbative lattice QCD determination of the scheme- and scale-independent low-energy constant $\ell_7$ using $N_f=2+1$ staggered fermion simulations, yielding a final result of $2.79(61) \times 10^{-3}$ that significantly improves upon previous estimates.

The Gist

Imagine the universe is built from a giant, complex Lego set. The most fundamental pieces are tiny particles called quarks, which stick together to form protons and neutrons (the building blocks of atoms). The "glue" that holds these quarks together is a force called Quantum Chromodynamics (QCD).

However, QCD is incredibly complicated. Trying to calculate exactly how these particles interact at every energy level is like trying to predict the exact weather pattern of every single raindrop in a hurricane. It's too messy and computationally impossible.

The Solution: A "Cheat Sheet" for the Low Energies

Physicists use a tool called Effective Field Theory. Think of this as a "cheat sheet" or a simplified map. Instead of tracking every single raindrop, the map just tells you the general direction of the wind and the average rainfall. This map works perfectly for the "low-energy" world (where things move slowly and aren't exploding), which is where our everyday matter exists.

This map has a few missing numbers on it, called Low-Energy Constants (LECs). These numbers are like the specific ingredients in a secret sauce recipe. If you don't know the exact amount of salt (one of these constants), your sauce (the theory) won't taste right, and your predictions about the universe will be off.

The Mystery Ingredient: $\ell_7$

One of these missing ingredients is a number called $\ell_7$ (ell-seven).

• What does it do? It controls a very subtle difference between two types of particles: the charged pion (a particle with a positive or negative electric charge) and the neutral pion (a particle with no charge).
• The Analogy: Imagine you have two identical twins. One wears a heavy winter coat (charged), and the other wears a light t-shirt (neutral). In a perfect world, they would weigh exactly the same. But in our universe, the "coat" makes the charged twin slightly heavier. The number $\ell_7$ tells us exactly how much heavier the charged twin is due to the "strong force" glue, ignoring electricity.
• Why do we care? This tiny difference is crucial for understanding Axions. Axions are hypothetical particles that scientists think might be the "Dark Matter" holding galaxies together. To figure out how heavy an axion is, or how it interacts with other particles, we need to know the exact value of $\ell_7$. If $\ell_7$ is wrong, our predictions for Dark Matter could be completely off.

The Problem: The Recipe Was Guesswork

For a long time, scientists had to guess the value of $\ell_7$ based on rough estimates or indirect clues. It was like trying to bake a cake by guessing how much sugar to add because you lost the recipe. The guesses varied wildly, and the uncertainty was huge. This made it hard to make precise predictions about the axion and Dark Matter.

The Experiment: Simulating the Universe on a Computer

In this paper, a team of physicists decided to stop guessing and start measuring. They used a massive supercomputer to simulate the universe from scratch. This is called Lattice QCD.

• The Grid: Imagine a 3D grid (like a giant Rubik's cube) representing space and time. They placed the quarks on the corners of this grid.
• The Simulation: They ran millions of calculations to see how the quarks interacted on this grid.
• The Challenge: The grid has a "pixel size" (called lattice spacing). If the pixels are too big, the image is blurry. If the pixels are too small, the computer takes forever to calculate. They had to run the simulation on grids of different sizes and with different "weights" for the quarks to see how the result changed.

The Breakthrough: Finding the Signal in the Noise

The team used a clever method to isolate $\ell_7$. They looked at the tiny mass difference between the charged and neutral pions.

• The Analogy: Imagine trying to hear a whisper in a rock concert. The "whisper" is the tiny mass difference caused by $\ell_7$. The "rock concert" is the massive amount of background noise from the computer simulation.
• The Innovation: Previous attempts using a specific type of computer code (Wilson fermions) worked but were slow. This team used a different, faster code called Staggered Fermions. However, this code had a quirk: it sometimes made the "whisper" disappear into the noise.
• The Fix: They discovered that by looking at the pions from a different "angle" (using a specific mathematical structure called the axial-vector channel), the whisper became loud and clear. It was like putting on noise-canceling headphones that only let the specific frequency of the whisper through.

The Result: A Precise Measurement

After running simulations on 12 different "universes" (computer grids) with varying sizes and quark weights, they extrapolated the results to the "perfect" universe (where the grid pixels are infinitely small and the quarks have their real physical weights).

Their Final Answer:
They determined the value of $\ell_7$ with much higher precision than ever before.

• Old Guess: "It's somewhere between 2 and 10."
• New Result: "It is 2.79, and we are very confident it's between 2.6 and 3.0."

Why This Matters

This is a huge step forward for physics.

1. Better Dark Matter Hunting: With this precise number, scientists can now calculate the properties of the axion (the leading Dark Matter candidate) much more accurately. This helps experimentalists know exactly what to look for in their detectors.
2. Proof of Concept: They proved that you can use this specific, fast computer code (Staggered Fermions) to solve these difficult problems, opening the door for faster and more precise calculations in the future.

In summary: The team built a digital microscope to zoom in on the fundamental forces of nature. They found a tiny, previously blurry number ($\ell_7$) that acts as a key to unlocking the secrets of Dark Matter, turning a rough guess into a precise measurement.

Technical Summary

1. Problem Statement

The paper addresses the precise determination of the low-energy constant (LEC) $\ell_7$ within the framework of $N_f=2$ Chiral Perturbation Theory ($\chi$PT) at Next-to-Leading Order (NLO).

• Physical Role: $\ell_7$ is the unique scheme- and scale-independent LEC that parametrizes strong isospin-breaking (IB) effects. It governs the mass splitting between charged and neutral pions ($M_{\pi^+}^2 - M_{\pi^0}^2$) arising solely from the up/down quark mass difference ($m_d - m_u$), distinct from electromagnetic effects.
• Importance: Accurate knowledge of $\ell_7$ is critical for:
  • Predicting axion properties (mass and self-coupling), where current uncertainties in $\ell_7$ dominate the error budget for axion physics.
  • Constraining Beyond Standard Model (BSM) phenomenology.
• Current Status: Previous determinations were either phenomenological estimates with large systematic errors (e.g., $\sim 7(4) \times 10^{-3}$) or indirect lattice calculations with large uncertainties (e.g., $6.5(3.8) \times 10^{-3}$). A direct, high-precision lattice calculation with controlled extrapolations was lacking.

2. Methodology

The authors perform a first-principles lattice QCD calculation using staggered fermions to generalize a method previously proposed for Wilson fermions.

A. Simulation Setup

• Fermions: $N_f = 2 + 1$ flavors of rooted stout-smeared staggered fermions.
• Gauge Action: Tree-level Symanzik-improved action.
• Ensembles: 12 distinct gauge ensembles generated using the Rational Hybrid Monte Carlo (RHMC) algorithm.
  • Parameters: 3 different Lines of Constant Physics (LCPs) corresponding to light-to-strange quark mass ratios $R = m_\ell/m_s^{phys} \approx 4, 6, 9$ times the physical ratio.
  • Lattice Spacing: 4 different lattice spacings ($a \approx 0.075 - 0.15$ fm) to enable continuum extrapolation.
  • Volume: Fixed physical volume $L \approx 3$ fm, ensuring $M_\pi L \geq 4$ to suppress finite-size effects.

B. The Mass-Splitting Method (RM123 Approach)

The calculation relies on the RM123 method, which treats isospin breaking as a perturbation around the isosymmetric point ($m_u = m_d$).

1. Theoretical Basis: The pion mass splitting is expanded in powers of the mass difference $\Delta m = (m_d - m_u)/2$. At leading order in $\Delta m$, the splitting is proportional to $\ell_7$.
2. Derivative Extraction: Instead of simulating non-degenerate masses directly, the method computes the second derivative of the pion mass difference with respect to $\Delta m$ at $\Delta m = 0$.
   $$\ell_7 \propto \frac{1}{(\Delta m)^2} (M_{\pi^+}^2 - M_{\pi^0}^2) \bigg|_{\Delta m \to 0}$$
3. Correlation Functions: This derivative is extracted from the difference between charged and neutral pion correlation functions expanded to $O((\Delta m)^2)$. This involves calculating connected and disconnected 4-point correlation functions.
4. Interpolating Operators: A critical technical innovation was the choice of interpolating operators. The standard pseudoscalar channel yielded a disconnected signal dominated by noise. The authors successfully utilized the axial-vector channel ($\xi_A$), where the disconnected signal is non-zero and statistically robust, allowing for a reliable extraction of $\ell_7$.

C. Extrapolation Strategy

• Continuum Limit: Data is extrapolated to $a \to 0$ using linear fits in $a^2$ and Symanzik Effective Theory (SyMFT) inspired fits involving the taste-symmetry breaking parameter $\Delta = M_{\pi, A}^2 - M_{\pi, P}^2$.
• Chiral Limit: Results at finite quark masses are extrapolated to the $SU(2)$ chiral limit ($m_\ell \to 0$) using NNLO $\chi$PT formulas.
  • The authors derived new NNLO expressions for the pion mass splitting, identifying the specific chiral logarithms and linear terms dependent on $\ell_7$ and other NNLO LECs.
  • Correction factors were applied to remove leading chiral logs before performing the linear extrapolation in $R$.

3. Key Contributions

1. First Direct Staggered Determination: This is the first lattice calculation of $\ell_7$ using staggered fermions, providing a crucial cross-check of universality against previous Wilson fermion results.
2. Methodological Generalization: Successfully adapted the RM123 mass-splitting method to staggered fermions, overcoming the specific challenge of noisy disconnected diagrams by switching to the axial-vector channel.
3. Controlled Systematics: The study features the first calculation of $\ell_7$ with fully controlled extrapolations to the continuum limit, infinite volume, and chiral limit.
4. NNLO $\chi$PT Calculation: The paper provides the first complete derivation of the NNLO contributions to the strong pion mass splitting, including isospin-breaking operators, which is essential for the chiral extrapolation.
5. Auxiliary Results: The study also provides updated determinations for $\bar{\ell}_3$ and $\bar{\ell}_4$ (other NLO LECs) and confirms the restoration of taste symmetry in the continuum limit.

4. Results

The authors report the final value for the low-energy constant $\ell_7$ in the chiral limit:

$$\ell_7 \times 10^3 = 2.79(58)_{\text{stat}}(19)_{\text{syst}} = 2.79(61)_{\text{tot}}$$

• Precision: The total uncertainty is approximately 2.2%, a significant improvement over previous lattice results (which had uncertainties of $\sim 60\%$ or more).
• Comparison:
  • Compatible with the indirect RBC-UKQCD result ($6.5(3.8)$) but with much higher precision.
  • Compatible with the previous direct Wilson fermion result ($2.5(1.4)$) but with a factor of $\sim 5$ reduction in uncertainty.
  • Compatible with phenomenological estimates but provides a rigorous first-principles constraint.
• Error Budget: The error is dominated by statistical uncertainties (approx. 85%), primarily arising from the chiral extrapolation. The signal-to-noise ratio deteriorates rapidly as the quark mass approaches the physical point, limiting further precision without significantly increased computational resources.

5. Significance

• Axion Physics: The result significantly reduces the theoretical uncertainty in axion mass predictions and self-couplings. Since the axion quartic self-coupling error was previously dominated by $\ell_7$ uncertainty (approx. 6%), this new precision tightens constraints on axion models and dark matter bounds.
• Isospin Breaking: It provides a robust, non-perturbative benchmark for strong isospin-breaking effects, essential for precision tests of the Standard Model (e.g., in hadron spectroscopy and weak decays).
• Lattice QCD Validation: The agreement between staggered and Wilson fermion results, despite their different discretization artifacts, strongly validates the control over systematic errors in modern lattice QCD simulations.

In conclusion, this paper represents a major step forward in the precision of QCD low-energy constants, delivering a high-precision determination of $\ell_7$ that resolves previous ambiguities and enables more accurate predictions in particle and astroparticle physics.