Electromagnetic pion mass splitting using a Pauli-Villars-regulated photon propagator

This paper presents a lattice QCD calculation of the charged-neutral pion mass splitting using a Pauli-Villars-regulated photon propagator to avoid power-law finite-volume effects, yielding a result of 4.56(22) MeV that agrees well with experimental measurements and validates the formalism for future electromagnetic corrections.

The Gist

Imagine the universe as a giant, complex machine built from tiny, invisible Lego bricks called quarks. When these quarks stick together, they form particles like protons, neutrons, and pions. For a long time, scientists thought that two specific types of pions—one positively charged ($\pi^+$) and one neutral ($\pi^0$)—should be identical twins, weighing exactly the same.

However, in reality, they aren't twins; they are cousins with a slight difference in weight. The charged pion is just a tiny bit heavier than the neutral one. This tiny difference is caused by electromagnetism (the force that makes magnets stick and lightning strike).

This paper is a report from a team of scientists who used a supercomputer to calculate exactly how big this weight difference is, and to prove that their new method for doing the math is reliable.

Here is a breakdown of what they did, using simple analogies:

1. The Problem: Simulating a "Perfect" World vs. Reality

To understand why the pions weigh different amounts, scientists use a technique called Lattice QCD. Imagine the universe is a giant 3D grid (like a chessboard, but in 4 dimensions). They place the quarks on this grid and simulate how they interact.

Usually, to make the math easier, scientists pretend the grid is infinite and the world is perfectly symmetrical. But in the real world:

• The grid is actually finite (it has edges).
• There is electromagnetism (photons) zipping around, which makes the math messy.

When you try to simulate electromagnetism on a finite grid, you get "echoes" or "ghosts" bouncing off the walls. In physics terms, these are called finite-volume effects. It's like trying to measure the sound of a whisper in a small room; the echo makes it hard to hear the true volume. Previous methods struggled with these echoes, making the calculations very difficult and prone to errors.

2. The Solution: The "Pauli-Villars" Filter

The authors of this paper used a clever new trick involving something called a Pauli-Villars (PV) regulated photon propagator.

Think of the photon (the particle of light) as a messenger running between the quarks.

• Old Method: The messenger runs forever. On a finite grid, the messenger hits the wall and bounces back, creating noise (the echoes).
• New Method (PV): The scientists put a "speed limit" or a "filter" on the messenger. They introduced a scale called $\Lambda$ (Lambda). This filter acts like a pair of noise-canceling headphones. It stops the messenger from running into the walls of the simulation box.

Because of this filter, the simulation behaves as if it were in an infinite universe, even though the computer grid is finite. This removes the "echoes" and makes the calculation much cleaner.

3. The Challenge: Removing the Filter

There is a catch. The filter ($\Lambda$) is an artificial tool. In the real world, there is no such speed limit on photons. So, the scientists had to do a two-step dance:

1. Run the simulation with the filter set to different strengths (different values of $\Lambda$).
2. Turn the filter off (let $\Lambda$ go to infinity) to see what the result looks like in the real world.

They found that the "charged" part of the pion's mass (the part that comes from the photon interacting with the pion itself) was the biggest player. They could calculate this part using a known formula (the Cottingham formula), which is like using a trusted map to verify their GPS.

4. The Result: A Perfect Match

After running thousands of simulations on different grid sizes and removing the artificial filter, they calculated the final weight difference:

The charged pion is heavier than the neutral pion by 4.56 MeV (with a tiny margin of error).

• Why this matters: This number matches the experimental measurement (what we see in real particle accelerators) almost perfectly.
• The "Connected" vs. "Disconnected" parts: The calculation had two parts. The main part (the "connected" diagram) was the heavy lifter. The second part (the "disconnected" diagram) was like a tiny, faint whisper in the background. They calculated it too, and found it was very small, confirming that the main calculation was doing the heavy lifting.

5. The Conclusion: A New Tool for the Toolbox

The paper doesn't just give us a number; it proves that their new "noise-canceling headphone" method (the PV-regulated propagator) works.

• Validation: They showed that this new way of handling electromagnetism on a computer grid is robust and accurate.
• Future Use: Because this method works so well for pions, the scientists are now confident they can use it to solve even harder puzzles, like calculating the mass difference between protons and neutrons or improving calculations for the "muon g-2" (a famous experiment looking for new physics).

In summary: The scientists built a new, quieter simulation room to study how light affects the weight of pions. They proved that by using a special filter to block out the "echoes" of the simulation box, they could calculate the weight difference with high precision, matching reality perfectly. This success means they now have a powerful new tool to study the fundamental forces of nature.

Technical Summary

Technical Summary: Electromagnetic Pion Mass Splitting using a Pauli-Villars-Regulated Photon Propagator

