Reconciling the kinematical constraint with the JIMWLK evolution equation: correlation functions non-local in rapidity

In the high-energy limit of DIS experiments the effective degrees of freedom of QCD are Wilson line operators. Their evolution in the rapidity variable is predicted by the set of Balitsky-JIMWLK evolution equations. We analyze a new class of two-point correlation functions of Wilson line operators where the Wilson lines are taken at different values of the rapidity variable. Such correlation functions can appear in the discussion of the kinematical constraint. We find that in the Langevin formulation of the JIMWLK equation, such correlation functions are affected by an infrared divergence. We discuss a possible regularization of this divergence and its consequences for the implementation of the kinematical constraint for the JIMWLK equation.