Geometric scaling and QCD evolution

We study the impact of the QCD Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution on the geometric scaling of gluon distributions that is expected to hold at small x within the saturation models. With this aim we solve the DGLAP evolution equations with the initial conditions provided along the critical line Q²=Q_s²(x) with Q_s²(x)∼x^-λ and satisfying geometric scaling. Both fixed and running coupling cases are studied. We show that in the fixed coupling case the geometric scaling at low x is stable against the DGLAP evolution for sufficiently large values of the parameter λ, and in the double logarithmic approximation of the DGLAP evolution this happens for λ≥4N_cα_s/π. In the running coupling case geometric scaling is found to be approximately preserved at very small x. The residual geometric scaling violation in this case can be approximately factored out and the corresponding form factor controlling this violation is found.