Improved calculation of radiative corrections to $\tau\to\pi\pi\nu_\tau$ decays

This Letter presents an improved, model-independent calculation of radiative corrections to $\tau\to\pi\pi\nu_\tau$ decays that incorporates structure-dependent effects beyond point-like pions and threshold enhancements, thereby refining the precision of the hadronic-vacuum-polarization contribution to the muon's anomalous magnetic moment derived from tau decay data.

The Gist

Here is an explanation of the paper using simple language and everyday analogies.

The Big Picture: The "Mystery of the Wobbly Muon"

Imagine the universe is a giant,精密 (precise) clockwork machine. Physicists have a manual (the Standard Model) that tells them exactly how every gear should turn. One of the most sensitive gears is a particle called the muon.

Scientists have measured how this muon wobbles (its "magnetic moment") with incredible precision. They found that the real-world wobble is slightly different from what the manual predicts. This suggests there might be a hidden gear in the machine that we haven't found yet—perhaps a new particle or a new force of nature.

However, to be sure, we need to make sure our calculation of the "known" gears is perfect. The biggest source of uncertainty in the manual comes from a messy, chaotic part of the machine called the Hadronic Vacuum Polarization. Think of this as a cloud of virtual particles (like pions) that pop in and out of existence, briefly interfering with the muon's wobble.

The Problem: Two Different Maps

To calculate the size of this "cloud," scientists usually look at data from electron-positron collisions ($e^+e^-$). But recently, a new experiment (CMD-3) measured this cloud and found a result that clashes with previous data. It's like two GPS apps giving you different routes to the same destination, and you don't know which one is right.

So, scientists decided to try a different map: Tau decays.
Instead of looking at electrons colliding, they looked at Tau particles (heavy cousins of the electron) decaying into pions. This is a complementary route. If both maps lead to the same destination, we can be confident in our result. If they differ, we know there's a problem.

The Obstacle: The "Static" in the Signal

Here is the catch: The Tau decay map is covered in "static."
When a Tau particle decays, it doesn't just spit out two pions cleanly. It also emits photons (light particles) due to electromagnetic effects. These photons mess up the measurement, just like static on a radio distorts a song.

To use the Tau data to fix the muon mystery, scientists must mathematically remove this static (radiative corrections) to see the pure signal underneath. Until now, the methods used to remove this static were a bit like using a blunt knife to cut a diamond—they worked, but they weren't precise enough, especially near the "resonance" (a specific energy level where the pions behave like a vibrating guitar string, the $\rho$-meson).

The Solution: A New, Sharper Knife

This paper presents a new, improved way to calculate and remove that static. The authors (Colangelo, Cottini, Hoferichter, and Holz) did three main things:

1. They stopped treating pions like billiard balls.

   • Old way: They assumed pions were simple, point-like dots.
   • New way: They realized pions have an internal structure (they are made of quarks). They used a "dispersive" method, which is like listening to the entire song of the pion's behavior rather than just a single note. This revealed that the "static" is actually much louder near the $\rho$-resonance than previously thought.

2. They fixed the "Threshold" problem.

   • Near the very bottom of the energy scale (where two pions just barely have enough energy to exist), the math gets very unstable, like trying to balance a pencil on its tip. The authors developed a new mathematical trick (changing the variables) to keep the calculation stable all the way down to the very edge.

3. They checked the data.

   • They took all the available experimental data from major labs (Belle, ALEPH, CLEO, OPAL) and fitted their new model to it. They found that their new model fits the data very well, but it also highlighted some small tensions between different experiments, suggesting we need even better data in the future (perhaps from the Belle II experiment).

The Result: A Shift in the Numbers

When they applied this new, sharper calculation to the muon mystery:

• The correction they needed to apply to the Tau data changed significantly.
• The "cloud" of virtual particles contributes less to the muon's wobble than some previous estimates suggested.
• This shifts the theoretical prediction for the muon's wobble by about 2.5 standard deviations.

What does this mean?
It doesn't solve the mystery yet, but it tightens the noose. The gap between the "Standard Model" prediction and the "Real World" measurement has shifted. The uncertainty in the calculation has been reduced by a factor of three.

