⚛️ phenomenology

Double-real corrections to color singlet decay in a parton-shower inspired scheme

This paper introduces a local infrared subtraction method for next-to-next-to-leading order QCD calculations in color singlet decays, utilizing scalar radiators and pure splitting functions to handle overlapping singularities and demonstrating the finiteness and numerical convergence of the double-real remainder for $e^+e^-\to q\bar{q}$ processes.

Original authors: John M. Campbell, Stefan Höche, Max Knobbe

Published 2026-06-26

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: John M. Campbell, Stefan Höche, Max Knobbe

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 you are trying to predict exactly how a billiard ball will break when it hits a rack of other balls. In the world of particle physics, this "break" is a collision that creates a shower of new particles. Scientists use a set of rules called Quantum Chromodynamics (QCD) to calculate these outcomes.

For decades, these calculations have been incredibly accurate. However, as our telescopes (particle colliders like the LHC) get sharper and our measurements become more precise, the old calculation methods are starting to show their limits. They are like trying to measure the width of a hair with a ruler meant for measuring a football field. To keep up with the new "Tera-Z" precision goals of future colliders, physicists need a new, ultra-fine ruler.

This paper introduces a new, smarter way to do these complex calculations, specifically for a type of event called a "color singlet decay" (think of it as a clean, isolated explosion of particles).

Here is the breakdown of their new method using simple analogies:

1. The Problem: Overlapping Messes

When particles collide, they sometimes spit out extra "gluons" (the glue holding particles together). Sometimes, they spit out two at once.

The Issue: In the math, these extra gluons create "singularities." Imagine trying to count the number of people in a room where some people are standing on top of each other, and others are invisible. The numbers blow up to infinity, making the calculation impossible.
The Old Way: Previous methods tried to fix this by cutting the room into different zones and counting people in each zone separately. It worked, but it was messy and hard to connect to how computers simulate particle showers (like a video game engine).
The New Way: The authors propose a "local subtraction" method. Instead of trying to fix the whole room at once, they identify exactly where the "people on top of each other" are and subtract a counter-amount right there.

2. The Solution: The "Scalar Radiator" and "Splitting"

The authors built their new tool using two main ingredients:

Scalar Radiators: Think of these as a generic "noise meter." They measure how much "gluon noise" is being created by the collision without worrying about the specific spin or direction of the particles yet. It's a simplified, robust way to measure the chaos.
Pure Splitting Functions: These are the specific rules for how one particle breaks into two.

The Magic Trick: Untangling the Knots
The hardest part of the math is when two singularities overlap (like two knots tied together). The authors use a technique called partial fractioning.

Analogy: Imagine you have a tangled ball of yarn with two knots. Instead of trying to pull the whole ball apart, you use a special tool to separate the two knots so you can untie them one by one.
In their math, they separate the overlapping "noise" so that every singularity can be assigned to a specific, unique cause. This ensures the math stays clean and doesn't double-count anything.

3. The Map: Following the Trail

To make sure the math matches reality, they need a map of where the particles go.

They use a mapping system that is essentially a "step-by-step" guide. It's like a GPS that tells you how to get from the starting point to the destination by taking one turn, then another.
Crucially, this map is designed to look exactly like the maps used in parton showers (the computer programs that simulate how particles evolve). This is a big deal because it means their high-precision math can be easily plugged into existing simulation software, bridging the gap between "perfect theory" and "practical simulation."

4. The Test: Does it Work?

The authors tested their new method on a specific scenario: an electron and a positron colliding to create a quark and an antiquark, which then emit extra gluons.

The Result: They showed that when they applied their subtraction method, the "infinite" numbers canceled out perfectly, leaving a finite, sensible result.
The Stability: They ran millions of simulations. Even when they pushed the system to the extreme limits (where particles are almost invisible or moving in the exact same direction), the method remained stable and accurate.
The Speed: They calculated that this method is fast enough to be used for the massive amounts of data expected from future colliders (the "Tera-Z" option), requiring only a modest amount of computing power to achieve extreme precision.

Summary

In short, this paper presents a new, cleaner, and more efficient way to calculate the messy details of particle collisions. By using a "noise meter" approach and a step-by-step mapping system that fits perfectly with existing simulation tools, they have created a method that is ready to handle the extreme precision demands of the next generation of particle physics experiments. They haven't just fixed the math; they've made it compatible with the tools scientists actually use to predict what will happen in the lab.

Technical Summary: Double-real corrections to color singlet decay in a parton-shower inspired scheme

Problem Statement
High-precision predictions for high-energy physics experiments, particularly for future facilities like the FCC-ee, require Next-to-Next-to-Leading Order (NNLO) Quantum Chromodynamics (QCD) calculations. A primary challenge in these calculations is the treatment of infrared (IR) singularities arising from soft-gluon divergences and collinear singularities. While the Bloch-Nordsieck and Kinoshita-Lee-Nauenberg theorems guarantee the cancellation of these singularities between real and virtual corrections, practical computation requires isolating them to compute finite remainders. Existing NNLO subtraction methods (e.g., sector decomposition, antenna subtraction, CoLoRFulNNLO) often struggle with the complexity of overlapping singularities in multi-parton final states and the difficulty of matching fixed-order results to parton showers with Next-to-Leading Logarithmic (NLL) accuracy. Specifically, there is a need for a subtraction scheme that utilizes simple kinematic mappings compatible with parton-shower recoil definitions and separates singular structures in a way that facilitates matching to resummed calculations.

