Search

We formulate the first-order dissipative anisotropic hydrodynamical theory for a relativistic conformal uncharged fluid, which generalizes the Bemfica-Disconzi-Noronha-Kovtun first-order viscous fluid framework. Our approach maintains causal behavior in the nonlinear regime with or without general relativity coupling, and we derive and analyze the constraints on transport coefficients imposed by causality. We demonstrate the causal and stable behavior of our theory in specific cases, including the discussion of nonlinear causality as well as stability for linearized perturbations. We apply our newly developed first-order anisotropic theory to the Bjorken flow and show how causality and stability impose constraints on the behavior of the early-time attractor., Comment: Matches the published version

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far-from-equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. K$\phi$MP$\phi$ST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding non-equilibrium Green's functions in the relaxation time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g \rightarrow q\bar{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. We show that this method alters the initial state charge geometries and is applicable in central to mid-central collisions., Comment: 66 pages, 18 figures

While it is well known that there is a significant amount of conserved charges in the initial state of nuclear collisions, the production of these due to gluon splitting has yet to be thoroughly investigated. The ICCING (Initial Conserved Charges in Nuclear Geometry) algorithm reconstructs these quark distributions, providing conserved strange, baryon, and electric charges, by sampling a given model for the $g \rightarrow q\bar{q}$ splitting function over the initial energy density, which is valid at top collider energies, even when $\mu_B=0$. The ICCING algorithm includes fluctuations in the gluon longitudinal momenta, a structure that supports the implementation of dynamical processes, and the c++ version is now open-source. A full analysis of parameter choices on the model has been done to quantify the effect these have on the underlying physics. We find there is a sustained difference across the different charges that indicates sensitivity to hot spot geometry.

In this work we investigate the impact of conformal symmetry breaking on hydrodynamization of a far-from-equilibrium fluid. We find a new kind of transseries solutions for the non-conformal hydrodynamic equations of a longitudinal boost invariant expanding plasma. The new transseries solutions unveil a rich physical structure which arises due to the interplay of different physical scales. In the perfect fluid case the non-conformal speed of sound slows down the cooling of the temperature due to the emergence of logarithmic corrections that depends on the mass of the particle. These terms propagate into the perturbative and non-perturbative sectors of the transseries once viscous corrections are included. The logarithmic mass contributions increase the asymptotic value of the Knudsen number while decreasing the damping rate of the transient non-hydrodynamic modes and thus, yielding to an extremely slow hydrodynamization process where flow lines merge to their forward attractor at extremely late times. The early time free streaming expansion is modified and receives logarithmic mass corrections induced by the shear-bulk couplings. The global flow structure and numerical analyses carried out in our work demonstrate the existence of the early and late-time attractors for the shear viscous tensor and bulk viscous pressure., Comment: 37 pages (main text) + 11 pages (appendices + references), 12 figures. v2: Improved discussion, added new figure

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Finally, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world., Comment: 103 pages

It is frequently supposed that quark-gluon plasma created in heavy-ion collisions undergoes free streaming at early times. We examine this issue based on the assumption that a universal attractor dominates the dynamics already at the earliest stages, which offers a way to connect the initial state with the start of the hydrodynamic expansion in an approximate but conceptually transparent fashion. We demonstrate that the centrality dependence of the measured particle multiplicities can be used to quantitatively constrain the pressure anisotropy and find that it strongly depends on the model of the initial energy deposition. As an illustration, we compare three initial state models and show that they predict rather different early-time values of the pressure anisotropy. This strongly suggests that assuming free streaming prior to hydrodynamization is not necessarily compatible with a generic initial state model and that features of the pre-hydrodynamic flow need to be matched with the model of the initial state., Comment: v2: significantly extended; v3: presentation improved

