Classical initial conditions for ultrarelativistic heavy ion collisions

We construct an analytical expression for the distribution of gluons in the state immediately following a heavy ion collision in the quasi-classical limit of QCD given by the McLerran-Venugopalan model. The resulting gluon number distribution function includes the effects of all multiple rescatterings of gluons with the nucleons of both colliding nuclei. The typical transverse momentum k _⊥ of the produced gluons is shown to be of the order of the saturation scale of the nuclei Q _s, as predicted by Mueller. We analyze the properties of the obtained distribution and demonstrate that due to multiple rescatterings it remains finite (up to logarithms of k _⊥) in the soft transverse momentum limit of k _⊥≪Q _s, unlike the usual perturbative initial conditions given by collinear factorization. We calculate the total number of produced gluons and show that it is proportional to the total number of gluons inside the nuclear wave function before the collision with the proportionality coefficient c≈2 ln2 .