Small x dynamics of the unpolarised color dipole gluon TMD PDFs for all transverse momenta

This paper derives a closed-form expression for the unpolarized color dipole gluon TMD distribution across all transverse momenta by Fourier-Bessel transforming a general solution of the Balitsky-Kovchegov equation, revealing a characteristic inversion of x-ordering at the saturation scale that serves as a model-independent signature of gluon saturation.

The Gist

Imagine a proton not as a solid marble, but as a bustling, chaotic city made of tiny particles called gluons. These gluons are the "glue" holding the proton together, but when the proton is moving incredibly fast, the behavior of these gluons changes dramatically.

This paper is like a team of physicists (Mariyah, Nahid, and Mushood) trying to draw a perfect, single map of this gluon city. Their goal is to understand how these gluons are distributed not just in how fast they move forward, but also in how they wiggle side-to-side (their "transverse momentum").

Here is the story of their discovery, broken down into simple concepts:

1. The Problem: A City with Two Different Rules

For a long time, scientists had to use two different rulebooks to describe the gluon city, depending on how "crowded" it was:

• The Empty City (Low Density): When gluons are spread out, they behave like independent travelers. Scientists had a good map for this.
• The Packed City (High Density/Saturation): When you zoom in or look at very high speeds, the gluons get so crowded they start bumping into each other and merging. This is called "saturation." In this zone, the old maps broke down, and scientists had to use a completely different, complicated set of rules.

The big problem was that no one had a single, smooth map that worked for the whole city, from the empty suburbs to the packed downtown. Previous attempts were like stitching two different maps together with a jagged seam in the middle.

2. The Solution: A Master Key (The BK Equation)

The authors found a "Master Key" called the Balitsky–Kovchegov (BK) equation. Think of this equation as a mathematical recipe that describes how the gluon city grows and changes as you speed up.

While others had only figured out parts of this recipe, these authors used a general solution (a complete version of the recipe) that works everywhere. They treated the gluons as a "color dipole" (a pair of particles acting like a tiny antenna) and asked: If we send this antenna through the proton city, how does it scatter?

3. The Magic Trick: Turning the Map Inside Out

To get their final map, they performed a mathematical "magic trick" called a Fourier-Bessel transform.

• Imagine you have a blurry photo of a city taken from far away (the "dipole size").
• This trick converts that blurry photo into a sharp, high-definition map of the traffic flow (the "gluon momentum").

They did the math and found something surprising: The messy, infinite numbers that usually pop up in these calculations (divergences) simply vanished. It was as if the universe itself cancelled out the errors, leaving behind a clean, perfect formula.

4. The Result: The "One-Size-Fits-All" Map

They produced a single, elegant equation (Equation 13 in the paper) that describes the gluons perfectly across the entire spectrum. Here is what the map shows:

• The Deep Downtown (Low Momentum): When gluons are very slow and the city is super-packed, the number of gluons drops off sharply. It's like a "Sudakov suppression"—a force that keeps the city from collapsing under its own weight.
• The Peak (The Saturation Boundary): As you move out from the center, the number of gluons rises to a distinct, smooth peak. This is the "busy hour" of the proton.
• The Suburbs (High Momentum): As you go further out, the number of gluons falls off smoothly, like a gentle hill.

5. The "Time Travel" Surprise (The x-Ordering Inversion)

The most fascinating part of their map is how it changes when you look at the proton at different speeds (represented by a variable called x).

• Before the Peak: If you look at the "slow" gluons, the proton looks "fuller" when you are moving slower (higher x).
• After the Peak: But once you pass the peak and look at the "fast" gluons, the rule flips! The proton looks "fuller" when you are moving faster (lower x).

The authors call this a "characteristic inversion." It's like walking through a crowd: from the front, the people look close together; but if you run past them, the people at the back suddenly look like they are rushing toward you faster than the people in front. This "crossing" behavior is a unique fingerprint of gluon saturation.

6. Why This Matters for the Future

The paper mentions that this new map is crucial for the Electron-Ion Collider (EIC), a massive new machine being built to take pictures of protons and nuclei.

• Because this map is smooth and unified, scientists won't have to guess where to switch between different rulebooks.
• It allows them to measure the "size" of the proton's gluon cloud with much higher precision.
• It confirms that the "inversion" effect is a real, universal feature of nature, not just a quirk of one specific model.

In summary: These physicists found a single, smooth mathematical formula that perfectly describes how gluons are packed inside a proton, from the densest core to the outer edges. They proved that the "rules" of the proton change in a specific, predictable way as you speed up, providing a clear guide for future experiments to explore the hidden structure of matter.

Technical Summary

Technical Summary: Small-x Dynamics of the Unpolarised Color Dipole Gluon TMD

