⚛️ phenomenology

Next-to-leading order FsQED corrections to radiative pion pair production

This paper computes next-to-leading order FsQED corrections to the radiative pion pair production process $e^+ e^- \to \pi^+ \pi^- \gamma$ , compares these results with alternative theoretical approaches and previous models, and implements the findings along with additional structure-dependent mechanisms into the BabaYaga@NLO Monte Carlo generator to improve the precision of radiative return measurements at flavor factories.

Original authors: Carlo M. Carloni Calame, Marco Ghilardi, Andrea Gurgone, Guido Montagna, Mauro Moretti, Oreste Nicrosini, Fulvio Piccinini, Francesco P. Ucci

Published 2026-07-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Carlo M. Carloni Calame, Marco Ghilardi, Andrea Gurgone, Guido Montagna, Mauro Moretti, Oreste Nicrosini, Fulvio Piccinini, Francesco P. Ucci

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, bustling dance floor where tiny particles are the dancers. Sometimes, these dancers, like electrons and positrons (the antimatter twins of electrons), crash into each other and create new pairs of dancers, such as pions. Physicists love watching these crashes because the way the dancers move tells them secrets about the invisible forces holding the universe together. One of the biggest mysteries they are trying to solve is why a specific particle, the muon (a heavier cousin of the electron), wobbles a little differently than our best theories predict. To solve this, they need to calculate a "correction" to the muon's wobble caused by the vacuum of space itself acting like a foggy mirror. This calculation depends heavily on knowing exactly how pions behave when they are created.

To get this information, scientists use "radiative return" experiments. Think of this like a game of billiards where a fast-moving cue ball hits a target, but instead of just hitting it, the cue ball also shoots off a tiny spark of light (a photon) right at the moment of impact. This spark slows the cue ball down just enough so that the energy of the crash matches the perfect speed to create a pair of pions. By measuring the energy of that spark and the pions, scientists can reconstruct the "dance steps" (the form factor) of the pions. However, calculating these steps is incredibly tricky because pions aren't simple, solid balls; they are fuzzy, composite clouds made of even smaller quarks. For a long time, physicists have used a simplified map to draw these clouds, assuming they were perfectly smooth. But this map might be missing some of the fuzzy details, and those missing details could be the key to solving the muon mystery.

This paper is like a team of cartographers who decided to redraw that map with much higher resolution. The authors, a group of theoretical physicists, have calculated the next level of detail—what they call "next-to-leading order" corrections—for the radiative return process. Instead of treating the pions as simple, featureless dots, they used a sophisticated method called FsQED (Form-factor Scalar Quantum Electrodynamics) to weave the actual, complex structure of the pion directly into the math of the calculation. They wanted to see if this new, more detailed map changed the results compared to the old, simplified one.

What they found is fascinating. When they looked at the total number of pions produced (the "invariant mass"), the new, detailed map changed the prediction by only a tiny amount—about a few parts in a thousand. It's like realizing that a mountain is actually a few meters taller than you thought; it's a precise correction, but it doesn't change the fact that it's a mountain. However, when they looked at the direction the pions were flying (the "scattering angle") and the difference between pions flying forward versus backward, the new map made a much bigger difference, changing the predictions by about one percent. This is a significant shift, like realizing that while the mountain is the same height, the wind blowing around it is actually much stronger than we guessed.

The team also compared their new, detailed map against two other ways of drawing the same picture: a very simple "smooth ball" model and another complex model called GVMD. Surprisingly, their new FsQED results agreed almost perfectly with the GVMD model, even though the two methods use completely different mathematical tricks. This gives scientists a lot of confidence that they are on the right track. They also checked if there were any other hidden effects, like "neutral pseudoscalar poles" (which are like ghostly, invisible intermediate steps in the dance), but found that these contribute less than a part in a thousand, meaning they are too small to worry about for now.

