# Flow in heavy-ion collisions - Theory Perspective

## Abstract

I review recent developments in the field of relativistic hydrodynamics and its application to the bulk dynamics in heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron Collider (LHC). In particular, I report on progress in going beyond second order relativistic viscous hydrodynamics for conformal fluids, including temperature dependent shear viscosity to entropy density ratios, as well as coupling hydrodynamic calculations to microscopic hadronic rescattering models. I describe event-by-event hydrodynamic simulations and their ability to compute higher harmonic flow coefficients. Combined comparisons of all harmonics to recent experimental data from both RHIC and LHC will potentially allow to determine the desired details of the initial state and the medium properties of the quark-gluon plasma produced in heavy-ion collisions.

## 1 Introduction

The observed large elliptic flow at RHIC and recently at the LHC is one of the most striking observations in heavy-ion collision experiments. This asymmetry of particle production in the transverse plane of the collision is interpreted as the hydrodynamic response to the initial geometry. The applicability of hydrodynamics demands a short mean free path with respect to the system size. Therefore it is concluded that the created quark-gluon plasma is strongly coupled and behaves like a nearly perfect fluid. Recently interest has cascaded to all higher harmonics, including odd ones such as the triangular flow coefficient , which are non-zero in single events. A lot of progress has been made in extracting medium properties from hydrodynamic calculations, with the largest uncertainty emerging from the limited knowledge of the initial conditions.

## 2 Theoretical framework

The current standard for the viscous hydrodynamic description of relativistic heavy-ion collisions has been established in [1], where all terms up to second order in gradients for a conformal fluid have been derived. Additional terms for non-conformal fluids with non-zero bulk viscosity have been derived in [2].

In the ideal case, the evolution of the system created in relativistic heavy-ion collisions is described by the following 5 conservation equations

(1) |

where is the energy-momentum tensor and is the net baryon current. These are usually re-expressed using the time-like flow four-vector as

(2) |

where is the energy density, is the pressure, is the baryon density and is the metric tensor. The equations are then closed by adding the equilibrium equation of state

(3) |

as a local constraint on the variables.

In the first-order, or Navier-Stokes formalism for viscous hydrodynamics, the stress-energy tensor is decomposed into where is given in Eq. (2) and the viscous part of the stress energy tensor is given by

(4) |

where is the local 3-metric and is the local spatial derivative. Note that is transverse with respect to the flow velocity since and . Hence, is also an eigenvector of the whole stress-energy tensor with the same eigenvalue . is the shear viscosity of the medium.

The Navier-Stokes form is conceptually simple but introduces unphysical super-luminal signals that lead to numerical instabilities.

The second-order Israel-Steward formalism [3, 4, 5] avoids this super-luminal propagation, as does the more recent approach in [6]. In the Israel-Stewart formalism for a conformal fluid, derived in [1], the stress-energy tensor is decomposed as The evolution equations are

(5) | |||||

(6) |

where we neglected vorticity and terms that turn out to be numerically irrelevant. For the role of vorticity in heavy-ion collisions when including fluctuations see [7].

Simulations of bulk dynamics in heavy-ion collisions using this formalism have been performed in 2+1 dimensions in [8, 9, 10, 11, 12, 13], and within the equivalent Öttinger-Grmela [14] formalism in [15]. Recently, 3+1 dimensional viscous calculations have also become available [16, 17].

Having established the standard theoretical framework, in the following I will discuss recent developments in the field of relativistic viscous hydrodynamics and its application to heavy-ion collisions.

## 3 Equation of state

Typically, the equation of state (3) used in hydrodynamic simulations of heavy-ion collisions is determined from lattice QCD calculations combined with a hadron gas model. A recent parametrization and a detailed comparison of different lattice equations of state has been performed in [18]. The different equations of state [19, 20, 21, 22, 13], which have different results for the trace anomaly and the speed of sound, lead to different evolution of the momentum anisotropy when used in an ideal hydrodynamic evolution. However, the difference in the final spectra and elliptic flow results turns out to be negligible. This leads to the conclusion that one cannot determine details about the equation of state from the comparison of hydrodynamic calculations to experimental data. On the other hand, using one of the latest lattice equations of state will do a good enough job when extracting medium properties such as transport coefficients from hadronic observables. Thermal photon production, however, is potentially more sensitive to the equation of state, because photons are produced throughout the evolution and hence probe the dynamics of the system more directly [23]. Another relevant aspect is the inclusion of chemical freeze-out in the equation of state (and the freeze-out procedure) to reproduce the correct final particle ratios [24, 25, 26].

