The three-loop hadronic vacuum polarization in chiral perturbation theory

This paper presents a next-to-next-to-next-to-leading order calculation of the three-loop hadronic vacuum polarization in chiral perturbation theory, achieving unprecedented precision for low-energy contributions and providing a framework to improve lattice QCD computations of the muon magnetic moment.

The Gist

The Big Picture: Why Are We Doing This?

Imagine you are trying to measure the weight of a feather (the muon, a tiny particle) with a scale that is incredibly sensitive. You want to know exactly how heavy it is to test if our current understanding of the universe (the Standard Model) is perfect.

However, the feather is sitting in a room full of invisible, bouncing balloons (virtual particles). These balloons bump into the feather, making it feel slightly heavier or lighter. In physics, this "bumping" is called Hadronic Vacuum Polarization (HVP).

For years, scientists have been arguing about exactly how heavy these balloons are. This uncertainty is the biggest problem in our current measurements of the muon's magnetic moment. If we get the weight of these balloons wrong, we might think we've discovered a new law of physics when we haven't, or vice versa.

The Problem: The "Small Room" Issue

To calculate the weight of these balloons, scientists use supercomputers (Lattice QCD) to simulate the universe. But there's a catch: you can't simulate an infinite universe on a computer. You have to put the simulation in a finite box (a small room).

• The Analogy: Imagine trying to study how sound waves travel in a concert hall, but you are forced to do it inside a tiny closet. The sound bounces off the walls differently than it would in the open air.
• The Physics: In the computer simulation, the "virtual pions" (the balloons) hit the walls of the box. This creates an error called a Finite Volume Effect (FVE). To get the real answer, we need to know exactly how much the "closet walls" are messing up the calculation so we can subtract that error.

The Solution: A Better Map (Chiral Perturbation Theory)

To fix the "closet wall" error, the authors used a tool called Chiral Perturbation Theory (ChPT).

• The Analogy: Think of ChPT as a simplified map of the universe. It's not as detailed as the supercomputer simulation, but it's perfect for understanding how things behave at very low energies (like the slow-moving pions near the walls).
• The Goal: The authors wanted to use this map to calculate exactly how the "closet walls" distort the physics, so the supercomputer results can be corrected.

The Challenge: The "Three-Loop" Mountain

In physics, calculations are done in "orders" of complexity.

• 1-loop: A simple path.
• 2-loop: A path with one detour.
• 3-loop: A path with two detours, getting very twisty.

Previous scientists had only managed to map the "1-loop" and "2-loop" paths. This paper is the first time anyone has successfully climbed the "3-loop" mountain.

Why is this hard?
Most of the paths on this mountain are straightforward. You can break them down into simple pieces (like Lego bricks). However, there are six specific paths (diagrams) that are twisted knots. You can't break them apart.

• The Metaphor: Imagine trying to untangle a knot of headphones. Most of the time, you can just pull one end out. But these six knots are so complex that standard math tools (logarithms) fail to describe them. You need a special, exotic tool called Elliptic Functions (think of these as a "super-math" language that can describe complex, wavy shapes).

The Breakthrough: Solving the Knots

The authors didn't just give up on the six twisted knots. They developed a new strategy:

1. Integration by Parts (IBP): They used a mathematical trick to simplify the equations, reducing thousands of possibilities down to just six "Master Integrals" (the six most important knots).
2. The Dimension Trick: They found a way to translate these complex 4-dimensional knots into simpler 2-dimensional versions, making them easier to handle.
3. The "Schouten" Secret: They discovered a hidden mathematical rule (a symmetry) that made the messy parts cancel each other out, ensuring the final answer was clean and consistent.

The Result: A Blueprint for Precision

The authors have now produced the infinite-volume part of this complex calculation.

• What does this mean? They have built the "perfect map" of the physics without the closet walls.
• Why does it matter? Now, when the supercomputer simulations (which are in the "closet") are compared to this new "perfect map," scientists can calculate the exact error caused by the walls and subtract it.
• The Impact: This allows for a much more precise calculation of the muon's magnetic moment. This is crucial because it helps us decide if the "bumps" in the data are just measurement errors or if they are signs of new physics (like undiscovered particles).

Summary in One Sentence

The authors have successfully solved a massive, three-layered mathematical puzzle involving complex "knots" (elliptic integrals) to create a high-precision map that helps scientists correct errors in supercomputer simulations, bringing us closer to understanding the fundamental secrets of the universe.

(Note: The paper ends with a beautiful image of these mathematical "knots" visualized as a colorful flower arrangement, showing that even the most abstract math can have a hidden beauty.)

Technical Summary

1. Problem Statement

The Hadronic Vacuum Polarization (HVP) is a critical component in low-energy Quantum Chromodynamics (QCD) and is the dominant source of theoretical uncertainty in the calculation of the muon anomalous magnetic moment ($g-2$).

• Lattice QCD Limitations: While Lattice QCD is the primary method for computing HVP from first principles, it suffers from systematic errors related to Finite Volume Effects (FVEs). These errors are particularly severe for low-energy pion modes, which dominate the long-distance part of the HVP.
• The Role of ChPT: Chiral Perturbation Theory (ChPT), the effective field theory of low-energy QCD, excels at describing these low-energy regimes. The difference between ChPT observables calculated in finite and infinite volume provides a precise estimate of the FVEs in Lattice QCD.
• The Gap: Previous FVE corrections utilized ChPT at Next-to-Leading Order (NLO, one-loop) and Next-to-Next-to-Leading Order (NNLO, two-loop). However, to achieve the precision required for future high-accuracy lattice results, corrections at Next-to-Next-to-Next-to-Leading Order (N3LO, three-loops) are necessary. Until this work, the N3LO calculation was considered computationally prohibitive due to the complexity of the required loop integrals.

