Electromagnetic pion mass splitting using PV-regulated photon propagator

This paper applies a Pauli-Villars regulated photon propagator defined in continuum infinite volume to lattice QCD calculations on CLS ensembles to determine the electromagnetic pion mass splitting, thereby avoiding power-law finite-size effects and enabling direct cross-checks with phenomenology.

The Gist

Imagine you are trying to weigh two very similar twins: one is a "charged" pion (let's call her Pion Plus) and the other is a "neutral" pion (Pion Zero). In the world of particle physics, these twins are made of the same basic ingredients (quarks), but Pion Plus has an electric charge, while Pion Zero does not.

Because of this electric charge, Pion Plus interacts with the electromagnetic force (light/photons), which makes her slightly heavier than her neutral sister. Physicists have known this difference exists for a long time, but calculating exactly how much heavier she is using the fundamental laws of the universe (Quantum Chromodynamics, or QCD) is incredibly difficult. It's like trying to weigh a feather while standing on a shaking, vibrating scale.

This paper by Alessandro De Santis and colleagues is a story about how they built a shock-absorbing, perfectly stable scale to measure this tiny difference with extreme precision.

Here is the breakdown of their journey, using some everyday analogies:

1. The Problem: The "Shaky Room" and the "Infinite Noise"

In the past, when physicists tried to simulate these particles on a computer, they had to put them in a "box" (a finite volume).

• The Shaky Room: Imagine trying to measure the weight of a feather in a room with walls. The air pressure from the walls pushes back on the feather, messing up your measurement. In physics, this is called "finite-volume effects." The smaller the box, the worse the error.
• The Infinite Noise: When you add electricity (photons) to the mix, the math gets messy. Photons can have infinite energy (ultraviolet divergences), which is like a microphone picking up a deafening static noise that drowns out the signal.

Usually, physicists have to use complex tricks to fix the "shaky room" and the "noise," but these tricks often make it hard to compare results between different research groups. It's like everyone using a different type of ruler; you can't easily check if your measurements match.

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

The authors used a clever new method involving something called a Pauli-Villars (PV) regulated photon propagator. Let's break that down:

• The Infinite Room: Instead of putting the particles in a small, shaky box, they simulated them in an infinite, empty universe. This immediately solves the "shaky room" problem. There are no walls to push back on the pions.
• The Noise Filter: To handle the "infinite noise" of the photons, they introduced a "cutoff scale" (called $\Lambda$). Think of this like a high-quality noise-canceling headphone.
  • Imagine the photons are a chaotic crowd shouting.
  • The "cutoff" is a filter that says, "We only care about the voices below a certain volume."
  • They do the math with this filter in place, and then they slowly turn the volume up (increasing the cutoff) until the filter disappears. Because the math is set up correctly, the result stabilizes, and the "noise" cancels itself out perfectly.

3. The Strategy: Short-Range vs. Long-Range

To get the final answer, they split the problem into two parts, like measuring a journey by looking at the first mile and the rest of the trip separately:

• The Short-Range (The Lattice): For the part of the journey close to the pions (where the physics is messy and complex), they used a supercomputer to simulate the actual particle interactions. This is the "lattice" part.
• The Long-Range (The Analytical Formula): For the part of the journey far away, where the particles are just drifting in the empty infinite space, they didn't need a supercomputer. They used a known mathematical formula (based on the "Cottingham formula") to calculate the effect analytically.

By combining the computer simulation for the messy middle and the math formula for the clean edges, they got a complete picture without the "shaky room" errors.

4. The Result: A Perfect Match

After running these simulations on many different "virtual universes" (using different grid sizes and particle masses) and then mathematically adjusting them to match our real world (where pions have their actual physical mass), they got a result:

• Their Calculation: The mass difference is 4.52 MeV (with a tiny margin of error).
• The Real World (Experiment): The measured difference is 4.59 MeV.

The Verdict: Their calculation matches the real-world experiment almost perfectly!

Why Does This Matter?

This paper is a "proof of concept." It shows that this new method (using the infinite volume and the PV filter) works beautifully for pions.

Why is this a big deal?

1. It's a New Standard: It provides a way for different research groups to compare apples to apples, because they aren't using different "boxes" or "rulers."
2. It's a Stepping Stone: If this works for pions, it can be used for much harder problems, like calculating the mass difference between a proton and a neutron (which makes up the bulk of our visible universe) or the magnetic moment of the muon (a particle that might hold the key to "New Physics" beyond our current understanding).

In a nutshell: The authors built a perfect, infinite, noise-canceling laboratory on a computer to weigh two subatomic twins. They found that their theoretical weight matches the real-world weight almost exactly, proving their new "shock-absorbing" method is ready to tackle even bigger mysteries in physics.

Technical Summary

Here is a detailed technical summary of the paper "Electromagnetic pion mass splitting using PV-regulated photon propagator" by Alessandro De Santis, Dominik Erb, and Harvey B. Meyer.

1. Problem Statement

Lattice QCD calculations have reached sub-percent precision, necessitating the inclusion of electromagnetic (QED) effects and strong isospin-breaking to accurately describe hadronic observables. However, standard methods for incorporating QED on the lattice face two significant challenges:

• Scheme Dependence: The prescription for including QED is formulation-dependent, making cross-checks between different lattice groups and comparisons with phenomenology difficult.
• Finite-Volume Effects: Common approaches (like $QED_L$) introduce power-law suppressed finite-volume effects due to the insertion of a photon propagator in a finite box. These effects complicate the extraction of physical results and hinder direct comparison with infinite-volume phenomenological predictions.

Additionally, the insertion of a photon propagator typically leads to Ultraviolet (UV) divergences if proper counterterms are not included, complicating the renormalization procedure.

2. Methodology

The authors propose and test a computational strategy that avoids these issues by utilizing a Pauli-Villars (PV) regulated photon propagator defined directly in the continuum and infinite volume.

• Theoretical Framework:

  • The method replaces the internal photon line with a PV-regulated propagator, $G_\Lambda(x)$, where $\Lambda$ is a photon virtuality cutoff scale.
  • The propagator is defined in coordinate space in infinite volume:
    $$G_\Lambda(x) = \frac{1}{4\pi^2|x|} - 2G_{\Lambda/\sqrt{2}}(x) + G_\Lambda(x)$$
    where $G_m(x)$ is the propagator of a massive scalar. This ensures the propagator is finite as $|x| \to 0$.
  • The leading-order electromagnetic contribution to the pion mass splitting ($\Delta M_\pi$) is calculated using the "mid-point method," expanding the QCD+QED action in powers of the electromagnetic coupling $\alpha_{em}$.

• Lattice Implementation:

  • Calculations were performed using CLS ($n_f=2+1$) ensembles with various lattice spacings ($a$) and pion masses ($M_\pi$).
  • The calculation involves evaluating Wick contractions for connected and disconnected diagrams. The authors currently neglect the disconnected contribution (estimated to be $\sim 1\%$ of the connected part) due to numerical challenges.
  • To handle the discrete nature of the lattice time integration, the mass splitting is split into two parts:
    1. Short-distance contribution: Calculated directly on the lattice for time separations $z_0 \in [-t_c, t_c]$.
    2. Long-distance contribution: Calculated analytically in the continuum/infinite volume for $z_0 > t_c$, assuming single-pion intermediate states (using a Vector Meson Dominance parametrization for the form factor).
  • This hybrid approach allows for the correction of residual finite-volume effects and avoids power-law artifacts.

• Extrapolation:

  • A chiral-continuum extrapolation is performed to reach the physical pion mass ($M_\pi^{phys} \approx 135$ MeV) and the continuum limit ($a \to 0$).
  • The squared mass splitting $\Delta M_\pi^2$ is fitted using an ansatz based on Chiral Perturbation Theory (ChPT), including terms for lattice artifacts ($a^2$) and pion mass dependence ($M_\pi^2 \log M_\pi^2$).

3. Key Contributions

• Validation of Infinite-Volume PV Propagator: The paper demonstrates that using a PV-regulated propagator defined in infinite volume successfully eliminates power-law finite-volume effects, a major advantage over standard lattice QED methods.
• Separation of Elastic and Inelastic Contributions: The introduction of the scale $\Lambda$ facilitates the separation of the pion mass splitting into:
  • Elastic contribution: Determined phenomenologically via the Cottingham formula using the modified PV propagator.
  • Inelastic contribution: Calculated from the lattice data by subtracting the elastic part.
• UV Finiteness Check: The authors show that while general hadron mass splittings diverge as $\Lambda \to \infty$ without counterterms, the pion mass splitting is naturally UV-finite. The inelastic contribution is shown to be independent of $\Lambda$ for sufficiently large cutoffs, validating the method.

4. Results

• Convergence: The inelastic contribution reaches a plateau quickly, stabilizing at $\Lambda \approx 20 m_\mu$ (where $m_\mu$ is the muon mass).
• Final Value:
  • The elastic contribution in the $\Lambda \to \infty$ limit is calculated as 4.33 MeV.
  • The inelastic contribution is found to be 0.19(15) MeV.
  • The total preliminary lattice result for the pion mass splitting is:
    $$\Delta M_\pi^{latt} = 4.52(15) \text{ MeV}$$
• Comparison with Experiment: This result is in excellent agreement with the experimental measurement:
  $$\Delta M_\pi^{exp} = 4.5936(5) \text{ MeV}$$

5. Significance

This work provides a robust, scheme-independent methodology for calculating electromagnetic corrections in lattice QCD. By avoiding power-law finite-volume effects and utilizing a continuum-defined propagator, the method offers a clearer path for cross-checking results between different lattice groups and comparing them with phenomenological predictions.

Although the pion mass splitting is UV-finite and does not strictly require the PV regulator for renormalization, this study validates the formalism (first proposed for the muon $g-2$). This paves the way for applying the same strategy to more complex observables where UV divergences are present, such as:

• The proton/neutron mass splitting.
• The hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment.

The successful application to the pion mass splitting confirms that the PV-regulated infinite-volume approach is a viable and powerful tool for high-precision hadronic physics.