Lattice Calculation of Light Meson Radiative Leptonic Decays

Using $N_f=2+1$ domain wall fermion ensembles and the infinite volume reconstruction method, this study calculates the branching ratios and form factors of radiative leptonic decays for pions and kaons, finding agreement with PIBETA and KLOE data for electron modes while confirming existing tensions between lattice results and specific experimental measurements for muon modes.

The Gist

Imagine the subatomic world as a bustling city where tiny particles called mesons (specifically pions and kaons) are like delivery trucks. Usually, these trucks deliver their cargo (a lepton and a neutrino) and disappear. But sometimes, in a rare event, the truck drops a package and accidentally sparks a tiny flash of light (a photon) as it goes. This is called a radiative leptonic decay.

Scientists want to understand exactly how these trucks are built inside. To do this, they need to measure how often these "spark-and-drop" events happen and what the light looks like. This paper is a report from a team of physicists who used a super-powerful digital simulation (called Lattice QCD) to calculate these events from first principles, essentially building the truck from scratch in a computer to see how it behaves.

Here is a breakdown of their journey, using simple analogies:

1. The Problem: The "Room Size" Limit

Imagine trying to study how a sound wave travels across a vast ocean, but you are forced to do it inside a small, tiled bathtub. In the bathtub, the waves bounce off the walls and create weird echoes that don't exist in the real ocean. This is the main problem with simulating particle physics on a computer: the "universe" of the simulation is a tiny box (the lattice).

The authors used a clever trick called Infinite-Volume Reconstruction (IVR). Think of this as a magical mirror that takes the data from the small bathtub and mathematically "unfolds" it to look like the vast ocean. This allowed them to remove the "echoes" (artifacts) caused by the small size of their computer simulation, giving them a clear picture of how the particles behave in the real, infinite world.

2. The "Electron vs. Muon" Difference

The team studied two types of delivery trucks:

• The Electron Truck: The electron is very light, like a feather.
• The Muon Truck: The muon is heavier, like a bowling ball.

The Feather Problem: When the light electron truck drops its package, it is so sensitive that it gets "jittery." It tends to emit extra, invisible sparks (photons) that are hard to see but change the math significantly. The paper explains that for the electron, these extra sparks create a massive "magnifying glass" effect (mathematically called a large logarithmic factor). If you ignore these extra sparks, your calculation is off by about 10%. That's a huge error in the world of particle physics.

The Bowling Ball: The muon is heavy and stable. It doesn't get jittery. For the muon truck, these extra sparks are negligible, so the math is much simpler.

3. The Results: Solving the Mystery

The team compared their computer-generated numbers with real-world experiments conducted by groups like PIBETA, KLOE, and E36.

• The Pion (π) Mystery: Previous computer simulations for the pion truck didn't match the real-world PIBETA experiment. The numbers were too high. However, once this team added the "jittery spark" corrections (the 10% fix mentioned above), their numbers perfectly matched the real experiment. It turns out the old simulations just forgot to account for the electron's jitter.
• The Kaon (K) Mystery: For the kaon truck, things are a bit more complicated.
  • KLOE vs. E36: Two different real-world experiments (KLOE and E36) got different results for the kaon. The authors suggest this is because the two experiments had different rules for what counts as a "spark." One experiment ignored extra sparks, while the other counted them. When the team applied the correct math for each experiment's specific rules, their results aligned with KLOE but showed a slight tension (a 1.7σ difference) with E36.
  • The Angle Issue: For the muon version of the kaon decay, the team confirmed a previous finding: when the muon and the photon fly off at wide angles, the computer predictions disagree with the ISTRA and OKA experiments. This suggests there might be something about the "truck's" internal structure that we still don't fully understand.

4. The "Blueprints" (Form Factors)

Beyond just counting how often the decay happens, the team mapped out the "blueprints" of the mesons. They calculated Form Factors, which are like a 3D map showing how the electric charge is distributed inside the meson.

