QCD Saturation in the Semi-classical Approach

In this paper the semi-classical approach to the solution of non-linear evolution equation is developed. We found the solution in the entire kinematic region to the non-linear evolution equation that governs the dynamics in the high parton density QCD. The large impact parameter ($b_t$) behavior of the solution is discussed as well as the way how to include the non-perturbative QCD corrections in this region of $b_t$. The geometrical scaling behavior and other properties of the solution in the saturation (Color Glass Condensate) kinematic domain are analyzed. We obtain the asymptotic behavior for the physical observables and found the unitarity bounds for them.