In an earlier work (arXiv:0808.0953) we established that causal Israel-Stewart viscous hydrodynamics is only accurate in RHIC applications at very low shear viscosities 4 pi eta_s / s < ~ 1.5-2. We show here that the region of applicability is significantly reduced if bulk viscosity plays a role in the dynamics.
The charged-particles final state spectrum is derived from an analytic perturbative solution for the relativistic viscous hydrodynamics. By taking into account the longitudinal acceleration effect in relativistic viscous hydrodynamics, the pseudorapidity spectrum describes well the nucleus-nucleus colliding systems at RHIC and LHC. Based on both the extracted longitudinal acceleration parameters $lambda^{*}$ and a phenomenological description of the $lambda^{*}$, the charged-particles pseudorapidity distributions for $sqrt{s_{NN}}$ = 5.44 TeV Xe+Xe collisions are computed from the final state expression in a limited space-time rapidity $eta_{s}$ region.
Hydrodynamics is a general theoretical framework for describing the long-time large-distance behaviors of various macroscopic physical systems, with its equations based on conservation laws such as energy-momentum conservation and charge conservation. Recently there has been significant interest in understanding the implications of angular momentum conservation for a corresponding hydrodynamic theory. In this work, we examine the key conceptual issues for such a theory in the relativistic regime where the orbital and spin components get entangled. We derive the equations for relativistic viscous hydrodynamics with angular momentum through Navier-Stokes type of gradient expansion analysis.
Correlations of different azimuthal flow harmonics $v_n$ and symmetry planes $Psi_n$ can add constraints to theoretical models, and probe aspects of the system that are independent of the traditional single-harmonic measurements. Using NeXSPheRIO, a hydrodynamical model which has accurately reproduced a large set of single-harmonic correlations at RHIC, we make predictions of these new observables for 200 A GeV Au+Au Collisions, providing an important baseline for comparison to correlations of flow harmonics, which contain non-trivial information about the initial state. We also point out how to properly compare theoretical calculations to measurements using wide centrality bins and non-trivial event weighting.
We utilize nonequilibrium covariant transport theory to determine the region of validity of causal Israel-Stewart dissipative hydrodynamics (IS) and Navier-Stokes theory (NS) for relativistic heavy ion physics applications. A massless ideal gas with 2->2 interactions is considered in a 0+1D Bjorken scenario, appropriate for the early longitudinal expansion stage of the collision. In the scale invariant case of a constant shear viscosity to entropy density ratio eta/s ~ const, we find that Israel-Stewart theory is 10% accurate in calculating dissipative effects if initially the expansion timescale exceeds half the transport mean free path tau0/lambda0 > ~2. The same accuracy with Navier-Stokes requires three times larger tau0/lambda0 > ~6. For dynamics driven by a constant cross section, on the other hand, about 50% larger tau0/lambda0 > ~3 (IS) and ~9 (NS) are needed. For typical applications at RHIC energies s_{NN}**(1/2) ~ 100-200 GeV, these limits imply that even the Israel-Stewart approach becomes marginal when eta/s > ~0.15. In addition, we find that the naive approximation to Israel-Stewart theory, which neglects products of gradients and dissipative quantities, has an even smaller range of applicability than Navier-Stokes. We also obtain analytic Israel-Stewart and Navier-Stokes solutions in 0+1D, and present further tests for numerical dissipative hydrodynamics codes in 1+1, 2+1, and 3+1D based on generalized conservation laws.
It is argued that the use of the initial Gaussian energy density profile for hydrodynamics leads to much better uniform description of the RHIC heavy-ion data than the use of the standard initial condition obtained from the Glauber model. With the modified Gaussian initial conditions we successfully reproduce the transverse-momentum spectra, v2, and the pionic HBT radii (including their azimuthal dependence). The emerging consistent picture of hadron production hints that a solution of the long standing RHIC HBT puzzle has been found.