Matching of perturbative and exponentiated initial state radiation corrections to $e^+e^-$-annihilation

This paper analyzes higher-order initial state radiation corrections in electron-positron annihilation by presenting numerical results for future collider energies, estimating uncertainties, and proposing a modified scheme that matches exponentiated photonic and non-singlet pair corrections with existing analytic calculations using a new DIS-like subtraction method.

The Gist

Imagine you are trying to measure the exact speed of a car (an electron) as it zooms toward a wall to crash into another car (a positron). In the world of particle physics, this crash is called annihilation, and it creates a burst of new particles. Scientists want to predict exactly how this crash will look to test their theories about the universe.

However, there's a problem. As the cars speed up, they don't just travel in a straight line; they constantly emit tiny sparks of light (photons) and occasionally spit out small pairs of new particles. These are called Initial State Radiation (ISR). If you ignore these sparks, your prediction of the crash will be wrong.

This paper is about how to calculate the effect of these "sparks" with extreme precision, especially for future, super-powerful particle colliders. Here is the breakdown of their solution using simple analogies:

1. The Two Ways to Count the Sparks

The authors discuss two different methods for counting these sparks, and they realized they needed to combine them.

• Method A: The "Step-by-Step" Calculator (Perturbation)
  Imagine trying to count every single spark by hand, one by one. You calculate the effect of one spark, then two sparks, then three. This is very accurate for the first few sparks, but as you try to count the 10th or 10th spark, the math gets incredibly messy and hard to finish. This is the "perturbative" approach. It's great for the obvious, big effects but struggles with the infinite number of tiny, faint sparks.

• Method B: The "Magic Formula" (Exponentiation)
  Imagine instead of counting every spark, you use a magic formula that assumes the sparks happen in a specific, predictable pattern (like a crowd of people leaving a stadium). This formula, called "exponentiation," is great at predicting the overall behavior of millions of tiny sparks all at once. However, it might miss some specific, weird details that only show up in the "Step-by-Step" method.

The Paper's Solution:
The authors created a "hybrid" system. They took the "Step-by-Step" results (which are known to be very accurate for the first few orders) and "matched" them with the "Magic Formula."

• They used the Magic Formula to handle the millions of tiny, soft sparks.
• They used the Step-by-Step math to handle the specific, hard-to-calculate details.
• Crucially, they made sure not to count the same sparks twice (avoiding "double counting").

2. The "Tail" and the "Residue"

When you mix these two methods, there's a leftover bit of math called the "tail."

• Think of the "Step-by-Step" method as a detailed map of a city.
• Think of the "Magic Formula" as a satellite view of the whole country.
• The authors figured out how to subtract the parts of the satellite view that are already on the detailed map, so they only add the new information the satellite provides. This ensures their final prediction is the most accurate possible version of both maps combined.

3. Changing the Rules of the Game (The Subtraction Scheme)

In physics, sometimes you have to choose a "ruler" or a "scheme" to measure things. The standard ruler (called the MS scheme) works well, but it makes the math for the "Magic Formula" very complicated because it includes some messy, extra terms that cancel out later but are annoying to carry around.

The authors invented a new ruler (a new subtraction scheme).

• Analogy: Imagine you are baking a cake. The standard recipe tells you to measure flour, then sift it, then measure sugar, then sift it. It works, but it's tedious.
• The authors' new recipe says, "Let's measure the flour and sugar together in a specific way so we don't have to sift them separately."
• This new method makes the math much cleaner and easier to handle, especially when the particles are moving very fast (near the speed of light).

4. How Precise Are They?

The authors ran the numbers for future colliders (like the FCC-ee and CEPC).

• They found that their new hybrid method reduces the "guesswork" (theoretical uncertainty) to a tiny fraction of a percent.
• Specifically, at the energy level where the famous "Z boson" is created, their uncertainty is about 0.0004%.
• To put that in perspective: If you were measuring the distance from the Earth to the Moon, their method would be accurate to within a few centimeters.

Summary

The paper doesn't claim to discover a new particle or cure a disease. Instead, it provides a better calculator for physicists.

1. It combines a detailed, step-by-step counting method with a powerful, all-encompassing formula.
2. It invents a new way to organize the math to make it less messy.
3. It proves that this combination allows scientists to predict the results of future particle collisions with unprecedented precision, ensuring that when they build these massive machines, they know exactly what to look for.

The authors conclude that while their method is a huge improvement, the work isn't done; they need to keep refining the math to include even more subtle effects, like particles interacting after the crash (Final State Radiation).

Technical Summary

Technical Summary: Matching of Perturbative and Exponentiated Initial State Radiation Corrections to $e^+e^-$-Annihilation

