⚛️ high-energy experiments

Factorization of the Energy-Energy Correlation in the two-jet limit in the massive case

This paper investigates non-logarithmic heavy-quark mass effects in the two-jet limit of the Energy-Energy Correlation by introducing a novel partial event fraction calculable via real emission diagrams in four dimensions and proposing an improved factorization scheme that ensures a smooth massless limit through an angle-dependent coefficient function.

Original authors: Ugo Giuseppe Aglietti, Giancarlo Ferrera, Lorenzo Rossi

Published 2026-06-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ugo Giuseppe Aglietti, Giancarlo Ferrera, Lorenzo Rossi

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 at a massive, high-energy particle collider, like a giant cosmic pinball machine. When two particles smash together, they shatter into a spray of smaller pieces. Physicists want to understand how these pieces fly apart. One specific way they measure this is called the Energy-Energy Correlation (EEC).

Think of the EEC as a way to measure the "angle of separation" between pairs of particles flying out of the crash. If two particles fly off in exactly opposite directions (back-to-back), that's a "two-jet" event. If they fly off in the same direction, that's a "forward" event. The EEC tells us how likely it is to find pairs at any given angle.

For a long time, physicists have been able to calculate this very precisely for massless particles (like light). But real heavy particles, like the "bottom" or "charm" quarks, have mass. This mass acts like a heavy backpack the particles are wearing, making their flight paths slightly different and much harder to calculate.

Here is what this paper does, broken down into simple concepts:

1. The Problem: The "Heavy Backpack" Mess

When physicists tried to calculate the EEC for these heavy particles, they ran into a mathematical wall.

The Old Way: To get the full picture, they had to calculate two things: the "real" events (where a particle actually flies out) and the "virtual" events (ghostly, invisible loops of energy that exist only in the math).
The Issue: When you add mass to the mix, the math for these "ghostly" loops becomes incredibly messy and requires complex, abstract math tricks (changing the number of dimensions) to solve. It's like trying to solve a puzzle where the pieces keep changing shape.

2. The New Trick: A "Partial" View

The authors of this paper came up with a clever shortcut. Instead of trying to measure every possible angle (from back-to-back to forward), they decided to only look at the back-to-back region (the two-jet limit) and ignore the forward region.

The Analogy: Imagine you are trying to count how many people are walking out of a stadium. The old method required you to count everyone leaving through every single gate, including the ones where the crowd is so dense you can't see anything.
The New Method: The authors say, "Let's just count people leaving through the back exit."
Why it helps: By ignoring the forward region, they found that they didn't need to calculate those messy "ghostly" virtual loops anymore. They could just look at the real particles flying out. It's like solving the puzzle using only the visible pieces, which is much easier and doesn't require changing the rules of the game (mathematical dimensions).

3. The "Smoothness" Problem

After doing the math for the heavy particles, they noticed something weird.

The Glitch: If you take a heavy particle and slowly make it lighter and lighter until it has zero mass, the math result should smoothly turn into the result for a massless particle.
The Reality: In their standard calculation, the result "jumped" or "snapped" when the mass hit zero. It was like a car driving smoothly down a road that suddenly had a giant, invisible pothole right at the finish line. This "discontinuity" made no physical sense.

4. The Solution: A New Map

To fix this "snapping" problem, the authors built a new factorization scheme (a new way of organizing the math).

The Fix: They realized that the "heavy backpack" effect wasn't just a constant number; it depended on the angle at which the particles were flying.
The Analogy: Imagine you are adjusting the suspension on a car. The old way was to have a fixed setting that worked for smooth roads but broke on bumpy ones. The new way is to have a suspension that automatically adjusts based on the specific bump you are hitting.
The Result: By making their "coefficient function" (the part of the math that handles the heavy mass) depend on the angle, they smoothed out the road. Now, as the particle gets lighter, the math flows perfectly into the massless result without any jumps.

5. What They Actually Did

They calculated the energy distribution for heavy quarks and gluons flying out of a collision.
They verified their numbers by comparing them to older, computer-heavy calculations from the 1980s and found they matched up very well (within about 2%).
They created a new, cleaner mathematical formula that works for heavy particles but behaves perfectly when those particles become light.

Summary

This paper is about cleaning up the math used to predict how heavy particles fly apart after a collision. The authors found a way to ignore the messy parts of the calculation by focusing on a specific angle, and then fixed a "glitch" in the math so that heavy particles and light particles fit together seamlessly. They didn't invent a new machine or a new drug; they just made the theoretical map of the particle world more accurate and easier to read.

Technical Summary: Factorization of the Energy-Energy Correlation in the Two-Jet Limit in the Massive Case

Problem Statement
The Energy-Energy Correlation (EEC) function is a precision observable in perturbative QCD (pQCD), describing the angular separation between hard particle pairs in $e^+e^-$ annihilation. While high-accuracy calculations exist for the massless case, the inclusion of heavy-quark mass effects is essential for high-precision phenomenology. Previous work by the authors addressed logarithmic mass effects (terms of the form $\alpha_S^n \log^k(m^2/Q^2)$ ), which are resummed into a massive Sudakov form factor. This paper focuses on the remaining non-logarithmic mass effects (terms of the form $\alpha_S^n m^k/Q^k$ ), which enter the coefficient function and the remainder function of the standard factorization formula.

