Search

In this letter, we demonstrate for the first time that the moment expansion for an ultrarelativistic gas undergoing Bjorken flow diverges. We then show how this series can be resummed using the Borel-Pad\'e method and use this to determine the single-particle distribution function of the gas. Finally, we compare the exact resummed solution of the single-particle distribution function with solutions of the Boltzmann equation in the hydrodynamic limit and verify that the system displays considerable deviations from local equilibrium., Comment: Main text (8 pages, 5 figures) + Supplemental Material (12 pages, 6 figures)

We investigate the first-order transport coefficients of a fluid made of quasiparticles with a temperature-dependent mass extracted from chiral models. We describe this system using an effective kinetic theory, given by the relativistic Boltzmann equation coupled to a temperature-dependent background field determined from the thermal masses. We then simplify the collision term using the relaxation time approximation and implement a Chapman-Enskog expansion to calculate all first-order transport coefficients. In particular, we compute the bulk and shear viscosities using thermal masses extracted from the linear sigma model coupled with constituent quarks and the NJL model., Comment: 11 pages, 17 figures

We argue that the ratio between the shear viscosity and the shear relaxation time, $\eta/\tau_\pi$, should be defined as a thermodynamic quantity obtained from the equal-time symmetric correlator of the shear-stress tensor. In kinetic theory, we show that this ratio does not depend on the type of interaction. Similarly, an exact expression for this ratio is obtained for holographic gauge theories. We also determine how stochastic fluctuations change $\eta/\tau_\pi$ in transient relativistic hydrodynamics and show that thermal fluctuations do not spoil causality and stability., Comment: 13 pages (including Appendix)

We propose a new observable derived from a centrality-dependent scaling of transverse particle spectra. By removing the global scales of total particle number and mean transverse momentum, we isolate the shape of the spectrum. In hydrodynamic simulations, while the multiplicity and mean transverse momentum fluctuate significantly, the scaled spectrum is found to be almost constant even at an event-by-event level and after resonance decays. This universality survives when averaging over events in each centrality bin before scaling. We then investigate the presence of this scaling in experimental data from the ALICE collaboration in Pb-Pb, Xe-Xe, and p-Pb collisions. We find a remarkable universality in the experimentally observed scaled spectra at low transverse momentum, compatible with hydrodynamic predictions. The data show a minor breaking of universality at large transverse momentum and hints of evolution with the system size that are not seen in simulations. Our results motivate further theoretical and experimental investigations of this new observable to bring to light the collective and non-collective behavior encoded in the transverse particle spectrum of different collision systems., Comment: V1: 6 pages, 3 figures

Using an analytical result for the eigensystem of the linearized collision term for a classical system of massless scalar particles with quartic self-interactions, we show that the shear-stress linear response function possesses a branch-cut singularity that covers the whole positive imaginary semi-axis. This is demonstrated in two ways: (1) by truncating the exact, infinite linear system of linear equations for the rank-two tensor modes, which reveals the cut touching the origin; and (2) by employing the Trotterization techniques to invert the linear response problem. The former shows that the first pole tends towards the origin and the average separation between consecutive poles tends towards zero as power laws in the dimension of the basis. The latter allows one to obtain the response function in closed form in terms of Tricomi hypergeometrical functions, which possess a branch-cut on the above-mentioned semi-axis. This suggests that the presence of a cut along the imaginary frequency axis of the shear stress correlator, inferred from previous numerical analyses of weakly coupled scalar $\lambda \varphi^4$ theories, does not arise due to quantum statistics but instead emerges from the fundamental properties of this system's interactions., Comment: 18 pages, 3 figures, v2: published version

We calculate all transport coefficients of second order transient hydrodynamics in two effective kinetic theory models: a hadron-resonance gas and a quasiparticle model with thermal masses tuned to reproduce QCD thermodynamics. We compare the corresponding results with calculations for an ultrarelativistic single-component gas, that are widely employed in hydrodynamic simulations of heavy ion collisions. We find that both of these effective models display a qualitatively different normalized bulk viscosity, when compared to the calculation for the single-component gas. Indeed, $\zeta/[\tau_{\Pi}(\varepsilon_{0} + P_{0})] \simeq 16.91(1/3-c_{s}^{2})^{2}$, for the hadron-resonance gas model, and $\zeta/[\tau_{\Pi}(\varepsilon_{0} + P_{0})] \simeq 5 (1/3-c_{s}^{2})$ for the quasiparticle model. Differences are also observed for many second-order transport coefficients, specially those related to the bulk viscous pressure. The transport coefficients derived are shown to be consistent with fundamental linear stability and causality conditions., Comment: 21 pages, 5 figures v2 - published version