The Bottom Line

Think of this paper as upgrading the noise-canceling headphones for the muon experiment.

• Before: The headphones removed most of the noise, but left a hum that made it hard to hear the music clearly.
• Now: The new headphones (the improved calculation) remove the noise much more effectively, especially the tricky, loud parts.

This allows physicists to hear the "music" of the muon more clearly. While the music still sounds slightly different from the manual's prediction (hinting at new physics), we are now much more confident that the difference is real and not just a glitch in our calculation. The paper calls for even more precise measurements in the future to finally solve the mystery.

Technical Summary

Here is a detailed technical summary of the paper "Improved calculation of radiative corrections to $\tau \to \pi\pi\nu_\tau$ decays" by Colangelo et al.

1. Problem Statement

The anomalous magnetic moment of the muon, $a_\mu$, is a critical test of the Standard Model (SM). Recent experimental results from Fermilab have reduced the experimental uncertainty significantly, creating a tension with the SM prediction. The dominant source of theoretical uncertainty in the SM prediction is the leading-order hadronic vacuum polarization (HVP) contribution, $a_\mu^{\text{HVP, LO}}$.

Traditionally, this contribution is evaluated using data from $e^+e^- \to \text{hadrons}$ cross-sections. However, recent tensions in the $e^+e^- \to \pi^+\pi^-$ data (specifically from the CMD-3 experiment) have cast doubt on the reliability of the $e^+e^-$-based HVP evaluation. Consequently, there is a renewed interest in using $\tau \to \pi\pi\nu_\tau$ decays as a complementary data source.

To utilize $\tau$ decay data for $a_\mu$, one must perform an isospin rotation from the charged channel ($\pi^\pm\pi^0$) to the neutral channel ($\pi^+\pi^-$). This process requires the removal of $\tau$-specific radiative corrections and the application of isospin-breaking (IB) corrections. The paper identifies that the current calculation of these radiative corrections is insufficiently precise, relying heavily on Chiral Perturbation Theory (ChPT) combined with resonance models, which makes robust uncertainty quantification difficult. Specifically, previous calculations neglected structure-dependent (SD) virtual corrections arising from the internal structure of the pions, which are significant near the $\rho(770)$ resonance.

2. Methodology

The authors present a model-independent, dispersive approach to calculate the electromagnetic correction factor $G_{\text{EM}}(s)$, which encapsulates the radiative corrections to the differential decay rate.

A. Virtual Corrections (Dispersive Approach)

• Beyond Point-Like Pions: Unlike previous ChPT-based calculations that treated pions as point-like particles at low orders, this work incorporates the internal structure of the pion via a dispersive representation of the pion form factor, $f_+(s)$.
• Box Diagram Reinterpretation: The virtual corrections are dominated by the "box diagram" (Fig. 2 in the paper). Instead of calculating this purely within ChPT, the authors reinterpret it as a unitarity diagram involving the pion form factor.
• Calculation: They insert a once-subtracted dispersion relation for $f_+(s)$ into the loop calculation. This reduces the problem to standard one-loop Passarino–Veltman functions ($B_0, C_0, D_0$).
  • This ensures the result is manifestly UV finite.
  • IR singularities and chiral logarithms are verified to match ChPT limits.
• Matching to ChPT: To account for ChPT diagrams not included in the dispersive box diagram (e.g., contact terms), the authors match the dispersive result to ChPT at $s=t=0$. This suppresses sensitivity to the high-energy asymptotic behavior of the form factor.
• Scheme Independence: The calculation addresses the scheme ambiguity in the electroweak correction factor $S_{\pi\pi}^{\text{EW}}$ by ensuring the product $S_{\pi\pi}^{\text{EW}} G_{\text{EM}}(s)$ is scale and scheme independent.

B. Real Corrections

