⚛️ lattice

Lattice QCD Calculation of the Pion Mass Splitting

Using the infinite volume reconstruction method on six $2+1$ -flavor domain wall fermion ensembles, this study calculates the charged-neutral pion mass splitting as $4.534(42)(43)$ MeV with exponentially suppressed finite-volume errors, achieving excellent agreement with experimental measurements.

Original authors: Xu Feng, Luchang Jin, Michael Joseph Riberdy

Published 2026-08-04

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Xu Feng, Luchang Jin, Michael Joseph Riberdy

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is built from a giant, invisible LEGO set. The smallest, most fundamental bricks are called quarks, and they snap together to form larger structures like protons and neutrons. But these bricks don't just stick together; they are glued by a super-strong force called the "strong interaction," which is so powerful that it creates almost all the mass of the things we see around us. However, there's a second, much weaker force at play: electromagnetism. You know this as the force behind magnets and lightning. While the strong force is the heavyweight champion, electromagnetism is the subtle artist that can still tweak the final shape of the LEGO creations.

One of the most famous creations in this cosmic LEGO set is the pion. Think of pions as the "glue" that holds atomic nuclei together, acting as messengers between protons and neutrons. Nature gives us two main types of pions: the neutral pion (which has no electric charge) and the charged pion (which carries a positive or negative charge). In a perfect world where electromagnetism didn't exist, these two twins would weigh exactly the same. But because the charged pion feels the electromagnetic force and the neutral one doesn't, they end up with slightly different weights. This tiny difference in weight is called "mass splitting." Physicists care deeply about this because measuring it with extreme precision is like checking the fine print of the universe's rulebook; if our calculations don't match the real world, it means we might be missing a piece of the puzzle about how the strong and electromagnetic forces dance together.

In this study, a team of scientists used a powerful digital microscope called "Lattice Quantum Chromodynamics" (Lattice QCD) to calculate this tiny weight difference from scratch. Instead of just looking at real pions in a lab, they built a virtual, 3D grid of space-time on a supercomputer and simulated the behavior of quarks and gluons. They used a clever trick called the "infinite volume reconstruction method," which is like taking a photo of a crowded room through a small window and mathematically reconstructing the entire room so you don't miss anything happening in the corners. By running these simulations on six different virtual grids of varying sizes and resolutions, they were able to filter out the "noise" caused by the small size of their digital world.

The researchers found that the charged pion is heavier than the neutral pion by exactly 4.534(42)stat(43)sys MeV. To put that in perspective, the numbers in the parentheses represent the margin of error, showing that their calculation is incredibly precise. When they compared their result to the actual weight difference measured in real-world experiments, which is 4.5936(5) MeV, the two numbers matched up beautifully. The team also discovered that they had to account for two different types of interactions: "connected" diagrams, where the particles talk directly to each other, and "disconnected" diagrams, which are like ghostly whispers where particles interact through a complex, indirect loop. They confirmed that even these ghostly interactions play a role, and by including them, their simulation became even more accurate.

The study also tested two different mathematical "languages" (called Feynman gauge and Coulomb gauge) to describe how the electromagnetic force works. Just as you can describe a landscape using a map or a 3D model, these two methods look different on paper but should lead to the same destination. The authors found that as they made their digital grid finer and finer (approaching the "continuum limit"), both methods agreed perfectly, giving them confidence that their result wasn't just a fluke of the math they chose. This work doesn't just confirm a number; it proves that we can now simulate the messy, complex interplay of the strong and electromagnetic forces with such high precision that we can trust our digital models to explore the deepest secrets of matter.

Technical Summary: Lattice QCD Calculation of the Pion Mass Splitting

Problem Statement
The primary objective of this work is to calculate the mass splitting between the charged ( $\pi^\pm$ ) and neutral ( $\pi^0$ ) pions ( $\Delta m_\pi$ ) using Lattice Quantum Chromodynamics (QCD) combined with Quantum Electrodynamics (QED). While the strong interaction binds quarks into hadrons, the mass difference between charged and neutral pions arises primarily from electromagnetic (EM) interactions, as isospin-breaking effects due to the up/down quark mass difference are suppressed by a factor of $\sim 10^{-4}$ . Previous lattice calculations of this splitting have struggled with finite-volume errors and varying levels of precision. This study aims to provide a high-precision calculation at the physical pion mass, explicitly resolving contributions from both connected and disconnected quark diagrams.

