Hadronic contributions to $\alpha(Q^{2})$ and $\sin^{2}\theta_{W}(Q^{2})$ from spectral reconstruction of lattice-QCD data

This paper presents preliminary lattice-QCD results using $N_f = 2+1+1$ HISQ ensembles to calculate the hadronic contributions to the running of the electromagnetic coupling and the electroweak mixing angle through a proposed spectral-reconstruction strategy.

The Gist

The Cosmic "Running" Problem: A Simple Guide to the Paper

Imagine you are trying to measure the strength of a magnet. If you hold it close, it feels incredibly strong. If you move it a few inches away, it feels much weaker. The "strength" of the magnet isn't a fixed number; it changes depending on how close you are to it.

In the world of particle physics, the fundamental forces of nature (like electromagnetism) work the same way. They have a "strength" that changes depending on how much energy you use to probe them. Physicists call this "running."

This paper is about scientists trying to calculate exactly how these strengths "run" (change) with extreme precision.

---

1. The Problem: The "Foggy" Middleman

To understand how these forces change, physicists need to look at what happens in the "empty" space between particles. According to quantum physics, "empty" space isn't actually empty; it’s bubbling with tiny, short-lived particles called hadrons.

Think of these hadrons as a thick, swirling fog that sits between two magnets. When you try to measure the magnetic force, you aren't just measuring the magnets; you are measuring the magnets through the fog.

The catch: We can’t see the fog directly. We can only see how it affects the force. Because the math governing this "fog" (called QCD) is incredibly messy and chaotic, we can't just use a simple pen-and-paper formula. We have to use supercomputers to simulate the entire universe on a tiny grid—this is called Lattice QCD.

2. The Challenge: The "Sticky" Data Problem

The researchers used a method called TMR (Time-Moment Representation). Imagine you are trying to take a video of a dancer to understand their movement.

The TMR method is like taking thousands of individual photos. However, because the photos are taken so quickly, each photo is almost identical to the one before it. In math terms, the data points are "strongly correlated."

If you try to combine all these "sticky," nearly identical photos to create a smooth, continuous movie of the dancer, the computer gets confused. It’s like trying to draw a smooth line through a series of dots that are all clumped together in a way that makes the math "explode" or become unstable.

3. The Solution: The "Spectral Reconstruction" (The Magic Lens)

To fix this, the authors propose a new strategy: Spectral Reconstruction (specifically using a method called HLT).

Instead of trying to stitch together those "sticky" individual photos to make a movie, they are doing something much smarter. They are trying to reconstruct the "DNA" of the dance itself.

Instead of looking at the dancer at every millisecond, they are looking at the underlying rhythm and energy that causes the dance. By reconstructing this "spectral density" (the underlying rhythm), they can mathematically "smooth out" the fog. This allows them to create a beautiful, continuous curve that shows exactly how the force changes, rather than just a series of disconnected, shaky dots.

4. Why does this matter? (The MUonE Connection)

Why go to all this trouble? Because we are currently in a "high-precision" era of science. We are looking for tiny cracks in our understanding of the universe—cracks that might reveal "New Physics" (things like Dark Matter).

The researchers mention a future experiment called MUonE. Think of MUonE as a high-powered flashlight that will shine through the "hadronic fog" to measure these forces experimentally.

The scientists in this paper are building the ultimate theoretical map. By having a perfect map from their supercomputer simulations, they can compare it to the "flashlight" data from the MUonE experiment. If the map and the flashlight don't match, it means we've discovered something brand new about how the universe works!

---

Summary in a Nutshell

• The Goal: Calculate how the strength of nature's forces changes at different energies.
• The Obstacle: A "fog" of particles makes the math messy, and the data points are too "sticky" to connect easily.
• The Innovation: Instead of connecting dots, they are reconstructing the "rhythm" of the particles to create a smooth, continuous prediction.
• The Big Picture: This provides a benchmark to help us find "New Physics" in upcoming global experiments.

Technical Summary

Technical Summary: Hadronic Contributions to $\alpha(Q^2)$ and $\sin^2 \theta_W(Q^2)$

