Employing the so-called hybrid formalism, we calculate the cross section of inclusive hadron production in proton-proton collisions at forward rapidity in small-$x$ formalism at one-loop order. For the case of hadron production at forward rapidity, we can uses collinear parton distributions for projectile proton and $k_\perp$ dependent gluon distribution for target proton. We show that collinear divergences associated with initial and final state parton radiations are renormalized into parton distributions and fragmentation functions in terms of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equation, respectively. Furthermore, rapidity divergence can be absorbed into the wave function of target proton which gives rise to the well-known Balitsky-Fadin-Kuraev-Lipatov equation. These divergences are completely separated from the short distance partonic hard parts, which is now finite at the next-to-leading order accuracy. The result presented in this paper can be reckoned as a baseline calculation without any non-linear QCD effects in small-$x$ formalism. As a consistency check, we compare our results with the previous calculation for non-linear proton-nucleus collisions in the small-$x$ formalism and find complete agreement in the dilute and large $N_c$ limit. In phenomenology, the direct comparison of the above two separate calculations can reveal the role and strength of the non-linear dynamics in high energy QCD, and thus help us reliably study the onset of gluon saturation when genuine non-linear interactions become important.
Read full abstract