Methodology
The authors employ the Infinite Volume Reconstruction (IVR) method, originally proposed in Ref. [28], to calculate the EM self-energy of stable hadrons. This approach allows for the computation of electromagnetic corrections with only exponentially suppressed finite-volume errors. The calculation proceeds through the following steps:

Formalism: The mass shift $\Delta M$ is expressed as an infinite volume Euclidean space-time integral of the hadronic tensor $H_{\mu,\nu}(x)$ convolved with an analytically known photon propagator $S^\gamma_{\mu,\nu}(x)$ .
$\Delta M = \frac{1}{2} \int d^4x H_{\mu,\nu}(x) S^\gamma_{\mu,\nu}(x)$
The hadronic tensor is defined via the time-ordered product of electromagnetic currents $J_\mu$ within the pion state.
IVR Implementation: The integral is split into a "short-distance" part ( $|x_t| < t_s$ ) and a "long-distance" part ( $|x_t| \ge t_s$ ).
- The short-distance part is calculated directly on the finite lattice volume using lattice-calculable matrix elements.
- The long-distance part is approximated using the lightest single-particle intermediate states, weighted by an infinite-volume QED function $L_{\mu,\nu}$ .
- The separation parameter is set to $t_s = L/2$ to ensure both excited-state contamination and finite-volume errors are exponentially suppressed.
Gauge Fixing: Calculations are performed in two distinct gauges to test consistency: Feynman gauge and Coulomb gauge. The photon propagator and weighting functions are derived analytically for both gauges. Special regularization is applied to the photon propagator at the origin ( $x=0$ ) to handle continuum divergences on the lattice.
Ensembles and Diagrams: The calculation utilizes six gauge ensembles generated by the RBC and UKQCD collaborations using $2+1$ -flavor domain wall fermions. Five ensembles are at the physical pion mass ( $m_\pi \approx 135$ MeV), and one is at $m_\pi \approx 340$ MeV.
- Connected and Disconnected Diagrams: The calculation includes all relevant quark-field contractions. Crucially, this work explicitly includes the quark disconnected diagram (related to $\pi^0-\eta-\eta'$ mixing), which had not been clearly resolved in previous lattice studies of this splitting.
- Techniques: To improve statistical precision, the authors utilize All-modes-averaging (AMA), $z$ -Möbius, and compressed eigenvector deflation.

Key Results
After performing continuum extrapolations and correcting for residual finite-volume effects, the authors obtain the following results for the pion mass splitting:

Feynman Gauge: $\Delta m_\pi = 4.534(42)_{\text{stat}}(43)_{\text{sys}}$ MeV.
Coulomb Gauge: $\Delta m_\pi = 4.560(46)_{\text{stat}}(41)_{\text{sys}}$ MeV.

The results from both gauges are consistent with each other and with the experimental value of $4.5936(5)$ MeV. The total uncertainty is at the percent level, representing an improvement of approximately 5–10 times over previous lattice QCD calculations.

The breakdown of contributions reveals:

The connected diagrams provide the dominant contribution ( $\sim 4.5$ MeV).
The disconnected diagrams contribute a smaller but non-negligible amount ( $\sim 0.05$ MeV).
The Coulomb potential (time-component of the photon propagator) accounts for approximately $1.86$ MeV of the splitting.

Significance and Claims
The paper claims several specific achievements:

First Application of IVR to Pion Splitting: This is the first application of the Infinite Volume Reconstruction method to the pion mass splitting, demonstrating its viability for calculating electromagnetic corrections to stable hadron masses.
Resolution of Disconnected Contributions: For the first time in the literature, the contribution from the quark disconnected diagram to the pion mass splitting is clearly resolved and included, addressing a component related to $\pi^0-\eta-\eta'$ mixing.
Gauge Consistency: The agreement between results obtained in Feynman and Coulomb gauges, particularly as the lattice spacing approaches zero, validates the methodology and the handling of gauge-dependent terms.
Precision: The calculation achieves a precision that matches experimental measurements, confirming the success of combining lattice QCD with QED corrections using the IVR method.

The authors conclude that this high-precision agreement demonstrates the success of both lattice QCD and the IVR method, establishing a foundation for future lattice QCD + QED spectroscopy calculations.

Technical Summary: Lattice QCD Calculation of the Pion Mass Splitting

More like this