The hadronic contribution to the running of the electroweak gauge couplings

Using $N_f=2+1$ CLS lattice ensembles with a refined analysis strategy involving telescopic series and new kernel functions, this paper presents an updated, high-precision ab initio determination of the hadronic vacuum polarization contribution to the running of the electromagnetic coupling and the electroweak mixing angle, aiming to meet the projected precision requirements of future electroweak measurements at next-generation colliders.

The Gist

The Big Picture: Why Do We Care?

Imagine the universe is a giant,精密 (precision) clockwork machine. Physicists call this the Standard Model. To make sure the clock is ticking correctly, they need to know the exact size of every gear and spring. One of the most important "springs" is the electromagnetic force (the force behind electricity, magnetism, and light).

However, this force isn't constant. It changes strength depending on how close you are to the action. Think of it like a flashlight: the light is dim when you are far away, but blindingly bright when you are right up close.

Physicists need to know exactly how bright the "flashlight" is at the specific moment two particles smash together at the highest energies possible (like at the Large Hadron Collider). This brightness is called the electromagnetic coupling.

The Problem: The "Dirty" Middle Ground

Calculating how this force changes is easy for some parts (like the light particles, or "leptons") because they behave like predictable billiard balls. But there is a messy middle ground involving quarks and gluons (the stuff inside protons and neutrons).

This messy part is called Hadronic Vacuum Polarization (HVP).

• The Analogy: Imagine trying to measure the temperature of a room, but the room is filled with a thick, churning fog. You can't just stick a thermometer in and get a reading; the fog itself distorts the measurement.
• The Issue: For years, scientists tried to measure this "fog" by looking at data from particle collisions (like watching the fog swirl from a distance). But different experiments gave slightly different answers, creating a "tension" or disagreement.

The Solution: Building a Virtual Room

Instead of watching the fog from the outside, this team of scientists decided to build a virtual room inside a supercomputer and simulate the fog from the ground up. This is called Lattice QCD.

• The Lattice: Imagine a giant 3D grid (like a fishing net) stretching through space. They place the particles on the intersections of this grid.
• The Simulation: They run a simulation where the particles interact on this grid, effectively recreating the "fog" of the vacuum.

What Did They Do Differently?

The authors (a team from Germany, Switzerland, Japan, and the US) didn't just run the simulation; they upgraded the tools they used to look at the results.

1. The "Telescopic" Lens:

   • The Problem: The "fog" behaves differently depending on how close you look. Up close, the grid lines of the computer simulation mess things up (discretization errors). Far away, the simulation gets noisy and fuzzy (statistical noise).
   • The Fix: They used a telescopic strategy. They broke the problem into three zones:
     • Short Distance: Looked at the immediate neighborhood where the grid lines matter most.
     • Medium Distance: The transition zone.
     • Long Distance: The far-away zone where the physics is smooth but noisy.
   • By treating each zone with a different mathematical "lens," they could clean up the errors in each specific area.

2. The "Split" Technique:

   • To get the final answer for the high-energy collisions (at the $Z$-boson pole), they had to connect their computer simulation (which works in "Euclidean" time, a weird mathematical space) to real-world time.
   • They used a split technique: They let the computer do the hard, messy part (the low-energy fog), and then used standard math (perturbation theory) to bridge the gap to the high-energy part. It's like using a boat to cross a swamp, then switching to a car for the highway.

The Results: A Sharper Picture

• Precision: They achieved a precision that is roughly twice as good as previous estimates.
• The Tension: Their results confirmed a long-standing mystery. The "fog" they calculated from first principles is slightly different from the "fog" measured by looking at experimental data (the $R$-ratio). This suggests either our understanding of the experimental data is slightly off, or there is new physics we don't understand yet.
• The $Z$-Boson: They calculated the value of the electromagnetic coupling at the energy level of the $Z$-boson (a heavy particle). This is a crucial number for future experiments.

Why Does This Matter for the Future?

The next generation of particle colliders (like the proposed FCC-ee) will be incredibly precise. They want to measure things with an accuracy of 1 part in 10,000.

• The Bottleneck: Currently, the "fog" (HVP) is the biggest source of uncertainty. It's like trying to aim a laser at a target, but your glasses are slightly blurry.
• The Goal: This paper shows that by improving the computer simulations (making the grid finer and the math smarter) and combining them with better theoretical math, we can clear up the blur.
• The Verdict: They showed that with about double the effort in computer power and a slightly higher "matching scale" (looking at higher energies in the simulation), we can reach the precision needed for these future machines.

Summary in One Sentence

This paper uses advanced supercomputer simulations to clear up the "fog" of quantum particles, providing a much sharper, more precise map of how the electromagnetic force behaves, which is essential for testing the laws of physics at the highest energy levels.

Technical Summary

1. Problem Statement

Precision tests of the Standard Model (SM) rely heavily on accurate determinations of fundamental electroweak parameters, specifically:

• The electromagnetic coupling at the Z-boson pole, $\alpha(M_Z^2)$: This is a critical input for global electroweak fits and high-energy collider physics.
• The running of the weak mixing angle, $\sin^2\theta_W(q^2)$: Its scale dependence provides complementary tests of the SM and sensitivity to Beyond Standard Model (BSM) physics.

The dominant theoretical uncertainty in determining these quantities arises from the Hadronic Vacuum Polarization (HVP) contribution to the running of the couplings, denoted as $\Delta\alpha^{(5)}_{\text{had}}$. Traditionally, this has been calculated using dispersive approaches based on experimental $e^+e^- \to \text{hadrons}$ cross-section data. However, tensions between different experimental datasets and the need for independent, first-principles calculations motivate the use of Lattice QCD.