2. Methodology

The authors performed a calculation of the infinite-volume part of the HVP at N3LO within the framework of $SU(2)$ ChPT (two isospin-symmetric pion flavors coupled to a photon).

Theoretical Framework

• Lagrangian: The calculation uses the standard ChPT Lagrangian expanded in powers of momenta and quark masses.
• Counterterms: The calculation involves a massive number of Low-Energy Constants (LECs). Specifically, there are 475 counterterms at N3LO.
• Power Counting: The amplitude is constructed using Feynman diagrams up to three loops. The expansion parameter is $\xi = M_\pi^2 / (16\pi^2 F_\pi^2) \approx 0.014$.

Computational Techniques

The primary challenge was the evaluation of three-loop integrals, specifically six non-factorizable diagrams highlighted in the paper.

1. Factorizable vs. Non-Factorizable: Most diagrams decompose into products of 1-loop tadpoles and bubbles. However, six specific 3-loop diagrams are non-factorizable and require elliptic functions for their solution, rather than standard polylogarithms.
2. Integration-by-Parts (IBP) Reduction: The authors used the Laporta algorithm (via the LiteRed 2 package) to reduce the thousands of integrals to a finite set of Master Integrals (MIs).
   • They identified six elliptic master integrals: $E_1, \dots, E_6$.
3. Dimensional Recursion (Tarasov's Trick): To handle divergences in 4 dimensions, the authors applied Tarasov's relation, which allows 4D integrals to be expressed as linear combinations of 2D integrals. This separates the issue of divergence from the integration itself.
4. Differential Equations: The master integrals satisfy systems of inhomogeneous differential equations with respect to the photon virtuality $t$.
   • For $E_1, E_2, E_3$, the equations were solved analytically in terms of elliptic polylogarithms.
   • For $E_5$ and $E_6$, the coefficients in the differential equations diverge as $d \to 4$, breaking standard renormalization logic.
5. Schouten Relations: To resolve the divergence and ensure renormalizability, the authors derived non-IBP relations (Schouten identities) based on the Gram matrix of loop momenta. These relations proved that the elliptic functions cancel out among themselves in the final renormalized amplitude, a non-trivial mathematical consistency check.
6. Numerical Solution: The remaining finite parts of the difficult master integrals ($\bar{E}_5, \bar{E}_6$) were solved numerically to high precision using the derived differential equations.

3. Key Contributions

• First N3LO HVP Calculation: This is the first computation of the HVP in ChPT at the three-loop level (N3LO).
• Elliptic Master Integrals: The work successfully evaluates a class of three-loop integrals with massive propagators that require elliptic functions, pushing the boundaries of multiloop integral technology.
• Renormalization Proof: The authors demonstrated that despite the presence of elliptic functions in intermediate steps, the final physical amplitude is renormalizable, relying on hidden Schouten relations to cancel the elliptic contributions.
• Software Infrastructure: The calculation utilized the new ChPTlib library within the FORM system, providing a reusable toolkit for future high-order ChPT calculations.

4. Results

The authors derived the transverse part of the photon vacuum polarization function, $\Pi_T(t)$, expanded in powers of $\xi$:
$$16\pi^2 [\Pi_T(t) - \Pi_T(0)] = \bar{\Pi}_{NLO} + \xi \bar{\Pi}_{NNLO} + \xi^2 \bar{\Pi}_{N3LO} + O(\xi^3)$$

• Analytic Structure: The N3LO result ($\bar{\Pi}_{N3LO}$) is expressed as a linear combination of:
  • Known NLO and NNLO Low-Energy Constants (LECs).
  • A few unknown N3LO LECs (specifically combinations denoted as $\tilde{c}_i$).
  • The function $E(t)$, which is a complex linear combination of the elliptic master integrals ($E_{1,2,3}$ and $\bar{E}_{5,6}$), the basic one-loop integral $B(t)$, and higher polylogarithms.
• Physical Impact: The authors note that terms not containing the elliptic functions or the basic loop integral $B(t)$ have limited physical impact on Finite Volume Effects. Therefore, the determination of specific unknown LECs (like $r_{V1}$, related to the pion charge radius) is the primary remaining task for phenomenological application.
• Validation: The result passed rigorous checks, including the Ward-Takahashi identity ($\Pi_L = 0$) and invariance under field redefinitions.

5. Significance

• Precision for Lattice QCD: This calculation provides the theoretical "blueprint" for estimating Finite Volume Effects in Lattice QCD simulations at the N3LO level. This is crucial for reducing the systematic error in the muon $g-2$ calculation to match the precision of experimental data.
• Toolbox Expansion: The techniques developed for handling elliptic master integrals and Schouten relations in massive theories extend the toolkit available for precision calculations in both QCD and other quantum field theories.
• Control over LECs: While ChPT at this order has limited predictivity due to many unknown constants, the result shows that the HVP depends on a surprisingly small number of physically relevant unknowns. This opens pathways to determine these constants through other phenomenological observables (e.g., pion charge radius).
• Future Outlook: This work lays the foundation for the subsequent calculation of the finite-volume counterpart, which will directly enable improved lattice QCD analyses.

In summary, this paper represents a major theoretical breakthrough in low-energy QCD, overcoming significant mathematical hurdles to provide the necessary precision tools for resolving the muon $g-2$ anomaly.