Numerical evaluation of the nonlinear Gribov-Levin-Ryskin-Mueller-Qiu evolution equations for nuclear parton distribution functions

Abstract

We numerically study for the first time the nonlinear GLR-MQ evolution equations for nuclear parton distribution function (nPDFs) to next-to-leading order accuracy and quantify the impact of gluon recombination at small x. Using the nCTEQ15 nPDFs as input, we confirm the importance of the nonlinear corrections for small x≲10−3, whose magnitude increases with a decrease of x and an increase of the atomic number A. We find that at x=10−5 and for heavy nuclei, after the upward evolution from Q0=2 GeV to Q=10 GeV, the quark singlet Ω(x,Q2) and the gluon G(x,Q2) distributions become reduced by 9%–15%, respectively. The relative effect is much stronger for the downward evolution from Q0=10 GeV to Q=2 GeV, where we find that Ω(x,Q2) is suppressed by 40%, while G(x,Q2) is enhanced by 140%. These trends propagate into the FA2(x,Q2) nuclear structure function and the FAL(x,Q2) longitudinal structure function, which after the downward evolution become reduced by 45% and enhanced by 80%, respectively. Our analysis indicates that the nonlinear effects are most pronounced in FAL(x,Q2) and are already quite sizable at x∼10−3 for heavy nuclei. We have checked that our conclusions very weakly depend on the choice of input nPDFs. In particular, using the EPPS21 nPDFs as input, we obtain quantitatively similar results.