We review progress in the hydrodynamic description of heavy-ion collisions, focusing on recent developments in modeling the fluctuating initial state and event-by-event viscous hydrodynamic simulations. We discuss how hydrodynamics can be used to extract information on fundamental properties of quantum-chromo-dynamics from experimental data, and review successes and challenges of the hydrodynamic framework. arXiv:1301.5893v1 [nucl-th]

We present results for the elliptic and triangular flow coefficients v(2) and v(3) in Au+Au collisions at √s=200  AGeV using event-by-event D=3+1 viscous hydrodynamic simulations. We study the effect of initial state fluctuations and finite viscosities on the flow coefficients v(2) and v(3) as functions of transverse momentum and pseudorapidity. Fluctuations are essential to reproduce the measured centrality dependence of elliptic flow. We argue that simultaneous measurements of v(2) and v(3) can determine η/s more precisely.

We compute initial conditions in heavy ion collisions within the color glass condensate framework by combining the impact parameter dependent saturation model with the classical Yang-Mills description of initial Glasma fields. In addition to fluctuations of nucleon positions, this impact parameter dependent Glasma description includes quantum fluctuations of color charges on the length scale determined by the inverse nuclear saturation scale Q(s). The model naturally produces initial energy fluctuations that are described by a negative binomial distribution. The ratio of triangularity to eccentricity ε(3)/ε(2) is close to that in a model tuned to reproduce experimental flow data. We compare transverse momentum spectra and v(2,3,4)(p(T)) of pions from different models of initial conditions using relativistic viscous hydrodynamic evolution.

Anisotropic flow coefficients v(1)-v(5) in heavy ion collisions are computed by combining a classical Yang-Mills description of the early time Glasma flow with the subsequent relativistic viscous hydrodynamic evolution of matter through the quark-gluon plasma and hadron gas phases. The Glasma dynamics, as realized in the impact parameter dependent Glasma (IP-Glasma) model, takes into account event-by-event geometric fluctuations in nucleon positions and intrinsic subnucleon scale color charge fluctuations; the preequilibrium flow of matter is then matched to the music algorithm describing viscous hydrodynamic flow and particle production at freeze-out. The IP-Glasma+MUSIC model describes well both transverse momentum dependent and integrated v(n) data measured at the Large Hadron Collider and the Relativistic Heavy Ion Collider. The model also reproduces the event-by-event distributions of v(2), v(3) and v(4) measured by the ATLAS Collaboration. The implications of our results for better understanding of the dynamics of the Glasma and for the extraction of transport properties of the quark-gluon plasma are outlined.

We present music, an implementation of the Kurganov-Tadmor algorithm for relativistic 3+1 dimensional fluid dynamics in heavy-ion collision scenarios. This Riemann-solver-free, second-order, high-resolution scheme is characterized by a very small numerical viscosity and its ability to treat shocks and discontinuities very well. We also incorporate a sophisticated algorithm for the determination of the freeze-out surface using a three dimensional triangulation of the hyper-surface. Implementing a recent lattice based equation of state, we compute pT -spectra and pseudorapidity distributions for Au+Au collisions at √ s = 200 GeV and present results for the anisotropic flow coefficients v2 and v4 as a function of both pT and pseudorapidity η. We were able to determine v4 with high numerical precision, finding that it does not strongly depend on the choice of initial condition or equation of state.