In this paper, we generalized the work of Ref. [1]}, exploring
high energy evolution with non-vanishing quark masses. We computed the hadronic light-cone wave function of a fast moving projectile, focusing on the effects originating from non-vanishing quark masses. As an application of the wave function, we derived the NLO evolution kernels for the
JIMWLK Hamiltonian, particularly the kernels \(K^m_{q\bar{q}}\)  and \(K^m_{JSJ}\) . Both enter
the dipole evolution at NLO [2]} through the form presented in [3]}.