• They found that for the pion, the map is fairly smooth and predictable.
• For the kaon, the map shows a slight "bump" or curve, suggesting the presence of internal resonances (like a hidden gear inside the truck) that makes it behave slightly differently than the simplest theories predicted.

Summary

In short, this paper is a high-precision engineering report. The team built a better "mathematical mirror" (IVR) to simulate particle decays without the distortion of a small computer box. They discovered that for the lightest particles (electrons), you must account for a specific type of "static electricity" (collinear radiation) to get the right answer. Once they did, their computer models finally agreed with the real-world data for pions, and provided a new, detailed explanation for the mixed results seen in kaon experiments. This work helps physicists refine the "Standard Model" of the universe, ensuring our understanding of how matter is built is as accurate as possible.

Technical Summary

1. Problem Statement

The paper addresses the theoretical calculation of radiative leptonic decays of light pseudoscalar mesons ($P \to \ell \nu_\ell \gamma$, where $P = \pi, K$ and $\ell = e, \mu$). These processes are critical for:

• Testing the Standard Model (SM): They probe hadronic structure and allow for constraints on new physics (e.g., tensor interactions).
• CKM Unitarity: Precise determination of CKM matrix elements $V_{ud}$ and $V_{us}$ is required to test the unitarity of the first row ($|V_{ud}|^2 + |V_{us}|^2 + |V_{ub}|^2 = 1$). Current determinations show a $\sim 2.4\sigma$ tension with unitarity.
• Radiative Corrections: To extract precise CKM elements from leptonic decays, one must account for $O(\alpha)$ radiative corrections. According to the Bloch-Nordsieck theorem, this requires summing virtual photon loops and real photon emissions.
• Existing Discrepancies: Previous lattice QCD calculations (notably by the Rome-Southampton collaboration) showed tensions with experimental data (PIBETA, KLOE, E36, ISTRA, OKA, E787) in specific kinematic regions. Additionally, for electron channels, $O(\alpha^2)$ collinear radiative corrections enhanced by large logarithms ($\ln(m_P^2/m_e^2)$) were often neglected, potentially explaining some discrepancies.

2. Methodology

The authors employ Lattice QCD using $N_f = 2+1$ domain wall fermion ensembles generated by the RBC and UKQCD collaborations at the physical pion mass. The study introduces several key methodological advancements:

• Infinite-Volume Reconstruction (IVR):

  • Standard lattice calculations suffer from finite-volume effects and discrete momentum constraints.
  • The authors use the IVR method to reconstruct infinite-volume hadronic matrix elements from finite-volume lattice data.
  • This involves separating the calculation into a short-range contribution (computed directly on the lattice) and a long-range contribution (reconstructed using ground-state dominance and exponential ansatz).
  • Finite-volume effects are corrected using a point-particle approximation and structure-dependent models, reducing power-law finite-volume errors to exponentially suppressed ones.

• Scalar Function Method:

  • To handle the hadronic matrix element $H_{\mu\nu}^M$, the authors decompose it into six coordinate-space scalar functions ($I_1$ to $I_6$).
  • These are transformed into momentum-space scalar functions ($\tilde{I}_i$) using weight functions involving spherical Bessel functions. This allows the extraction of form factors at arbitrary momenta, bypassing the need for twisted boundary conditions.

• Radiative Corrections ($O(\alpha^2)$):

  • A major focus is the inclusion of collinear radiative corrections for electron channels ($P \to e \nu_e \gamma$).
  • These corrections are enhanced by large logarithmic factors: $\ln(m_P^2/m_e^2)$ and $\ln(2E_{\gamma, \max}^*/m_P)$.
  • The authors implement corrections for both fully inclusive processes (PIBETA, E36) and processes with specific experimental cuts (KLOE), accounting for the angle-independent laboratory-frame energy cutoff on the second photon.