Finally, the authors looked at a specific energy range around a particle called the $\phi$ -meson, where things get a bit chaotic. They found that in this specific zone, other processes—like the $\phi$ -meson decaying directly into pions and a photon—become important and can't be ignored. These extra effects can change the shape of the data by several percent, acting like a loud drumbeat that drowns out the subtle rhythm of the main dance.

All of these new calculations and insights have been baked into a powerful computer program called BabaYaga@NLO. This tool is now ready for experimentalists at "flavor factories" (massive particle colliders) to use. By using this updated software, scientists can simulate these particle crashes with much higher precision, ensuring that when they finally measure the muon's wobble, they aren't tripping over a tiny, fuzzy detail in the pion's dance steps. The paper doesn't claim to have solved the muon mystery yet, but it has provided a much sharper pair of glasses for the scientists trying to solve it.

Technical Summary: Next-to-leading order FsQED corrections to radiative pion pair production

Problem Statement
The hadronic cross section in $e^+e^-$ annihilation at low energies is a critical input for the data-driven computation of the hadronic vacuum polarisation (HVP) contribution to the muon anomalous magnetic moment ( $a_\mu$ ). The dominant channel, $e^+e^- \to \pi^+\pi^-$ , currently represents the primary source of theoretical uncertainty. While recent tensions between data-driven predictions and experimental world averages for $a_\mu$ have spurred high-precision measurements via energy scans and radiative return methods, theoretical tools must evolve to control systematic uncertainties at the sub-percent level.

A specific challenge arises in Next-to-Leading Order (NLO) calculations involving pions. The standard "naive factorised scalar QED" (F $\times$ sQED) approach, which multiplies point-like scalar QED amplitudes by an overall pion form factor, fails to fully capture the composite nature of the pion in loop integrals. Previous work established that for the energy scan process ( $e^+e^- \to \pi^+\pi^-$ ), structure-dependent corrections—where the pion form factor is inserted analytically into loop momenta—are essential for agreement with data, particularly for forward-backward asymmetries. However, a corresponding calculation for the radiative return process ( $e^+e^- \to \pi^+\pi^-\gamma$ ) using the rigorous FsQED (Form-factor scalar QED) approach was missing from the literature.

Methodology
This work computes the NLO structure-dependent corrections for $e^+e^- \to \pi^+\pi^-\gamma$ using the FsQED approach. The methodology involves:

