We study BFKL evolution and, in particular, the energy dependence of the saturation momentum in the presence of saturation boundaries limiting the region of linear BFKL evolution. In the case of fixed coupling evolution we confirm the previously found exponential term in Q s (Y ) and determine the prefactor Y and α dependences. In the running coupling case we find Y 1/6 corrections to the Y 1/2 exponential behavior previously known. Geometrical scaling of the scattering amplitude is valid in a wide range of momenta for fixed coupling evolution and in a more restricted region for running coupling evolution.

We show that the Balitsky-JIMWLK equations proposed to describe non-linear evolution in QCD at high energy fail to include the effects of fluctuations in the gluon number, and thus to correctly describe both the low density regime and the approach towards saturation. On the other hand, these fluctuations are correctly encoded (in the limit where the number of colors is large) in Mueller's color dipole picture, which however neglects saturation. By combining the dipole picture at low density with the JIMWLK evolution at high density, we construct a generalization of the Balitsky hierarchy which includes the particle number fluctuations, and thus the pomeron loops. After an additional coarse-graining in impact parameter space, this hierarchy is shown to reduce to a Langevin equation in the universality class of the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov (sFKPP) equation. This equation implies that the non-linear effects in the evolution become important already in the high momentum regime where the average density is small, which signals the breakdown of the BFKL approximation.