## 4 Bulk viscosity and limits of second order viscous hydrodynamics

QCD is a non-conformal theory and contains a finite bulk viscosity. In [27] the effect of bulk viscosity on elliptic flow and in [28, 29] the combined effect of shear and bulk viscosity has been studied. Bulk viscosity is expected to peak (possibly along with its relaxation time) around the critical temperature , where the system may develop large correlation lengths [30, 31, 32, 33]. In [29], within the range of applicability of second order viscous hydrodynamics, corrections from bulk viscosity are found to be small compared to those from shear viscosity if the bulk relaxation time peaks around . This means that if bulk and shear viscous correction do not become so large as to render the hydrodynamic expansion invalid for the relevant part of a heavy-ion collision, the extraction of the shear viscosity should be possible to reasonable accuracy when neglecting the effect of bulk viscosity.

Apart from viscous corrections to the thermal equilibrium distribution functions becoming large, bulk and shear viscous corrections can lead to a negative longitudinal pressure, and hence the breakup of the system into droplets [34, 35, 36]. Second order viscous hydrodynamics can hint at when such cavitation happens but is not suited to describe the process as it happens outside the range of its applicability.

## 5 Progress beyond second order viscous hydrodynamics

The problem with the viscous correction to the stress energy tensor becoming larger than the equilibrium part, leading to negative pressure, arises especially in the early stage of the evolution, when the local momentum distribution is not yet equilibrated but highly anisotropic due to rapid longitudinal expansion. Recent progress has been made in describing this early time evolution and late time hydrodynamics within the same framework by performing the hydrodynamic expansion around an anisotropic distribution [37, 38, 39, 40]. This results in new equations of motion, including one for the degree of anisotropy of the distribution function. The procedure can reproduce both the limits of free streaming and ideal hydrodynamics and results in second order viscous hydrodynamics when expanding around a small anisotropy parameter, which so far has been shown in the one dimensional case [37].

Another logical step is to expand to third order in gradients as done in [41]. Numerical differences to second order Israel Stewart theory are completely negligible for and become significant for . It has also been pointed out that the hydrodynamic equations depend on the details of their derivation. While Israel and Stewart used the second moment of the Boltzmann equation to derive hydrodynamic equations for the dissipative currents, in [42] the definition of the latter was used directly. This leads to equations of motion of the same form but with different coefficients. Microscopic transport calculations [43, 44] show very good agreement with the new equations of motion up to , while the Israel-Stewart equations show large differences for . This demonstrates that the details of the derivation are relevant in particular when hydrodynamics is being matched to kinetic theory at late times.

## 6 Temperature dependent

What is extracted from experimental data by comparison with viscous hydrodynamic calculations using constant is at best an effective, or average, . In reality should depend on the local temperature of the medium, dropping from large values at high temperatures to a minimum at , and rising with decreasing temperature in the hadronic phase [45, 46, 47, 48].

It has been shown that when including such modeled temperature dependence in calculations at RHIC energies, the details of in the quark-gluon plasma phase have little influence on the final elliptic flow result, while hadronic modifies strongly [48]. At the highest LHC energies the conclusion is the opposite: weak dependence on the hadronic but strong dependence on in the QGP phase. Interestingly, at RHIC energies, a significant dependence has been found on the minimum value of around [49]. This might indicate that we are determining such a minimal value, rather than , when using a constant .

In addition, it was pointed out in [50] that there is a strong dependence on the initial value of when starting with a large . So, particularly when including a temperature dependent , it is essential to gain a better understanding of the pre-equilibrium stage in heavy-ion collisions to determine the initial conditions for viscous hydrodynamics.

## 7 Viscous corrections to particle distribution functions

When translating the dissipative to particles in the Cooper-Frye formalism [51], corrections to the distribution function have to be taken into account:

