Search

We obtain a formal integral solution to the 3+1 D Boltzmann Equation in relaxation time approximation. The gradient series obtained from this integral solution contains exponentially decaying non-hydrodynamic terms. It is shown that this gradient expansion can have a finite radius of convergence under certain assumptions of analyticity. We then argue that, in the relaxation time model, proximity to local thermal equilibrium is not necessary for the system to be described by hydrodynamic equations., Comment: 12 pages, 0 figures

We investigate the impact of baryon stopping on the temporal evolution of electromagnetic fields in vacuum at low-energy Au+Au collisions with $\sqrt{s_{NN}} = 4$-$20$ GeV. Baryon stopping is incorporated into the Monte-Carlo Glauber model by employing a parameterized velocity profile of participant nucleons with non-zero deceleration. The presence of these decelerating participants leads to noticeable changes in the centrality and $\sqrt{s_{NN}}$ dependence of electromagnetic fields compared to scenarios with vanishing deceleration. The influence of baryon stopping differs for electric and magnetic fields, also exhibiting variations across their components. We observe slight alteration in the approximate linear dependency of field strengths with $\sqrt{s_{NN}}$ in the presence of deceleration. Additionally, the longitudinal component of the electric field at late times becomes significant in the presence of baryon stopping., Comment: 10 pages and 10 figures

We study the effect of external proper time-dependent longitudinal forces on the evolution of the distribution function using the Boltzmann Equation with a relaxation time collision kernel under Bjorken flow. We derive an exact solution and study the early time attractor behaviour of different components of energy-momentum tensor under Bjorken symmetry. We show that the arbitrary initial data approaches the attractor solution but is influenced by the external force with slightly broken Bjorken symmetry., Comment: 13 pages,8 figures

Using the blast-wave model, we explore the effect of electric fields on spectra and flow harmonics (especially the elliptic flow) for charged pions and protons. We incorporate the first-order correction to the single-particle distribution function due to the electric fields and the dissipative effect while calculating the invariant yields of hadron in the Cooper-Frey prescription at the freezeout hypersurface. We find a noticeable correction to the directed and elliptic flow of pions and protons for unidirectional and azimuthal asymmetric electric fields in the transverse plane of magnitude $\sim m_{\pi}^{2}$. Further, we observe mass dependency of the directed flow generated due to the electric fields. The splitting of particle and antiparticle's elliptic flow is also discussed., Comment: 10 pages, 8 figures

The recently proposed connection between the Lorentz invariance of stability and the speed of signal propagation has been tested for a first-order relativistic dissipative hydrodynamic theory. The fact that the stability situation in different reference frames agrees with each other only as long as the signal propagation respects causality, has been explicitly established for the theory, which is microscopically derived from the covariant kinetic equation in general hydrodynamic frames with arbitrary momentum-dependent interactions., Comment: Accepted for the publication in Physics Letters B

