Abstract
Quark mass effects in the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution appear at the next to leading order (NLO) of QCD, i.e., at 𝒪(α2s) of the QCD coupling for the JIMWLK evolution kernels. Based on the light-cone hadronic wave function approach, we obtain the analytic expressions for the JIMWLK evolution kernels, which are expressed in terms of one dimensional integrals that in principle can be numerically useful for future phenomenological applications. In comparison to the situation with massless quarks, a new divergence arises. We illustrate the origin of this divergence and how the divergence is canceled within the renormalization procedures of light-cone quantization.