We derive the equations of motion of relativistic magnetohydrodynamics, as well as microscopic expressions for all of its transport coefficients, from the Boltzmann equation using the method of moments. In contrast to reference Phys. Rev. D 98(7) 2018, where a single component gas was considered, we perform our derivation for a locally neutral fluid composed of two massless particle species with opposite charges. We demonstrate that the magnetohydrodynamical equations of motion become dramatically different for this more realistic system. The shear-stress tensor no longer obeys a single differential equation; it breaks into three non-degenerate components with respect to the magnetic field, each evolving according to different dynamical equations. For large magnetic fields, we further show that the solution of this theory displays oscillatory behaviour that can no longer be described by an Israel-Stewart-like theory. Finally, we investigate the derived equations in a Bjorken flow scenario., Comment: 18 pages, 12 figures

Relativistic dissipative fluid dynamics finds widespread applications in high-energy nuclear physics and astrophysics. However, formulating a causal and stable theory of relativistic dissipative fluid dynamics is far from trivial; efforts to accomplish this reach back more than 50 years. In this review, we give an overview of the field and attempt a comparative assessment of (at least most of) the theories for relativistic dissipative fluid dynamics proposed until today and used in applications., Comment: 73 pages, no figures, invited review for journal "Entropy"

We investigate the causality and stability of three different relativistic dissipative fluid-dynamical formulations emerging from a system of classical, ultra-relativistic scalar particles self-interacting via a quartic potential. For this particular interaction, all transport coefficients of Navier-Stokes, Bemfica-Disconzi-Noronha-Kovtun and second-order transient theories can be computed in analytical form. We first show that Navier-Stokes theory is acausal and unstable regardless of the matching conditions. On the other hand, BDNK theory can be linearly causal and stable for a particular set of matching choices that does not contain the so-called exotic Eckart prescription. In particular, using the Li\'enard-Chipart criterion, we obtain a set of sufficient conditions that guarantee the stability of the theory. Last, second-order transient hydrodynamic theory in Landau matching is shown to be linearly causal and stable., Comment: 21 pages, 5 figures

Heavy-ion collisions, such as Pb-Pb or p-Pb, produce extreme conditions in temperature and density that make the hadronic matter transition to a new state, called quark-gluon plasma (QGP). Simulations of heavy-ion collisions provide a way to improve our understanding of the QGP's properties. These simulations are composed of a hybrid description that results in final observables in agreement with accelerators like LHC and RHIC. However, recent works pointed out that these hydrodynamic simulations can display acausal behavior during the evolution in certain regions, indicating a deviation from a faithful representation of the underlying QCD dynamics. To pursue a better understanding of this problem and its consequences, this work simulated two different collision systems, Pb-Pb and p-Pb at $\sqrt{s_{NN}} = 5.02$ TeV. In this context, our results show that causality violation, even though always present, typically occurs on a small part of the system, quantified by the total energy fraction residing in the acausal region. In addition, the acausal behavior can be reduced with changes in the pre-hydrodynamic factors and the definition of the bulk-viscous relaxation time. Since these aspects are fairly arbitrary in current simulation models, without solid guidance from the underlying theory, it is reasonable to use the disturbing presence of acausal behavior in current simulations to guide improvements towards more realistic modeling. While this work does not solve the acausality problem, it sheds more light on this issue and also proposes a way to solve this problem in simulations of heavy-ion collisions., Comment: 15 pages, 17 figures

