Heavy-ion collisions at RHIC

1. Heavy-ion collisions at RHIC Jean-Yves Ollitrault, IPhT Saclay IPN Orsay, December 16, 2009

2. Experiments: The relativistic heavy-ion programme RHIC at Brookhaven : Two beams of atomic Au or Cu nuclei, accelerated at energies up to 100 GeV per nucleon (since 2000) LHC at CERN: Two beams of atomic Pb nuclei will be accelerated at energies up to 2.7 TeV per nucleon (2010?)

3. Theory:The phase diagram of strong interactions Can we learn anything about hot gauge theories by smashing heavy ions together? Do heavy-ion collisions have anything to do with temperature and thermodynamics?

4. Phases of a heavy-ion collision at RHIC • Lorentz contraction: colliding pancakes • Strongly-interacting matter is formed: this is the nonperturbative regime of QCD. • If interactions are strong enough, it may reach local thermal equilibrium. Then, it expands like a fluid. • Eventually, one measures yields, single-particle spectra and correlations between identified hadrons. • Can we make sense out of all this?

5. Outline • Heavy-ion chemistry and the phase diagram • Microscopic modeling of heavy-ion collisions: how much nuclear structure do we need? • Elliptic flow and the « perfect fluid » estimating the viscosity of hot QCD from RHIC data • Higher harmonics and flow fluctuations open issues

6. Heavy-ion chemistry Hadron yields in heavy-ion collisions are well described by statistical models : grand-canonical ensemble for an ideal gas of hadrons at temperature T, baryon chemical potential μ, volume V: 3 parameters only • Technical subtleties: • particle identification is over a limited range in momentum: must be extrapolated using a reasonable spectrum • only stable particles (with respect to strong interactions) are detected : include « feeddown » from resonance decays Andronic et al nucl-th/0511071

7. Heavy-ion chemistry Using data at different energies, the statistical model gives a line in the (T,μ) plane. The temperature at RHIC is compatible with the deconfinement temperature predicted by lattice QCD calculations (can be done at μ=0 only) The « hadronization curve » might tell us about the phases of QCD Andronic at al 0911.4806

8. Modeling heavy-ion collisions • Inputs: • nuclear density profile(from electron scattering experiments) • nucleon-nucleon inelastic cross sectionσ(pp data + smooth interpolation) Miller et al nucl-ex/0701025

9. Monte-Carlo Glauber calculations • Assume independent nucleons (interactions and Pauli blocking neglected except for possible hard-core interactions) • Draw randomly positions of nucleons according to the nuclear density • Treat nucleons as black disks of surface σ. A collision occurs between projectile and target nucleons if their disks overlap in transverse space • . Nucleons which collide at least once (filled disks) are participants Miller et al nucl-ex/0701025

10. Why we need Glauber calculations They are used to estimate experimentally the centrality of the collision from the multiplicity Nch of particles seen in detectors, assuming Nch=a * (number of participants) +b * (number of nucleon-nucleon collisions) And also less trivial information about the size and shape of the fireball… Miller et al nucl-ex/0701025

11. RHIC matter as an expanding fluid Fluid dynamics = only well-controlled description of an expanding, strongly-interacting system. Macroscopic description: systematic gradient expansion ∂μ((є+P)uμuν-Pgμν) + η∂μ∂ρ…+ ∂μ∂ρ ∂σ … = 0 ideal viscous 2nd order terms ~1/R ~1/R2where R = system size A large system expands as an ideal fluid Is a Au (or Cu) nucleus large enough ? Not sure… Or, equivalently, is the viscosity ηsmall enough ?

12. Viscosity and strong coupling The stronger the interactions, the smaller the viscosity η: in transport theory, ηscales like the mean free path λ So η may well be small in strongly-coupled QCD… but lattice QCD is not yet able to give a reliable estimate of η An exact result for N=4 supersymmetricgauge theories at strong coupling, obtained using the gauge/string duality: η/s=1/4π(s=entropy density) Kovtun Son Starinets hep-th/0405231 How close is the QCD ηto the string theory prediction?