The discovery and characterization of hot and dense QCD matter, known as Quark Gluon Plasma (QGP), remains the most international collaborative effort and synergy between theorists and experimentalists in modern nuclear physics to date. The experimentalists around the world not only collect an unprecedented amount of data in heavy-ion collisions, at Relativistic Heavy Ion Collider (RHIC), at Brookhaven National Laboratory (BNL) in New York, USA, and the Large Hadron Collider (LHC), at CERN in Geneva, Switzerland but also analyze these data to unravel the mystery of this new phase of matter that filled a few microseconds old universe, just after the Big Bang. In the meantime, advancements in theoretical works and computing capability extend our wisdom about the hot-dense QCD matter and its dynamics through mathematical equations. The exchange of ideas between experimentalists and theoreticians is crucial for the progress of our knowledge. The motivation of this first conference named "HOT QCD Matter 2022" is to bring the community together to have a discourse on this topic. In this article, there are 36 sections discussing various topics in the field of relativistic heavy-ion collisions and related phenomena that cover a snapshot of the current experimental observations and theoretical progress. This article begins with the theoretical overview of relativistic spin-hydrodynamics in the presence of the external magnetic field, followed by the Lattice QCD results on heavy quarks in QGP, and finally, it ends with an overview of experiment results., Comment: Compilation of the contributions (148 pages) as presented in the `Hot QCD Matter 2022 conference', held from May 12 to 14, 2022, jointly organized by IIT Goa & Goa University, Goa, India

We derive the first and second-order expressions for the shear, the bulk viscosity, and the thermal conductivity of a relativistic hot boson gas in a magnetic field using the relativistic kinetic theory within the Chapman-Enskog method. The order-by-order off-equilibrium distribution function is obtained in terms of the associate Laguerre polynomial with magnetic field-dependent coefficients using the relativistic Boltzmann-Uehling-Uhlenbeck transport equation. The order-by-order anisotropic transport coefficients are evaluated in powers of the dimensionless ratio of kinetic energy to the fluid temperature for finite magnetic fields. In a magnetic field, the shear viscosity (in all order) splits into five different coefficients. Four of them show a magnetic field dependence as seen in a previous study \cite{Ashutosh1} using the relaxation time approximation for the collision kernel. On the other hand, bulk viscosity, which splits into three components (in all order), is independent of the magnetic field. The thermal conductivity shows a similar splitting but is field-dependent. The difference in the first and second-order results are prominent for the thermal conductivities than the shear viscosity; moreover, the difference in the two results is most evident at low temperatures. The first and second-order results seem to converge rapidly for high temperatures., Comment: 35 pages, 2 figures

Here we study the wave propagation and stability of general relativistic non-resistive dissipative second-order magnetohydrodynamic equations in curved space-time. We solve the Boltzmann equation for a system of particles and antiparticles using the relaxation time approximation and the Chapman-Enskog-like gradient expansion for the off-equilibrium distribution function, truncating beyond second-order in curved space-time in electromagnetic fields. Unlike holographic calculation~\cite{Baier:2007ix}, we show that the viscous evolution equations do not explicitly depend on the curvature of space-time. Also, we have tested the causality and stability of the second-order theory in curved space-time in the presence of linearised metric perturbation and derived dispersion relations for various modes. Interestingly, we found the coupling of gravitational modes with the usual magneto-sonic modes in the small wave-number limit. Also, we show additional non-hydrodynamical modes arise due to gravity for a bulk-viscous fluid., Comment: 17 pages, 2 figures

We derive a first-order, stable and causal, relativistic hydrodynamic theory from the microscopic kinetic equation using the gradient expansion technique in a general frame. The general frame is introduced from the arbitrary matching conditions for hydrodynamic fields. The interaction is introduced in the relativistic Boltzmann equation through the momentum-dependent relaxation time approximation (MDRTA) with the proposed collision operator that preserves the conservation laws. We demonstrate here for the first time that not only the general frame choice, but also the momentum dependence of microscopic interaction rate, captured through MDRTA, is imperative for producing the essential field corrections that give rise to a causal and stable first-order relativistic theory., Comment: 6 pages, 3 figures, Accepted for publication in Physical Review D Letter

Intense transient electric ({\bf E}) and magnetic ({\bf B}) fields are produced in the high energy heavy-ion collisions. The electromagnetic fields produced in such high-energy heavy-ion collisions are proposed to give rise to a multitude of exciting phenomenon including the Chiral Magnetic Effect. We use a Monte Carlo (MC) Glauber model to calculate the electric and magnetic fields, more specifically their scalar product $\bf{E}\cdot\bf{B}$, as a function of space-time on an event-by-event basis for the Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV for different centrality classes. We also calculate the same for the isobars Ruthenium and Zirconium at $\sqrt{s_{NN}}=200$ GeV. In the QED sector $\bf{E}\cdot\bf{B}$ acts as a source of Chiral Separation Effect, Chiral Magnetic Wave, etc., which are associated phenomena to the Chiral Magnetic Effect. We also study the relationships between the electromagnetic symmetry plane angle defined by $\bf{E}\cdot\bf{B}$ ($\psi_{E.B}$) and the participant plane angle $\psi_{P}$ defined from the participating nucleons for the second-fifth order harmonics.

Here we derive the relativistic resistive dissipative second-order magnetohydrodynamic evolution equations using the Boltzmann equation, thus extending our work from the previous paper \href{https://link.springer.com/article/10.1007/JHEP03(2021)216}{JHEP 03 (2021) 216} where we considered the non-resistive limit. We solve the Boltzmann equation for a system of particles and antiparticles using the relaxation time approximation and the Chapman-Enskog like gradient expansion for the off-equilibrium distribution function, truncating beyond second-order. In the first order, the bulk and shear stress are independent of the electromagnetic field, however, the diffusion current, shows a dependence on the electric field. In the first order, the transport coefficients~(shear and bulk stress) are shown to be independent of the electromagnetic field. The diffusion current, however, shows a dependence on the electric field. In the second-order, the new transport coefficients that couple electromagnetic field with the dissipative quantities appear, which are different from those obtained in the 14-moment approximation~\cite{Denicol:2019iyh} in the presence of the electromagnetic field. Also we found out the various components of conductivity in this case., Comment: 11 pages, 0 figures

We derive the relativistic non-resistive, viscous second-order magnetohydrodynamic equations for the dissipative quantities using the relaxation time approximation. The Boltzmann equation is solved for a system of particles and antiparticles using Chapman-Enskog like gradient expansion of the single-particle distribution function truncated at second order. In the first order, the transport coefficients are independent of the magnetic field. In the second-order, new transport coefficients that couple magnetic field and the dissipative quantities appear which are different from those obtained in the 14-moment approximation \cite{Denicol:2018rbw} in the presence of a magnetic field. However, in the limit of the weak magnetic field, the form of these equations are identical to the 14-moment approximation albeit with a different values of these coefficients. We also derive the anisotropic transport coefficients in the Navier-Stokes limit., Comment: 20 pages, 1 figure

In this article, there are 18 sections discussing various current topics in the field of relativistic heavy-ion collisions and related phenomena, which will serve as a snapshot of the current state of the art. Section 1 reviews experimental results of some recent light-flavored particle production data from ALICE collaboration. Other sections are mostly theoretical in nature. Very strong but transient magnetic field created in relativistic heavy-ion collisions could have important observational consequences. This has generated a lot of theoretical activity in the last decade. Sections 2, 7, 9, 10 and 11 deal with the effects of the magnetic field on the properties of the QCD matter. There are several unanswered questions about the QCD phase diagram. Sections 3, 11 and 18 discuss various aspects of the QCD phase diagram and phase transitions. Recent years have witnessed interesting developments in foundational aspects of hydrodynamics and their application to heavy-ion collisions. Sections 12, 15, 16 and 17 of this article probe some aspects of this exciting field. Transport coefficients together with their temperature- and density-dependence, are essential inputs in hydrodynamical calculations. Sections 5, 8 and 14 deal with calculation/estimation of various transport coefficients (shear and bulk viscosity, thermal conductivity, relaxation times, etc.) of quark matter and hadronic matter. Sections 4, 6 and 13 deals with interesting new developments in the field. Section 4 discusses color dipole gluon distribution function at small transverse momentum in the form of a series of Bells polynomials. Section 6 discusses the properties of Higgs boson in the quark gluon plasma using Higgs-quark interaction. Section 13 discusses modification of coalescence model to incorporate viscous corrections and application of this model., Comment: 109 pages, 49 captioned figures, compilation of the contributions as presented in the `Workshop on Dynamics of QCD Matter', 15th to 17th August 2019, NISER Bhubaneswar, India, published version

