Beyond scalar QED radiative corrections: the $\rho^{\pm}-\rho^0$ width difference, FSR corrections and their impact on $\Delta a_{\mu}^{\rm HVP, LO}[\tau]$

This paper re-evaluates the $\rho^{\pm}-\rho^0$ width difference and calculates structure-dependent final state radiation corrections by incorporating the full electromagnetic structure of charged mesons and Lorentz vertices, thereby providing crucial inputs for isospin breaking corrections in the muon $g-2$ hadronic vacuum polarization derived from $\tau$ lepton data.

The Gist

Imagine you are trying to weigh a very specific, tiny object to solve a massive mystery about the universe: Why does the muon (a heavy cousin of the electron) wiggle differently than our current laws of physics predict?

This "wiggle" is called the muon's magnetic moment (or $g-2$). Scientists have measured this wiggle with incredible precision, and it doesn't quite match the math. This suggests there are invisible particles or forces messing with the muon. To find the culprit, physicists need to calculate exactly how much "noise" the known particles create.

The biggest source of this noise comes from a cloud of virtual particles popping in and out of existence, specifically involving pions (tiny particles made of quarks). To calculate this noise, physicists usually look at data from two different sources:

1. Electron-Positron Collisions ($e^+e^-$): Smashing electrons and positrons together to create pions.
2. Tau Decays: Watching a heavy particle called a "tau" decay into pions.

Ideally, these two sources should tell the same story because of a rule called Isospin Symmetry (think of it as a cosmic rule saying "positive and negative charges should behave exactly the same way"). But in reality, they don't. There are tiny differences, like a slight weight difference between a positive pion and a neutral pion.

The Problem: The "Rho" Meson and the "Fuzzy" Math

The main character in this story is the Rho meson ($\rho$). It's a short-lived particle that decays into pions. It's the "bridge" between the tau decay data and the electron collision data.

In the past, scientists tried to calculate the tiny differences between the charged Rho ($\rho^\pm$) and the neutral Rho ($\rho^0$). They used a simplified model called Scalar QED.

• The Analogy: Imagine trying to calculate the drag on a car by pretending the car is a perfect, smooth, featureless sphere. It's easy to do the math, but it ignores the spoilers, the tires, and the engine.
• The Flaw: This "smooth sphere" model ignored the internal structure of the particles. It assumed the particles were just point-like dots with no internal parts. This led to a "fuzzy" uncertainty in the final answer, which was too big to solve the muon mystery.

The Solution: Looking Under the Hood

This paper is like the team that decided to stop pretending the car is a smooth sphere. They built a much more detailed model that accounts for the internal structure of the particles.

Here is what they did, step-by-step:

1. The "Structure-Dependent" Upgrade
Instead of treating the particles as simple dots, they used a model called Vector Meson Dominance (VMD).

• The Analogy: Imagine you are studying how a ball bounces off a wall. The old model said, "The ball is a point, so it bounces perfectly." The new model says, "Wait, the ball is made of rubber, filled with air, and has a specific texture. When it hits the wall, it squishes, vibrates, and bounces differently."
• By including this "squishiness" (the electromagnetic structure), they could calculate the radiative corrections (tiny energy adjustments due to light/photons) much more accurately.

2. The "Width" Difference
Every unstable particle has a "width," which is basically a measure of how fast it decays. The charged Rho and the neutral Rho decay at slightly different rates.

• The Old Result: Previous calculations (using the "smooth sphere" model) said the charged Rho decayed slower than the neutral one by a certain amount.
• The New Result: With the new "detailed structure" model, the authors found that the difference is actually negative. The charged Rho decays faster than the neutral one (or rather, the difference flips sign).
• Why it matters: This tiny flip changes the final calculation of the muon's wiggle significantly.

3. The "Final State Radiation" (FSR)
When particles decay, they sometimes spit out a photon (a particle of light). This is called "Final State Radiation."