Problem Statement
The internal structure of hadrons at small Bjorken-$x$ is governed by a rapid increase in gluon density, leading to non-linear QCD effects and the transition from a dilute to a dense gluon regime (saturation). While Transverse-Momentum-Dependent distributions (TMDs) are crucial for describing parton correlations in momentum space, a complete, unified analytical description of the unpolarized color dipole gluon TMD, $xF^{(1)}(x, k_\perp)$, valid across the entire transverse momentum spectrum ($k_\perp$) has remained an open problem. Existing analytical results are strictly piecewise: they describe the perturbative tail ($k_\perp \gg Q_s$), the saturation boundary ($k_\perp \sim Q_s$), and the deep saturation region ($k_\perp \ll Q_s$) using different approximations (e.g., McLerran-Venugopalan, Levin-Tuchin). There is a lack of a single closed-form expression that smoothly interpolates between these regimes without artificial matching, which is essential for precision phenomenology at the Electron-Ion Collider (EIC).

Methodology
The authors address this by employing a general analytic solution to the Balitsky-Kovchegov (BK) equation, derived in Ref. [33], which is valid across the full kinematic range of saturation dynamics.

1. S-Matrix Ansatz: The study utilizes a dipole S-matrix ansatz expressed in terms of the dilogarithm function ($Li_2$):
   $$S = S_0 \exp\left(\tau Li_2\left(-\lambda r_\perp^2 Q_s^2\right)\right)$$
   This form is constructed to reproduce the Gaussian behavior of the McLerran-Venugopalan model in the dilute limit and the quadratic logarithmic dependence of the Levin-Tuchin solution in the black-disc limit.
2. Fourier-Bessel Transform: The unpolarized gluon TMD is defined via the Fourier transform of the dipole S-matrix (Eq. 1). The authors substitute the general S-matrix ansatz into this integral.
3. Analytical Evaluation: The resulting two-dimensional Fourier integral is evaluated using methods involving Bell polynomials and the expansion of Gamma function ratios.
4. Regularization and Expansion: The authors decompose the result into divergent and finite parts. They demonstrate that the potential ultraviolet divergence (a $1/\eta$ pole arising from the Taylor series expansion) vanishes identically due to a Gaussian suppression factor. The finite part is then systematically expanded in the leading-logarithm (LL) approximation.

Key Contributions and Results
The primary contribution of the work is the derivation of a unified, closed-form analytical expression for the unpolarized color dipole gluon TMD (Eq. 13):
$$xF^{(1)}(x, k_\perp) = \frac{S_\perp N_c}{2\pi^3 \alpha_s} \exp\left[-\frac{\tau}{2} \ln^2\left(\frac{k_\perp^2}{4\lambda Q_s^2}\right)\right] \left[1 + \frac{\tau \gamma_E k_\perp^2}{2\lambda Q_s^2}\right]$$
where $\tau \approx 0.2$, $\lambda = 7.22$, and $\gamma_E$ is the Euler-Mascheroni constant.

The paper establishes the following properties of this result:

• Unified Validity: The expression is valid for all $k_\perp$, smoothly interpolating between the deep saturation, saturation boundary, and perturbative regimes without piecewise matching.
• Reproduction of Limits:
  • In the deep saturation region ($k_\perp \ll Q_s$), the linear correction is negligible, and the distribution reduces to a pure log-Gaussian (Sudakov-like) form, recovering the Levin-Tuchin result.
  • In the perturbative tail ($k_\perp \gg Q_s$), the linear correction dominates, yielding a behavior consistent with the MV and GBW models and JIMWLK evolution.
• Geometric Scaling: The distribution depends on $k_\perp$ and the saturation scale $Q_s$ exclusively through the ratio $k_\perp^2/Q_s^2$, confirming the geometric scaling property.
• Numerical Behavior: Numerical evaluations for EIC-relevant $x$ values ($0.01, 0.001, 0.0001$) show a smooth, peak-shaped distribution. Crucially, the results exhibit a characteristic "inversion of x-ordering" at the saturation scale: while higher $x$ distributions are larger in the pre-saturation region, lower $x$ distributions overtake them beyond the saturation peak and grow more steeply.

Significance and Claims
The authors claim that this result provides a model-independent signature of gluon saturation. The derived expression:

1. Resolves the Piecewise Problem: It eliminates the systematic uncertainties inherent in matching different analytical solutions across kinematic boundaries.
2. Connects Evolution Schemes: The double-logarithmic suppression structure in the deep saturation region parallels the Sudakov factor in the Collins-Soper-Sterman (CSS) framework, suggesting a deep connection between non-linear small-$x$ evolution and soft-gluon resummation physics.
3. Phenomenological Utility: The unified expression is directly applicable to processes probing the full $k_\perp$ spectrum, including inclusive and diffractive Deep Inelastic Scattering (DIS), direct photon-jet correlations, and forward dijet production in $pA$ and $eA$ collisions, as well as Semi-Inclusive DIS (SIDIS) at the EIC.
4. Validation: The numerical results, particularly the $x$-ordering inversion, align closely with the unintegrated color dipole gluon distribution computed within the Bartels–Golec-Biernat–Kowalski (BGK) model, reinforcing the physical reliability of the derived analytic form.

The paper concludes that this unified analytic expression constitutes a timely contribution to the precision program of small-$x$ gluon tomography at the EIC, offering a robust tool for extracting the saturation scale and non-perturbative parameters from future experimental data.