We investigate the causality and the stability of the relativistic viscous magneto-hydrodynamics in the framework of the Israel-Stewart (IS) second-order theory, and also within a modified IS theory which incorporates the effect of magnetic fields in the relaxation equations of the viscous stress. We compute the dispersion relation by perturbing the fluid variables around their equilibrium values. In the ideal magnetohydrodynamics limit, the linear dispersion relation yields the well-known propagating modes: the Alfv\'en and the magneto-sonic modes.In the presence of bulk viscous pressure, the causality bound is found to be independent of the magnitude of the magnetic field. The same bound also remains true, when we take the full non-linear form of the equation using the method of characteristics. In the presence of shear viscous pressure, the causality bound is independent of the magnitude of the magnetic field for the two magneto-sonic modes. The causality bound for the shear-Alfv\'en modes, however, depends both on the magnitude and the direction of the propagation. For modified IS theory in the presence of shear viscosity, new non-hydrodynamic modes emerge but the asymptotic causality condition is the same as that of IS. In summary, although the magnetic field does influence the wave propagation in the fluid, the study of the stability and asymptotic causality conditions in the fluid rest frame shows that the fluid remains stable and causal given that they obey certain asymptotic causality condition., Comment: 38 pages, 7 figures, published in JHEP

An intense transient magnetic field is produced in high energy heavy-ion collisions mostly due to the spectator protons inside the two colliding nucleus. The magnetic field introduces anisotropy in the medium and hence the isotropic scalar transport coefficients become anisotropic and split into multiple components. Here we calculate the anisotropic transport coefficients shear, bulk viscosity, electrical conductivity, and the thermal diffusion coefficients for a multicomponent Hadron- Resonance-Gas (HRG) model for a non-zero magnetic field by using the Boltzmann transport equation in a relaxation time approximation (RTA). The anisotropic transport coefficient component along the magnetic field remains unaffected by the magnetic field, while perpendicular dissipation is governed by the interplay of the collisional relaxation time and the magnetic time scale, which is inverse of the cyclotron frequency. We calculate the anisotropic transport coefficients as a function of temperature and magnetic field using the HRG model. The neutral hadrons are unaffected by the Lorentz force and do not contribute to the anisotropic transports, we estimate within the HRG model the relative contribution of isotropic and anisotropic transports as a function of magnetic field and temperature. We also give an estimation of these anisotropic transport coefficients for the hadronic gas at finite baryon chemical potential., Comment: 15 pages, 7 figures