(7) |

for an component system with the degeneracy of species . The ansatz

(8) |

where is a unit vector in the direction, leaves some freedom and the usual procedure is to assume that all coefficients are equal, even though they should depend on the individual particle species’ interaction rate (see [52]), and use , which is derived within a relaxation time plus Boltzmann approximation:

(9) |

It has been shown in [53] that , where can take on values from 1 to 2, which is the case for example for a system with radiative and elastic energy loss that has . The exact form of the correction has a large effect on the differential elliptic flow for . It should be noted, however, that the analysis of experimental data in [54] indicates . A more general problem is that corrections can become large compared to . For hadrons this problem can be reduced when using a hadronic afterburner and switching at intermediate temperatures of , but for photons that are produced throughout the whole evolution this becomes a serious concern [23].

## 8 Hadronic afterburner

The use of hybrid models coupling early hydrodynamic evolution to a microscopic hadronic cascade in the later stage has a long history [55, 56, 57, 58, 59, 60, 61]. Recently 2+1 dimensional viscous hydrodynamics has been coupled to UrQMD [62, 63, 64, 65, 66]. When using a chemically frozen equation of state, transverse momentum spectra of produced particles show very little dependence on the temperature at which the switching between hydro and UrQMD is performed. However, does show a dependence unless an increasing is introduced in the hadronic stage of the hydro simulation to emulate the effect of UrQMD’s larger dissipation [62]. This rises to , for which the applicability of second order viscous hydrodynamics becomes questionable. Also, it turns out that the dynamics in UrQMD cannot be described by viscous hydrodynamics with due to different relaxation times [62], underlining the need for the more realistic microscopic rescattering.

## 9 Initial conditions and event-by-event hydrodynamics

Fluctuating initial conditions for hydrodynamic simulations of heavy-ion collisions have been argued to be very important for the exact determination of collective flow observables and to describe features of multi-particle correlation measurements in heavy-ion collisions [67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 16, 17, 84]. Real event-by-event hydrodynamic simulations have been performed and show modifications to spectra and flow from “single-shot” hydrodynamics with averaged initial conditions [81, 16, 17, 84, 85]. An important advantage of event-by-event hydrodynamic calculations is the possibility to study higher flow harmonics such as , which are entirely due to fluctuations in the initial conditions. Different depend differently on and the details of the initial condition, which is determined by the dynamics and fluctuations of partons in the incoming nuclear wave functions. This observation can be used to determine these long sought after details of the initial state and medium properties in heavy-ion collisions by performing a systematic analysis of all harmonics , up to e.g. as a function of and the initial state properties and compare to experimental data. First predictions of [16] agree extremely well with experimental data from RHIC [86]. Furthermore, it has now been shown that at low (and (ALICE), (ATLAS)), the main features of dihadron correlations in the angular difference between the hadron momenta can be described by flow, i.e. the sum of to only [87, 88]. The double-peak structure on the away-side is hence described mostly by (triangular) flow as predicted in [77, 79].

## 10 Summary of LHC predictions, conclusions and outlook

Calculations using ideal and viscous hydrodynamics at LHC energies hint at little or no increase of [89, 90, 60, 61, 17, 65, 91]. However, to be consistent a temperature dependent should be employed to extend RHIC calculations to LHC (see e.g. [65]). Also, essentially all models underestimate for , which might be explained by contributions from non-thermalized particles [92]. The issue is however not settled yet. Nevertheless, the overall early success of hydrodynamics in describing the main features of flow of charged hadrons at the LHC indicates that the QGP at LHC is also a nearly perfect fluid. Higher harmonics are also being computed for the LHC and comparison of experimental data to first predictions from event-by-event simulations [17] show that hydrodynamic results are in the right ball park. As mentioned earlier, a systematic analysis will be needed to potentially determine and the details of the initial conditions at LHC with the use of all .

Better understanding of viscous corrections to the distribution functions and dependence on the model for coupling to a hadronic afterburner, detailed studies of 3+1 dimensional hydrodynamics with viscosity, and above all a better understanding of the pre-equilibrium stage and its transition to hydrodynamics are the next important steps in developing a reliable viscous hydrodynamic description of the bulk dynamics in heavy-ion collisions.

