
Single Inclusive Hadron Production in DIS at Small x: Next to Leading Order
Corrections
Filip Bergabo1, 2, ∗and Jamal Jalilian-Marian1, 2, †
1Department of Natural Sciences, Baruch College, CUNY,
17 Lexington Avenue, New York, NY 10010, USA
2City University of New York Graduate Center, 365 Fifth Avenue, New York, NY 10016, USA
We calculate the one-loop corrections to single inclusive hadron production in Deep Inelastic
Scattering (DIS) at small xin the forward rapidity region using the Color Glass Condensate for-
malism. We show that the divergent parts of the next to leading order (NLO) corrections either
cancel among each other or lead to x(rapidity) evolution of the leading order (LO) dipole cross
section according to the JIMWLK evolution equation and DGLAP evolution of the parton-hadron
fragmentation function. The remaining finite parts constitute the NLO (αs) corrections to the LO
single inclusive hadron production cross section in DIS at small x.
I. INTRODUCTION
Gluon saturation [1, 2] at small xas encoded in the Color Glass Condensate (CGC) formalism [3–7] has been the
subject of intense theoretical studies and experimental searches. Theoretical work based on leading order (LO) or
leading log (LL) approximations to gluon saturation have successfully described structure functions, suppression of
the single inclusive hadron transverse momentum spectrum and disappearance of the away side peak in dihadron
angular correlations in high energy proton(deuteron)-gold/lead collisions at RHIC and the LHC [8–54]. Nevertheless
firmly establishing gluon saturation as the QCD dynamics responsible for these experimental observations requires
more precise theoretical calculations. The ongoing work on improving the accuracy of leading order CGC calculations
can be broadly put into three categories; higher order in αscorrections to leading order results [55–80], sub-eikonal
corrections which aim to relax the infinite energy assumption inherent to eikonal approximation [81–92], and inclusion
of intermediate/large xdynamics into CGC in order to generalize CGC to include DGLAP evolution and collinear
factorization and high ptphysics [93–101]. Here we will focus on next to leading order corrections to single inclusive
hadron production in Deep Inelastic Scattering (DIS) at small xin the forward rapidity region [102] (virtual photon
going direction) for the case when the virtual photon is longitudinal. We note that leading order results for single
inclusive hadron production in DIS in the midrapidity region were obtained in [103].
The ideal environment in which to investigate gluon saturation and CGC is DIS experiments at high energy as the
incoming virtual photon does not interact strongly. Single inclusive hadron production in DIS (SIDIS) at small xis
one of the most attractive channels for gluon saturation studies as it is not sensitive to Sudakov effects which can
obscure saturation dynamics in dihadron production and angular correlations. Furthermore, it is more discriminatory
than the total cross section (structure functions) so that it contains more information about the QCD dynamics of
the target. While there exists leading order calculations of single inclusive hadron production in DIS at small xin
the CGC framework [102, 103] it is highly desirable and in fact urgently needed to perform a next to leading order
calculation which can then be used for quantitative studies of the transverse momentum spectra of produced hadrons
in DIS with proton and nuclear targets at the proposed Electron Ion Collider (EIC).
Here we calculate the next to leading order corrections to single inclusive hadron production in DIS at small xin
the forward rapidity region using the Color Glass Condensate formalism. To do so we use our recent results for next
to leading order corrections to dihadron production [80] in DIS and integrate out one of the final state partons. As
expected we encounter various divergences which appear when we integrate over the phase space of the final state
parton. We show that UV and soft divergences cancel among each other while the collinear divergences associated
with radiation of a massless parton are absorbed into the parton-hadron fragmentation function. We show that all
quadrupole terms appearing in the intermediate steps of the calculation cancel among various terms and one is left
with dipoles (and squared dipoles) only. The rapidity divergences arising from integrating over longitudinal phase
space of the final state parton are absorbed into evolution of the dipoles describing the target dynamics and lead
to JIMWLK evolution of the leading order cross section. The remaining terms are finite and constitute the O(αs)
corrections to leading order single inclusive hadron production in DIS at small x.
∗fbergabo@gradcenter.cuny.edu
†jamal.jalilian-marian@baruch.cuny.edu
arXiv:2210.03208v3 [hep-ph] 20 Dec 2022