The Boltzmann equation is solved in the relaxation time approximation using a hierarchy of angular moments of the distribution function. Our solution is obtained for an azimuthally symmetric radially expanding boost-invariant conformal system that is undergoing Gubser flow. The solution of moments that we get after truncating the infinite set of equations at various orders is compared to the exact kinetic solution. The dynamics of transition is described by the presence of fixed points which describes the evolution of the system from an early time collisionless free streaming to the hydrodynamic regime at intermediate times and back to free streaming at late times. The attractor solution is found for various orders of moments as an interpolation between these fixed points. The relation of moments to various approximations of relativistic viscous hydrodynamics is investigated., Comment: 7 pages, 4 figures v2: major rewrite of hydrodynamic attractors section, new figures added, accepted for publication in Physics Letters B

Finding the existence and the location of the QCD critical point is one of the main goals of the RHIC beam energy scan program. To make theoretical predictions and corroborate with the experimental data requires modeling the space-time evolution of the matter created in heavy-ion collisions by dynamical models such as the relativistic hydrodynamics with an appropriate Equation of State (EoS). In the present exploratory study, we use a viscous 2+1 dimensional event-by-event (e-by-e) hydrodynamic code at finite baryon densities with two different EoSs (i) Lattice QCD + HRG with a crossover transition and (ii) EoS with a first-order phase transition to studying the normalized symmetric cumulants of charged pions $v_n$ $(n=2-4)$. We show that the normalized symmetric cumulants can differentiate the two EoSs while all other conditions remain the same. The conclusion does not change for various initial conditions and shear viscosity. This indicates that these observables can be used to gain information about the QCD EoS from experimental data and can be used as an EoS meter., Comment: Substantial changes from v2. Results with finite viscosity has been added

We calculate the $\delta f$ correction to the one particle distribution function in presence of magnetic field and non-zero shear viscosity within the relaxation time approximation. The $\delta f$ correction is found to be electric charge dependent. Subsequently, we also calculate one longitudinal and four transverse shear viscous coefficients as a function of dimensionless Hall parameter $\chi_{H}$ in presence of the magnetic field. We find that a proper linear combination of the shear viscous coefficients calculated in this work scales with the result obtained from Grad's moment method in \cite{Denicol:2018rbw}. Calculation of invariant yield of $\pi^{-}$ in a simple Bjorken expansion with cylindrical symmetry shows no noticeable change in spectra due to the $\delta f$ correction for realistic values of the magnetic field and relaxation time. However, when transverse expansion is taken into account using a blast wave type flow field we found noticeable change in spectra and elliptic flow coefficients due to the $\delta f$ correction. The $\delta f$ is also found to be very sensitive on the magnitude of magnetic field. Hence we think it is important to take into account the $\delta f$ correction in more realistic numerical magnetohydrodynamics simulations., Comment: 14 pages, 6 figures, revised version, new section added, new figures added, published in EPJA

