Radiative Corrections in Bound States: Recent Results

Technical Summary: Radiative Corrections in Bound States

Problem Statement
This work addresses two distinct problems regarding radiative corrections in bound-state quantum systems where previous theoretical estimates or calculations contained significant discrepancies:

1. Muon Decay in Orbit (DIO): A long-standing disagreement exists between analytical perturbative calculations (accurate to order $(\alpha Z)^2$) and fully relativistic numerical results for the decay rate of a muon bound to light nuclei. Specifically, for oxygen ($Z=8$), a fully relativistic partial-wave expansion reported a decay rate ratio $\Gamma/\Gamma_0 \approx 0.994$, whereas the analytical prediction suggested $\Gamma/\Gamma_0 \approx 0.9981$. The discrepancy ($\sim 4 \times 10^{-3}$) was too large to be explained by the expected next-to-leading order corrections of $(\alpha Z)^3$.
2. Parapositronium Decay: The decay of parapositronium (pPs, the spin-0 ground state of positronium) into three photons ($pPs \to 3\gamma$) is forbidden in pure Quantum Electrodynamics (QED) by charge-conjugation ($C$) symmetry but is enabled by the $C$-violating weak interaction. Previous estimates of this decay rate were significantly higher than recent explicit calculations including $Z$-boson contributions, suggesting a need to understand the suppression mechanisms, particularly the vanishing of specific diagram classes.

Methodology
The authors employ two distinct theoretical frameworks corresponding to the two problems:

• Bound-Muon Decay:

  • The decay process is factorized into a cascade: $(Z\mu^-) \to [Ze^-]A$, followed by the decay of a fictitious massive neutral vector boson $A$ (carrying the neutrino flavor) into $\bar{\nu}_e \nu_\mu$.
  • The calculation utilizes a fully relativistic approach based on the Dirac equation. The initial muon state is described by the closed-form Dirac–Coulomb ground state, while the final-state continuum electron is obtained by numerically integrating the Dirac equation using the RADIAL code.
  • The decay rate is expressed as a sum over partial waves (characterized by the quantum number $\kappa_e$). The authors perform high-precision numerical integration of radial integrals using the adaptive DQAGS routine from QUADPACK with a tolerance of $10^{-9}$.
  • Convergence is achieved by extending the partial-wave sum until term-to-term changes reach $10^{-5}$ to $10^{-6}$, with the remaining tail estimated geometrically.

• Parapositronium Decay ($pPs \to 3\gamma$):

  • The analysis focuses on the one-loop amplitude involving $W$ and $Z$ bosons.
  • Calculations are performed in the unitary gauge to explicitly demonstrate the vanishing of specific contributions.
  • The authors utilize Ward identities derived from the Feynman rules for $ZWW$ and $ZAW$ vertices to show that the contraction of the $Z$-boson momentum with the amplitude leads to a telescopic cancellation of terms around the closed $W$-loop.

Key Contributions and Results

1. Resolution of the DIO Discrepancy:

   • The authors calculated the decay rate ratio $\Gamma/\Gamma_0$ for light nuclei with atomic numbers $4 \le Z \le 9$.
   • Result: The converged numerical results (e.g., for $Z=8$, $\Gamma/\Gamma_0 = 0.9976$) align with the analytical perturbative prediction of Ref. [7] ($0.9981$) when higher-order terms are accounted for, rather than contradicting it.
   • Diagnosis: The discrepancy with the previous numerical result (Ref. [8]) was traced to a premature truncation of the slowly convergent partial-wave series. Ref. [8] terminated the expansion at $L=31$ (yielding $0.994$), which corresponded to an unconverged partial sum in the current work. Only by extending the summation to $\kappa_{max} = 59$ and extrapolating the tail was the correct, converged value obtained.
   • Nuclear Effects: Comparisons with Ref. [19], which included finite nuclear charge distributions, confirm that for light nuclei ($Z \le 9$), the effect of the finite nuclear size on the decay rate is negligible.

2. Suppression of Parapositronium Decay:

   • The paper confirms that the rate for $pPs \to 3\gamma$ is many orders of magnitude smaller than early estimates.
   • Mechanism: The authors provide a detailed proof of why certain classes of diagrams vanish exactly. Specifically, they demonstrate that the amplitude for the process mediated by a virtual $Z$ boson decaying into three photons via a $W$-loop vanishes.
   • Proof: By analyzing the amplitude structure $\mathcal{M} \propto q_\mu \mathcal{M}^\mu(Z^* \to 3\gamma)$, where $q$ is the $Z$-boson momentum, they show that the contraction with $q$ results in a telescopic sum of differences between adjacent propagator insertions. Due to the closed-loop nature of the diagrams, these differences sum to zero, proving the vanishing of the complete amplitude for this topology.

Significance
The paper claims significance in two primary areas:

• Precision Physics for Muonic Atoms: By resolving the discrepancy between analytical and numerical DIO rates, the work establishes a reliable theoretical baseline for experiments searching for charged-lepton-flavor-violating (CLFV) processes (e.g., $\mu^- N \to e^- N$). Accurate knowledge of the irreducible background from ordinary muon decay is critical for these searches.
• Understanding Weak Corrections in Bound States: The detailed analysis of the $pPs \to 3\gamma$ decay clarifies the suppression mechanisms of $C$-violating weak interactions in bound states. It corrects previous overestimates of the decay rate by demonstrating that specific diagram topologies vanish exactly due to gauge invariance and symmetry properties, rather than merely being suppressed by small mass ratios.

The work does not propose new experimental setups but serves to refine the theoretical inputs required for interpreting existing and future precision measurements in muonic atoms and positronium physics.