• Threshold Behavior: The real photon emission corrections exhibit a strong threshold enhancement ($\propto 1/(s-4M_\pi^2)$), making numerical evaluation unstable near the two-pion threshold.
• Stable Implementation: The authors developed a method using specific angular integration variables to yield stable expressions for the amplitude times the Jacobian, allowing for robust numerical evaluation down to the threshold.
• Resonance Contributions: The calculation includes contributions from vector ($\rho$) and axial-vector ($a_1$) meson exchanges, as well as the Wess–Zumino–Witten (WZW) anomaly.
• Coupling Constants: Due to uncertainties in the axial-vector coupling $F_A$, the authors define their central values as the average of fits using two different strategies for determining couplings: Short-Distance (SD) constraints and data-driven determinations. The difference is assigned as a $1\sigma$ uncertainty.

C. Global Fits to Data

• Iterative Procedure: To ensure self-consistency, the authors employ an iterative procedure:
  1. Start with an Omnès function for the pion form factor $f_+(s)$.
  2. Calculate a first approximation of $G_{\text{EM}}(s)$.
  3. Fit $f_+(s)$ to experimental $\tau$ spectral function data (Belle, ALEPH, CLEO, OPAL) using the calculated $G_{\text{EM}}(s)$.
  4. Recalculate $G_{\text{EM}}(s)$ with the improved form factor and repeat until convergence.
• Fit Strategy: Fits were performed with varying numbers of parameters ($N=3, 4, 5$) to describe inelastic effects and excited resonances ($\rho', \rho''$). $N=4$ was chosen as the central result to balance fit quality and avoid overfitting.

3. Key Contributions

1. First Inclusion of Structure-Dependent Virtual Corrections: This is the first calculation of radiative corrections to $\tau \to \pi\pi\nu_\tau$ that includes the effects of the pion's internal structure (via the dispersive form factor) in the loop diagrams.
2. Dispersive Framework: The shift from purely model-dependent resonance models to a dispersive representation allows for a more rigorous handling of unitarity and analyticity, particularly near the $\rho(770)$ resonance.
3. Stable Threshold Implementation: The development of a stable numerical method for real corrections near the two-pion threshold resolves previous issues with threshold singularities.
4. Quantification of Isospin-Breaking: The work provides a precise, data-driven assessment of the isospin-breaking corrections required to convert $\tau$ data into an $a_\mu$ contribution.

4. Results

• Correction Factor $G_{\text{EM}}(s)$: The new calculation shows significant deviations from previous results (e.g., Flores-Baéz et al., Miranda & Roig) in the vicinity of the $\rho(770)$ resonance due to the structure-dependent virtual corrections.
• Impact on $a_\mu$:
  • The total radiative correction to the $\tau$-based HVP contribution is calculated as:
    $$\Delta a_\mu^{\text{HVP, LO}}[\pi\pi, \tau]_{\text{Full}} = -24.8(1.4) \times 10^{-10}$$
  • This represents a shift of approximately 2.5$\sigma$ (about 3 units in $10^{-10}$) compared to previous evaluations (e.g., Ref. [9]).
  • The uncertainty due to $G_{\text{EM}}(s)$ has been reduced by nearly a factor of three compared to previous estimates.
• Data Tensions: The global fits reveal tensions between different experimental data sets (Belle vs. ALEPH/CLEO/OPAL) and a moderate tension between the low-energy spectrum and the $\rho$ region, driven by constraints from analyticity and unitarity.

5. Significance

• Resolving the $a_\mu$ Tension: By providing a more reliable and precise calculation of radiative corrections, this work strengthens the case for using $\tau$ decay data as a cross-check for the $e^+e^-$-based HVP evaluation. This is crucial given the current discrepancies in $e^+e^-$ data.
• Reduced Theoretical Uncertainty: The reduction in the uncertainty of the radiative corrections (by a factor of three) makes the $\tau$-based determination of $a_\mu^{\text{HVP}}$ competitive with lattice QCD results and $e^+e^-$ data.
• Future Directions: The authors highlight that while their theoretical framework is improved, the current experimental data sets show tensions. They strongly motivate new, high-statistics measurements of the $\tau \to \pi\pi\nu_\tau$ spectral function at Belle II to resolve these tensions and further reduce the uncertainty in the Standard Model prediction for the muon $g-2$.

In summary, this paper represents a major step forward in the theoretical precision of $\tau$-decay based HVP calculations, utilizing modern dispersive techniques to account for pion structure and providing a robust framework for future comparisons with muon $g-2$ experiments.