We investigate the effect of large magnetic fields on the $2+1$ dimensional reduced-magnetohydrodynamical expansion of hot and dense nuclear matter produced in $\sqrt{s_{\rm NN}}$ = 200 GeV Au+Au collisions. For the sake of simplicity, we consider the case where the magnetic field points in the direction perpendicular to the reaction plane. We also consider this field to be external, with energy density parametrized as a two-dimensional Gaussian. The width of the Gaussian along the directions orthogonal to the beam axis varies with the centrality of the collision. The dependence of the magnetic field on proper time ($\tau$) for the case of zero electrical conductivity of the QGP is parametrized following [Deng 2012], and for finite electrical conductivity following [Tuchin 2013]. We solve the equations of motion of ideal hydrodynamics for such an external magnetic field. For collisions with non-zero impact parameter we observe considerable changes in the evolution of the momentum eccentricities of the fireball when comparing the case when the magnetic field decays in a conducting QGP medium and when no magnetic field is present. The elliptic-flow coefficient $v_2$ of $\pi^{-}$ is shown to increase in the presence of an external magnetic field and the increment in $v_2$ is found to depend on the evolution and the initial magnitude of the magnetic field., Comment: 12 pages, 6 figures, new discussion added, published in Phys.Rev.C

The interplay of magnetic field and thermal vorticity in a relativistic ideal fluid might generate fluid vorticity during the fluid evolution provided the flow fields and the entropy density of the fluid is inhomogeneous \cite{Mahajan:2010}. Exploiting this fact and assuming large magnetic Reynolds number we study the evolution of generalised magnetic field ($\hat{B}$) which is defined as a combination of the usual magnetic field ($\vec{B}$) and relativistic thermal vorticity ($\omega^{\mu\nu}$), in a 2(space)+1(time) dimensional isentropic evolution of Quark Gluon Plasma (QGP) with longitudinal boost invariance. The temporal evolution of $\hat{B}$ is found to be different than $\vec{B}$ , and the $\hat{B}$ evolution also depends on the position of the fluid along the beam direction (taken along the z axis) with respect to the mid-plane $z=0$. Further it is observed that the transverse components ($\hat{B_{x}}$, $\hat{B_{y}}$) evolve differently around the mid-plane., Comment: 8 pages, 2 figures

Relativistic hydrodynamics has been quite successful in explaining the collective behaviour of the QCD matter produced in high energy heavy-ion collisions at RHIC and LHC. We briefly review the latest developments in the hydrodynamical modeling of relativistic heavy-ion collisions. Essential ingredients of the model such as the hydrodynamic evolution equations, dissipation, initial conditions, equation of state, and freeze-out process are reviewed. We discuss observable quantities such as particle spectra and anisotropic flow and effect of viscosity on these observables. Recent developments such as event-by-event fluctuations, flow in small systems (proton-proton and proton-nucleus collisions), flow in ultra central collisions, longitudinal fluctuations and correlations and flow in intense magnetic field are also discussed., Comment: 36 pages, 16 figures, invited review, published version

