⚛️ phenomenology

Simultaneous Color Glass Condensate fit to deep inelastic scattering and forward hadron production at HERA, RHIC, and the LHC

This paper presents the first simultaneous Color Glass Condensate fit to HERA deep inelastic scattering and RHIC/LHC forward hadron production data using LO BK evolution with running coupling, demonstrating that a global fit with collider-specific K-factors achieves excellent agreement with data while revealing complementary constraints on dipole evolution and identifying fragmentation functions and factorization scales as dominant systematic uncertainties.

Original authors: Piotr Korcyl, Truong My Hau Le, Farid Salazar, Tomasz Stebel

Published 2026-07-28

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Piotr Korcyl, Truong My Hau Le, Farid Salazar, Tomasz Stebel

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 the universe is built from tiny, invisible Lego bricks called quarks and gluons. These aren't just sitting still; they are zooming around at nearly the speed of light, glued together by a force so strong it's like trying to pull apart two magnets that have been fused with superglue. This is the world of Quantum Chromodynamics (QCD), the rulebook for how matter holds itself together. But here's the twist: when you smash these particles together at incredibly high speeds, they don't just bounce off each other. Instead, the gluons multiply so wildly that they start to crowd each other out, creating a dense, saturated "soup" where the rules of the game change. Scientists call this state the "Color Glass Condensate." It's like a traffic jam where the cars (gluons) are so packed that they can't move freely anymore. Understanding this jam is a huge goal for physicists because it reveals how the universe behaves under extreme pressure, something we can only see in giant particle accelerators like the ones at CERN or Brookhaven.

Now, picture a team of scientists acting like cosmic detectives. They have two different crime scenes to investigate: one is a deep-sea dive into the proton's structure using electron beams (from the HERA collider), and the other is a high-speed crash test where protons smash into each other head-on (from the RHIC and LHC colliders). The paper you're about to read is the report from these detectives. They tried to solve a puzzle: Can one single, universal "map" of the proton's gluon soup explain what happens in both of these very different crash scenarios?

The team used a sophisticated mathematical tool called the Balitsky–Kovchegov (BK) equation to evolve their map. Think of this equation as a time-machine that predicts how the gluon crowd grows and changes as you zoom in or speed up. They tested two versions of this time-machine: one that accounts for the changing strength of the glue (running coupling) and another that adds extra safety rules to prevent the math from going crazy (kinematical constraints). To make their predictions match the real-world data, they had to introduce "K-factors." You can think of these as "fudge factors" or volume knobs. Since their math was a simplified version of reality (leading order), they needed to turn the volume up or down to match the loudness of the actual experimental signals.

Here is what they found: The team successfully built a single map that fits both the deep-sea electron data and the high-speed proton crash data simultaneously. It's a rare win in physics to get two different types of experiments to agree on the same underlying picture. They discovered that the "volume knob" (the K-factor) needed to be turned up much higher for the RHIC collisions than for the LHC collisions—about twice as high. This isn't a mistake; it actually confirms other theories that suggest lower-energy crashes need more "correction" because they are closer to the edge of what's physically possible.

The most exciting part is that the two types of data didn't fight each other. Instead, they helped each other. The proton crash data (SIHP) acted like a tightrope walker, forcing the map to evolve at just the right speed, while the electron data (DIS) kept the map's shape steady. Together, they narrowed down the uncertainties, making the picture of the gluon soup sharper than ever before. While there are still some small bumps in the road—like needing to figure out exactly how the "fudge factors" change with different settings—the study suggests that the Color Glass Condensate theory is on the right track. It proves that the same fundamental rules govern the proton whether it's being probed by an electron or smashed into another proton, bringing us one step closer to understanding the dense, chaotic heart of matter.

Technical Summary: Simultaneous Color Glass Condensate Fit to DIS and Forward Hadron Production

Problem and Motivation
The Color Glass Condensate (CGC) effective theory provides a framework for describing the high-energy limit of Quantum Chromodynamics (QCD), specifically the regime of gluon saturation characterized by a dynamically generated saturation scale $Q_s(x)$ . While the universality of the CGC suggests that a single dipole amplitude, evolved via the Balitsky–Kovchegov (BK) equation, should describe diverse observables across different collision systems, previous phenomenological efforts have largely been restricted. Traditionally, non-perturbative initial conditions are constrained by fitting HERA Deep Inelastic Scattering (DIS) data, with predictions subsequently made for other observables. Conversely, simultaneous fits have been limited to single collision systems (e.g., only DIS or only ultra-peripheral collisions). There has been a lack of attempts to describe two distinct observables—reduced DIS cross-sections and forward single inclusive hadron production (SIHP)—across vastly different colliding systems (electron-proton vs. proton-proton) and energy scales (HERA, RHIC, LHC) using a unified dipole amplitude.