Methodology
The authors introduce a local IR subtraction method for NNLO double-real corrections in color singlet decays (specifically $e^+e^- \to q\bar{q}$ ). The method is built upon a recently developed algorithm for assembling splitting functions from spin-independent and spin-dependent components.

Structure of Counterterms: The subtracted double-real squared matrix element ( $M_{RR}$ ) is constructed by combining NLO ( $S_1$ ) and NNLO ( $S_2$ ) subtraction terms. The scheme separates the calculation into strongly ordered terms ( $S_1 \otimes S_1$ ), pure NNLO remainders ( $S_2 - \sum S_1 \otimes S_1$ ), and mixed terms involving finite $n+1$ -parton contributions.
Kinematics and Gauge Choice: The method employs a momentum mapping based on identified particles and recoil momenta, iterating NLO mappings to handle double emissions. Crucially, the authors utilize a time-like gauge vector (identified with the recoil momentum) rather than a light-like axial gauge. This choice aligns with the recoil definitions of NLL-safe parton-shower algorithms, avoiding the impracticality of negative splitting functions and azimuthal averaging associated with axial gauges in additive matching schemes.
Disentangling Overlapping Singularities: To address overlapping singularities in scalar radiators (where gluons are radiated coherently off QCD charge multipoles), the authors employ partial fractioning techniques. This disentangles the soft-enhanced components, attributing singularities to unique collinear configurations.
- NLO Case: Partial fractioning is applied to eikonal factors using angular variables.
- NNLO Case (Abelian & Non-Abelian): The scalar radiators for double-gluon emission and $q\bar{q}$ pair emission are decomposed. The non-abelian contributions, which contain entangled triple-collinear poles, are partitioned using angular weights to map them to specific collinear sectors.
Spin Correlations: The authors define an improved method for constructing polarization tensors using partial polarization vectors. This ensures consistency between on-shell counterterms and off-shell terms, avoiding gauge-dependent poles that arise from mismatched momentum assignments in the subtraction terms.

Key Contributions and Results

Analytic Integration: The authors computed the phase-space integrals of the scalar counterterms analytically for the back-to-back kinematic configuration. Using reverse unitarity and the Kira package for integral reduction, they derived the pole structures ( $\epsilon^{-4}$ to $\epsilon^0$ ) for abelian, non-abelian, and quark-antiquark emission contributions.
Numerical Validation:
- Infrared Limits: The method was validated using the Pepper and Sherpa event generators. The subtracted remainder ( $\Delta$ ) was shown to vanish linearly with the infrared scaling parameter $\lambda$ in single-unresolved limits and quadratically in double-unresolved limits, confirming the correct cancellation of leading and sub-leading singularities.
- Numerical Integration: The authors verified that the phase-space integrals of the scalar counterterms can be computed numerically using sector decomposition (pySecDec) combined with the specific phase-space parametrization. The numerical results for the integrated radiators matched the analytic back-to-back results with high precision.
- Performance: A preliminary integration of the subtracted double-real contribution to the $e^+e^- \to q\bar{q}$ cross-section at the Z pole yielded a statistical uncertainty of $654$ fb ( $< 2 \cdot 10^{-5}$ relative to the Born cross-section) with $5.1 \cdot 10^8$ phase-space points.
Jet Rates: Differential jet rates ( $2 \to 3$ and $3 \to 4$ ) were computed, demonstrating the algorithm's performance in NLO-like situations and confirming sufficient numerical convergence for precision physics.

Significance and Claims
The paper claims that this new scheme provides a straightforward separation of singular structures and establishes a direct correspondence to resummed calculations in the form of parton showers. By using simple, iterated kinematic mappings that match the recoil definitions of NLL-safe parton showers, the method facilitates the combination of fixed-order NNLO calculations with numerical resummation. The authors posit that this approach will enable the use of double-real counterterms as local higher-order corrections to parton-shower splitting kernels, thereby aiding the development of fully differential parton showers at higher precision. The work represents a step toward closing the precision gap between theory and experiment anticipated at future electron-positron colliders, specifically noting the requirements of a Tera-Z option at the FCC-ee. The authors remain modest, noting that the current results are unphysical without the inclusion of real-virtual and double-virtual corrections, and that further work is required to integrate one-loop scalar radiator functions and derive corresponding counterterms to complete the full NNLO calculation for $e^+e^- \to \text{hadrons}$ .

1. The Problem: Overlapping Messes

2. The Solution: The "Scalar Radiator" and "Splitting"

3. The Map: Following the Trail

4. The Test: Does it Work?

Summary

More like this