We derive and investigate several hydrodynamic formalisms that emerge from a system of classical, ultra-relativistic scalar particles self-interacting via a quartic potential. The specific form of the total cross-section enables the analytical computation of all transport coefficients that appear in Navier-Stokes (NS), Bemfica-Disconzi-Noronha-Kovtun (BDNK), and second-order transient hydrodynamic theories. We solve all these formalisms in a Bjorken flow scenario and show that NS and BDNK theories display unphysical features when gradients become sufficiently large. This implies that these hydrodynamic approaches may not be suitable to describe the early stages of heavy ion collisions., Comment: 30 pages, 8 figures

We derive a linearly causal and stable third-order relativistic fluid-dynamical theory from the Boltzmann equation using the method of moments. For this purpose, we demonstrate that such theory must include novel degrees of freedom, corresponding to irreducible tensors of rank 3 and 4. The equations of motion derived in this work are compared with numerical solutions of the Boltzmann equation, considering an ultrarelativistic, classical gas in the highly symmetric Bjorken flow scenario. These solutions are shown to be in good agreement for a wide range of values of shear viscosity and initial temperatures., Comment: 16 pages, 4 figures

We extend our previous investigation of the effects of pre-hydrodynamic evolution on final-state observables in heavy-ion collisions to smaller systems. We use a state-of-the-art hybrid model for the numerical simulations with optimal parameters obtained from a previous Bayesian study. By studying p-Pb collisions, we find that the effects due to the assumption of a conformal evolution in the pre-hydrodynamical stage are even more important in small systems. We also show that this effect depends on the time duration of the pre-equilibrium stage, which is further enhanced in small systems. Finally, we show that the recent proposal of a free-streaming with subluminal velocity for the pre-equilibrium stage, thus effectively breaking conformal invariance, can alleviate the contamination of final state observables. Our study further reinforces the need for moving beyond conformal approaches in pre-equilibrium dynamics modeling, especially when extracting transport coefficients from hybrid models in the high-precision era of heavy-ion collisions., Comment: 15 pages, 14 figures

We analytically determine all the eigenvalues and eigenfunctions of the linearized Boltzmann collision operator in massless scalar $\lambda \phi^4$ theory in the high-temperature (classical) regime. This is used to exactly compute the shear viscosity and particle diffusion transport coefficients of this system. The corresponding relaxation time approximation for this linearized Boltzmann equation is also derived., Comment: 10 pages including supplemental material

We use a novel formulation of the relaxation time approximation to consistently calculate the bulk and shear viscosity coefficients using QCD-inspired energy-dependent relaxation times and phenomenological thermal masses obtained from fits to lattice QCD thermodynamics. The matching conditions are conveniently chosen to simplify the computations., Comment: Proceedings for the XXIXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions (Quark Matter 2022). 6 pages, 1 figure

We investigate the far-from-equilibrium behavior of the Boltzmann equation for a gas of massless scalar field particles with quartic (tree level) self-interactions ($\lambda \phi^4$) in Friedmann-Lemaitre-Robertson-Walker spacetime. Using a new covariant generating function for the moments of the Boltzmann distribution function, we analytically determine a subset of the spectrum and the corresponding eigenfunctions of the linearized Boltzmann collision operator. We show how the covariant generating function can be also used to find the exact equations for the moments in the full nonlinear regime. Different than the case of a ultrarelativistic gas of hard spheres (where the total cross section is constant), for $\lambda \phi^4$ the fact that the cross section decreases with energy implies that moments of arbitrarily high order directly couple to low order moments. Numerical solutions for the scalar field case are presented and compared to those found for a gas of hard spheres., Comment: 19 pages, 5 figures

Hydrodynamics can be formulated in terms of a perturbative series in derivatives of the temperature, chemical potential, and flow velocity around an equilibrium state. Different formulations for this series have been proposed over the years, which consequently led to the development of various hydrodynamic theories. In this work, we discuss the relativistic generalizations of the perturbative expansions put forward by Chapman and Enskog, and Hilbert, using general matching conditions in kinetic theory. This allows us to describe, in a comprehensive way, how different out-of-equilibrium definitions for the hydrodynamic fields affect the development of the hydrodynamic perturbative series. We provide a perturbative method for systematically deriving the hydrodynamic formulation recently proposed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK) from relativistic kinetic theory. The various transport coefficients that appear in BDNK (at first-order) are explicitly computed using a new formulation of the relaxation time approximation for the Boltzmann equation. Assuming Bjorken flow, we also determine the hydrodynamic attractors of BDNK theory and compare the overall hydrodynamic evolution obtained using this formulation with that generated by the Israel-Stewart equations of motion and also kinetic theory., Comment: 36 pages, 8 figures