We study the temporal expansion of an ultracold Bose gas in two-dimensional, square optical lattices. The gas is described by the Bose-Hubbard model deep in the superfluid regime, with initially all bosons condensed in the central site of the lattice. We use the previously developed nonequilibrium propagator method for capturing the time evolution of an interacting bosonic system, where the many-body Hamiltonian is represented in an appropriate local basis and the corresponding field operators are separated into the classical [Bose-Einstein condensate (BEC)] part and quantum mechanical fluctuations. After a quench, i.e. after a sudden switch of the lattice nearest-neighbor hopping, the expanding, bosonic cloud separates spatially into a fast, ballistic forerunner and a slowly expanding central part controlled by selftrapping. We show that the forerunner expansion is driven by the coherent dynamics of the BEC and that its velocity is consistent with the Lieb-Robinson bound. For smaller lattices we analyze how quasiparticle collisions lead to enhanced condensate depletion and oscillation damping., Comment: 13 pages, 8 figures; published version, extended discussion, improved figure presentation

Non-equilibrium Green's functions provide an efficient way to describe the evolution of the energy-momentum tensor during the early time pre-equilibrium stage of high-energy heavy ion collisions. Besides their practical relevance they also provide a meaningful way to address the question when and to what extent a hydrodynamic description of the system becomes applicable. Within the kinetic theory framework we derive a new method to calculate time dependent non-equilibrium Green's functions describing the evolution of energy and momentum perturbations on top of an evolving far-from-equilibrium background. We discuss the approach towards viscous hydrodynamics along with the emergence of various scaling phenomena for conformal systems. By comparing our results obtained in the relaxation time approximation to previous calculations in Yang-Mills kinetic theory, we further address the question which macroscopic features of the energy momentum tensor are sensitive to the underlying microscopic dynamics., Comment: 18 pages, 8 figures

At top collider energies where baryon stopping is negligible, the initial state of heavy ion collisions is overall charge neutral and predominantly composed of gluons. Nevertheless, there can also be significant local fluctuations of the baryon number, strangeness, and electric charge densities about zero, perturbatively corresponding to the production of quark/antiquark pairs. These previously ignored local charge fluctuations can permit the study of charge diffusion in the quark-gluon plasma (QGP), even at top collider energies. In this paper we present a new model denoted ICCING (Initial Conserved Charges in Nuclear Geometry) which can reconstruct the initial conditions of conserved charges in the QGP by sampling a ($g \rightarrow q\bar{q}$) splitting probability over the initial energy density. We find that the new charge distributions generally differ from the bulk energy density; in particular, the strangeness distribution is significantly more eccentric than standard bulk observables and appears to be associated with the geometry of hot spots in the initial state. The new information provided by these conserved charges opens the door to studying a wealth of new charge- and flavor-dependent correlations in the initial state and ultimately the charge transport parameters of the QGP., Comment: 60 pages, 22 figures

We initialize the Quantum Chromodynamic conserved charges of baryon number, strangeness, and electric charge arising from gluon splitting into quark-antiquark pairs for the initial conditions of relativistic heavy-ion collisions. A new Monte Carlo procedure that can sample from a generic energy density profile is presented, called Initial Conserved Charges in Nuclear Geometry (ICCING), based on quark and gluon multiplicities derived within the color glass condensate (CGC) effective theory. We find that while baryon number and electric charge have nearly identical geometries to the energy density profile, the initial strangeness distribution is considerable more eccentric and is produced primarily at the hot spots corresponding to temperatures of $T\gtrsim 400$ MeV for PbPb collisions at $\sqrt{s_{NN}}=5.02$ TeV., Comment: 7 pages, 8 figures

