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(perpendicular to) 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(perpendicular to)) in the soft transverse momentum limit of k(perpendicular to) much less than 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 approximate to 2 ln 2. (C) 2001 Elsevier Science B.V. All rights reserved.