We investigate the transport properties of a kinetic theory model that is tuned to describe the thermodynamic properties of QCD at zero chemical potential using a new formulation of the relaxation time approximation. In contrast to previous approaches, the latter is constructed to preserve the fundamental properties of the collision term of the Boltzmann equation for any energy-dependence of the relaxation time. A novel choice of matching conditions is implemented to ensure that the background mean-field depends only on the temperature even when the system is out of equilibrium. We provide a consistent analysis of how the transport coefficients of relativistic Navier-Stokes theory vary with the energy dependence of the relaxation time. We also show that the entropy production of this theory is consistent with the second law of thermodynamics and verify that it is independent of the matching conditions employed. We used this fact to calculate the matching independent combination of transport coefficients., Comment: 15 pages, 8 figures

Recent works have revealed that matching conditions play a major role on general consistency properties of relativistic fluid dynamics such as causality, stability and wellposedness of the equations of motion. In this paper we derive transient fluid dynamics from kinetic theory, using the method of moments as proposed by Israel and Stewart, without imposing an specific matching condition. We then investigate how the equations of motion and their corresponding transport coefficients are affected by the choice of matching condition., Comment: 20 pages, 7 figures

We demonstrate that a Bjorken expanding strongly coupled $\mathcal{N}=4$ Supersymmetric Yang-Mills plasma can display dynamically-driven violations of the dominant and also the weak energy condition during hydrodynamization. In addition, we find that a period of vanishing entropy production in far-from-equilibrium stages induces later violations of the dominant energy condition in the strongly coupled plasma. Such violations cannot occur in a classical description of hydrodynamization and suggest that the inclusion of quantum effects in transport can lead to new phenomena in these regimes, even for systems without anomalies or spin., Comment: 9 pages, 1 figure; v3: to appear in Physical Review D

We show that the widely used relaxation time approximation to the relativistic Boltzmann equation contains basic flaws, being incompatible with microscopic and macroscopic conservation laws. We propose a new approximation that fixes such fundamental issues and maintains the basic properties of the linearized Boltzmann collision operator. We show how this correction affects transport coefficients, such as the bulk viscosity and particle diffusion., Comment: 7 pages, 4 figures

We investigate the effects of pre-hydrodynamic evolution on final-state observables in heavy-ion collisions using state-of-the art event simulations coupled to different pre-hydrodynamic scenarios, which include the recently-developed effective kinetic transport theory evolution model KoMPoST. Flow observables are found to be insensitive to the details of pre-hydrodynamic evolution. The main effect we observe is in the $p_T$ spectra, particularly the mean transverse momentum. However, at least part of this effect is a consequence of the underlying conformal invariance assumption currently present in such approaches, which is known to be violated in the temperature regime probed in heavy-ion collisions. This assumption of early time conformal invariance leads to an artificially large out-of-equilibrium bulk pressure when switching from (conformal) pre-hydrodynamic evolution to hydrodynamics (using the non-conformal QCD equation of state), which in turn increases the transverse momentum. Our study indicates that a consistent treatment of pre-hydrodynamic evolution in heavy-ion collisions requires the use of non-conformal models of early time dynamics., Comment: 17 pages, 15 figures, 1 appendix

We investigate whether hydrodynamic attractors are present in simulations of the quark-gluon plasma formed in heavy-ion collisions. We argue that Lagrangian schemes to solve the relativistic viscous fluid equations can be particularly useful to characterize the properties of attractors in heavy ion collisions. Preliminary results are shown in Israel-Stewart theory to support such a claim. We also discuss how to perform the slow-roll expansion in Israel-Stewart theory undergoing a general flow, which naturally leads to the definition of a resummed shear viscosity coefficient that depends on the gradients of the hydrodynamic fields., Comment: 4 pages, 1 figure, contribution to the Quark Matter 2019 proceedings

