Small- x physics 2- The saturation regime of QCD and the Color Glass Condensate

1. Small-x physics2- The saturation regime of QCD and the Color Glass Condensate Cyrille Marquet Columbia University

2. the saturation regime: for with this regime is non-linear, yet weakly coupled magnitude of Qs x dependence the idea of the CGC is to take into account this effect via strong classical fields The saturation regime of QCD • the BFKL approximation summing terms leads to a growth of the gluon density in the hadronic wave function but high-density effects are missing when the gluon density is large enough, gluon recombination is important • gluon recombination in the hadronic wave function gluon density per unit area recombination cross-section recombinations important when gluon kinematics

3. The Color Glass Condensate • effective description of a high-energy hadron/nucleus the CGC is moving in the direction and the gauge is the current is solving Yang-Mills equations gives  short-lived fluctuations separation between the long-lived high-x partons and the short-lived low-x gluons CGC wave function in practice we deal with CGC averages such as we are interested in the evolution of with x, which sums both the BFKL equation is recovered in the low-density regime

4. Outline of the second lecture • The evolution of the CGC wave functionthe JIMWLK equation and the Balitsky hierarchy • A mean-field approximation: the BK equationsolutions: QCD traveling waves the saturation scale and geometric scaling • Beyond the mean field approximationstochastic evolution and diffusive scaling • Computing observables in the CGC frameworksolving evolution equation vs using dipole models

5. The evolution of theCGC wave function:the B-JIMWLK equations

6. the Wilson lines sum powers of adjoint representation the JIMWLK equation gives the evolution of the wave function for small enough x the equivalent Balitsky equations are obtained by considering the scattering of simple test projectiles (dipoles) off the CGC The JIMWLK equation Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov, Kovner • a functional equation for the x evolution of : with

7. eigenvalues scattering amplitude off the (dense) target CGC the interaction with the CGC conserves transverse position Scattering off the CGC • eigenstates ? their interaction should conserve their spin, polarisation, color, momentum … in the high-energy limit, the eigenstates are simple when the partons transverse momenta are Fourier transformed to transverse coordinates in the large-Nc limit, further simplification: the eigenstates are only made of dipoles x, y: transverse coordinates other quantum numbers aren’t explicitely written, the eigenvalues depend only on the transverse coordinates

8. Dipoles as test projectiles • dipoles are ideal projectiles to probe small distances the dipole: u: quark space transverse coordinate v: antiquark space transverse coordinate scattering of the quark: dependence kept implicit in the following JIMWLK equation → evolution equation for the dipole correlators

9. in the large Nc limit, the hierarchy is restricted to dipoles general structure: solving the B-JIMWLK equations gives , … which can then be used to compute observables BFKL saturation we will now derive the first equation of the hierarchy The Balitsky hierarchy • equations for dipoles scattering of the CGC Balitsky (1996) the BFKL kernel an hierarchy of equation involving correlators with more and more dipoles

10. the dipole dressed with one gluon color index from QCD at order gS the dipole rapidity: Recall the dipole wave function • the valence component (using the mixed space): x: quark space transverse coordinate y: antiquark space transverse coordinate such that only perturbative sizes r << 1/ QCD are included other degrees of freedom are not explicitely written

11. scattering of the antiquark scattering of the quark let’s compute the elastic scattering amplitude Ael(Y ) in the dipole-CGC collision where is the total rapidity Y = YCGC + Yd • in the frame in which and Elastic scattering of the dipole • the dipole scattering is described by Wilson lines for the qqg component:

12. = 1/2 -1/(2Nc) • frame independence: our two expressions for Ael(Y ) should be identical, this requirement in the limit Y - YCGC dY  0 gives the first Balitsky equation Frame independence of Ael(Y) • in the frame in which and using the following identity to get rid of the adjoint Wilson line one obtains

13. A mean field approximation: the Balitsky-Kovchegov equation

14. at small Y, is small, and the quadratic term can be neglected, the equation reduces then to the linear BFKL equation and rises exponentially with Y as gets close to 1 (the stable fixed point of the equation), the non-linear term becomes important, and , saturates at 1 with increasing Y, the unitarization scale get bigger The BK equation • obtain by neglecting correlations the BK equation is a closed equation for obtained by assuming • solutions: qualitative behavior r = dipole size let’s consider impact-parameter independent solutions

15. Coordinate vs momentum space • let’s go to momentum space  due to conformal invariance, the linear part of the equation is the same for and coordinate space momentum space dipole scattering amplitude unintegrated gluon distribution linear BFKL equation linear BFKL equation genuine saturation: no real saturation: both equation are in the same universality class as the F-KPP equation

16. The F-KPP equation Fisher, Kolmogorov, Petrovsky, Piscounov • same features as the BK equation time derivative space derivative linear term meaning exponential growth non linear term meaning saturation • dictionary F-KPP → BK when expanding to second order, the equations are the same (in momentum space) in spite of small difference, same universality class the precise forms of the space derivatives and of the non-linear term doesn’t matter these equations belong to the same universality class F-KPP and BK asympotic (in t or Y) solutions are the same: traveling waves

