A global analysis of Energy-Energy Correlation data: determination of $α_S$ and non-perturbative QCD parameters

This paper presents a comprehensive global analysis of Energy-Energy Correlation data across a wide range of center-of-mass energies, utilizing N³LL resummation matched to $\mathcal{O}(\alpha_S^3)$ fixed-order calculations and including non-perturbative power corrections to precisely determine the strong coupling constant $\alpha_S(m_Z^2) = 0.119 \pm 0.002$ and non-perturbative QCD parameters.

The Gist

Imagine you are a detective trying to understand the rules of a chaotic crime scene, but instead of fingerprints, you are looking at the paths of tiny, invisible particles flying apart after a high-speed collision. This is exactly what the physicists in this paper did, but their "crime scene" is the subatomic world of Quantum Chromodynamics (QCD), the theory that explains how the universe's building blocks stick together.

Here is the story of their investigation, broken down into simple concepts.

The Crime Scene: Electron-Positron Annihilation

Picture two tiny particles, an electron and a positron, smashing into each other at nearly the speed of light. When they collide, they vanish in a flash of energy and instantly transform into a shower of new particles (mostly quarks and gluons) that fly out in all directions.

The physicists are interested in a specific pattern called the Energy-Energy Correlation (EEC). Think of this as measuring the "angle of separation" between pairs of particles in the resulting shower.

• The Back-to-Back Case: Most of the time, the energy shoots out in two opposite directions (like a firecracker popping). This is the "two-jet" region.
• The Messy Case: Sometimes, the energy spreads out more randomly.

The Problem: The "Mathematical Fog"

For decades, physicists have tried to predict exactly how these particles spread out using math. However, there's a catch.

• The "Fog": When particles fly almost directly opposite each other (the back-to-back case), the standard math equations get clogged with "infinite" terms (logarithms). It's like trying to drive a car through a thick fog; you can't see the road ahead, so your predictions become unreliable.
• The "Hidden Rules": There are also invisible forces at play that the standard math can't see. These are non-perturbative effects—think of them as the "stickiness" of the particles as they form into real matter. Standard math ignores this stickiness because it's too complex to calculate directly.

The Solution: A New Toolkit

The authors of this paper built a super-powered toolkit to clear the fog and measure the stickiness. They did three main things:

1. Resumming the Fog (N3LL): They developed a new way to "resum" (re-summarize) the infinite mathematical terms. Imagine taking a blurry, foggy photo and using a super-algorithm to sharpen it until you can see every detail clearly. They did this to the highest level of precision ever attempted (called N3LL).
2. The "Dispersive" Model (The Stickiness): To handle the invisible "stickiness" (non-perturbative effects), they didn't just guess. They used an analytic dispersive model.
   • Analogy: Imagine trying to predict how a rubber band snaps. You can calculate the physics of the snap (perturbative), but you also need to know how stretchy the rubber is (non-perturbative). Instead of using a generic "rubber band" from a toy store, they created a custom mathematical model for the specific "rubber" of the universe, fitting it directly to the data.
3. The Global Fit (The Master Puzzle): Instead of looking at just one experiment, they gathered data from 691 different measurements spanning 30 years.
   • They used data from the highest energy collisions (91 GeV, like at the Z-boson peak) down to very low energies (7.7 GeV).
   • Analogy: It's like a detective solving a mystery by looking at crime scenes from 1990, 2000, and 2020 simultaneously. By seeing how the "stickiness" changes as the energy changes, they could separate the rules of the game from the specific conditions of the day.

The Big Discoveries

1. The "Glue" Strength ($\alpha_S$)
The most famous result of this paper is a new, incredibly precise measurement of the Strong Coupling Constant ($\alpha_S$).

• Analogy: If the universe were a giant Lego set, this number tells you how sticky the glue is between the bricks.
• The Result: They found the glue strength to be 0.119. This matches perfectly with other global averages, confirming that our understanding of the universe's glue is solid.

2. The "Evolution" of the Stickiness (Collins-Soper Kernel)
This is the most novel part of the paper. They discovered that the "stickiness" of the universe isn't constant; it changes depending on how hard you smash the particles.

