Nonlinear corrections to the nuclear heavy flavor structure functions
F. Abdi, B. Rezaei
TL;DR
The paper investigates nonlinear corrections to nuclear gluon and heavy-quark structure functions at small $x$ by numerically solving the GLR-MQ-ZRS evolution equation for $G^A(x,Q^2)$. It extends the collinear high-energy DAS framework to nuclei, computing $F_2^{c\bar c-A}$ and $F_L^{c\bar c-A}$ via convolution with nonlinear gluon distributions and includes LO/NLO coefficient functions, with results benchmarked against $nCTEQ15$, RGK, and H1/ZEUS data. Key findings show nonlinear recombination effects become significant at $x\lesssim 10^{-4}$ and low $Q^2$, growing with the nucleus size $A$ and depending on the correlation length $R$, while remaining compatible with existing nuclear PDF analyses and providing robust predictions for LHeC/EIC kinematics. The work advances saturation phenomenology in nuclei, informs global nPDF fits, and offers concrete predictions for future electron-ion collider measurements where nonlinear QCD dynamics are expected to be most pronounced.
Abstract
We study numerically the small-$x$ behaviour of the nuclear gluon distribution function $ G^A(x,Q^2)$ at next-to-leading order (NLO) approximation of the Gribov-Levin-Ryskin-Mueller-Qiu, Zhu-Ruan-Shen (GLR-MQ-ZRS) nonlinear equation and quantify the impact of gluon recombination in the kinematic range of $x \le 10^{-2}$ and $Q^2 \ge 5 \text{GeV}^2$ respectively. The results are comparable to the Rausch-Guzey-Klasen (RGK) [J.Rausch, V.Guzey and M.Klasen, Phys. Rev. D 107, 054003 (2023)] nuclear gluon distributions and the nCTEQ15 parametrization at the corresponding $Q^2$ values. Using the solutions of $ G^A(x,Q^2)$ in the framework of the nonlinear GLR-MQ-ZRS evolution equation, the linear and nonlinear behavior of the charm structure functions of nuclei per nucleon $F^{c\bar{c}-A}_2(x,Q^2)$ and $F^{c\bar{c}-A}_L(x,Q^2)$ are considered. The results reveal that nonlinear corrections play an important role in charm nuclear reduced cross sections at small-$ x$ and low $Q^2$ values. The computed results are compared with experimental data from the H1 and ZEUS Collaborations.