Problem Statement
Precise theoretical predictions for radiative corrections in electron-positron annihilation are critical for future high-energy colliders (e.g., FCC-ee, CEPC) aiming for precision tests of the Standard Model. A primary source of uncertainty in these predictions is the theoretical description of Initial State Radiation (ISR). While multi-loop calculation methods have advanced, computing higher-order corrections for realistic differential observables remains complex. Existing methods face a dichotomy: direct perturbative calculations are limited by computational order, while exponentiation methods (resummation) capture dominant logarithmic terms but may miss specific finite-order contributions. The paper addresses the need to reconcile these approaches to minimize theoretical uncertainties and avoid double-counting of corrections.

Methodology
The authors employ the QED structure function (parton distribution function) approach, which systematically calculates corrections enhanced by large logarithms ($L = \ln(\mu_F^2/\mu_R^2)$). The methodology involves three main technical components:

1. Perturbative Expansion: The differential cross-section is expressed as a convolution of electron parton distribution functions (PDFs) and partonic cross-sections, expanded up to Next-to-Next-to-Leading Logarithmic (NNLL) accuracy. This includes contributions up to $O(\alpha^2 L^0)$ and specific higher-order logarithmic terms ($c_{kl}$ coefficients).
2. Exponentiation Scheme: The authors utilize the Gribov-Lipatov solution for the pure photonic part of the electron PDF and propose a modified scheme for the simultaneous exponentiation of pure photonic and non-singlet (NS) pair corrections. This involves defining a modified parameter $\tilde{\beta}$ to account for pair corrections within the exponent.
3. Matching Procedure: To combine the advantages of both methods, the authors construct an improved cross-section ($\sigma_{imp}$) by subtracting the series expansion of the exponentiated result (up to the orders present in the perturbative calculation) from the full exponentiated result, and then adding the known perturbative result. This ensures that higher-order terms in the exponent are retained while lower-order terms are replaced by exact perturbative calculations, thereby avoiding double counting.
4. Scheme Modification: A new subtraction scheme is proposed to replace the standard $\overline{\text{MS}}$ scheme. In the $\overline{\text{MS}}$ scheme, the evolution kernel generates higher powers of $\ln(1-z)$ that do not correspond to physical singularities and complicate the matching with exponentiated forms. The new scheme modifies the initial condition $d^{(1)}_{ee}(x)$ to eliminate these unphysical terms, simplifying convolutions and improving behavior near $z \to 1$.

Key Contributions

• Unified Matching Framework: The paper presents a formalism to match existing analytic higher-order perturbative calculations (up to $O(\alpha^2)$ and specific logarithmic orders) with exponentiated results for both pure photonic and non-singlet pair corrections.
• Modified Exponentiation: A specific expression is derived for the joint exponentiation of photonic and non-singlet pair corrections, adjusting the exponent parameter $\beta$ to $\tilde{\beta}$ to ensure consistency with perturbative expansions at $O(\alpha^2 L^1)$.
• New Subtraction Scheme: The authors introduce a modified factorization scheme that removes unphysical logarithmic terms ($\ln(1-z)$) from the splitting functions, facilitating easier convolution calculations and better compatibility with exponentiation.
• Uncertainty Estimation: A systematic method is proposed to estimate theoretical uncertainties by comparing the magnitude of uncalculated higher-order terms (e.g., estimating $h_{66}$ based on $h_{55}$ and $h_{44}$).

Results

• Numerical Predictions: Numerical results for relative corrections ($h_{kl}$) are provided for various center-of-mass energies ($M_Z \pm 1/2$, 160 GeV, 240 GeV, and 3 TeV) and different kinematic cuts ($z_{min}$).
• Radiative Return Effects: The analysis highlights that at energies of 160 GeV and 240 GeV, radiative return to the $Z$-resonance significantly enhances corrections (up to ten times larger than at the $Z$-peak) for small $z_{min}$ values.
• Theoretical Uncertainty: At the $Z$-peak ($\sqrt{s} = M_Z$), the estimated total theoretical uncertainty in the description of ISR for the total cross-section is approximately $4 \times 10^{-4}\%$. The dominant contribution to this uncertainty is identified as the uncalculated NNLL term ($h_{31}$).
• Scheme Comparison: The new subtraction scheme is shown to simplify the calculation of convolutions, particularly in the region $z \to 1$, compared to the standard $\overline{\text{MS}}$ scheme.

Significance and Claims
The paper claims that the proposed matching scheme allows for the combination of the precision of perturbative calculations with the resummation power of exponentiation, thereby improving the overall precision of theoretical predictions for $e^+e^-$ annihilation. The authors state that this framework is designed to be flexible, allowing future perturbative results to be easily incorporated. The work is intended for implementation in the ZFITTER program and serves as a foundation for differential treatments in Monte Carlo generators (specifically ReneSANCe). The authors modestly note that while the current uncertainty is small, further improvements require the inclusion of new perturbative calculations, specifically regarding final state radiation (FSR) and ISR-FSR interference, which are not fully addressed in this study.