• Analogy: Imagine stretching a piece of taffy. If you pull it slowly, it stretches one way. If you yank it fast, it behaves differently.
• They measured how this "taffy" changes over different energy levels. This allowed them to extract a fundamental rule called the Collins-Soper kernel.
• Why it matters: This kernel is a universal rule that applies not just to particle collisions, but also to how particles behave inside protons and in other high-energy experiments. It's like finding a universal law of physics that connects different types of experiments.

3. Heavy Quarks (The Bottom Mass)
They also checked if the heavy "bottom" quark (a heavy type of particle) changed the results.

• The Finding: Including the heavy mass didn't change the final answer much; it just shifted the "stickiness" numbers slightly. This tells them that while heavy particles matter, the main rules of the game are dominated by the lighter, more common particles.

The Bottom Line

This paper is a triumph of "global thinking." By combining the sharpest mathematical tools (to clear the fog) with a clever model for the invisible forces (the stickiness), and by looking at data from every era of particle physics, the authors:

1. Confirmed the strength of the universe's glue with high precision.
2. Mapped out how the "stickiness" of the universe evolves with energy.
3. Proved that their mathematical framework works perfectly across a massive range of energies, from the very small to the very large.

In short, they took a messy, foggy, 30-year-old puzzle and solved it with a single, elegant picture that fits every piece perfectly.

Technical Summary

1. Problem Statement

The Energy-Energy Correlation (EEC) function in electron-positron ($e^+e^-$) annihilation is a fundamental observable for testing Quantum Chromodynamics (QCD). However, a comprehensive global analysis faces several theoretical and phenomenological challenges:

• Perturbative Breakdown: In the back-to-back (two-jet) region ($\chi \to \pi$), fixed-order perturbative QCD expansions fail due to large Sudakov logarithms (double logarithms of infrared origin).
• Non-Perturbative (NP) Effects: Realistic descriptions require modeling hadronization (NP) effects, which are power-suppressed but significant, especially at lower energies.
• Energy Evolution: Previous analyses were often limited to the $Z$-boson resonance ($\sqrt{s} \approx 91.2$ GeV). A global fit across a wide energy range ($7.7 \text{ GeV} \leq \sqrt{s} \leq 91.2 \text{ GeV}$) is necessary to disentangle energy-dependent NP effects from the strong coupling constant, $\alpha_S$.
• Heavy Quark Masses: The impact of heavy quark (specifically bottom) mass effects on the Sudakov form factor and the overall fit quality had not been systematically included in a global EEC analysis.

2. Methodology

The authors developed a theoretical framework combining high-order perturbative calculations with an analytic dispersive model for non-perturbative corrections.

Theoretical Framework

• Resummation: In the back-to-back region, logarithmically enhanced contributions are resummed to Next-to-Next-to-Next-to-Leading Logarithmic (N3LL) accuracy.
• Fixed-Order Matching: The resummed results are consistently matched to fixed-order calculations up to Next-to-Next-to-Leading Order (NNLO, $\mathcal{O}(\alpha_S^3)$).
• Formalism:
  • The calculation is performed in impact-parameter ($b$) space to handle transverse momentum conservation and Sudakov logarithms.
  • The distribution is factorized into a hard function $H(\alpha_S)$ and a Sudakov form factor $S(b, Q)$.
  • The $b_*$ prescription is used to regulate the Landau singularity at low $b$ (large impact parameters), freezing the impact parameter at $b_{max}$.
• Non-Perturbative Model:
  • An analytic dispersive approach (based on the Dokshitzer-Marchesini-Webber model) is used to parameterize NP effects.
  • The NP Sudakov form factor $S_{NP}$ includes a term $f(b)$ (energy-independent hadronization) and a term $g_K(b)$ encoding the Collins-Soper evolution kernel (energy-dependent evolution).
  • The functional forms are:
    $$f(b) = f_1 b + f_2 b^2$$
    $$g_K(b) = g_0 \left\{ 1 - \exp\left[ -\frac{C_F \alpha_S(b_0^2/b_*^2)}{\pi g_0} \frac{b^2}{b_{max}^2} \right] \right\}$$
• Heavy Quark Masses: Finite mass effects for the bottom quark are incorporated into the Sudakov form factor, distinguishing between regions where the bottom mass is negligible and where it is relevant.

Data and Fitting Procedure

