  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.
