Third-order mixed electroweak-QCD corrections to the W-boson mass prediction from the muon lifetime

This paper presents the calculation of the previously missing ${\cal O}(\alpha^2\alpha_\mathrm{s})$ corrections to the quantity $\Delta r$ involving one closed fermion loop, which significantly increases the Standard Model prediction for the W-boson mass by over 3 MeV.

The Gist

Imagine the Standard Model of particle physics as the ultimate "Rulebook of the Universe." It tells us how the tiniest building blocks of nature interact. One of the most important characters in this book is the W-boson, a heavy particle that acts like a cosmic delivery truck, carrying the "weak force" that allows particles to change their identity (like a neutron turning into a proton).

For decades, physicists have been trying to measure the exact "weight" (mass) of this W-boson with incredible precision. Why? Because if the measured weight doesn't match the weight predicted by the Rulebook, it means there's a hidden chapter in the book we haven't found yet—perhaps a new particle or a new force.

The Problem: The "Rounding Error"

Think of calculating the W-boson's mass like trying to balance a very delicate scale. You have a known weight (the muon, a heavy cousin of the electron) and you use it to predict the W-boson's weight.

However, the universe is messy. Particles don't just sit still; they constantly pop in and out of existence as "virtual" particles, creating a chaotic foam of activity around them. To get the right answer, physicists have to add up the effects of all these virtual particles.

