No Arabic abstract
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.
The nonextensive one-dimensional version of a hydrodynamical model for multiparticle production processes is proposed and discussed. It is based on nonextensive statistics assumed in the form proposed by Tsallis and characterized by a nonextensivity parameter $q$. In this formulation the parameter $q$ characterizes some specific form of local equilibrium which is characteristic for the nonextensive thermodynamics and which replaces the usual local thermal equilibrium assumption of the usual hydrodynamical models. We argue that there is correspondence between the perfect nonextensive hydrodynamics and the usual dissipative hydrodynamics. It leads to simple expression for dissipative entropy current and allows for predictions for the ratio of bulk and shear viscosities to entropy density, $zeta/s$ and $eta/s$, to be made.
We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in textit{A}+textit{A} collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time ($tau_{text{rel}}$) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small $tau_{text{rel}}$ it also allows one to catch the viscous effects in hadronic component - hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion $m_{T}$ spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher $p_{T}$ particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.
We review integrated dynamical approaches to describe heavy ion reaction as a whole at ultrarelativistic energies. Since final observables result from all the history of the reaction, it is important to describe all the stages of the reaction to obtain the properties of the quark gluon plasma from experimental data. As an example of these approaches, we develop an integrated dynamical model, which is composed of a fully (3+1) dimensional ideal hydrodynamic model with the state-of-the-art equation of state based on lattice QCD, and subsequent hadronic cascade in the late stage. Initial conditions are obtained employing Monte Car
We study charm production in ultra-relativistic heavy-ion collisions by using the Parton-Hadron-String Dynamics (PHSD) transport approach. The initial charm quarks are produced by the PYTHIA event generator tuned to fit the transverse momentum spectrum and rapidity distribution of charm quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) calculations. The produced charm quarks scatter in the quark-gluon plasma (QGP) with the off-shell partons whose masses and widths are given by the Dynamical Quasi-Particle Model (DQPM), which reproduces the lattice QCD equation-of-state in thermal equilibrium. The relevant cross sections are calculated in a consistent way by employing the effective propagators and couplings from the DQPM. Close to the critical energy density of the phase transition, the charm quarks are hadronized into $D$ mesons through coalescence and/or fragmentation. The hadronized $D$ mesons then interact with the various hadrons in the hadronic phase with cross sections calculated in an effective lagrangian approach with heavy-quark spin symmetry. The nuclear modification factor $R_{AA}$ and the elliptic flow $v_2$ of $D^0$ mesons from PHSD are compared with the experimental data from the STAR Collaboration for Au+Au collisions at $sqrt{s_{NN}}$ =200 GeV and to the ALICE data for Pb+Pb collisions at $sqrt{s_{NN}}$ =2.76 TeV. We find that in the PHSD the energy loss of $D$ mesons at high $p_T$ can be dominantly attributed to partonic scattering while the actual shape of $R_{AA}$ versus $p_T$ reflects the heavy-quark hadronization scenario, i.e. coalescence versus fragmentation. Also the hadronic rescattering is important for the $R_{AA}$ at low $p_T$ and enhances the $D$-meson elliptic flow $v_2$.
A key ingredient of hydrodynamical modeling of relativistic heavy ion collisions is thermal initial conditions, an input that is the consequence of a pre-thermal dynamics which is not completely understood yet. In the paper we employ a recently developed energy-momentum transport model of the pre-thermal stage to study influence of the alternative initial states in nucleus-nucleus collisions on flow and energy density distributions of the matter at the starting time of hydrodynamics. In particular, the dependence of the results on isotropic and anisotropic initial states is analyzed. It is found that at the thermalization time the transverse flow is larger and the maximal energy density is higher for the longitudinally squeezed initial momentum distributions. The results are also sensitive to the relaxation time parameter, equation of state at the thermalization time, and transverse profile of initial energy density distribution: Gaussian approximation, Glauber Monte Carlo profiles, etc. Also, test results ensure that the numerical code based on the energy-momentum transport model is capable of providing both averaged and fluctuating initial conditions for the hydrodynamic simulations of relativistic nuclear collisions.