• Lattice Setup:

  • Two ensembles (48I and 64I) with physical pion masses but different lattice spacings ($a^{-1} \approx 1.73$ GeV and $2.36$ GeV) are used for continuum extrapolation.
  • A partially quenched ensemble (64Ipq) is utilized to correct for kaon mass mismatches between ensembles.

3. Key Contributions

1. First IVR Application to Radiative Decays: This work initiates the application of the IVR method to $P \to \ell \nu_\ell \gamma$, serving as a stepping stone toward a complete $O(\alpha)$ radiative correction calculation (including virtual loops).
2. Resolution of $\pi \to e \nu_e \gamma$ Discrepancy: The authors demonstrate that the previously observed tension between lattice results and PIBETA experimental data is resolved once $O(\alpha^2)$ collinear radiative corrections are included.
3. Clarification of $K \to e \nu_e \gamma$ Tensions: The study highlights that the $4\sigma$ discrepancy between KLOE and E36 measurements is largely due to their different treatments of additional inner-bremsstrahlung photons. The lattice results, when corrected for these specific experimental cuts, align with KLOE but show a $1.7\sigma$ tension with E36.
4. Confirmation of $K \to \mu \nu_\mu \gamma$ Discrepancies: For muon channels (where collinear corrections are negligible), the results confirm existing tensions between lattice predictions and ISTRA/OKA (at large photon energies) and E787 (at large muon-photon angles).
5. Form Factor Determination: The paper provides high-precision determinations of the vector ($F_V$) and axial-vector ($F_A$) form factors across the full kinematic phase space, confirming their linear behavior (with mild deviations for $K$ $F_V$) and consistency with Chiral Perturbation Theory ($O(p^4)$).

4. Key Results

• Branching Ratios (Normalized):
  • $\pi \to e \nu_e \gamma$: Lattice results with radiative corrections agree with PIBETA data within statistical errors. Without corrections, a discrepancy existed.
  • $K \to e \nu_e \gamma$:
    • Inclusive result (Region 1-5): $16.9(1.3) \times 10^{-6}$.
    • With KLOE-like cut ($E_{\gamma2}^{lab} < 20$ MeV): $16.1(1.3) \times 10^{-6}$.
    • Comparison: The KLOE-like result is consistent with KLOE data ($14.83(67) \times 10^{-6}$) within $1\sigma$. The inclusive result shows a $1.7\sigma$ tension with the E36 measurement ($19.8(1.1) \times 10^{-6}$).
  • $K \to \mu \nu_\mu \gamma$: Results confirm previous lattice findings and experimental tensions at large angles and photon energies.
• Form Factors:
  • $F_V$ and $F_A$ were extracted for 128 uniformly spaced values of $x_\gamma$ (photon energy fraction).
  • Results are statistically consistent with previous Rome-Southampton lattice calculations.
  • $F_V$ for the kaon shows a mild deviation from linearity, suggesting pole-like behavior from low-lying resonances.
• Finite Volume Effects:
  • Corrections were found to be $\sim 1\%$ for pion decays and negligible for kaon decays in the studied regions, validating the IVR approach.

5. Significance

• Precision Physics: The work significantly improves the theoretical precision required for extracting $V_{us}$ from leptonic decays, which is crucial for resolving the CKM unitarity tension.
• Methodological Validation: It validates the IVR method as a robust tool for handling real-photon emissions in lattice QCD, overcoming the limitations of finite volume and discrete momenta.
• Experimental Guidance: By explicitly modeling different experimental cuts and radiative correction schemes, the paper provides a framework for interpreting current and future data from experiments like NA62 and KLOE-2.
• Future Outlook: The study sets the stage for the complete calculation of radiative corrections (including virtual loops) using the IVR method, which will further reduce theoretical uncertainties in the Standard Model tests.

In summary, this paper resolves long-standing discrepancies in pion radiative decays through rigorous inclusion of higher-order radiative corrections and validates the lattice approach for kaon decays, while highlighting remaining tensions that may point to new physics or require further experimental/theoretical refinement.