2. Methodology

The authors present an updated Lattice QCD calculation using $N_f = 2+1$ CLS ensembles (non-perturbatively $O(a)$-improved Wilson fermions and Lüscher-Weisz gauge action). The study covers five lattice spacings ($a \approx 0.039\text{--}0.085$ fm) and pion masses ranging from unphysical values down to the physical point.

Key methodological innovations include:

• Telescopic Window Decomposition: To manage distinct systematic errors across different Euclidean distance scales, the subtracted HVP function $\bar{\Pi}(-Q^2)$ is decomposed into three virtuality regimes:
  $$\bar{\Pi}(-Q^2) = \tilde{\Pi}(-Q^2) + \tilde{\Pi}(-Q^2/4) + \bar{\Pi}(-Q^2/16)$$
  These correspond to High-Virtuality (HV), Mid-Virtuality (MV), and Low-Virtuality (LV) terms. This separation allows for tailored extrapolation strategies:
  • HV: Dominated by short-distance physics; sensitive to discretization (cutoff) effects.
  • LV: Dominated by long-distance physics; sensitive to chiral dependence and statistical noise.
• Kernel Subtraction Strategy: To further suppress short-distance cutoff effects, the Time-Momentum Representation (TMR) kernels are modified by subtracting a term calculable in continuum perturbation theory. This cancels log-enhanced cutoff effects in the very short Euclidean-time region.
• Noise Reduction: The calculation employs Low-Mode Averaging (LMA) for light-quark connected contributions, bounding methods for long-distance tails, and spectral reconstruction techniques to improve the signal-to-noise ratio.
• Finite Volume & Isospin Breaking: Finite-volume effects are corrected using a hybrid strategy (Hansen-Patella for short distances, Meyer-Lellouch-Lüscher for long distances). Isospin-breaking (IB) effects are estimated using lattice QCD+QED and phenomenological models, found to be subleading.
• Euclidean Split Technique: To convert the space-like lattice results ($Q^2 > 0$) to the time-like Z-pole ($M_Z^2$), the authors use the Euclidean split technique. This splits the total running into:
  1. Non-perturbative lattice contribution at a matching scale $Q_0^2$.
  2. Perturbative QCD (pQCD) evolution from $Q_0^2$ to $M_Z^2$.
  3. A small perturbative continuation from $-M_Z^2$ to $+M_Z^2$.

3. Key Contributions

• Enhanced Precision: The study achieves roughly twice the precision of previous phenomenological estimates and significantly improves upon the authors' own 2022 calculation.
• Refined Analysis Strategy: The introduction of the telescopic series and new kernel functions enables a clean separation of short-distance cutoff effects from long-distance chiral dependence, allowing for reliable extrapolation to higher Euclidean momenta ($Q^2$ up to $12 \text{ GeV}^2$).
• Ab Initio Estimate: Provides a first-principles determination of $\Delta\alpha^{(5)}_{\text{had}}(M_Z^2)$ with a total relative uncertainty of approximately 1.7‰.
• Future Roadmap: The paper analyzes improvement scenarios required to meet the precision targets of next-generation colliders (e.g., FCC-ee).

4. Results

• Space-like Region ($0.25 \le Q^2 \le 12 \text{ GeV}^2$):
  • The lattice results for $\Delta\alpha^{(5)}_{\text{had}}(-Q^2)$ are consistent with the authors' previous work but show mild discrepancies with results from the BMW collaboration.
  • There is a persistent tension (systematic deviation) between these lattice determinations and dispersive evaluations based on the $R$-ratio (experimental data), with lattice results lying systematically higher.
• Z-Pole ($M_Z^2$):
  • Using the Euclidean split technique with the AdlerPy package for pQCD evolution, the authors determine $\Delta\alpha^{(5)}_{\text{had}}(M_Z^2)$.
  • The tension observed in the space-like region is significantly reduced at the Z-pole once perturbative running is included.
  • The final result is slightly larger (at the $1\text{--}2\sigma$ level) than standard dispersive evaluations but remains compatible with determinations based on CMD-3 data.
• Weak Mixing Angle: The paper also provides updated determinations for the hadronic contribution to the running of the weak mixing angle, $(\Delta \sin^2\theta_W)_{\text{had}}$.

5. Significance and Future Outlook

• Impact on Electroweak Precision: The reduction of the non-perturbative uncertainty by a factor of two is crucial for interpreting future high-precision electroweak measurements.
• Collider Requirements: Future colliders like FCC-ee aim for an absolute precision of $(3\text{--}5) \times 10^{-5}$ on $\alpha(M_Z^2)$. The current work is approaching this, but further improvements are needed.
• Optimization Strategy: The authors model the total uncertainty as a function of the matching scale $Q_0^2$ and lattice precision. They conclude that the most efficient path to the required precision ($\sim 3 \times 10^{-5}$) is not simply increasing $Q_0^2$ or improving lattice precision alone, but a combined approach:
  • Reducing lattice uncertainty by a factor of $\sim 2$.
  • Increasing the matching scale to $Q_0^2 \sim 20 \text{ GeV}^2$.
• Next Steps: Future work must focus on extending lattice calculations to higher virtualities, incorporating full isospin-breaking effects, and developing strategies to further reduce discretization effects at large momenta.

In summary, this paper represents a major step forward in the lattice determination of electroweak running couplings, offering a robust, first-principles alternative to data-driven methods and providing a clear roadmap for meeting the stringent precision requirements of future particle physics experiments.