Quark and gluon distributions at the earliest stage of heavy ion collision

Using the general framework of quantum field kinetics we consider new principles to compute initial distribution of quarks and gluons after the first hard interaction of heavy ions. We start by rewriting the integral equations of QCD in the form which is generalizations of the familiar QCD evolution equations. These equations describe both space-time-- and $(x,Q^2)$–evolution before the collision, and allow one to use the $ep$ DIS data without reference to parton phenomenology. New technique generate perturbation theory that avoid double count of the processes, does not contain an artificial factorization scale, and does not require low-momentum cut-offs since infrared behavior is controlled by the DIS data.