In general, the constituents of the bulk matter produced in heavy-ion collisions carry, besides electric charge, multiple other conserved quantum numbers like baryon number and strangeness. Therefore, an electric field will not only generate an electric current but, at the same time, also currents in baryon number and strangeness. We propose that the impact of the electric field on these conserved currents should be characterized by additional transport coefficients, which we call cross-conductivities. In this paper, we introduce and present a calculation of these cross-conductivities from the Green-Kubo formalism within the transport code SMASH for different chemical compositions of hadron resonance gases. We find that the coefficients underlie an ordering in the active degrees of freedom and that thus the chemical composition of the system plays a crucial role. Further, we argue that in future comparisons of lattice QCD calculations with these findings, one could constrain which degrees of freedom and their corresponding charge properties are relevant for the QCD dynamics of the system., Comment: 6 pages, 5 figures

The bulk nuclear matter produced in heavy ion collisions carries a multitude of conserved quantum numbers: electric charge, baryon number, and strangeness. Therefore, the diffusion processes associated to these conserved charges cannot occur independently and must be described in terms of a set of coupled diffusion equations. This physics is implemented by replacing the traditional diffusion coefficients for each conserved charge by a diffusion coefficient matrix, which quantifies the coupling between the conserved quantum numbers. The diagonal coefficients of this matrix are the usual charge diffusion coefficients, while the off-diagonal entries describe the diffusive coupling of the charge currents. In this paper, we show how to calculate this diffusion coefficient matrix from kinetic theory and provide results for a hadron resonance gas and a gas of partons. We further find that the off-diagonal entries can reach similar magnitudes compared to the diagonal entries. In order to provide some insight on the influence that the coupling between the net charge diffusion currents can have on heavy ion observables, we present first results for the diffusive evolution of a hadronic system in a simple (1+1)D-fluid dynamics approach, and study different configurations of the diffusion matrix., Comment: 27 pages, 10 Figures

We derive the general analytical solution of the viscous hydrodynamic equations for an ultrarelativistic gas of hard spheres undergoing Bjorken expansion, taking into account effects from particle number conservation, and use it to analytically determine its attractor at late times. Differently than all the cases considered before involving rapidly expanding fluids, in this example the gradient expansion converges. We exactly determine the hydrodynamic attractor of this system when its microscopic dynamics is modeled by the Boltzmann equation with a fully nonlinear collision kernel. The exact late time attractor of this system can be reasonably described by hydrodynamics even when the gradients are large., Comment: 7 pages, 2 figures; version accepted for publication in Physical Review Letters

The collective behavior of hadrons and of electromagnetic radiation in heavy-ion collisions has been widely used to study the properties of the the Quark-Gluon Plasma (QGP). Indeed this collectivity, as measured by anisotropic flow coefficients, can be used to constrain the transport properties of QGP. The goal of this contribution is to investigate the influence of the specific bulk viscosity ($\zeta/s$) on dilepton production, both at Relativistic Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) energies. We explore the sensitivity of dileptons to dynamical features that bulk viscosity induces on the evolution of a strongly-interacting medium, and highlight what makes them a valuable probe in the pursuit to also constrain $\zeta/s$., Comment: 25 pages, 16 figures

We derive the equations of motion of relativistic, resistive, second-order dissipative magnetohydrodynamics from the Boltzmann-Vlasov equation using the method of moments. We thus extend our previous work [Phys. Rev. D 98, 076009 (2018)], where we only considered the non-resistive limit, to the case of finite electric conductivity. This requires keeping terms proportional to the electric field $E^\mu$ in the equations of motions and leads to new transport coefficients due to the coupling of the electric field to dissipative quantities. We also show that the Navier-Stokes limit of the charge-diffusion current corresponds to Ohm's law, while the coefficients of electrical conductivity and charge diffusion are related by a type of Wiedemann-Franz law., Comment: 11 pages

A hybrid (hydrodynamics + hadronic transport) theoretical framework is assembled to model the bulk dynamics of relativistic heavy-ion collisions at energies accessible in the Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC) and the NA61/SHINE experiment at CERN. The system's energy-momentum tensor and net baryon current are evolved according to relativistic hydrodynamics with finite shear viscosity and non-zero net baryon diffusion. Our hydrodynamic description is matched to a hadronic transport model in the dilute region. With this fully integrated theoretical framework, we present a pilot study of the hadronic chemistry, particle spectra, and anisotropic flow. Phenomenological effects of a non-zero net-baryon current and its diffusion on hadronic observables are presented for the first time. The importance of the hadronic transport phase is also investigated., Comment: 22 pages, 22 figures