We study the one-dimensional, longitudinally boost-invariant motion of an ideal fluid with infinite conductivity in the presence of a transverse magnetic field, i.e., in the ideal transverse magnetohydrodynamical limit. In an extension of our previous work Roy et al., [Phys. Lett. B 750, 45 (2015)], we consider the fluid to have a non-zero magnetization. First, we assume a constant magnetic susceptibility $\chi_{m}$ and consider an ultrarelativistic ideal gas equation of state. For a paramagnetic fluid (i.e., with $\chi_{m}>0$), the decay of the energy density slows down since the fluid gains energy from the magnetic field. For a diamagnetic fluid (i.e., with $\chi_{m}<0$), the energy density decays faster because it feeds energy into the magnetic field. Furthermore, when the magnetic field is taken to be external and to decay in proper time $\tau$ with a power law $\sim\tau^{-a}$, two distinct solutions can be found depending on the values of $a$ and $\chi_m$. Finally, we also solve the ideal magnetohydrodynamical equations for one-dimensional Bjorken flow with a temperature-dependent magnetic susceptibility and a realistic equation of state given by lattice-QCD data. We find that the temperature and energy density decay more slowly because of the non-vanishing magnetization. For values of the magnetic field typical for heavy-ion collisions, this effect is, however, rather small. Only for magnetic fields which are about an order of magnitude larger than expected for heavy-ion collisions, the system is substantially reheated and the lifetime of the quark phase might be extended., Comment: 12 pages, 9 figures; minor revisions, two references added, journal version accepted by PRD

The initial energy density distribution and fluctuation in the transverse direction lead to anisotropic flows of final hadrons through collective expansion in high-energy heavy-ion collisions. Fluctuations along the longitudinal direction, on the other hand, can result in decorrelation of anisotropic flows in different regions of pseudo rapidity ($\eta$). Decorrelation of the $2$nd and $3$rd order anisotropic flows with different $\eta$ gaps for final charged hadrons in high-energy heavy-ion collisions is studied in an event-by-event (3+1)D ideal hydrodynamic model with fully fluctuating initial conditions from A Multi-Phase Transport (AMPT) model. The decorrelation of anisotropic flows of final hadrons with large $\eta$ gaps are found to originate from the spatial decorrelation along the longitudinal direction in the AMPT initial conditions through hydrodynamic evolution. The decorrelation is found to consist of both a linear twist and random fluctuation of the event-plane angles. The agreement between our results and recent CMS data in most centralities suggests that the string-like mechanism of initial parton production in AMPT model captures the initial longitudinal fluctuation that is responsible for the measured decorrelation of anisotropic flows in Pb+Pb collisions at LHC. Our predictions for Au+Au collisions at the highest RHIC energy show stronger longitudinal decorrelation, indicating larger longitudinal fluctuations at lower beam energies. Our study also calls into question some of the current experimental methods for measuring anisotropic flows and extraction of transport coefficients through comparisons to hydrodynamic simulations that do not include longitudinal fluctuations.

We estimate the event-by-event (e-by-e) distribution of the ratio ($\sigma$) of the magnetic field energy to the fluid energy density in the transverse plane of Au-Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. A Monte-Carlo (MC) Glauber model is used to calculate the $\sigma$ in the transverse plane for impact parameter b=0, 12 fm at time $\tau_i\sim$0.5 fm. The fluid energy density is obtained by using Gaussian smoothing with two different smoothing parameter $\sigma_g$=0.25 , 0.5 fm. For $b=0~\rm fm$ collisions $\sigma$ is found to be $\ll$ 1 in the central region of the fireball and $\sigma\gtrsim$ 1 at the periphery. For b=12 fm collisions $\sigma\gtrsim$ 1. The e-by-e correlation between $\sigma$ and the fluid energy density ($\varepsilon$) is studied. We did not find strong correlation between $\sigma$ and $\varepsilon$ at the centre of the fireball, whereas they are mostly anti-correlated at the periphery of the fireball., Comment: 9 pages, 12 figures

We have calculated the temperature dependence of shear $\eta$ and bulk $\zeta$ viscosities of quark matter due to quark-meson fluctuations. The quark thermal width originating from quantum fluctuations of quark-$\pi$ and quark-$\sigma$ loops at finite temperature is calculated with the formalism of real-time thermal field theory. Temperature-dependent constituent-quark and meson masses, and quark-meson couplings are obtained in the Nambu--Jona-Lasinio model. We found a non-trivial influence of the temperature-dependent masses and couplings on the Landau-cut structure of the quark self-energy. Our results for the ratios $\eta/s$ and $\zeta/s$, where $s$ is the entropy density (also determined in the Nambu--Jona-Lasinio model in the quasi-particle approximation), are in fair agreement with results of the literature obtained from different models and techniques. In particular, our result for $\eta/s$ has a minimum very close to the conjectured AdS/CFT lower bound, $\eta/s = 1/4\pi$., Comment: 10 pages, 11 figures