A specific technical challenge arises in the two-jet limit ( $\chi \to 0$ ). The standard "total" event fraction, which integrates over the entire angular range including the forward region ( $\chi \to \pi$ ), requires the calculation of one-loop virtual corrections at the endpoints to satisfy unitarity and boundary conditions. In the massive case, these one-loop virtual corrections are currently unknown. Furthermore, the standard factorization scheme, when naively extended to the massive case, exhibits a discontinuity in the massless limit ( $\eta = 2m/Q \to 0$ ), failing to reproduce the known massless remainder function continuously.

Methodology
To address these issues, the authors employ the following methodology:

Definition of a "Partial" Event Fraction:
The authors define a new observable, a "partial" event fraction, restricted to the two-jet region ( $0 \le \chi \le \chi_M$ ) and explicitly excluding the forward region ( $\chi_M < \pi$ , typically $\chi_M \sim \pi/2$ ).
- Advantage: In this restricted region, the first-order ( $O(\alpha_S)$ ) calculation requires only real emission diagrams ( $e^+e^- \to Q\bar{Q}g$ ) in $D=4$ space-time dimensions. Virtual corrections and dimensional regularization ( $D \neq 4$ ) are not required, as the infrared singularities associated with the forward region are excluded by the definition of the observable.
Kinematic and Matrix Element Evaluation:
- The authors perform a detailed kinematic analysis of the three-body final state ( $Q\bar{Q}g$ ) with a massive quark-antiquark pair. They derive the allowed region in the quark-antiquark energy plane (Dalitz plot) and the functional dependence of the correlation angle $\chi$ on the parton energies.
- They evaluate the squared matrix elements for both vector and axial currents, generalizing known massless formulas to include the mass parameter $\eta$ .
- The calculation involves integrating over the phase space. The authors reduce the problem to a two-dimensional integral over quark and antiquark energies. The first integration is performed analytically, while the second is evaluated numerically.
- Special attention is paid to the quark-antiquark correlation, which is more complex than the quark-gluon correlation due to the presence of two massive partons. This leads to a kinematic constraint equation that, when squared to eliminate radicals, introduces spurious solutions. The authors rigorously identify the physical domains for the two distinct solutions ( $x^*$ and $x^{**}$ ) of the kinematic constraint.
Factorization and Resummation:
- The authors extend the standard massless factorization scheme to the massive case. They numerically evaluate the massive remainder function and the coefficient function by subtracting the first-order expansion of the massive Sudakov form factor (using an improved form factor that is smooth in $\chi$ and valid for arbitrary $m/Q$ ) from the computed EEC spectrum.
- They analyze the behavior of these functions as $\eta \to 0$ .

Key Contributions and Results

Numerical Evaluation of First-Order Corrections: The paper provides the first-order massive EEC spectrum for the partial event fraction. The results are expressed as one-dimensional integrals depending on the correlation angle $\chi$ and the mass parameter $\eta$ .
Validation: The numerical results for the EEC spectrum are compared with previous numerical calculations from the 1980s [39, 40]. The authors find good agreement, with discrepancies attributed to normalization uncertainties in the older works (estimated at $\sim 2\%$ ).
Discovery of a Discontinuity: The authors demonstrate that the standard massive factorization scheme is not smooth in the massless limit. Specifically, the massive remainder function $Rem^{(1)}(\chi; \eta)$ does not converge to the ab initio massless remainder function $Rem^{(1)}_0(\chi)$ as $\eta \to 0$ . Instead, it is shifted by a large, non-vanishing amount $a_\eta \approx 6.5$ . This discontinuity arises because the limits $\eta \to 0$ and $\chi \to 0$ do not commute in the standard scheme.
Proposed Improved Factorization Scheme: To resolve the discontinuity, the authors propose a new factorization scheme.
- They "move" the mass-dependent term from the remainder function to the coefficient function.
- This results in a new coefficient function, $\tilde{C}^{(1)}$ , which explicitly depends on the correlation angle $\chi$ .
- In this new scheme, the remainder function becomes identical to the massless one ( $\tilde{Rem}^{(1)} = Rem^{(1)}_0$ ), ensuring a smooth limit as $\eta \to 0$ .
- The new coefficient function provides a positive correction in the two-jet region, enhancing the predicted rate compared to the massless case.

Significance
The paper claims that this work provides the first complete theoretical control over non-logarithmic heavy-quark mass effects in the EEC function within the two-jet limit. By defining the partial event fraction, the authors bypass the need for unknown virtual corrections, allowing for a purely real-emission calculation in four dimensions.

The primary theoretical significance lies in the identification and resolution of the non-smooth massless limit in the standard factorization scheme. The proposed improved scheme, which introduces a $\chi$ -dependent coefficient function, restores continuity and physical consistency in the limit of vanishing quark mass. This framework enables precision phenomenological analyses of heavy-quark mass effects in EEC data, complementing previous work on logarithmic mass effects. The authors note that while the current work is at next-to-leading order (NLO), the established framework paves the way for future extensions to higher orders ( $O(\alpha_S^2)$ ), though such calculations would require significant numerical effort (e.g., Monte Carlo integration).