Hydrodynamics flow and flow fluctuations JeanYves Ollitrault IPh

Hydrodynamics, flow, and flow fluctuations Jean-Yves Ollitrault IPh. T-Saclay Hirschegg 2010: Strongly Interacting Matter under Extreme Conditions International Workshop XXXVIII on Gross Properties of Nuclei and Nuclear Excitations January 19, 2010 In collaboration with Clément Gombeaud and Matt Luzum

Outline • Elliptic flow and v 4 • A simple, universal prediction from ideal hydrodynamics • Comparison with data from RHIC • v 4 in ideal and viscous hydrodynamics • Flow fluctuations, their effect on v 4 • Conclusion Gombeaud, JYO, ar. Xiv: 0907. 4664, Phys. Rev. C 81, 014901 (2010) Luzum, Gombeaud, JYO, in preparation

Initial azimuthal distribution of particles Random parton -parton collisions occurring on scales << nuclear radius. No preferred direction in the production process. Isotropic azimuthal distribution

Final azimuthal distribution of particles Bar length = number of particles in the direction = Azimuthal (φ) distribution plotted in polar coordinates φ (for pions with transverse momentum ~ 2 Ge. V/c) Strong elliptic flow created by pressure gradients in the overlap area

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

Azimuthal distribution without v 4 A small effect: Average value 0. 3%, maximum value 3% Should we care? The beauty is in the details!

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

The simplicity of v 4 • Within the approximation that particle momentum p and fluid velocity v are parallel (valid for large momentum pt and low freeze-out temperature T) d. N/dφ=exp(2ε pt cos(2φ)/T) • Expanding to order ε, the cos(2φ) term is v 2=ε pt/T • Expanding to order ε 2, the cos(4φ) term is v 4=½ (v 2)2 Hydrodynamics has a universal prediction for v 4/(v 2)2 ! Should be independent of equation of state, initial conditions, centrality, particle momentum and rapidity, particle type Borghini JYO nucl-th/0506045

PHENIX results for v 4 PHENIX data for charged pions Au-Au collisions at 100+100 Ge. V 20 -60% most central The ratio is independent of p. T, as predicted by hydro. The ratio is also independent of particle species within errors. But… the value is significantly larger than 0. 5. Can detailed (ideal or viscous) hydro calculations explain this?

v 4 and coalescence • The pt range where we have data for v 4 is the range where quark coalescence is thought to be important • Coalescence predicts, with n_q=2 or 3 constituent quarks (d. N/dφ)hadron(pt)=(d. N/dφ)qn_q (pt/n_q) • Our simple hydrodynamical picture (d. N/dφ)=exp(2ε pt cos(2φ)/T) is stable under quark coalescence • In particular, v 4/(v 2)2 is closer to ½ for hadrons than for parent quarks • Discrepancy data/hydro is worse for the parent quarks: coalescence does not help here.

v 4 in viscous hydrodynamics Viscosity changes the fluid evolution and distorts the momentum distribution of particles emitted at freeze-out. fviscous(p)=e-E/T(1+δf(p)) where δf=Cχ(p) (pipj/p 2 -δij/3)∂iuj and χ(p) depends on the microscopic interactions and C is a normalization fixed by matching with the fluid Tμν Most calculations use χ(p)=p 2 (quadratic ansatz) but it has been recently pointed out that other choices are possible such as χ(p)=p (linear ansatz) Dusling Moore Teaney ar. Xiv: 0909. 0754

Results from ideal and viscous hydro Hydro parameters tuned to fit spectra and v 2 at RHIC (Luzum and Romatschke) In particular, Tf=140 Me. V Ideal hydro predicts a flat ratio as expected Viscous hydro with a linear ansatz is also OK Viscous hydro with (usual) quadratic ansatz fails badly at large pt. M. Luzum, work in progress Hydro is unable to explain a ratio larger than 0. 5. We need something more

More data : centrality dependence Au-Au collision 100+100 Ge. V per nucleon STAR: Yuting Bai, Ph. D thesis Utrecht PHENIX: Roy Lacey, private communication Data > hydro Small discrepancy between STAR and PHENIX data

Estimating experimental errors v 2 and v 4 are not measured directly but inferred from azimuthal correlations (more later on this). There are many sources of correlations (jets, resonance decays, …): this is the « nonflow » error which we can estimate (order of magnitude only) Difference between STAR and PHENIX data compatible with non-flow error How do we understand the discrepancy with hydrodynamics? ?

Eccentricity scaling We understand elliptic flow as the consequence of the almond shape of the overlap area It is therefore natural to expect that v 2 scales like the eccentricity ε of the initial density profile, defined as : (this is confirmed by numerical hydro calculations) y x

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 v 2 scales like the eccentricity, eccentricity fluctuations translate into v 2 (elliptic flow) fluctuations

We need fluctuations to understand v 2 results (see next talk by Raimond Snellings) Results using various methods (STAR) After correcting for fluctuations and nonflow JYO Poskanzer Voloshin, PRC 2009

Eccentricity fluctuations in central collisions Central collisions are azimuthally symmetric, except for fluctuations: In the most central bin, v 2 and v 4 are all fluctuations! Eccentricity fluctuations are to a good approximation Gaussian in the transverse plane (2 -dimensional Gaussian distribution). This implies <ε 4>=2 <ε 2>2 The value of v 4 for central collisions rather suggests <(v 2)4>~3 <(v 2)2>2. It is therefore unlikely that elliptic flow fluctuations are solely due to fluctuations in the initial eccentricity. We need ideas.

Conclusions • The fourth harmonic, v 4, of the azimuthal distribution gives a further, independent indication that the matter produced at RHIC expands like a relativistic fluid • v 4 is mostly induced by v 2 as a second order effect. • v 4 may help us constrain models based on viscous hydrodynamics, in particular viscous corrections at freeze-out : standard quadratic ansatz ruled out? • v 4 is a sensitive probe of elliptic flow fluctuations. The standard model of eccentricity fluctuations fails for central collisions. We need a better understanding of fluctuations.

Backup slides

More results from viscous hydro Glauber initial conditions and smaller viscosity also reproduces the measured v 2.