CGC and saturation approach: Impact-parameter dependent model of perturbative QCD and combined HERA data

This paper presents an impact-parameter dependent perturbative QCD model based on the Color Glass Condensate framework, utilizing an analytical solution to the Balitsky-Kovchegov equation and a Froissart-consistent exponential saturation momentum, which successfully describes a wide range of combined HERA data at small-$x$ and offers a reliable foundation for future high-energy experiments like the Electron-Ion Collider.

The Gist

Imagine the proton (a tiny particle inside an atom) not as a solid marble, but as a bustling, chaotic city filled with invisible messengers called gluons. When you smash a high-speed electron into this proton city, you are essentially throwing a probe into a storm. The goal of this paper is to build a better weather map to predict exactly how that storm behaves.

Here is a simple breakdown of what the authors did, using everyday analogies:

The Problem: The "Too Fast" Storm

Physicists have a theory called the Color Glass Condensate (CGC). Think of this as a rulebook for how the proton's "gluon storm" behaves when you hit it with high energy.

• The Old Rulebook: Previous versions of this theory were like a weather model that predicted the storm would get infinitely bigger and stronger the faster you drove your car. But when they checked the real data from the HERA particle collider (a giant microscope that smashed electrons into protons), the storm didn't behave that way. It was too wild.
• The Missing Piece: The old models also treated the proton like a perfect circle with a fuzzy edge that faded out like a Gaussian bell curve (a smooth, symmetrical hill). The authors realized this was wrong. In reality, the "edge" of the proton's influence should fade away more sharply, like a light dimming exponentially, to obey the laws of physics (specifically a rule called the Froissart theorem).

The Solution: A New, Smarter Map

The authors took the existing CGC rulebook and fixed two major issues to make it match the real-world data from HERA:

1. The "Non-Linear" Fix: They used a more advanced mathematical solution (the Balitsky-Kovchegov equation) that accounts for the fact that the gluons in the proton aren't just flying around; they are bumping into each other and merging. It's like realizing that in a crowded dance floor, dancers don't just move in straight lines; they bump, merge, and change the flow of the crowd. This "non-linear" correction stops the storm from growing too fast.
2. The "Shape" Fix: Instead of using the old "smooth hill" shape for the proton's edge, they used a new shape that fades away exponentially (like a light turning off quickly). This ensures the model respects the fundamental laws of physics regarding how big a collision can get.

The Experiment: Tuning the Radio

To make their new map work, the authors had to "tune the radio." They had four "dials" (parameters) they could adjust:

• How strong the interaction is at the start.
• How fast the storm grows as energy increases.
• The "size" of the proton's core.
• A mass scale related to how the gluons are confined.

They turned these dials while looking at the combined data from the H1 and ZEUS experiments (the two teams that ran the HERA collider). They kept turning the dials until their mathematical predictions matched the experimental data as closely as possible.

The Results: A Perfect Match

Once they found the right settings, they tested their new map against a wide variety of "weather patterns" (different types of particle collisions):

• Standard Collisions: They predicted how the proton's structure changes (the $F_2$ and $F_L$ functions).
• Heavy Quark Collisions: They predicted what happens when the collision creates heavy "charm" particles.
• Exclusive Collisions: They predicted what happens when the proton stays intact but creates a new particle (like a $J/\psi$ meson) or bounces a photon off it (DVCS).

The Outcome: The new model fit the data incredibly well across a huge range of energies. It was like their weather forecast predicted the rain, wind, and temperature perfectly, even for storms they hadn't seen before.

Why It Matters (According to the Paper)

The authors claim this is a major step forward because:

• It is based on solid theoretical ground (first principles) rather than just guessing shapes.
• It respects the "Froissart theorem," a fundamental law of physics that the old models violated.
• It provides a reliable tool for predicting what will happen in future, even more powerful colliders like the Electron-Ion Collider and the LHeC.

In short: The authors took a theoretical model of how protons behave at high speeds, fixed its shape and its internal logic, tuned it to match real-world data, and found that it now predicts particle collisions with high accuracy. They believe this is the best path forward for understanding the fundamental forces of nature at the highest energies.

Technical Summary

Technical Summary: CGC/Saturation Approach with Impact-Parameter Dependence

Problem Statement
The primary challenge addressed in this work is the reliable description of non-linear evolution in high-energy scattering processes within Quantum Chromodynamics (QCD). While the leading-order Color Glass Condensate (CGC) approach provides a framework, it predicts an energy dependence for the scattering amplitude and saturation momentum that grows unchecked at high energies, contradicting experimental data and the Froissart theorem. Furthermore, standard saturation models often fail to reproduce the correct behavior of the scattering amplitude at large impact parameters ($b$), typically assuming a Gaussian dependence ($Q_s^2 \propto \exp(-b^2/B)$) which violates theoretical constraints regarding analyticity and unitarity at high energies. The authors aim to confront a specific CGC/saturation approach (Ref. [1]) with combined HERA experimental data to determine its phenomenological parameters and validate its ability to describe deep inelastic scattering (DIS) and exclusive diffractive processes across a wide kinematic range.

Methodology
The authors employ a dipole model based on the CGC/saturation effective theory for high-energy QCD. The methodology involves the following key components:

1. Analytical Solution to the BK Equation: The model utilizes an analytical solution to the non-linear Balitsky-Kovchegov (BK) evolution equation at Next-to-Leading Order (NLO) in the BFKL kernel. This solution incorporates re-summation procedures to fix the BFKL kernel and introduces non-linear evolution to ensure the correct high-energy asymptotic behavior of the scattering amplitude.
2. Impact Parameter Dependence: A critical feature of the model is the treatment of the impact parameter ($b$) dependence. Instead of the common Gaussian parametrization, the saturation momentum $Q_s$ is parametrized to exhibit an exponential decay at large $b$:
   $$Q_s^2(b, Y) \propto (S(b, m))^{1/\bar{\gamma}} \quad \text{where} \quad S(b) \sim (mb K_1(mb))$$
   This leads to $Q_s^2 \propto \exp(-mb)$ for $b \gg 1/m$, ensuring the scattering amplitude falls off exponentially at large distances, consistent with the Froissart theorem. This behavior also yields a power-like decrease in the scattering amplitude at high momentum transfers ($Q_T$), aligning with perturbative QCD.
3. Fitting Procedure: The model is fitted to the combined reduced inclusive DIS cross-section ($\sigma_r$) data from the H1 and ZEUS collaborations. The fit is restricted to the kinematic range $0.85 \text{ GeV}^2 < Q^2 < 30 \text{ GeV}^2$ and $x \leq 10^{-2}$.
4. Parameters: Four phenomenological parameters are determined via $\chi^2$ minimization:
   • $\bar{\alpha}_S$: The effective strong coupling constant.
   • $N_0$: The value of the scattering amplitude at the saturation scale ($\tau=1$).
   • $Q_0^2$: The initial saturation scale at $Y=0, b=0$.
   • $m$: The mass scale governing the impact parameter dependence (related to the gluon mass scale).
   • Light quark masses ($m_{u,d,s}$) and charm quark mass ($m_c$) are also varied to assess stability.

Key Contributions

• Theoretical Consistency: The paper demonstrates that using an exponential impact parameter dependence ($Q_s^2 \propto \exp(-mb)$) rather than a Gaussian one resolves the violation of the Froissart theorem and ensures consistency with perturbative QCD at large momentum transfers.
• Analytical Solution: The work applies a specific analytical solution to the non-linear BK equation that incorporates NLO corrections, moving beyond simple ansatzes.
• Comprehensive Data Comparison: The model is tested not only against inclusive structure functions ($F_2, F_L, F_2^{c\bar{c}}$) but also against exclusive diffractive processes, including vector meson production ($J/\psi, \phi, \rho$) and Deeply Virtual Compton Scattering (DVCS).

Results

• Inclusive Processes: The model achieves good agreement with HERA data for the proton structure function $F_2$, the longitudinal structure function $F_L$, and the charm structure function $F_2^{c\bar{c}}$ over a wide kinematic range ($Q^2 \in [0.85, 120] \text{ GeV}^2$ and $x < 10^{-2}$). The $\chi^2/d.o.f$ values are generally close to unity (e.g., 1.239 for $F_2$ vs $x_B$ with $m_c=1.4$ GeV).
• Extrapolation: The model successfully predicts data outside the fitting range, including high $Q^2$ values where $F_2^{c\bar{c}}$ was not included in the fit, and extends to $x > 10^{-2}$ with robust agreement until $x \approx 10^{-1}$.
• Exclusive Processes: The model provides a reasonable description of exclusive vector meson production and DVCS. However, the fit for the $\phi$ meson production shows higher $\chi^2$ values (e.g., 2.48 for $\sigma$ vs $Q^2$), indicating a potential discrepancy. The slope parameter $B_D$ of the $t$-distribution is reproduced, though the model faces challenges in the DVCS sector where the predicted energy dependence of the slope deviates slightly from data.
• Parameter Sensitivity: The fit is more sensitive to light quark masses than to the charm mass. The extracted value for $\bar{\alpha}_S$ is approximately $0.10 - 0.11$, and the parameter $m$ is found to be around $0.54$ GeV, consistent with the gluon mass scale.
• Limitations: The agreement is less satisfactory for $Q^2 < 0.85 \text{ GeV}^2$, likely due to the lack of a solid theoretical foundation for photon wave functions in the non-perturbative regime. Additionally, the model exhibits a high correlation among its four free parameters, which complicates the extraction of unique values.

Significance and Claims
The authors claim that their results provide a "strong guide" for developing a Color Glass Condensate/saturation effective theory approach capable of making reliable predictions from first principles. They emphasize that the model's ability to describe a wide range of observables (inclusive and exclusive) using a single set of parameters, while adhering to fundamental theoretical constraints like the Froissart theorem, validates the underlying approach.

The paper positions this model as a necessary step toward understanding high-energy QCD dynamics for future experiments, specifically citing the Electron-Ion Collider (EIC) and the LHeC. It also suggests applications for predicting dijet production in $p-p$ and $p-Pb$ collisions and understanding the soft Pomeron structure. The authors conclude that this approach, grounded in the non-linear BK evolution and correct impact parameter behavior, represents a viable path for extending QCD predictions to higher energies, including those accessible at the LHC, without relying on purely phenomenological ansatzes that violate theoretical bounds.