Problem Statement
Lattice QCD calculations have achieved sub-percent precision for many phenomenological observables, yet standard approximations—specifically the isosymmetric limit (degenerate light quark masses) and the neglect of electromagnetic effects ($\alpha_{em} = 0$)—are no longer sufficient for meaningful comparison with experimental data. Including Quantum Electrodynamics (QED) and isospin-breaking effects introduces significant technical challenges. While perturbative expansions around the isosymmetric point are widely adopted, the implementation of the photon propagator on a finite lattice is non-unique. Various regularization schemes (e.g., QED$_L$, QED$_\infty$) exist, but they often suffer from power-law finite-volume effects or require complex renormalization counterterms to remove ultraviolet (UV) divergences. This complicates cross-checks between collaborations and comparisons with phenomenology. The specific observable addressed here is the charged-neutral pion mass splitting, $M_{\pi^+} - M_{\pi^0}$, at order $O(\alpha_{em})$. While this splitting is a leading-order electromagnetic effect, its precise calculation requires a robust framework to handle finite-volume artifacts and the continuum limit.

Methodology
The authors employ a recently proposed framework [14] that utilizes a Pauli-Villars (PV) regulated photon propagator defined directly in the continuum and infinite-volume limit. The key features of this methodology include:

1. PV-Regulated Propagator: The photon propagator is replaced by a PV-regulated version, $G_\Lambda(x)$, characterized by a scale $\Lambda$. This propagator acts as a UV-finite weight function for QCD correlation functions. The physical result is obtained by taking the limit $\Lambda \to \infty$ after the continuum limit ($a \to 0$).
2. Perturbative Expansion: The calculation treats $\alpha_{em}$ and the quark mass difference $\delta m_\ell$ as perturbative corrections to pure QCD simulations performed at the isosymmetric point. The mass splitting is computed via the "midpoint method," inserting the photon propagator between two vector currents in two-point correlation functions.
3. Diagram Topologies: The calculation includes both quark-connected and quark-disconnected diagrams. The connected contribution dominates the splitting, while the disconnected contribution (which vanishes in the $SU(2)$ chiral limit) is computed on a subset of ensembles.
4. Finite-Volume Treatment: By defining the propagator in infinite volume, the method avoids the power-law finite-volume effects typical of schemes like QED$_L$. The long-distance contribution to the mass splitting is reconstructed analytically in the continuum using the elastic contribution (single-pion intermediate states) and the Cottingham formula, while the short-distance part is calculated on the lattice.
5. Ensembles and Extrapolation: Data were generated using 15 CLS ensembles with $N_f=2+1$ dynamical quarks, covering four lattice spacings ($a \in [0.049, 0.085]$ fm) and pion masses up to $\sim 360$ MeV. The results are extrapolated to the physical pion mass and the continuum limit. The dependence on the cutoff $\Lambda$ is studied for values ranging from $3 m_\mu$ to $80 m_\mu$.

Key Contributions

• Validation of PV Formalism: This work provides a complete, first-principles validation of the PV-regulated photon propagator formalism in a theoretically clean setting. It demonstrates that results obtained at fixed $\Lambda$ are universal and can be compared across different lattice actions.
• Separation of Elastic and Inelastic Contributions: The authors decompose the mass splitting into elastic and inelastic parts. They show that the elastic contribution, calculable via the Cottingham formula, dominates the $\Lambda$ dependence, while the inelastic contribution saturates quickly as $\Lambda$ increases.
• Handling of Disconnected Diagrams: The paper includes a dedicated analysis of the disconnected contribution, confirming its small magnitude ($\sim 1\%$ of the total splitting) and providing a numerical estimate.
• Systematic Error Control: The study rigorously addresses systematic uncertainties, including finite-volume effects, discretization errors, and the extrapolation to the physical point and $\Lambda \to \infty$. An alternative analysis strategy (extrapolating $\Lambda \to \infty$ before the continuum limit) is presented in Appendix B as a consistency check.

Results
The final result for the electromagnetic pion mass splitting, after removing the cutoff scale $\Lambda$ and including both connected and disconnected contributions, is:
$$M_{\pi^+} - M_{\pi^0} = 4.56(22) \text{ MeV}$$
where the uncertainty combines statistical and systematic errors.

• Comparison with Experiment: This result is in excellent agreement with the experimental value of $4.5936(5)$ MeV.
• Comparison with Other Lattice Calculations: The result is consistent with previous lattice determinations using different regularization schemes (e.g., $4.622(95)$ MeV [17] and $4.534(60)$ MeV [18]).
• Elastic vs. Inelastic: The elastic contribution accounts for more than 90% of the mass splitting. The inelastic contribution is determined to be $0.29(20)$ MeV.
• Cutoff Dependence: The study confirms that the pion mass splitting is UV-finite in the continuum limit, allowing for large values of $\Lambda$ (up to $8.5$ GeV) without requiring counterterms, unlike generic hadron mass splittings.

Significance
The paper claims that this calculation serves as a robust validation of the PV-regulated photon propagator formalism. By successfully reproducing the experimental pion mass splitting with high precision, the authors demonstrate the consistency and practicality of this approach. The methodology offers a pathway to compute electromagnetic corrections for more complex observables, such as the hadronic vacuum polarization contribution to the muon $g-2$ and the kaon mass splitting, where UV finiteness is not guaranteed and renormalization is more challenging. The ability to define the photon propagator in infinite volume and separate elastic/inelastic contributions provides a powerful tool for controlling finite-volume effects and isolating specific physical mechanisms in lattice QCD+QED calculations.