In this work we introduce the generic conditions for the existence of a non-equilibrium attractor that is an invariant manifold determined by the long-wavelength modes of the physical system. We investigate the topological properties of the global flow structure of the Gubser flow for the Israel-Stewart theory and a kinetic model for the Boltzmann equation by employing Morse-Smale theory. We present a complete classification of the invariant submanifolds of the flow and determine all the possible flow lines connecting any pair of UV/IR fixed points. The formal transseries solutions to the Gubser dynamical system around the early-time (UV) and late-time (IR) fixed points are constructed and analyzed. It is proven that these solutions are purely perturbative (or power-law asymptotic) series with a finite radius of convergence. Based on these analyses, we find that Gubser-like expanding kinetic systems do not hydrodynamize owing to the failure of the hydrodynamization process which heavily relies on the classification of (non)hydrodynamic modes in the IR regime. This is in contrast to longitudinal boost-invariant plasmas where the asymptotic dynamics is described by a few terms of the hydrodynamic gradient expansion. We finally compare our results for both Bjorken and Gubser conformal kinetic models., Comment: 46 pages, 5 figures, 1 table. Minor changes in the text and modifications in Fig.3

We study the behavior of the bulk viscosity $\zeta$ in QCD near a possible critical endpoint in the phase diagram. We verify the expectation that $(\zeta/s)\sim a(\xi/\xi_0)^{x_\zeta}$, where $s$ is the entropy density, $\xi$ is the correlation length, $\xi_0$ is the non-critical correlation length, $a$ is a constant and $x_\zeta\simeq 3$. Using a recently developed equation of state that includes a critical point in the universality class of the Ising model we estimate the constant of proportionality $a$. We find that $a$ is typically quite small, $a\sim O(10^{-4})$. We observe, however, that the result is sensitive to the commonly made assumption that the Ising temperature axis is approximately aligned with the QCD baryon chemical potential axis. If this is not the case, then the critical $\zeta/s$ can approach the non-critical value of $\eta/s$, where $\eta$ is the shear viscosity, even if the enhancement of the correlation length is modest, $\xi/\xi_0\sim 2$., Comment: 18 pages, 3 figures; version to appear in PRD

In this paper, we compute a number of cross sections for the production of multiple particles at mid-rapidity in the semi-dilute / dense regime of the color-glass condensate (CGC) effective field theory. In particular, we present new results for the production of two quark-antiquark pairs (whether the same or different flavors) and for the production of one quark-antiquark pair and a gluon. We also demonstrate the existence of a simple mapping which transforms the cross section to produce a quark-antiquark pair into the corresponding cross section to produce a gluon, which we use to obtain various results and to cross-check them against the literature. We also discuss hadronization effects in the heavy flavor sector, writing explicit expressions for the production of various combinations of $D$ and $\bar D$ mesons, $J/\psi$ mesons, and light hadrons. The various multiparticle cross sections presented here contain a wealth of information and can be used to study heavy flavor production, charge-dependent correlations, and "collective" flow phenomena arising from initial-state dynamics., Comment: 51 pages, 8 figures

We discuss the non-equilibrium attractors of systems undergoing Gubser flow within kinetic theory by means of nonlinear dynamical systems. We obtain the attractors of anisotropic hydrodynamics, Israel-Stewart (IS) and transient fluid (DNMR) theories. These attractors are non-planar and the basin of attraction is three dimensional. We compare the asymptotic attractors of each hydrodynamic model with the one obtained from the exact Gubser solution of the Boltzmann equation within the relaxation time approximation. Anisotropic hydrodynamics matches, up to high numerical accuracy, the attractor of the exact theory while the other hydrodynamic theories fail to do so. Thus, anisotropic hydrodynamics is an effective theory for far-from-equilibrium fluids, which consists of the dissipative (nonperturbative) contributions at any order in the gradient expansion., Comment: 3 pages, 1 figure, proceedings from XXVIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2018)

We study the small frequency behavior of the bulk viscosity spectral function using stochastic fluid dynamics. We obtain a number of model independent results, including the long-time tail of the bulk stress correlation function, and the leading non-analyticity of the spectral function at small frequency. We also establish a lower bound on the bulk viscosity which is weakly dependent on assumptions regarding the range of applicability of fluid dynamics. The bound on the bulk viscosity $\zeta$ scales as $\zeta_{\it min} \sim (P-\frac{2}{3}{\cal E})^2 \sum_i D_i^{-2}$, where $D_i$ are the diffusion constants for energy and momentum, and $P-\frac{2}{3}{\cal E}$, where $P$ is the pressure and ${\cal E}$ is the energy density, is a measure of scale breaking. Applied to the cold Fermi gas near unitarity, $|\lambda/a_s|\geq 1$ where $\lambda$ is the thermal de Broglie wave length and $a_s$ is the $s$-wave scattering length, this bound implies that the ratio of bulk viscosity to entropy density satisfies $\zeta/s \geq 0.1\hbar/k_B$. Here, $\hbar$ is Planck's constant and $k_B$ is Boltzmann's constant., Comment: 23 pages, accepted for publication in Physical Review A

