Classical initial conditions for ultrarelativistic heavy ion collisions Original Research Article

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 View the MathML source of the produced gluons is shown to be of the order of the saturation scale of the nuclei View the MathML source, 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 View the MathML source) in the soft transverse momentum limit of View the MathML source, 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 View the MathML source.