We derive the equations of motion of relativistic, non-resistive, second-order dissipative magnetohydrodynamics from the Boltzmann equation using the method of moments. We assume the fluid to be composed of a single type of point-like particles with vanishing dipole moment or spin, so that the fluid has vanishing magnetization and polarization. In a first approximation, we assume the fluid to be non-resistive, which allows to express the electric field in terms of the magnetic field. We derive equations of motion for the irreducible moments of the deviation of the single-particle distribution function from local thermodynamical equilibrium. We analyze the Navier-Stokes limit of these equations, reproducing previous results for the structure of the first-order transport coefficients. Finally, we truncate the system of equations for the irreducible moments using the 14-moment approximation, deriving the equations of motion of relativistic, non-resistive, second-order dissipative magnetohydrodynamics. We also give expressions for the new transport coefficients appearing due to the coupling of the magnetic field to the dissipative quantities., Comment: 18 pages, 2 figures; Discussions, references and corrections added. Fig. 1a corrected, typos fixed

Perturbative expansions, such as the well-known gradient series and the recently proposed slow-roll expansion, have been recently used to investigate the emergence of hydrodynamic behavior in systems undergoing Bjorken flow. In this paper we determine for the first time the large order behavior of these perturbative expansions in relativistic hydrodynamics in the case of Gubser flow. While both series diverge, the slow-roll series can provide a much better overall description of the system's dynamics than the gradient expansion when both series are truncated at low orders. The truncated slow-roll series can also describe the attractor solution of Gubser flow as long as the system is sufficiently close to equilibrium near the origin (i.e., $\rho=0$) in $dS_3 \otimes \mathbb{R}$. Differently than the case of Bjorken flow, here we show that the Gubser flow attractor solution is not solely a function of the effective Knudsen number $\tau_R \sqrt{\sigma_{\mu\nu}\sigma^{\mu\nu}} \sim \tau_R\, \tanh\rho$. Our results give further support to the idea that new \emph{resummed} constitutiv relations between dissipative currents and the gradients of conserved quantities can emerge in systems far from equilibrium that are beyond the regime of validity of the usual gradient expansion., Comment: 24 pages, 19 figures; minor modifications

In order to characterize nuclear matter under extreme conditions, we calculate all diffusion transport coefficients related to baryon, electric and strangeness charge for both a hadron resonance gas and a simplified kinetic model of the quark-gluon plasma. We demonstrate that the diffusion currents do not depend only on gradients of their corresponding charge density. Instead, we show that there exists coupling between the different charge currents, in such a way that it is possible for density gradients of a given charge to generate dissipative currents of another charge. Within this scheme, the charge diffusion coefficient is best viewed as a matrix, in which the diagonal terms correspond to the usual charge diffusion coefficients, while the off-diagonal terms describe the coupling between the different currents. In this letter, we calculate for the first time the complete diffusion matrix including the three charges listed above. We find that the baryon diffusion current is strongly affected by baryon charge gradients, but also by its coupling to gradients in strangeness. The electric charge diffusion current is found to be strongly affected by electric and strangeness gradients, whereas strangeness currents depend mostly on strange and baryon gradients., Comment: 6 pages, 2 figures

We find the general analytical solution of the viscous relativistic hydrodynamic equations (in the absence of bulk viscosity and chemical potential) for a Bjorken expanding fluid with a constant shear viscosity relaxation time. We analytically determine the hydrodynamic attractor of this fluid and discuss its properties. We show for the first time that the slow-roll expansion, a commonly used approach to characterize the attractor, diverges. This is shown to hold also in a conformal plasma. The gradient expansion is found to converge in an example where causality and stability are violated., Comment: 22 pages, 7 figures