This paper presents preliminary results from a lattice QCD study conducted by the Fermilab Lattice and MILC Collaborations. The research focuses on calculating the hadronic vacuum polarization (HVP) contributions to the running of two fundamental electroweak parameters: the electromagnetic coupling $\alpha(Q^2)$ and the electroweak mixing angle $\sin^2 \theta_W(Q^2)$.

1. The Problem: Precision and Correlation

The precision of Standard Model (SM) electroweak tests is currently limited by the uncertainty in the running of fundamental couplings. While leptonic contributions are well-understood via perturbation theory, the hadronic contributions ($\Delta\alpha_{\text{had}}$ and $\Delta\sin^2\theta_W$) involve low-energy, non-perturbative QCD dynamics.

Traditionally, these are estimated using experimental $e^+e^- \to \text{hadrons}$ data and dispersion relations. While lattice QCD provides a first-principles alternative, the standard Time-Momentum Representation (TMR) method faces a significant numerical hurdle: the vector correlators used in TMR are highly correlated across different momentum scales ($Q^2$). This leads to a nearly singular correlation matrix, making it mathematically unstable to perform a simultaneous continuum extrapolation for the entire running curve.

2. Methodology

The study utilizes $N_f = 2+1+1$ Highly Improved Staggered Quark (HISQ) ensembles at physical quark masses. The methodology is divided into two complementary approaches:

• Point-by-Point Continuum Extrapolation (TMR):
  The authors compute the subtracted vacuum polarization functions $\hat{\Pi}_{XX}(Q^2)$ for individual flavor components (light, strange, and charm) and disconnected terms. To reach the continuum limit, they perform independent extrapolations for each $Q^2$ scale using a Bayesian model. This model accounts for discretization effects ($a^2, a^4$ terms) and potential sea-quark mass mistuning ($\delta_s$).

• Spectral-Density Reconstruction (HLT Method):
  To overcome the correlation issues of TMR and obtain a continuous curve, the authors implement the Hansen-Lupo-Tantalo (HLT) method. Instead of calculating $\Pi(Q^2)$ directly, they reconstruct "smeared" spectral densities $\rho_\epsilon(E)$. This involves:

  1. Relating the correlator $G(t)$ to the spectral density via a Laplace transform.
  2. Approximating the target energy kernel using a linear combination of basis functions.
  3. Using a smearing width ($\epsilon$) to balance the tradeoff between approximation bias and statistical noise.
  4. Extrapolating to the continuum at fixed $\epsilon$, and subsequently taking the unsmeared limit ($\epsilon \to 0$).

3. Key Contributions

• Flavor Decomposition: The work provides a detailed breakdown of contributions from light ($l$), strange ($s$), and charm ($c$) quarks, as well as flavor-singlet disconnected terms ($D_{ls}$), allowing for the isolation of systematic effects per flavor.
• Dual-Method Framework: The paper establishes a workflow that combines the stability of point-by-point TMR with the continuity of the HLT spectral reconstruction.
• Advanced Statistical Modeling: The use of Bayesian inference (via the lsqfit package) to handle discretization artifacts and mass mistuning provides a rigorous way to quantify systematic uncertainties.

4. Preliminary Results

• Continuum Extrapolation: Figure 1 demonstrates successful point-by-point extrapolations for light and strange quark contributions, showing stability when testing different fit forms (including $a^4$ terms).
• Running Curves: Figure 2 presents the preliminary, blinded running of $\Delta\alpha(Q^2)$ and $\Delta\sin^2\theta_W(Q^2)$ for $0 < Q^2 < 7 \text{ GeV}^2$. The results show high statistical precision, with lattice data points closely following the continuum-extrapolated curves.
• Spectral Reconstruction: Preliminary HLT reconstruction (Figure 3) confirms the ability to extract smeared spectral densities, providing a path toward a continuous, stable determination of the vacuum polarization.

5. Significance and Outlook

This research is highly significant for the high-precision frontier of particle physics. By providing a first-principles lattice QCD determination of these couplings, the work offers a vital independent check on phenomenological estimates.

Furthermore, the continuous determination of $\Delta\alpha(Q^2)$ via spectral reconstruction serves as a strategic benchmark for the upcoming MUonE experiment, which aims to measure these values in the space-like region. Once the full systematic error budget (finite-volume, taste-breaking, and disconnected contributions) is finalized, these results will provide a critical test of the consistency of the Standard Model at the $M_Z$ scale.