If and how an isolated quantum system thermalizes despite its unitary time evolution is a long-standing, open problem of many-body physics. The eigenstate thermalization hypothesis (ETH) postulates that thermalization happens at the level of individual eigenstates of a system's Hamiltonian. However, the ETH requires stringent conditions to be validated, and it does not address how the thermal state is reached dynamically from an inital non-equilibrium state. We consider a Bose-Einstein condensate (BEC) trapped in a double-well potential with an initial population imbalance. We find that the system thermalizes although the initial conditions violate the ETH requirements. We identify three dynamical regimes. After an initial regime of undamped Josephson oscillations, the subsystem of incoherent excitations or quasiparticles (QP) becomes strongly coupled to the BEC subsystem by means of a dynamically generated, parametric resonance. When the energy stored in the QP system reaches its maximum, the number of QPs becomes effectively constant, and the system enters a quasi-hydrodynamic regime where the two subsystems are weakly coupled. In this final regime the BEC acts as a grand-canonical heat reservoir for the QP system (and vice versa), resulting in thermalization. We term this mechanism dynamical bath generation (DBG)., Comment: Published version, 20 pages, 19 figures, minor grammar and typos corrected

A Bose gas in a double-well potential, exhibiting a true Bose-Einstein condensate (BEC) amplitude and initially performing Josephson oscillations, is a prototype of an isolated, non-equilibrium many-body system. We investigate the quasiparticle (QP) creation and thermalization dynamics of this system by solving the time-dependent Keldysh-Bogoliubov equations. We find avalanche-like QP creation due to a parametric resonance between BEC and QP oscillations, followed by slow, exponential relaxation to a thermal state at an elevated temperature, controlled by the initial excitation energy of the oscillating BEC above its ground state. The crossover between the two regimes occurs because of an effective decoupling of the QP and BEC oscillations. This dynamics is analogous to elementary particle creation in models of the early universe. The thermalization in our set-up occurs because the BEC acts as a grand canonical reservoir for the quasiparticle system., Comment: 5 pages, 4 Figures, final version accepted for publication in Phys. Rev. Lett

A new transient effective theory of the relativistic Boltzmann equation is derived for locally momentum-anisotropic systems. In the expansion of the distribution function around a local "quasi-equilibrium" state a non-hydrodynamic dynamical degree of freedom is introduced at leading order that breaks local momentum isotropy. By replacing the deviation of the distribution function from this quasi-equilibrium state in terms of moments of the leading-order distribution and applying a systematic power counting scheme that orders the non-hydrodynamic modes by their microscopic time scales, a closed set of equations for the dynamical degrees of freedom is obtained. Truncating this set at the level of the slowest non-hydroynamic mode we find that it exhibits transient oscillatory behavior -- a phenomenon previously found only in strongly coupled theories, where it appears to be generic. In weakly coupled systems described by the Boltzmann equation, these transient oscillations depend on the breaking of local momentum isotropy being treated non-perturbatively at leading order in the expansion of the distribution function., Comment: 5 pages, 1 figure; minor changes, published version