In the initial stage of relativistic heavy-ion collisions, strong magnetic fields appear due to the large velocity of the colliding charges. The evolution of these fields appears as a novel and intriguing feature in the fluid-dynamical description of heavy-ion collisions. In this work, we study analytically the one-dimensional, longitudinally boost-invariant motion of an ideal fluid in the presence of a transverse magnetic field. Interestingly, we find that, in the limit of ideal magnetohydrodynamics, i.e., for infinite conductivity, and irrespective of the strength of the initial magnetization, the decay of the fluid energy density $e$ with proper time $\tau$ is the same as for the time-honored "Bjorken flow" without magnetic field. Furthermore, when the magnetic field is assumed to decay $\sim \tau^{-a}$, where $a$ is an arbitrary number, two classes of analytic solutions can be found depending on whether $a$ is larger or smaller than one. In summary, the analytic solutions presented here highlight that the Bjorken flow is far more general than formerly thought. These solutions can serve both to gain insight on the dynamics of heavy-ion collisions in the presence of strong magnetic fields and as testbeds for numerical codes., Comment: 9 pages, 4 figures, corrected typos and symbols, accepted for publication in Phys.Lett.B

Fluctuations in the initial transverse energy-density distribution lead to anisotropic flows as observed in central high-energy heavy-ion collisions. Studies of longitudinal fluctuations of the anisotropic flows can shed further light on the initial conditions and dynamical evolution of the hot quark-gluon matter in these collisions. Correlations between anisotropic flows with varying pseudorapidity gaps in Pb+Pb collisions at the CERN Large Hadron Collider are investigated using both an event-by-event (3+1)-D ideal hydrodynamical model with fluctuating initial conditions and the a multiphase transport (AMPT) Monte Carlo model for high-energy heavy-ion collisions. Anisotropic flows at different pseudorapidities are found to become significantly decorrelated with increasing pseudo-rapidity gaps due to longitudinal fluctuations in the initial states of heavy-ion collisions. The longitudinal correlation of the elliptic flow shows a strong centrality dependence while the correlation of the triangular flow is independent of the centrality. Longitudinal fluctuations as a source of the decorrelation are further shown to consist of a twist or gradual rotation in flow angles between the forward and backward direction and additional fluctuations on top of the twist. Within the AMPT model, longitudinal correlations of anisotropic flows are also found to depend on the value of partonic cross sections. The implicatiosn of constraining the initial conditions and shear viscosity to entropy density ratio of the partonic matter in high-energy heavy-ion collisions are also discussed., Comment: 12 pages and 11 figures, final published version

Centrality dependence of charged particles multiplicity, transverse momentum spectra, integrated and differential elliptic flow, in $\sqrt{s}_{NN}$=200 GeV Au+Au collisions are analyzed using event by event ideal hydrodynamics. Monte-Carlo Glauber model of initial condition, constrained to reproduce experimental charged particle's multiplicity in 0-10% Au+Au collisions, reasonably well reproduces all the experimental observables, e.g. centrality dependence of charged particles multiplicity, integrated and differential elliptic flow. Model predictions for higher flow harmonics, $v_3$, $v_4$ however overestimate the experimental data, more in the peripheral collisions than in the central collisions. Correlation between initial (spatial) asymmetry measures and flow coefficients are also studied. With exception of the elliptic flow, for all the higher flow coefficients ($v_n$,n=3-5), correlation is reduced with collision centrality. In peripheral collisions, higher flow coefficients are only weakly correlated to the asymmetry measures. Elliptic flow however, remains strongly correlated with initial eccentricity in all the collision centralities., Comment: 9 pages, 10 figures. arXiv admin note: text overlap with arXiv:1210.2249