We determine the dynamical attractors associated with anisotropic hydrodynamics (aHydro) and the DNMR equations for a 0+1d conformal system using kinetic theory in the relaxation time approximation. We compare our results to the non-equilibrium attractor obtained from exact solution of the 0+1d conformal Boltzmann equation, Navier-Stokes theory, and second-order Mueller-Israel-Stewart theory. We demonstrate that the aHydro attractor equation resums an infinite number of terms in the inverse Reynolds number. The resulting resummed aHydro attractor possesses a positive longitudinal to transverse pressure ratio and is virtually indistinguishable from the exact attractor. This suggests that kinetic theory involves not only a resummation in gradients (Knudsen number) but also a novel resummation in inverse Reynolds number. We also demonstrate that the DNMR result provides a better approximation to the exact kinetic theory attractor than Mueller-Israel-Stewart theory. Finally, we introduce a new method for obtaining approximate aHydro equations which relies solely on an expansion in inverse Reynolds number, carry out this expansion to third order, and compare these third-order results to the exact kinetic theory solution., Comment: 25 pages, 7 figures

In this proceedings we discuss recent findings regarding the large order behavior of the Chapman-Enskog expansion in relativistic kinetic theory. It is shown that this series in powers of the Knudsen number has zero radius of convergence in the case of a Bjorken expanding fluid described by the Boltzmann equation in the relaxation time approximation. This divergence stems from the presence of non-hydrodynamic modes, which give non-perturbative contributions to the Knudsen series., Comment: 4 pages, 1 figure, proceedings for Quark Matter 2017

Starting from IP-Glasma initial conditions, we investigate the effects of bulk pressure on low mass dilepton production at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) energies. Though thermal dilepton $v_2$ is affected by the presence of both bulk and shear viscosity, whether or not these effects can be measured depends on the dilepton "cocktail" contribution to the the low mass dilepton $v_2$. Combining the thermal and "cocktail" dileptons, the effects of bulk viscosity on total dilepton $v_2$ is investigated., Comment: Proceedings for the 26th International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2017), February 5-11 2017, Chicago, Illinois, USA

Using a hybrid (viscous hydrodynamics + hadronic cascade) framework, we model the bulk dynamical evolution of relativistic heavy-ion collisions at Relativistic Heavy Ion Collider (RHIC) Beam Energy Scan (BES) collision energies, including the effects from non-zero net baryon current and its dissipative diffusion. This framework is in full (3+1)D, which allows us to study the non-trivial longitudinal structure and dynamics of the collision systems, for example baryon stopping and transport, as well as longitudinal fluctuations. For the first time, the quantitative effect of net-baryon diffusion on hadronic observables is addressed. Finally, we propose a dynamical initialization scheme to study the importance of the pre-equilibrium stage at the RHIC BES energies., Comment: 4 pages, 2 figures, proceedings for Quark Matter 2017

We describe ultra-relativistic heavy ion collisions at RHIC and the LHC with a hybrid model using the IP-Glasma model for the earliest stage and viscous hydrodynamics and microscopic transport for the later stages of the collision. We demonstrate that within this framework the bulk viscosity of the plasma plays an important role in describing the experimentally observed radial flow and azimuthal anisotropy simultaneously. We further investigate the dependence of observables on the temperature below which we employ the microscopic transport description.

Starting from IP-Glasma initial conditions, we investigate the effects of bulk pressure on thermal dilepton production at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) energies. Though results of the thermal dilepton $v_2$ under the influence of both bulk and shear viscosity is presented for top RHIC energy, more emphasis is put on LHC energy where such a calculation is computed for the first time. The effects of the bulk pressure on thermal dilepton $v_2$ at the LHC are explored through bulk-induced modifications on the dilepton yield., Comment: Talk given at the 8th International Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions (Hard Probes 2016), September 23-27 2016, Wuhan, China; 4 pages, LaTeX, 2 PDF figures

This work reports on investigations of the effects on the evolution of viscous hydrodynamics and on the flow coefficients of thermal dileptons, originating from a temperature-dependent specific shear viscosity $\eta/s (T)$ at temperatures beyond 180 MeV formed at the Relativistic Heavy-Ion Collider (RHIC). We show that the elliptic flow of thermal dileptons can resolve the magnitude of $\eta/s$ at the high temperatures, where partonic degrees of freedom become relevant, whereas discriminating between different specific functional forms will likely not be possible at RHIC using this observable., Comment: 25 pages, 20 figures