Baryon number density perturbations offer a possible route to experimentally measure baryon number susceptibilities and heat conductivity of the quark gluon plasma. We study the fluid dynamical evolution of local and event-by-event fluctuations of baryon number density, flow velocity and energy density on top of a (generalized) Bjorken expansion. To that end we use a background-fluctuation splitting and a Bessel-Fourier decomposition for the fluctuating part of the fluid dynamical fields with respect to the azimuthal angle, the radius in the transverse plane and rapidity. We examine how the time evolution of linear perturbations depends on the equation of state as well as on shear viscosity, bulk viscosity and heat conductivity for modes with different azimuthal, radial and rapidity wave numbers. Finally we discuss how this information is accessible to experiments in terms of the transverse and rapidity dependence of correlation functions for baryonic particles in high energy nuclear collisions., Comment: 36 pages, 11 figures; v3: change in the format of the latex file. Minor changes in the text. Typos corrected and updated references. Tweak in Fig. 1-11. Accepted for publication in Phys. Rev. C

We study the evolution of the one particle distribution function that solves exactly the relativistic Boltzmann equation within the relaxation time approximation for a conformal system undergoing simultaneously azimuthally symmetric transverse and boost-invariant longitudinal expansion. We show, for arbitrary values of the shear viscosity to entropy density ratio, that the distribution function can become negative in certain kinematic regions of the available phase space depending on the boundary conditions. For thermal equilibrium initial conditions, we determine numerically the physical boundary in phase space where the distribution function is always positive definite. The requirement of positivity of this particular exact solution restricts its domain of validity, and it imposes physical constraints on its applicability., Comment: 17 pages, 6 figures. v2: Minor corrections in the text. Updated references. Published version

We use leading-order anisotropic hydrodynamics to study an azimuthally-symmetric boost-invariant quark-gluon plasma. We impose a realistic lattice-based equation of state and perform self-consistent anisotropic freeze-out to hadronic degrees of freedom. We then compare our results for the full spatiotemporal evolution of the quark-gluon plasma and its subsequent freeze-out to results obtained using 1+1d Israel-Stewart second-order viscous hydrodynamics. We find that for small shear viscosities, 4 pi eta/s ~ 1, the two methods agree well for nucleus-nucleus collisions, however, for large shear viscosity to entropy density ratios or proton-nucleus collisions we find important corrections to the Israel-Stewart results for the final particle spectra and the total number of charged particles. Finally, we demonstrate that the total number of charged particles produced is a monotonically increasing function of 4 pi eta/s in Israel-Stewart viscous hydrodynamics whereas in anisotropic hydrodynamics it has a maximum at 4 pi eta/s ~ 10. For all 4 pi eta/s > 0, we find that for Pb-Pb collisions Israel-Stewart viscous hydrodynamics predicts more dissipative particle production than anisotropic hydrodynamics., Comment: 40 pages, 16 figures

We generalize the derivation of viscous anisotropic hydrodynamics from kinetic theory to allow for non-zero particle masses. The macroscopic theory is obtained by taking moments of the Boltzmann equation after expanding the distribution function around a spheroidally deformed local momentum distribution whose form has been generalized by the addition of a scalar field that accounts non-perturbatively (i.e. already at leading order) for bulk viscous effects. Hydrodynamic equations for the parameters of the leading-order distribution function and for the residual (next-to-leading order) dissipative flows are obtained from the three lowest moments of the Boltzmann equation. The approach is tested for a system undergoing (0+1)-dimensional boost-invariant expansion for which the exact solution of the Boltzmann equation in relaxation time approximation is known. Nonconformal viscous anisotropic hydrodynamics is shown to approximate this exact solution more accurately than any other known hydrodynamic approximation., Comment: 23 pages, 3 figures. arXiv admin note: text overlap with arXiv:1311.6720

We present new analytic solutions to the relativistic Boltzmann equation within the relaxation time approximation. We first obtain spherically expanding solutions which are the kinetic counterparts of the exact solutions of the Israel-Stewart equation in the literature. This allows us to compare the solutions of the kinetic and hydrodynamic equations at an analytical level. We then derive a novel boost-invariant solution of the Boltzmann equation which has an unconventional dependence on the proper time. The existence of such a solution is also suggested in second order hydrodynamics and fluid-gravity correspondence., Comment: 18 pages.v3: update to match the published version