13. A specific observable: elliptic flow φ Non-central collision seen in the transverse plane: the overlap area, where particles are produced, is not a circle. A particle moving at φ=π/2 from the x-axis is more likely to be deflected than a particle moving at φ=0, which escapes more easily. Initially, particle momenta are distributed isotropically in φ. Collisions result in more particles near φ=0 or π(elliptic flow)

14. Anisotropic flow Fourier series expansion of the azimuthal distribution: Using the φ→-φ and φ→φ+π symmetries of overlap area: dN/dφ=1+2v2cos(2φ)+2v4cos(4φ)+… v2=<cos(2φ)> (<…> means average value) is elliptic flow v4=<cos(4φ)> is a (much smaller) « higher harmonic » higher harmonics v6, etc are 0 within experimental errors.

15. Elliptic flow at RHIC One of the first observables measured at RHIC (2nd physics paper): elliptic flow is as large as predicted by ideal hydro

16. Eccentricity scaling y v2 scales by construction like the eccentricityε of the initial density profile, defined as : x ε depends on the collision centrality • Ideal fluid dynamics ∂μTμν=0 is scale invariant : xμ→λxμ • This implies: v2/εindependent of system size R • Viscous corrections are expected to scale like -η/R.

17. Viscous hydro results for elliptic flow Lines: fits using (1+α/R)-1, α= fit parameter ~ system size R Masui et al arxiv:0908.0403 Deviations to v2/ε=constant clearly scale like η/R A direct way of estimating the viscosity η experimentally !

18. Testing the scaling on RHIC data Masui et al arxiv:0908.0403 η/s~0.15-0.3 at RHIC Caveat: v2 is measured, ε comes from a model!

19. x Nucleus 1 y Nucleus 2 b Participant Region Estimating the initial eccentricity Until 2005, this was thought to be the easy part. But puzzling results came: 1. v2 was larger than predicted by hydro in central Au-Au collisions. 2. v2 was much larger than expected in Cu-Cu collisions. This was interpreted by the PHOBOS collaboration as an effect of fluctuations in initial conditions [Miller & Snellings nucl-ex/0312008] In 2005, it was also shown that the eccentricity depends significantly on the model chosen for initial particle production. We compare two such models, Glauber and Color Glass Condensate.

20. Eccentricity fluctuations Depending on where the participant nucleons are located within the nucleus at the time of the collision, the actual shape of the overlap area may vary: the orientation and eccentricity of the ellipse defined by participants fluctuates. Assuming that v2 scales like the eccentricity, eccentricity fluctuations translate into v2 (elliptic flow) fluctuations

21. ε in Monte-Carlo Glauber calculations The eccentricity of the ellipse defined by the participants, εpart, is much larger than the standard eccentricity εstd computed from a smooth density profile (esp. for smaller systems) Miller et al nucl-ex/0701025

22. Anisotropic flow in hydrodynamics • Ideal gas (weakly-coupled particles) in global thermal equilibrium. The phase-space distribution is (Boltzmann) dN/d3pd3x = exp(-E/T) Isotropic! • A fluid moving with velocity v is in (local) thermal equilibrium in its rest frame: dN/d3pd3x = exp(-(E-p.v)/T) Not isotropic: Momenta parallel to v preferred • At RHIC, the fluid velocity depends on φ: typically v(φ)=v0+2ε cos(2φ)

23. The simplicity of v4 • Within the approximation that particle momentum p and fluid velocity v are parallel (good for large momenta) dN/dφ=exp(2ε p cos(2φ)/T) • Expanding to order ε, the cos(2φ) term is v2=ε p/T • Expanding to order ε2, the cos(4φ) term is v4=½ (v2)2 Hydrodynamics has a universal prediction for v4/(v2)2 ! Should be independent of equation of state, initial conditions, centrality, particle momentum and rapidity, particle type

