New results on small-x resummation for splitting functions

This paper presents new analytical all-order results for small-$x$ resummation, most notably the first properly resummed $qg$ splitting kernel, which form the basis of a more robust and numerically stable implementation in the upcoming HELL 4.0 framework.

The Gist

Imagine you are trying to predict the outcome of a massive, high-speed collision between two particles, like smashing two cars together at 99.9% the speed of light. To do this, physicists use a "recipe" called Parton Distribution Functions (PDFs). Think of these PDFs as the "ingredient lists" for the proton or electron. They tell you how much "stuff" (quarks, gluons, photons) is inside the particle and how fast that stuff is moving.

Usually, this recipe works great. But there's a tricky corner case: Small-x.

The Problem: The "Tiny Fraction" Nightmare

In the world of particle physics, "x" represents how much of the total speed a tiny piece of the particle is carrying.

• Large x: A piece carrying a big chunk of the speed (like a heavy truck in traffic). Easy to predict.
• Small x: A piece carrying a tiny, almost invisible fraction of the speed (like a speck of dust).

When you smash particles together at the extreme energies of future colliders (like a Muon Collider), you start looking at these "specks of dust." The problem is that the standard math used to predict their behavior starts to break down. It's like trying to use a ruler to measure the width of a single atom; the ruler is too coarse, and the numbers start to go haywire, becoming unstable and infinite.

This happens because of "logarithms"—mathematical terms that get huge when x gets tiny. If you ignore them, your prediction is wrong. If you try to calculate them one by one, the math explodes. You need to resum them, which means summing up an infinite number of these terms at once to get a stable answer.

The Old Solution: A "Best Guess" Approximation

For years, physicists had a way to handle this small-x problem, implemented in a computer code called HELL. However, this method was like using a crude sketch to draw a masterpiece.

• They knew the first few "brushstrokes" (mathematical coefficients) of the solution.
• They used a clever trick called Borel-Padé to guess what the rest of the painting looked like based on those few strokes.
• It worked okay for normal conditions, but it was fragile. If you pushed the physics to extreme limits (like very high energies where the strong force gets very strong), the sketch would start to warp, giving nonsensical results.

The New Breakthrough: The "Master Blueprint"

This paper, by Marco Bonvini, Stefano Frixione, and Giovanni Stagnitto, is about throwing away the sketch and drawing the exact, all-order blueprint.

Here is what they achieved, using some analogies:

1. Solving the "Impossible" Equation
The core of the problem is a complex equation that describes how these tiny particles interact. Previously, no one could solve this equation perfectly; they could only approximate it.

• The Analogy: Imagine trying to find the exact shape of a shadow cast by a complex 3D object. Before, people were guessing the shape based on a few light sources.
• The New Result: The authors found a way to solve the equation exactly for all orders of complexity. They didn't just guess the next few terms; they found the mathematical "formula" that generates the entire infinite series.

2. The "Green Function" (The Master Key)
In their math, there's a specific function called the "Green function" (specifically for the transition from quarks to gluons, or $P_{qg}$).

• The Analogy: Think of this as the master key to the whole lock. Previously, they had a key that was cut roughly (the approximation). It opened the door most of the time, but sometimes it got stuck or broke the lock when you turned it too hard (high energy).
• The New Result: They forged a perfectly cut, all-metal master key. They derived an exact, closed-form expression for this function. This means they can now calculate the behavior of these particles with perfect precision, without relying on shaky guesses.

3. Handling the "Strong Force" at Extreme Levels
The paper was motivated by the need to study Muon Colliders. Muons are heavy cousins of electrons. When you smash them together, the energy is so high that the "strong force" (the glue holding quarks together) becomes incredibly intense.

• The Analogy: Imagine driving a car. At normal speeds (low energy), the engine runs smoothly. But if you floor it to 300 mph (high energy/strong coupling), the engine might overheat or behave unpredictably. The old math (the sketch) couldn't handle the engine overheating.
• The New Result: Their new blueprint is robust enough to handle the engine running at 300 mph. They improved the underlying math so that even when the "strong force" gets huge, the predictions remain stable and physical.

Why Does This Matter?

1. Future Colliders: If we build a Muon Collider (a machine that could reach energies far beyond what we have today), we need these new, precise recipes to know what we will see. Without them, the data would be a mess.
2. Better Predictions for Today: Even for current colliders like the Large Hadron Collider (LHC), this new math makes the "ingredient lists" (PDFs) more accurate, especially for rare events where tiny particles play a role.
3. Stability: The old method was like a house of cards; a small change in the input could make the whole thing collapse. The new method is like a steel structure; it stands firm even under extreme pressure.

Summary

The authors took a messy, approximate mathematical problem that physicists have struggled with for decades. They didn't just tweak the numbers; they found the exact, underlying mathematical structure that governs these tiny particles. They replaced a "best guess" sketch with a perfect blueprint, allowing us to predict the behavior of matter at the most extreme energies imaginable with confidence. This new code (HELL version 4.0) will be the new standard for physicists trying to understand the universe at its smallest scales.

Technical Summary

Here is a detailed technical summary of the paper "New results on small-x resummation for splitting functions" by Bonvini, Frixione, and Stagnitto.

1. Problem Statement

The paper addresses two interconnected challenges in Quantum Chromodynamics (QCD) and Parton Distribution Function (PDF) evolution:

• Small-$x$ Resummation Instability: At high energies (small momentum fraction $x$), fixed-order DGLAP evolution becomes unstable due to large logarithmic terms ($\ln x$). While Next-to-Leading Logarithmic (NLL) resummation exists, previous implementations (specifically in the HELL code) relied on approximate methods. These approximations, particularly for the $qg$ (quark-gluon) splitting kernel, suffered from numerical instabilities and a lack of rigorous all-order analytical control.
• Large-$\alpha_s$ Regime: The motivation for this work stems from the physics of future muon colliders (and high-energy lepton colliders). Unlike hadron colliders, muon PDF evolution starts at a scale $\mu \sim m_\mu$, which is below $\Lambda_{QCD}$. To evolve PDFs to high scales, the strong coupling $\alpha_s$ must be extrapolated to large values (up to $\alpha_s \sim 0.8$). Previous resummation frameworks were optimized for small $\alpha_s$ and exhibited unphysical behaviors (e.g., splitting functions decreasing with decreasing $x$) when $\alpha_s$ became large.

2. Methodology

The authors revisited the theoretical foundations of small-$x$ resummation within the ABF (Altarelli-Ball-Forte) framework, implemented in the HELL code. Their approach involved:

• Analytical Derivation of All-Order Functions: Instead of relying on numerical extraction of coefficients or Borel-Padé approximations, the authors derived exact, closed-form analytical expressions for key all-order functions.
  • They solved the duality relation between the DGLAP anomalous dimension $\gamma_+$ and the BFKL kernel $\chi$ to obtain new representations for the leading-logarithmic (LL) eigenvalue $\gamma_s$.
  • They utilized the collinear factorization of the quark Green function in the high-energy limit to derive an exact all-order expression for the function $h_{qg}$, which governs the $qg$ anomalous dimension.
• Resummation of Subleading Effects: They implemented a rigorous resummation of subleading running-coupling effects (terms proportional to $\beta_0$) using an integral representation involving the evolution operator $U$ and the function $h_{qg}$.
• Numerical Improvements for Large $\alpha_s$:
  • They modified the approximation of the fixed-order anomalous dimension used to resum BFKL collinear poles to ensure stability at large $\alpha_s$.
  • They introduced a shift parameter $\eta$ in the running-coupling approximation to eliminate a "spurious pole" that caused unphysical behavior at large $\alpha_s$.
  • They developed a procedure to handle branch-cut changes in the complex $N$-plane (Mellin space) that occur at large $\alpha_s$, using Chebyshev polynomial extrapolation to maintain momentum conservation and numerical stability.

3. Key Contributions

The paper presents several novel analytical and numerical results:

• Exact All-Order $h_{qg}$ Function: The most significant contribution is the derivation of a closed-form, all-order expression for the function $h_{qg}(z)$ (Eq. 3.24), which relates the $qg$ anomalous dimension to the LL eigenvalue $\gamma_s$.
  • Previous status: $h_{qg}$ was known only via a finite number of perturbative coefficients, requiring Borel-Padé resummation which was numerically fragile.
  • New status: The exact form allows for a stable, parameter-independent calculation of the $P_{qg}$ splitting kernel.
• New Representations for $\gamma_s$: The authors provided multiple new analytical representations for the LL anomalous dimension $\gamma_s$, including:
  • A perturbative series with closed-form coefficients.
  • Integral representations (Borel and Hankel integrals).
  • A fixed-point iterative solution.
• Finite Green Functions: They derived the first all-order expressions for the finite ($\epsilon \to 0$) $qg$ and $gg$ Green functions, establishing a direct relationship between them (Eq. 5.5).
• HELL 4.0 Implementation: These results form the core of the upcoming HELL 4.0 version, which is designed to be robust for both hadronic and leptonic (muon) PDF evolution.

4. Results

The authors validated their new framework through extensive numerical comparisons:

• Stability of $P_{qg}$: Comparisons between the new exact $h_{qg}$ implementation and the old Borel-Padé approximation (used in HELL 3) showed that while the old method worked well for moderate $\alpha_s$, it was highly sensitive to the choice of the Mellin inversion path. The new exact result is path-independent and numerically stable.
• Large $\alpha_s$ Behavior:
  • The old HELL 3 implementation exhibited unphysical behavior (splitting functions decreasing as $x$ decreases) when $\alpha_s > 0.35$ due to a spurious pole.
  • The new implementation, utilizing the modified running-coupling approximation (with parameter $\eta$), correctly reproduces the expected growth of splitting functions at small $x$ even for $\alpha_s$ up to 0.5.
• Singlet Matrix: The paper presents the full singlet splitting function matrix ($P_{gg}, P_{qg}, P_{gq}, P_{qq}$) matched to NNLO with NLL' resummation. The results show the expected "dip" in $P_{gg}$ at small $x$ and a milder rise for $P_{qg}$ compared to fixed-order predictions.
• Uncertainty Bands: The uncertainty bands for the resummed results widen significantly at large $\alpha_s$, indicating that while the method is stable, higher-order corrections become increasingly important in the large-$\alpha_s$ regime.

5. Significance

• Muon Collider Physics: This work is a prerequisite for precise predictions at future muon colliders. It provides the first robust framework to handle the simultaneous challenges of small-$x$ resummation and large-$\alpha_s$ evolution required for muon PDFs.
• Hadronic Applications: Although motivated by muon colliders, the improvements (exact $h_{qg}$, stable large-$\alpha_s$ handling) significantly enhance the reliability of small-$x$ resummation for standard hadron collider physics (e.g., LHC, FCC-hh), particularly in regimes where $\alpha_s$ is large or precision requires control over subleading terms.
• Theoretical Advancement: The derivation of exact all-order expressions for $h_{qg}$ and the finite Green functions resolves long-standing issues regarding the convergence and stability of small-$x$ resummation, moving the field from approximate numerical extrapolations to rigorous analytical control.

In summary, the paper transforms the small-$x$ resummation framework from a set of approximate, numerically fragile tools into a robust, analytically exact, and numerically stable engine capable of handling the extreme kinematic regimes of future high-energy lepton colliders.