Probing Saturon-like Limits in QCD Systems

High-occupancy QCD matter enters a saturated regime when its entropy or occupancy approaches the unitarity bound $\sim 1/α$, the ``saturon" criterion. We test this criterion for protons and nuclei at small $x$ using analytic and numerical solutions of the BK equation. From these solutions we construct the gluon occupancy $N_g(x)$ and a thermodynamic entropy $S(x)$ via an Unruh-like temperature $T = Q_s/(2π)$ and an emergent gluon mass $M_g \sim Q_s$. For protons, both $N_g$ and $S$ rise toward small $x$ yet stay below $1/α_s$ in our baseline setup. For nuclei, by contrast, the nuclear entropy $S_A$ attains the $1/α_s$ benchmark in a small-$x$ window where the proton does not. This singles out nuclei as the natural environment to search for saturon-like behavior and motivates precision small-$x$ measurements and high-occupancy $pA$ and $AA$ collisions.

💡 Research Summary

**
The manuscript investigates whether high‑occupancy QCD matter can reach the “saturon” limit, a theoretical bound in which the entropy of a system saturates at S ≈ 1/α (the inverse of the coupling). The authors focus on the small‑x regime of protons and nuclei, where gluon densities become large and non‑linear effects (gluon recombination) lead to saturation.

The central tool is the Balitsky‑Kovchegov (BK) equation, which governs the rapidity evolution of the dipole scattering amplitude N(k,Y). Both the integral‑differential form (solved numerically with the BKsolver algorithm) and the differential form (treated analytically via a mapping to the Fisher‑Kolmogorov‑Petrovsky‑Piscounov equation) are employed. Two coupling scenarios are considered: fixed coupling and running coupling. In the running‑coupling case the authors adopt a phenomenological “saturation‑scale prescription” (SSP) that identifies the scale of ᾱ_s with the saturation momentum Q_s, thereby incorporating the dynamical weakening of the coupling at higher Q_s.

From the solved amplitudes the unintegrated gluon distribution (UGD) is constructed, and the gluon occupation number N_g(x) is obtained by integrating over transverse momentum. An entropy measure is defined by introducing an Unruh‑like temperature T = Q_s/(2π) and an effective gluon mass M_g ≈ Q_s, leading to S(x) = (M_g/T) N_g = 2π N_g. This definition mirrors earlier proposals (e.g., Kutak) and ensures that when N_g reaches the bound 1/α_s, the entropy automatically saturates the saturon limit.

Numerical results show distinct behavior for protons and nuclei. For protons, both N_g and S increase as x decreases from 10⁻¹ to 10⁻⁷, but they remain well below the benchmark 1/α_s (≈ 5 for α_s ≈ 0.2). The proton’s saturation scale grows only modestly (Q_s ≈ 1 GeV), and gluon recombination prevents the system from reaching the maximal entropy.

In contrast, nuclei (especially heavy ones with A ≈ 200) start with a larger initial Q_s (scaled roughly as A^{1/3}) and therefore develop a much higher occupation. The calculated nuclear entropy S_A(x) approaches or slightly exceeds 1/α_s in the window x ≈ 10⁻⁵–10⁻⁶. This indicates that heavy nuclei can achieve the saturon bound in a realistic small‑x region, making them prime candidates for observing saturon‑like physics.

Sensitivity checks varying the momentum‑space saturation threshold (κ = 0.5–1.0) and the running‑coupling prescription show that the qualitative conclusion—nuclei reaching the saturon limit while protons do not—is robust against reasonable model variations.

The authors discuss experimental implications. The predicted entropy saturation in nuclei should manifest in observables accessible at future electron‑ion colliders (EIC) and in high‑multiplicity pA/AA collisions at the LHC or RHIC. Possible signatures include deviations from geometric scaling in the nuclear structure function F₂^A, enhanced diffractive cross sections, and specific patterns in multi‑gluon correlations that reflect a maximally entropic gluon ensemble.

Finally, the paper situates its findings within a broader theoretical context: the saturon bound mirrors the Bekenstein‑Hawking entropy of black holes, suggesting a universal “critical packing” principle that transcends gravity. While QCD dynamics differ from gravitation, both systems exhibit a maximal entropy proportional to the inverse coupling at the point of saturation. The work therefore provides a concrete, QCD‑based test of this universality and outlines clear avenues for future theoretical refinements (higher‑order BK corrections, impact‑parameter dependence) and experimental verification.