## Acknowledgments

I thank Charles Gale, Sangyong Jeon, Derek Teaney, and Raju Venugopalan for helpful comments on the manuscript and discussions. This work was supported in part by the US Department of Energy under DOE Contract No.DE-AC02-98CH10886, and by a Lab Directed Research and Development Grant from Brookhaven Science Associates.

## References

### References

- Baier R, Romatschke P, Son D T, Starinets A O and Stephanov M A 2008 JHEP 04 100
- Betz B, Henkel D and Rischke D 2009 J.Phys.G G36 064029
- Israel W 1976 Ann. Phys. 100 310–331
- Stewart J 1977 Proc.Roy.Soc.Lond. A357 59
- Israel W and Stewart J M 1979 Ann. Phys. 118 341–372
- Muronga A 2002 Phys. Rev. Lett. 88 062302
- Floerchinger S and Wiedemann U A 2011 (Preprint arXiv:1106.5227)
- Romatschke P 2007 Eur. Phys. J. C52 203–209
- Romatschke P and Romatschke U 2007 Phys. Rev. Lett. 99 172301
- Luzum M and Romatschke P 2008 Phys. Rev. C78 034915
- Heinz U W, Song H and Chaudhuri A K 2006 Phys. Rev. C73 034904
- Song H and Heinz U W 2008 Phys. Lett. B658 279–283
- Song H and Heinz U W 2008 Phys. Rev. C78 024902
- Grmela M and Ottinger H C 1997 Phys. Rev. E56 6620–6632
- Dusling K and Teaney D 2008 Phys. Rev. C77 034905
- Schenke B, Jeon S and Gale C 2011 Phys. Rev. Lett. 106 042301
- Schenke B, Jeon S and Gale C 2011 Physics Letters B , in press (Preprint arXiv:1102.0575)
- Huovinen P and Petreczky P 2010 Nucl. Phys. A837 26–53
- Laine M and Schroder Y 2006 Phys. Rev. D73 085009
- Chojnacki M and Florkowski W 2007 Acta Phys.Polon. B38 3249–3262
- Chojnacki M, Florkowski W, Broniowski W and Kisiel A 2008 Phys. Rev. C78 014905
- Bazavov A, Bhattacharya T, Cheng M, Christ N, DeTar C et al. 2009 Phys. Rev. D80 014504
- Gale C 2011 these proceedings
- Hirano T and Tsuda K 2002 Phys. Rev. C66 054905
- Kolb P F and Rapp R 2003 Phys. Rev. C67 044903
- Shen C, Heinz U, Huovinen P and Song H 2010 Phys. Rev. C82 054904
- Denicol G, Kodama T, Koide T and Mota P 2009 Phys. Rev. C80 064901
- Monnai A and Hirano T 2009 Phys. Rev. C80 054906
- Song H and Heinz U W 2010 Phys. Rev. C81 024905
- Paech K and Pratt S 2006 Phys. Rev. C74 014901
- Kharzeev D and Tuchin K 2008 JHEP 0809 093
- Meyer H B 2009 Nucl.Phys. A830 641C–648C
- Buchel A 2009 Phys.Lett. B681 200–203 (Preprint arXiv:0908.0108)
- Torrieri G and Mishustin I 2008 Phys. Rev. C78 021901
- Rajagopal K and Tripuraneni N 2010 JHEP 1003 018
- Bhatt J R, Mishra H and Sreekanth V 2011 (Preprint arXiv:1103.4333)
- Martinez M and Strickland M 2010 Nucl.Phys. A848 183–197
- Martinez M and Strickland M 2011 Nucl.Phys. A856 68–87
- Florkowski W and Ryblewski R 2011 Phys. Rev. C83 034907
- Ryblewski R and Florkowski W 2011 (Preprint arXiv:1103.1260)
- El A, Xu Z and Greiner C 2010 Phys. Rev. C81 041901
- Denicol G, Koide T and Rischke D 2010 Phys. Rev. Lett. 105 162501
- Xu Z and Greiner C 2005 Phys. Rev. C71 064901
- Xu Z and Greiner C 2008 Phys. Rev. Lett. 100 172301
- Nakamura A and Sakai S 2005 Phys. Rev. Lett. 94 072305
- Noronha-Hostler J, Noronha J and Greiner C 2009 Phys. Rev. Lett. 103 172302
- Denicol G, Kodama T and Koide T 2010 J.Phys.G G37 094040
- Niemi H, Denicol G S, Huovinen P, Molnar E and Rischke D H 2011 (Preprint arXiv:1101.2442)
- Niemi H 2011 these proceedings
- Shen C, Heinz U W, Huovinen P and Song H 2011 (Preprint arXiv:1105.3226)
- Cooper F and Frye G 1974 Phys. Rev. D10 186
- Molnar D 2011 these proceedings
- Dusling K, Moore G D and Teaney D 2010 Phys. Rev. C81 034907
- Lacey R A, Taranenko A, Wei R, Ajitanand N, Alexander J et al. 2010 Phys. Rev. C82 034910
- Bass S and Dumitru A 2000 Phys. Rev. C61 064909
- Teaney D, Lauret J and Shuryak E V 2001 Phys. Rev. Lett. 86 4783–4786
- Hirano T, Heinz U W, Kharzeev D, Lacey R and Nara Y 2006 Phys. Lett. B636 299–304
- Nonaka C and Bass S A 2007 Phys. Rev. C75
- Nonaka C 2010 AIP Conf.Proc. 1235 165–171
- Hirano T, Huovinen P and Nara Y 2011 Phys. Rev. C83 021902
- Hirano T, Huovinen P and Nara Y 2010 (Preprint arXiv:1012.3955)
- Song H, Bass S A and Heinz U 2011 Phys. Rev. C83 024912
- Song H, Bass S A, Heinz U, Hirano T and Shen C 2011 Phys. Rev. Lett. 106 192301
- Song H, Bass S A, Heinz U, Hirano T and Shen C 2011 Phys. Rev. C83 054910
- Song H, Bass S A and Heinz U W 2011 (Preprint arXiv:1103.2380)
- Heinz U 2011 these proceedings
- Andrade R et al. 2006 Phys. Rev. Lett. 97 202302
- Adare A et al. (PHENIX) 2008 Phys. Rev. C78 014901
- Abelev B I et al. (STAR) 2009 Phys. Rev. Lett. 102 052302
- Alver B et al. (PHOBOS) 2010 Phys. Rev. C81 024904
- Alver B et al. (PHOBOS) 2010 Phys. Rev. Lett. 104 062301
- Abelev B I et al. (STAR) 2009 Phys. Rev. C80 064912
- Miller M and Snellings R 2003 (Preprint nucl-ex/0312008)
- Broniowski W et al. 2007 Phys. Rev. C76 054905
- Andrade R et al. 2008 Phys. Rev. Lett. 101 112301
- Hirano T and Nara Y 2009 Nucl. Phys. A830 191c–194c
- Takahashi J et al. 2009 Phys. Rev. Lett. 103 242301
- Andrade R et al. 2010 J. Phys. G37 094043
- Alver B and Roland G 2010 Phys. Rev. C81 054905
- Werner K, Karpenko I, Pierog T, Bleicher M and Mikhailov K 2010 (Preprint arXiv:1004.0805)
- Holopainen H, Niemi H and Eskola K J 2010 (Preprint arXiv:1007.0368)
- Alver B et al. 2010 (Preprint arXiv:1007.5469)
- Petersen H et al. 2010 (Preprint arXiv:1008.0625)
- Qiu Z and Heinz U W 2011 (Preprint arXiv:1104.0650)
- Jeon S 2011 these proceedings
- Adare A et al. (PHENIX Collaboration) 2011 (Preprint arXiv:1105.3928)
- Jia J 2011 these proceedings
- Aamodt K et al. (ALICE Collaboration) 2011 (Preprint arXiv:1105.3865)
- Luzum M and Romatschke P 2009 Phys. Rev. Lett. 103 262302
- Luzum M 2011 Phys. Rev. C83 044911
- Roy V and Chaudhuri A 2011 (Preprint arXiv:1103.2870)
- Bozek P 2011 Phys.Lett. B699 283–286