• The Analogy: Imagine two runners (the pions) finishing a race. Sometimes, they trip and drop a shoe (a photon) as they cross the finish line.
• The old model calculated how often they dropped the shoe based on a simple rule. The new model realized that because the runners have "structure" (they aren't just points), the way they drop the shoe is different.
• The authors found that this "structure" reduces the amount of radiation correction by about 30% in the high-energy range.

The Big Picture Impact

Why should you care?

The authors combined their new, more accurate calculations of the Rho meson's decay width and the radiation corrections. They plugged these numbers into the formula for the Muon's Magnetic Moment.

• The Result: Their new calculation shifts the predicted value of the muon's wiggle.
• The Comparison: This new value is now closer to the results obtained from the most recent electron-positron collision experiments (specifically the CMD-3 collaboration).
• The Takeaway: By fixing the "fuzzy math" and looking at the internal structure of the particles, the gap between the "Tau data" and the "Electron data" is narrowing. This brings us closer to understanding if the discrepancy in the muon's wiggle is due to new physics (like undiscovered particles) or just imperfect calculations.

Summary in a Nutshell

Think of the universe as a giant, complex clock. The "Muon Wiggle" is a gear that isn't turning quite right.

• Old Scientists: Tried to fix the gear by guessing its shape (assuming it was a perfect circle).
• These Authors: Said, "No, let's look at the actual metal, the grooves, and the wear and tear."
• The Outcome: By measuring the actual shape, they adjusted the gear's position. It now fits much better with the other parts of the clock, suggesting that the mystery might be solved by better math rather than new, invisible gears.

This paper is a triumph of precision engineering in theoretical physics, showing that even the tiniest details of particle structure matter when you are trying to solve the universe's biggest puzzles.

Technical Summary

1. Problem Statement

The paper addresses a critical uncertainty in the calculation of the Hadronic Vacuum Polarization (HVP) contribution to the muon anomalous magnetic moment ($a_\mu$).

• Context: The dominant contribution to HVP below 2 GeV comes from the $e^+e^- \to \pi^+\pi^-$ channel. This can be estimated using $\tau$-decay data ($\tau^- \to \pi^-\pi^0\nu_\tau$) via the Conserved Vector Current (CVC) hypothesis, provided Isospin Breaking (IB) corrections are applied.
• The Gap: Current IB corrections rely heavily on the radiative corrections to the $\rho \to \pi\pi$ decay width. Previous calculations (e.g., Ref. [1]) used the Scalar QED (sQED) approximation, treating pions and rho mesons as point-like particles and considering only "convection terms" for the charged rho vertex to ensure ultraviolet (UV) finiteness.
• The Issue: The sQED approximation ignores the internal electromagnetic structure of the hadrons. This introduces a significant uncertainty (estimated at 10% in previous works) in the IB corrections, specifically regarding the width difference ($\Delta\Gamma_\rho = \Gamma_{\rho^0} - \Gamma_{\rho^\pm}$) and the Final State Radiation (FSR) factor. These uncertainties limit the precision of $a_\mu$ predictions derived from $\tau$ data.

2. Methodology

The authors extend their previous work by moving beyond the point-like sQED approximation to include the electromagnetic structure of the mesons.

• Generalized Vector Meson Dominance (GVMD):

  • Instead of point-like interactions, the photon-hadron vertices are modeled using a GVMD framework.
  • The photon propagator in loop diagrams is modified by a form factor $F(k^2)$:
    $$\frac{1}{k^2} \to \frac{1}{k^2}[F(k^2)]^2$$
    where $F(k^2) = \sum_V \frac{a_V m_V^2}{m_V^2 - k^2}$, summing over vector resonances ($\rho, \rho', \rho''$).
  • This modification renders the loop integrals for the charged rho decay ($\rho^+ \to \pi^+\pi^0$) finite without relying on the truncated convection approximation.

• Calculation of Radiative Corrections ($\delta$):

  • The total radiative correction $\delta$ is split into Virtual (one-loop) and Real (photon emission) parts.
  • Virtual Corrections: The amplitude is decomposed into an sQED part (point-like) and a structure-dependent (VMD) part. The VMD part accounts for the full electromagnetic vertex of the spin-1 $\rho$ meson.
  • Real Corrections: Calculated using Low's theorem for soft photons and including hard photon contributions. The separation scale $\omega_0$ is introduced to cancel infrared divergences between virtual and real parts.

• Width Difference Calculation:

  • The width difference is calculated using the formula:
    $$\Delta\Gamma_\rho[\pi\pi(\gamma)] = \Gamma(\rho^0 \to \pi^+\pi^-) \left[ \delta_0 - \delta_+ - 2\frac{\delta g}{g_0} - \dots \right]$$
  • The authors include contributions from phase space differences (pion mass splitting) and the calculated radiative corrections $\delta_0$ and $\delta_+$.

• Impact on $a_\mu$:

  • The new corrections are fed into the dispersive integral for the HVP contribution using $\tau$ data, specifically evaluating the shift in the isospin-breaking factor $R_{IB}(s)$.

3. Key Contributions

1. Full Electromagnetic Vertex: The paper provides the first calculation of radiative corrections for $\rho^+ \to \pi^+\pi^0$ using the full Lorentz structure of the $\rho\rho\gamma$ vertex, rather than the truncated convection term.
2. Structure-Dependent Effects: It quantifies the impact of hadronic structure (via GVMD) on both the decay width difference and the FSR corrections.
3. Validation: The sQED limit of their calculation is shown to reproduce existing results (Schwinger, Refs [18, 22]), serving as a cross-check.
4. Re-evaluation of $\Delta\Gamma_\rho$: A rigorous recalculation of the width difference between neutral and charged rho mesons, incorporating structure-dependent loop corrections.

4. Key Results

A. Radiative Corrections ($\delta$)

• Neutral Rho ($\rho^0 \to \pi^+\pi^-$): The inclusion of structure-dependent effects reduces the radiative correction by approximately 25% compared to the sQED result.
• Charged Rho ($\rho^+ \to \pi^+\pi^0$): The corrections are smaller in magnitude and flip sign compared to the sQED convection approximation.
• FSR Factor: The structure-dependent effects reduce the FSR correction ($FSR(s) - 1$) by about 30% in the region above the $\rho$ resonance ($\sqrt{s} > 1.2$ GeV).

B. Rho Width Difference ($\Delta\Gamma_\rho$)

• Previous Estimate (sQED): $\Delta\Gamma_\rho[\pi\pi(\gamma)] \approx +1.82$ MeV.
• New Estimate (GVMD): The radiative contribution drops to $\approx +0.69$ MeV.
• Total Width Difference: Including pion mass differences and other decay channels:
  $$\Delta\Gamma_\rho = (-0.525 \pm 0.220) \text{ MeV}$$
  • This result is negative, contrasting with the positive value ($+0.76$ MeV) used in previous analyses.
  • It aligns well with recent data-driven determinations (e.g., Ref. [15] giving $-0.58 \pm 1.04$ MeV).

C. Impact on Muon $g-2$ ($a_\mu$)

• The shift in the isospin-breaking corrections leads to a change in the HVP contribution derived from $\tau$ data.
• The total shift relative to previous calculations (using the same form factor model) is $+2.9 \times 10^{-10}$.
• Final Result for $a_\mu^{\text{HVP,LO}}[\tau, 2\pi]$:
  $$a_\mu^{\text{HVP,LO}}[\tau, 2\pi] = (520.2 \pm 1.9 \pm 2.2 \pm 1.7) \times 10^{-10}$$
  • This value is closer to the recent CMD-3 experimental result ($526.0 \pm 4.2$) and larger than the average of previous $e^+e^-$ based evaluations.
  • The uncertainty due to IB corrections is slightly reduced, but the central value is significantly modified.

5. Significance

• Resolution of Discrepancies: The negative width difference and reduced radiative corrections help reconcile the tension between $a_\mu$ values derived from $\tau$ data and those from $e^+e^-$ data, bringing the $\tau$-based result closer to the CMD-3 measurement.
• Reduction of Model Dependence: By moving beyond the sQED approximation, the authors reduce the theoretical uncertainty associated with "missing structure-dependent effects," which was previously a major source of error in IB corrections.
• Precision Physics: As experimental precision on $a_\mu$ improves (e.g., Fermilab, J-PARC), theoretical inputs must match. This work demonstrates that hadronic structure effects in radiative corrections are non-negligible and must be included to achieve the required precision for testing the Standard Model.
• Methodological Advance: The successful application of GVMD to render the charged rho loop finite without ad-hoc truncations provides a robust framework for future calculations of radiative corrections in vector meson decays.