Methodology
The authors perform the first simultaneous global fit of HERA reduced DIS cross-sections and forward SIHP data from RHIC (BRAHMS, STAR) and LHCb (5 TeV and 13 TeV). The analysis is conducted within the Leading Order (LO) CGC framework with the following key components:

Dipole Evolution: The dipole amplitude is evolved using the running-coupling BK (rcBK) equation. The authors investigate two schemes: rcBK with and without kinematical constraints (kcBK).
Initial Conditions: The initial dipole amplitude at $x_0 = 0.015$ is parameterized using the MV $\gamma$ model. Parameters fitted include the saturation scale $Q_{s0}^2$ , the anomalous dimension $\gamma$ , the proton transverse area $\sigma_0/2$ , and the coupling scale parameter $C_2$ .
Cross-Section Formulation:
- DIS: Calculated via the optical theorem using the forward scattering amplitude of a virtual photon on a proton, convolved with the dipole amplitude.
- SIHP: Calculated in the dilute-dense approximation (appropriate for forward production), where the cross-section is a convolution of parton distribution functions (PDFs), the Fourier-transformed dipole amplitude, and fragmentation functions (FFs).
Higher-Order Corrections (K-factors): Since the analysis is LO, constant K-factors are introduced to account for missing higher-order corrections. Crucially, separate K-factors are applied for RHIC ( $K_{\text{RHIC}}$ ) and LHC ( $K_{\text{LHC}}$ ) to account for the energy dependence of threshold resummation effects, which are expected to be larger at lower energies.
Statistical Analysis: Uncertainties are determined using the Hessian method, validated against a Monte Carlo Bayesian inference study. Theoretical uncertainties from fragmentation functions are propagated using the NNFF1.0 replica ensemble, and scale uncertainties are estimated via three-point scale variations.

Key Contributions

First Simultaneous Fit: This work establishes the first global fit of the BK-evolved dipole amplitude to both DIS and forward SIHP data across HERA, RHIC, and LHC energies.
Universality Test: It tests the universality of the dipole amplitude's shape and evolution speed across electron-proton and proton-proton collisions.
K-factor Analysis: The study quantifies the necessity and magnitude of K-factors, finding them to be approximately a factor of two larger at RHIC than at the LHC.
Complementarity of Observables: The analysis demonstrates that SIHP data provides complementary constraints to DIS data, specifically tightening the constraints on the evolution speed of the dipole amplitude without introducing tension with the DIS description.

Results

Global Fit Quality: Both evolution schemes (rcBK and kcBK) achieve a good simultaneous description of the data. The rcBK scheme yields a superior global $\chi^2/\text{d.o.f.} = 0.943$ , compared to $1.139$ for kcBK. The preference for rcBK is driven by better fits to the DIS sector and LHCb spectra.
K-factor Magnitude: The extracted K-factors are $K_{\text{RHIC}} \approx 5.6\text{--}6.1$ and $K_{\text{LHC}} \approx 3.0\text{--}3.3$ . The ratio $K_{\text{RHIC}}/K_{\text{LHC}} \approx 1.9$ is stable across schemes and aligns with threshold-resummed one-loop studies, suggesting that constant per-collider K-factors effectively capture the dominant energy-dependent corrections.
Parameter Constraints:
- The global fit shifts parameters slightly toward a faster evolution (smaller $C_2$ ) and a steeper initial amplitude (larger $\gamma$ and $Q_{s0}^2$ ) compared to DIS-only fits.
- The inclusion of SIHP data reduces the uncertainty on the evolution speed parameter $C_2$ by roughly 25%.
- No significant tension is observed between the two observables; the DIS description quality remains essentially unchanged ( $\chi^2/\text{d.o.f.}$ moves from 0.97 to 0.98) when SIHP data is included.
Systematic Uncertainties: The dominant systematic uncertainties in the SIHP sector are identified as the choice of fragmentation function sets and the factorization scale prescription. The K-factors are found to be largely uncorrelated with the dipole amplitude parameters, as the overall normalization is absorbed by the K-factors while the shape is constrained by the data.

Significance and Claims
The paper claims that a consistent, simultaneous description of DIS and forward hadron production is achievable within the LO CGC framework supplemented by per-collider K-factors. This result supports the universality of the BK-evolved dipole amplitude across different collision systems. The authors emphasize that while the SIHP data modestly improves constraints on the evolution speed, the current precision is limited by theoretical uncertainties in fragmentation functions and scale choices. Consequently, the constraining power of SIHP at this order should not be overstated. The work serves as a baseline for future extensions, including proton-nucleus collisions, the inclusion of charm production in DIS, and ultimately, a full Next-to-Leading Order (NLO) global analysis to further stress-test the CGC effective theory.

Technical Summary: Simultaneous Color Glass Condensate Fit to DIS and Forward Hadron Production

More like this