For a long time, the prediction was missing a specific piece of the puzzle: a very complex, three-layer interaction involving both the electroweak force (the W-boson's job) and QCD (the strong force that holds quarks together).

Imagine you are baking a cake (the W-boson mass). You've added the flour, sugar, and eggs (the known forces). But you forgot the pinch of salt that comes from a specific interaction between the strong force and the weak force. Without it, the cake tastes slightly off.

The Solution: The "Three-Layer Cake" Calculation

This paper is about finding that missing pinch of salt. The authors, a team of theoretical physicists, performed a massive calculation to figure out the O(α²αs) correction.

In plain English, this means they calculated the effects of three loops of virtual particles interacting at the same time.

• The Challenge: Calculating three loops is like trying to solve a Rubik's Cube while blindfolded, while the cube is on fire, and while someone is shouting math equations at you. The number of possible paths the particles can take is astronomical.
• The Method: They used a combination of "analytical" math (solving equations on paper) and "numerical" math (using supercomputers to crunch numbers). They had to be incredibly careful because if you miss a tiny detail, the whole result explodes into nonsense.

They also had to deal with a tricky mathematical ghost called γ5 (gamma-five). In the world of quantum mechanics, this represents a "handedness" or chirality. Usually, math breaks when you try to calculate this in certain ways. The authors had to use a special "shield" (a technique called Pauli-Villars regularization) to keep the math from breaking while they figured out how these particles twist and turn.

The Result: A Shift in the Universe

When they finally added this missing three-loop correction to their calculation, the result was significant.

• The Shift: The predicted mass of the W-boson went up by about 3.14 MeV (Mega-electronvolts).
• Why it matters: 3 MeV sounds tiny, but in the world of particle physics, it's a massive jump. It's the difference between a car traveling at 100 mph and 100.003 mph. But when you are measuring the speed of light, that difference is huge.

Previously, scientists thought the uncertainty in their prediction was about 4 MeV. This new calculation effectively "fills in the gap," reducing the uncertainty and making the prediction much sharper.

The Analogy: The GPS Update

Think of the Standard Model as a GPS navigation system.

• Old GPS: It got you to the right city, but maybe you arrived 5 minutes late or took a slightly wrong turn because the map was missing a new road.
• The New GPS (This Paper): The authors just uploaded a massive software update that accounts for a complex, winding road that was previously ignored.
• The Outcome: Now, the GPS predicts your arrival time with much higher precision. If the real-world arrival time still doesn't match the new, ultra-precise GPS prediction, we know for sure that there is a new road (new physics) that we haven't even discovered yet.

Conclusion

This paper is a triumph of human calculation. It fills a critical gap in our understanding of the universe's rulebook. By refining the prediction for the W-boson's mass, the authors have sharpened our "telescope" for looking at the universe.

Now, when experiments at the Large Hadron Collider (LHC) or future colliders measure the W-boson's mass, they will have a much clearer target. If the measurement still disagrees with this new, highly precise prediction, it will be a smoking gun for New Physics—a discovery that could rewrite our understanding of reality.

In short: They fixed a tiny, three-dimensional math error in the universe's rulebook, and that small fix makes our view of the cosmos much clearer.

Technical Summary

Here is a detailed technical summary of the paper "Third-order mixed electroweak-QCD corrections to the W-boson mass prediction from the muon lifetime."

1. Problem Statement

The W-boson mass ($M_W$) is a fundamental parameter of the Standard Model (SM) and a critical test of electroweak symmetry breaking. Its value is indirectly determined with high precision via the muon lifetime ($\tau_\mu$), which relates the Fermi constant ($G_\mu$) to $M_W$ through the radiative correction parameter $\Delta r$:
$$M_W^2 \left(1 - \frac{M_W^2}{M_Z^2}\right) = \frac{\pi \alpha}{\sqrt{2}G_\mu} (1 + \Delta r)$$
While experimental precision for $M_W$ is rapidly improving (targeting $\sim 0.18$ MeV at future lepton colliders like FCC-ee), the theoretical uncertainty in $\Delta r$ has been limited by missing higher-order corrections. Specifically, the mixed electroweak-QCD corrections of order $O(N_f \alpha^2 \alpha_s)$ (where $N_f$ is the number of closed fermion loops) were the last missing piece of the third-order (NNNLO) corrections for the subset with a single closed fermion loop. Previous calculations covered two-loop fermion loops ($O(N_f^2 \alpha^2 \alpha_s)$) and purely electroweak terms, but the genuine three-loop amplitudes with one fermion loop remained uncalculated.

2. Methodology

The authors performed a complete calculation of the $O(N_f \alpha^2 \alpha_s)$ corrections to muon decay ($\mu^- \to \nu_\mu e^- \bar{\nu}_e$) within the full Standard Model. The calculation involved several sophisticated technical steps:

• Diagram Generation & Algebra:

  • Generated Feynman diagrams using FeynArts, DiaGen, and custom electroweak model files.
  • Performed Lorentz and Dirac algebra using FeynCalc and FormCalc, cross-checked with in-house implementations.
  • $\gamma_5$ Treatment: A subset of diagrams involves fermion triangles with gluonic corrections (PCAC-violating). Since standard dimensional regularization is incompatible with the trace identity involving $\gamma_5$ in $d$ dimensions, the authors employed Pauli-Villars (PV) regularization for these specific vertex diagrams. They verified gauge invariance and consistency by testing large regulator masses ($\Lambda = 5 \times 10^7$ and $5 \times 10^8$ GeV).

• Integration and Reduction:

  • Reduced the resulting Feynman integrals to Master Integrals (MIs) using Integration-by-Parts (IBP) identities via Kira 3 and FIRE 5.
  • Evaluated the MIs numerically using AMFlow (based on auxiliary mass flow) and TVID2.
  • Cross-checked results against analytical expansions in $\epsilon = (4-d)/2$ using HypExp 2 where available.

• Renormalization Scheme:

  • Utilized the On-Shell (OS) renormalization scheme, defining masses via the pole of the radiatively corrected propagators.
  • Accounted for the complex poles of the $W$, $Z$, and top quark propagators (due to non-negligible decay widths) in the mass counterterms.
  • Calculated three-loop self-energies at non-zero momentum, a computationally demanding task requiring high numerical precision (8–10 digits) to handle cancellations.

• Matching to Fermi Theory:

  • To isolate the finite $\Delta r$, the SM amplitude was matched to the Fermi theory.
  • Infrared (IR) divergences arising from photon exchange between external leptons were canceled by subtracting the corresponding QED corrections calculated within the Fermi model (using a Fierz transformation to simplify the Dirac structure).

3. Key Contributions

• First Complete Calculation: This work provides the first complete calculation of the genuine three-loop $O(N_f \alpha^2 \alpha_s)$ corrections to muon decay, involving diagrams with one closed fermion loop (specifically top-bottom loops).
• Technical Validation: The calculation was validated by achieving cancellations of UV-divergent terms ($\epsilon^{-3}, \epsilon^{-2}, \epsilon^{-1}$) at 19, 15, and 8 digits of precision, respectively, across two independent calculation chains.
• Handling of $\gamma_5$: The successful application of Pauli-Villars regularization to handle the non-vanishing contributions from PCAC-violating diagrams (where top and bottom quarks have different masses, violating the Adler-Bardeen theorem conditions) is a significant methodological achievement.
• Breakdown of Contributions: The authors provided a detailed breakdown of the correction, separating contributions from IR-finite boxes, QED boxes (matched to Fermi theory), $\gamma_5$ traces, 1-particle reducible diagrams, and the weak mixing angle counterterm.

4. Results

Using standard input parameters ($M_Z = 91.1876$ GeV, $M_W = 80.385$ GeV, $M_H = 125.1$ GeV, $m_t = 173.2$ GeV), the authors derived the following:

• Correction to $\Delta r$:
  The total $O(N_f \alpha^2 \alpha_s)$ contribution is:
  $$\Delta r^{(N_f \alpha^2 \alpha_s)} = -1294.9 \, C_F \left(\frac{\alpha}{\pi}\right)^2 \frac{\alpha_s}{\pi}$$
  Isolating the new contribution beyond the previously known $\Delta \rho$ parameter (which captures the leading $m_t^4$ terms):
  $$\Delta r^{(N_f \alpha^2 \alpha_s)}_{\text{beyond } \Delta \rho} = -745.0 \, C_F \left(\frac{\alpha}{\pi}\right)^2 \frac{\alpha_s}{\pi} \approx -2.0 \times 10^{-4}$$

• Impact on $M_W$:
  This correction induces a shift in the predicted W-boson mass:
  $$\Delta M_W = +3.14 \text{ MeV}$$

• Validity of Large-$m_t$ Expansion:
  The authors compared their full numerical result with the leading $O(\alpha_t^2 \alpha_s)$ term (derived from a large-$m_t$ expansion). They found that the leading term is a poor approximation of the full result (ratio $\approx 1.81$), indicating that sub-leading terms in the $1/m_t$expansion are numerically significant. This suggests limitations in relying solely on large-$m_t$ expansions for precision electroweak observables at the three-loop level.

5. Significance

• Theoretical Precision: The new correction shifts the SM prediction for $M_W$ by 3.14 MeV, a magnitude comparable to the previous theoretical uncertainty estimate (3 MeV) attributed to these unknown corrections.
• Future Experiments: This result is crucial for the High-Luminosity LHC (HL-LHC) and future lepton colliders (FCC-ee, CEPC), which aim for experimental uncertainties on $M_W$ of $\sim 5$ MeV and $\sim 0.18$ MeV, respectively. Without this calculation, the theoretical error would dominate the experimental precision.
• Completing the Picture: With this calculation, the only missing third-order correction to the W-boson mass prediction from the muon lifetime is the purely electroweak $O(\alpha^3)$ term.
• Benchmark for New Physics: By reducing the SM theoretical uncertainty, this work enhances the sensitivity of $M_W$ measurements to potential deviations caused by physics Beyond the Standard Model (BSM).

In conclusion, this paper represents a major milestone in precision electroweak physics, resolving a long-standing theoretical gap and enabling the next generation of high-precision tests of the Standard Model.