The experimentally measured elliptic ($v_{2}$) and hexadecapole ($v_{4}$) flow of charged particles as a function of transverse momentum ($p_{T}$) at midrapidity in Pb-Pb collisions at $\sqrt{s_{\mathrm NN}}$ = 2.76 TeV are compared with the relativistic viscous hydrodynamic model simulations. The simulations are carried out for two different initial energy density profiles obtained from (i) Glauber model, and (ii) Color Glass Condensate (CGC) model. Comparison to experimental data for 10-20% to 40-50% centrality, shows that a centrality dependent shear viscosity to entropy density ($\eta/s$) ratio with values ranging between 0.0 to 0.12 are needed to explain the $v_{2}$ data for simulations with the Glauber based initial condition. Whereas for the CGC based initial conditions a slightly higher value of $\eta/s$ is preferred, around 0.08 to 0.16. From the comparison of the $v_{4}$ simulated results to the corresponding experimental measurements we observe that for the centralities 20-30% to 40-50% the $\eta/s$ values lies between 0.0 to 0.12 for both the initial conditions studied. The $\eta/s$ values obtained from our studies for Pb-Pb collisions at $\sqrt{s_{\mathrm NN}}$ = 2.76 TeV are compared to other studies which uses both transport and hydrodynamic approaches., Comment: 14 pages, 6 figures. One new figure added, text modified. arXiv admin note: text overlap with arXiv:1204.2347

In heavy ion collisions, event-by-event fluctuations in participating nucleon positions can lead to triangular flow. Generally, one uses Monte-Carlo Glauber model to obtain the participating nucleon positions. To use in a hydrodynamic model, the positions needs to be smoothened. We study the effect of smoothening of Glauber Monte-Carlo initial conditions on elliptic and triangular flow. It is shown that integrated as well as differential elliptic and triangular flow remain largely unaltered, irrespective of functional form of the smoothening function, or the smoothening parameter, Comment: 4 pages, 4 figures

Simulated results from a 2+1D relativistic viscous hydrodynamic model have been compared to the experimental data on the centrality dependence of invariant yield, elliptic flow ($v_{2}$), and hexadecapole flow ($v_{4}$) as a function of transverse momentum ($p_{T}$) of charged hadrons in Au-Au collisions at $\sqrt{s_{\rm {NN}}}$ = 200 GeV. Results from two types of initial transverse energy density profile, one based on the Glauber model and other based on Color-Glass-Condensate (CGC) are presented. We observe no difference in the simulated results on the invariant yield of charged hadrons for the calculations with different initial conditions. The comparison to the experimental data on invariant yield of charged hadrons supports a shear viscosity to entropy density ratio ($\eta/s$) between 0 to 0.12 for the 0-10% to 40-50% collision centralities. The simulated $v_{2}(p_{T})$ is found to be higher for a fluid with CGC based initial condition compared to Glauber based initial condition for a given collision centrality. Consequently the Glauber based calculations when compared to the experimental data requires a lower value of $\eta/s$ relative to CGC based calculations. In addition, a centrality dependence of the estimated $\eta/s$ is observed from the $v_{2}(p_{T})$ study. The $v_{4}(p_{T})$ for the collision centralities 0-10% to 40-50% supports a $\eta/s$ value between 0 - 0.08 for a CGC based initial condition. While simulated results using the Glauber based initial condition for the ideal fluid evolution under estimates the $v_{4}(p_{T})$ for collision centralities 0-10% to 30-40%., Comment: 11 pages, 11 figures, Accepted for publication in Physical Review C

The effect of a temperature dependent bulk viscosity to entropy density ratio~($\zeta/s$) along with a constant shear viscosity to entropy density ratio~($\eta/s$) on the space time evolution of the fluid produced in high energy heavy ion collisions have been studied in a relativistic viscous hydrodynamics model. The boost invariant Israel-Stewart theory of causal relativistic viscous hydrodynamics is used to simulate the evolution of the fluid in 2 spatial and 1 temporal dimension. The dissipative correction to the freezeout distribution for bulk viscosity is calculated using Grad's fourteen moment method. From our simulation we show that the method is applicable only for $\zeta/s<0.004$., Comment: 4 pages, 4 figures, Contribution to the Proceeding of the 7th International Conference on Critical Point and Onset of Deconfinement, 7-11 November,2011, Wuhan, China