FsQED Framework: The authors embed the pion form factor $F_\pi(q^2)$ into one-loop virtual amplitudes associated with Final-State Radiation (FSR) and Initial-Final State Interference (IFI). Unlike F $\times$ sQED, FsQED utilizes a dispersive representation of the form factor:
$F_\pi(q^2) = 1 + \frac{q^2}{\pi} \int_{4m_\pi^2}^\infty ds' \frac{\text{Im} F_\pi(s')}{s'(s' - q^2 - i\epsilon)}$
This allows the form factor to vary with the loop momentum.
One-Loop Calculation: The calculation separates virtual corrections into subsets based on the number of virtual photons ( $1\gamma^*$ and $2\gamma^*$ ) and the current emitting the real photon (ISR or FSR).
- Infrared (IR) Regularization: The authors address IR divergences by introducing a fictitious photon mass $\lambda$ . For the $2\gamma^*$ diagrams, which involve more complex topologies including genuine five-point functions, a careful "add-and-subtract" procedure is employed. This isolates the IR-sensitive real and imaginary parts, ensuring that dispersive integrals remain finite.
- Reduction: Five-point scalar functions are reduced to combinations of four-point and lower-point functions using the hexagon.m package. End-point singularities in the dispersive integration are treated analytically.
Beyond FsQED Estimates: Following recent literature, the authors estimate contributions beyond the charged pion pole (the leading FsQED term) by including neutral pseudoscalar poles ( $P = \pi^0, \eta, \eta'$ ) in the Hadronic Light-by-Light (HLbL) scattering kernel. These are modeled using Vector Meson Dominance (VMD) for the transition form factors.
$\phi$ -Resonance Mechanisms: For center-of-mass energies around the $\phi$ -meson resonance, the paper incorporates additional real radiation mechanisms beyond F $\times$ sQED. These include direct radiative $\phi$ decays via scalar resonances ( $f_0, \sigma$ ), double resonance processes, and bremsstrahlung involving vector and axial-vector mesons, modeled within Chiral Perturbation Theory ( $\chi$ PT).
Implementation: All calculations are implemented in the Monte Carlo event generator BabaYaga@NLO, matched with a Parton Shower (PS) algorithm to achieve NLOPS accuracy.

Key Contributions and Results
The paper presents numerical results for differential cross sections ( $d\sigma/dM_{\pi\pi}$ and $d\sigma/d\theta_+$ ) and the forward-backward asymmetry ( $A_{FB}$ ) across four experimental scenarios: KLOE (large and small angles), BESIII, and B-factories.

FsQED vs. F $\times$ sQED:
- Invariant Mass ( $M_{\pi\pi}$ ): The complete set of FsQED corrections to the invariant mass distribution is of the order of a few permille ( $\sim 0.1\% - 0.3\%$ ) across all scenarios. The $(1\gamma^*)$ subset dominates the correction in most cases, driven by collinear logarithms in KLOEII and suppressed by the small value of $F_\pi(s)$ in BESIII and B-factories.
- Angular Distribution ( $\theta_+$ ): Structure-dependent corrections to the angular distribution are larger, reaching the percent level. These are dominated by the $(2\gamma^*, \text{ISR})$ topology, which interferes asymmetrically with the leading-order ISR.
- Forward-Backward Asymmetry ( $A_{FB}$ ): Corrections to $A_{FB}$ are approximately 1% near the $\rho$ -peak and a few permille below it.
Model Comparison (FsQED vs. GVMD):
The authors compare their FsQED results with predictions obtained using the Generalised Vector Meson Dominance (GVMD) model (previously implemented in BabaYaga@NLO). The two formalisms, despite using different methodologies and form factor parameterizations, show striking agreement (at the $10^{-4}$ level) for both the invariant mass and the forward-backward asymmetry. This validates the robustness of the structure-dependent corrections for radiative return.
Beyond FsQED Contributions:
The estimated contributions from neutral pseudoscalar poles in the $(2\gamma^*, \text{FSR})$ subset are found to be below the permille level in all scenarios. The authors note that a complete quantitative estimate of other beyond-FsQED effects requires further investigation.
$\phi$ -Resonance Effects:
In the $\phi$ -resonance region ( $\sqrt{s} \approx m_\phi$ ), additional FSR mechanisms (direct radiative decays, double resonances, $\chi$ PT bremsstrahlung) significantly impact observables.
- The direct radiative decay via the $f_0$ resonance induces effects of several percent in the invariant mass spectrum and forward-backward asymmetry around 980 MeV.
- While these effects are negligible for the integrated cross section in energy scans (except in a narrow window around the $\phi$ ), they are crucial for accurate modeling of the radiative return process at the $\phi$ -factory.

Significance
The paper claims that the structure-dependent corrections to $e^+e^- \to \pi^+\pi^-\gamma$ are non-negligible for radiative return measurements aiming at sub-percent precision. The work completes the treatment of these corrections in the BabaYaga@NLO generator, allowing for the evaluation of modeling uncertainties associated with pion-photon interactions in both energy scan and radiative return experiments.

By demonstrating the agreement between the FsQED and GVMD approaches, the authors provide a cross-validation of theoretical tools essential for reducing the systematic uncertainty in the determination of the pion form factor and, consequently, the HVP contribution to the muon $g-2$ . The inclusion of $\phi$ -resonance specific mechanisms further refines the description of data from the KLOE experiment. The authors conclude that these novel features are now available in a new version of BabaYaga@NLO to support ongoing and future precision measurements at flavor factories.

More like this