• Datasets: The analysis includes 691 data points spanning $\sqrt{s}$ from 7.7 GeV to 91.2 GeV.
  • Novel Inclusions: For the first time in a global EEC fit, data from the ALEPH collaboration (reanalyzed archived data at 91.2 GeV) and the AMY collaboration (58.0 GeV) are included.
  • Legacy Data: Includes data from OPAL, DELPHI, L3, SLD, TOPAZ, TASSO, PLUTO, JADE, CELLO, MARKII, and MAC.
• Fitting Strategy:
  • A simultaneous fit is performed for $\alpha_S(m_Z^2)$ and NP parameters ($f_1, f_2, g_0$).
  • The fit covers the peak and intermediate regions ($1.8 \leq \chi \leq \pi$) to ensure numerical stability.
  • Uncertainty Estimation:
    • Experimental/Theoretical correlations are handled via a nuisance parameter method.
    • Bootstrap method (1000 replicas) is used to estimate uncertainties from experimental data.
    • Perturbative uncertainties are estimated by varying renormalization ($\mu_R$) and resummation ($\mu_Q$) scales within $[Q/2, 2Q]$.

3. Key Contributions

1. First Global N3LL+NNLO Fit: This is the first global analysis of EEC data spanning over one order of magnitude in energy using N3LL resummation matched to NNLO.
2. Inclusion of New Data: The integration of high-precision, reanalyzed ALEPH data and previously unused AMY data provides unique constraints on the energy evolution of the EEC distribution.
3. Extraction of the Collins-Soper Kernel: By analyzing data at different energies simultaneously, the authors successfully disentangle the energy-dependent component of the NP form factor, allowing for a novel extraction of the Collins-Soper evolution kernel from $e^+e^-$ data.
4. Analytic Dispersive Approach: The paper demonstrates that an analytic dispersive model, augmented with an explicit $Q$-dependent evolution term, provides a superior description of data across all energies compared to methods relying solely on Monte Carlo tuning.

4. Results

• Strong Coupling Constant: The analysis yields a precise determination of the strong coupling constant at the $Z$-boson scale:
  $$\alpha_S(m_Z^2) = 0.119 \pm 0.002$$
  This value is fully consistent with the current PDG global average.
• Non-Perturbative Parameters: The best-fit values for the NP parameters are:
  • $f_1 = 1.0 \pm 0.2 \text{ GeV}$
  • $f_2 = 0.5 \pm 0.2 \text{ GeV}^2$
  • $g_0 = 0.9 \pm 0.3$
    These parameters characterize the hadronization effects and the Collins-Soper kernel evolution.
• Fit Quality: The theoretical framework describes the entire dataset (691 points) with a $\chi^2/\text{d.o.f.} = 1.2$, indicating excellent agreement across all energies, from the $Z$-peak down to 7.7 GeV.
• Collins-Soper Kernel: The extracted kernel $K(b, Q)$ at $Q=2$ GeV shows good agreement with determinations from Drell-Yan processes and recent Lattice QCD calculations, validating the universality of the kernel.
• Heavy Quark Mass Effects: Including bottom mass effects shifts the central values of the NP parameters but does not significantly improve the $\chi^2$. This suggests that perturbative mass effects can be effectively absorbed into the NP parameters, though a full treatment including mass effects in the fixed-order remainder is recommended for future precision.

5. Significance

• Precision QCD: The result $\alpha_S(m_Z^2) = 0.119 \pm 0.002$ represents a high-precision extraction of the strong coupling constant, competitive with other major determinations, achieved through a rigorous treatment of resummation and NP effects.
• Universality Test: The successful extraction of the Collins-Soper kernel from $e^+e^-$ annihilation (a process free from initial-state hadron PDFs) provides a clean, complementary test of the kernel's universality compared to SIDIS and Drell-Yan processes.
• Methodological Advancement: The study validates the necessity of including energy-dependent NP terms (via the Collins-Soper kernel) for global fits. It demonstrates that analytic dispersive models are robust tools for separating perturbative and non-perturbative dynamics without relying on the biases of general-purpose Monte Carlo event generators.
• Future Projections: The paper provides theoretical predictions for EEC at higher energies (LEP2 range: 133–206 GeV), suggesting that future measurements could further constrain the energy dependence of non-perturbative effects.

In conclusion, this work establishes a new benchmark for the analysis of EEC data, combining state-of-the-art perturbative QCD with a sophisticated analytic NP model to extract fundamental QCD parameters and test the universality of evolution kernels across different energy scales.