17. the formation of the traveling wave position: BK solutions: same rapidity evolution quantitative features can be derived saturation region wave traveling at speed v, independently of the initial condition provided it is steep enough there the initial condition has not been erased yet Traveling wave solutions • what is a traveling wave the speed of the wave is determined only by the linear term of the equation

18. Recall the BFKL solutions • a superposition of waves with speeds initial condition in Mellin space • the minimal speed is obtained for and its value is this will be the speed of the asymptotic traveling wave • the saturation exponent for the traveling wave to form, the initial condition must feature

19. same features for except in the saturation region Y >Y0 Y0 called geometric scaling • the saturation scale   with then QCD traveling waves Munier and Peschanski (2004) • the initial condition is steep enough in QCD • asymptotic solutions of BK

20. Numerical solutions • numerical simulations confirm the results for the scattering amplitude for the saturation scale sub-asymptotic corrections (≡ geometric scaling violations) are also known:

21. consequences similar to those first obtained by the simpler substitution running coupling corrections slow down the increase of Qs with energy also confirmed by numerical simulations, however this asymptotic regime is reached for larger rapidities Running coupling corrections • running coupling corrections to the BK equation taken into account by the substitution Kovchegov Weigert Balitsky

22. Beyond the mean fieldapproximation

23. still there is something missing describing properly the dilute regime is also important the B-JIMWLK equations do not describe it properly the equation for is fine, its dilute limit is the BFKL equation, however the other equations are incomplete, they do not agree with Mueller’s dipole model in the limit incorrect results for have consequences for Recall the Balitsky hierarchy • equations for dipole correlators it describes the so-called dense regime where the BK equation is a good approximation general structure at large Nc

24. reformulation into a Langevin problem the correlators can all be obtained from a single stochastic quantity which obeys a Langevin equation average over the CGC wave function average over the realisations of the stochastic process A new hierarchy of equations • the hierarchy must be completed Iancu and Triantafyllopoulos (2005) fluctuations (included to recover Mueller’s dipole model in the dilute limit) important when BFKL saturation important when this QCD evolution is equivalent to a stochastic process

25. the stochastic F-KPP equation noise r = dipole size the reduction to one dimension introduces the noise strength parameter  The stochastic F-KPP equation • the QCD Langevin equation deterministic part of the equation: BK ! noise high-energy QCD evolution = stochastic process in the universality class of reaction-diffusion processes, of the sF-KPP equation Iancu, Mueller and Munier (2005)

26. the Langevin evolution generates, at each rapidity, a stochastic ensemble of solutions for the average speed of the waves: analytic expressions for v and Dare known only for very small • the probability distribution of Qs Qsis a stochastic variable distributed according to a Gaussian probability law: (for ) average saturation scale: corrections to the Gaussian law for improbable fluctuations also known, picture confirmed by numerical simulations C.M., Soyez and Xiao A stochastic saturation scale • the stochastic evolution the dispersion coefficient: D equivalent to a stochastic saturation scale Qs

27. A new scaling law • the average dipole scattering amplitude : the diffusion is negligible and one recovers geometric scaling : the diffusion is important and new regime: diffusive scaling very high energies

28. The diffusive scaling regime • important properties in the diffusive scaling regime: the amplitudes are dominated by events that feature the hardest fluctuation of in average the scattering is weak, yet saturation is the relevant physics • running coupling corrections Dumitru et al. push the diffusive scaling regime up to asymptotically high energies • conclusion at the moment using the mean field approximation is good enough for phenomenology

29. Computing observablesin the CGC framework

30. the 2-point function or dipole amplitude this is the most common average, (for instance it determines deep inelastic scattering) its Fourier transform is the unintegrated gluon distribution, it obeys the BK equation • more complicated correlators for 2-particle correlations for instance The CGC averages • the Wilson line correlators in the CGC framework, any cross-section is determined by colorless combinations of Wilson lines, averaged over the CGC wave function the energy evolution of cross-sections is encoded in the evolution of this wave function is mainly non-perturbative, but its evolution is known  Balitsky/BK equation for the correlators

31. what is done in practice solving evolution equations with running coupling is not the favored approach instead one uses parametrizations with features inspired from theoretical results it makes it easier to deal with impact parameter dependence of the dipole amplitude famous dipole models: GBW, IIM, IPsat many processes have been studied in the CGC framework using dipole models The strategy • what should be done in principle - express the cross-sections in terms of correlators of Wilson lines the same correlators enter in the formulation of very different processes the more exclusive the final state is, the more complicated the correlators are - use B-JIMWLK or BK equation to compute the correlators use appropriate initial conditions include non-perturbative effects such as the impact parameter dependence running coupling effects should also be included

32. Outline of the third lecture • The hadron wave functionsummary of what we have learned • The saturation modelsfrom GBW to the latest ones • Deep inelastic scattering (DIS)the cleanest way to probe the CGC/saturationallows to fix the model parameters • Diffractive DIS and other DIS processesthese observables are predicted • Forward particle production in pA collisionsand the success of the CGC picture at RHIC