We present the numerical solution of the nonlinear evolution equation for DIS on nuclei for x=10-2-10-7. We demonstrate that the solution to the nonlinear evolution equation is quite different from the Glauber-Mueller formula which was used as the initial condition for the equation. We illustrate the energy profit for performing DIS experiments on nuclei. However, it turns out that the gain is quite modest: xAu≃5xproton for the same parton density. We find that the saturation scale Q2s∝A1/3. For gold the saturation scale Qs,Au≃1.5 GeV at x=10-3. Such a large value leads to considerable contribution of the high-density QCD phase to RHIC data and reveals itself in essential damping for both xGA and F2A.
- Deep inelastic structure functions
- Nonlinear evolution equation
- Perturbative QCD
- Saturation scale