We investigate the effects of finite baryon density and temperature on the bulk properties of matter formed in relativistic heavy ion collisions within second-order dissipative hydrodynamics. The relativistic fluid evolution equations for heat flow and shear stress tensor are derived from kinetic theory by using Grads 14-moment approximation for the single-particle phase-space distribution function. The new equations provide a number of additional terms associated with heat-shear couplings as compared to the existing derivations based on entropy principle. The dissipative equations are encoded in non-boost-invariant hydrodynamic model simulation and studied for the evolution of high baryon density matter encountered at the beam energy scan program at RHIC. We find that thermal dissipation dominates shear pressure in defining the bulk observables at the low energy but its effect diminishes at ultra-relativistic energies.
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$.
We review studies of vortical motion and the resulting global polarization of $Lambda$ and $bar{Lambda}$ hyperons in heavy-ion collisions, in particular, within 3FD model. 3FD predictions for the global midrapidity polarization in the FAIR-NICA energy range are presented. The 3FD simulations indicate that energy dependence of the observed global polarization of hyperons in the midrapidity region is a consequence of the decrease of the vorticity in the central region with the collision energy rise because of pushing out the vorticity field into the fragmentation regions. At high collision energies this pushing-out results in a peculiar vortical structure consisting of two vortex rings: one ring in the target fragmentation region and another one in the projectile fragmentation region with matter rotation being opposite in these two rings.
The hot and dense matter generated in heavy-ion collisions contains intricate vortical structure in which the local fluid vorticity can be very large. Such vorticity can polarize the spin of the produced particles. We study the event-by-event generation of the so-called thermal vorticity in Au + Au collisions at energy region $sqrt{s}=7.7-200$ GeV and calculate its time evolution, spatial distribution, etc., in a multiphase transport (AMPT) model. We then compute the spin polarization of the $Lambda$ and $bar{Lambda}$ hyperons as a function of $sqrt{s}$, transverse momentum $p_T$, rapidity, and azimuthal angle. Furthermore, we study the harmonic flow of the spin, in a manner analogous to the harmonic flow of the particle number. The measurement of the spin harmonic flow may provide a way to probe the vortical structure in heavy-ion collisions. We also discuss the spin polarization of $Xi^0$ and $Omega^-$ hyperons which may provide further information about the spin polarization mechanism of hadrons.
We develop for charmed hadron production in relativistic heavy-ion collisions a comprehensive coalescence model that includes an extensive set of $s$ and $p$-wave hadronic states as well as the strict energy-momentum conservation, which ensures the boost invariance of the coalescence probability and the thermal limit of the produced hadron spectrum. By combining our hadronization scheme with an advanced Langevin-hydrodynamics model that incorporates both elastic and inelastic energy loss of heavy quarks inside the dynamical quark-gluon plasma, we obtain a successful description of the $p_mathrm{T}$-integrated and differential $Lambda_c/D^0$ and $D_s/D^0$ ratios measured at RHIC and the LHC. We find that including the effect of radial flow of the medium is essential for describing the enhanced $Lambda_c/D^0$ ratio observed in relativistic heavy-ion collisions. We also find that the puzzling larger $Lambda_c/D^0$ ratio observed in Au+Au collisions at RHIC than in Pb+Pb collisions at the LHC is due to the interplay between the effects of the QGP radial flow and the charm quark transverse momentum spectrum at hadronization. Our study further suggests that charmed hadrons have larger sizes in medium than in vacuum.
Based on a generalized side-jump formalism for massless chiral fermions, which naturally takes into account the spin-orbit coupling in the scattering of two chiral fermions and the chiral vortical effect in a rotating chiral fermion matter, we have developed a covariant and total angular momentum conserved chiral transport model to study both the global and local polarizations of this matter. For a system of massless quarks of random spin orientations and finite vorticity in a box, we have demonstrated that the model can exactly conserve the total angular momentum of the system and dynamically generate the quark spin polarization expected from a thermally equilibrated quark matter. Using this model to study the spin polarization in relativistic heavy-ion collision, we have found that the local quark spin polarizations depend strongly on the reference frame where they are evaluated as a result of the nontrivial axial charge distribution caused by the chiral vortical effect. We have further shown that because of the anomalous orbital or side-jump contribution to the quark spin polarization, the local quark polarizations calculated in the medium rest frame are qualitatively consistent with the local polarizations of Lambda hyperons measured in experiments.