24. PHENIX results for v4 PHENIX data for charged pions Au-Au collisions at 100+100 GeV 20-60% most central The ratio is independent of pT, as predicted by hydro. But… the value is significantly larger than 0.5 5

25. More data on v4: centrality dependence Au-Au collision 100+100 GeV per nucleon Data > hydro Small discrepancy between STAR and PHENIX data

26. Why fluctuations change v4/(v2)2 The reference directions x and y are not known experimentally: thus v2=<cos(2φ)> and v4=<cos(4φ)> are not measured directly v2 from 2-particle correlations: <cos(2φ1-2φ2))>=<(v2)2> v4 from 3-particle correlations: <cos(4φ1-2φ2-2φ3))>=<v4 (v2)2> If v2 and v4 fluctuate, the measured v4/(v2)2 is really <v4 (v2)2>/<(v2)2>2. Inserting the prediction from hydrodynamics, [v4/(v2)2]exp=½ <(v2)4>/ <(v2)2>2

27. Data versus eccentricity fluctuations Eccentricity fluctuations can be modelled using a Monte-Carlo program provided by the PHOBOS collaboration: Throw the dice for the positions of nucleons, with probability given by the nuclear density: No free parameter! Gombeaud JYO arxiv:0907.4664 Fluctuations explain most of the discrepancy between data and hydro, except for central collisions which suggest <(v2)4>/ <(v2)2>2=3 By symmetry, v2=0 for central collisions, except for fluctuations!

28. Summary • Fluid dynamics is the most promising approach to describe the bulk of particle production in heavy-ion collisions. • In 2005, the « perfect liquid » picture was popular. Now, we have learned that viscous corrections are important : elliptic flow is 25 % below the «ideal hydro », even for central Au-Au collisions ! • The phenomenology of heavy-ion collisions has made a lot of progress in the last few years: we now have quantitative relations between observables and properties of hot QCD. • It is crucial to have a better control on the initial stage of the collision: although thermalization loses memory of initial conditions, so that the detailed structure of the initial state is lost, but nontrivial information remains through fluctuations. Fluctuations tell us about the (partonic or nucleonic) structure of incoming nuclei, and are not yet understood.

29. Backup slides

30. Dimensionless numbers in fluid dynamics They involve intrinsic properties of the fluid (mean free path λ, thermal/sound velocity cs, shear viscosity η, mass density ρ)as well as quantities specific to the flow pattern under study (characteristic size R, flow velocity v) Knudsen number K=λ/R K «1 : local equilibrium (fluid dynamics applies) Mach number Ma= v/cs Ma«1 : incompressible flow Reynolds number R= Rv/(η/ρ) R»1 : non-viscous flow (ideal fluid) They are related ! Transport theory: η/ρ~λcsimpliesR * K ~ Ma Remember: compressible+viscous = departures from local eq.

31. A simple (non-standard)test of ideal-fluid behaviour Ideal hydrodynamics∂μTμν=0is scale invariant : xμ→λxμ is still a solution. Scale invariance can be tested by varying the size of the colliding nuclei (Cu-Cu versus Au-Au) or the centrality of the collision (peripheral<central) Pt spectra are centrality independent to a good approximation, and this is also true in ideal hydro.

32. How well do we know the initial density profile? This is the initial eccentricity obtained assuming a smooth density profile. It is strongly model-dependent! We need a good control of initial conditions Drescher Dumitru Hayashigaki Nara, nucl-th/0605012

33. PHOBOS data for v2 • Phobos data for v2 • εobtained using Glauber or CGC initial conditions +fluctuations • Fit with • v2=αε/(1+1.4 K) • assuming • 1/K=(σ/S)(dN/dy) • with the fit parameters σ and α. K~0.3 for central Au-Au collisions v2 : 30% below ideal hydro!