Search

In this review article, we discuss the current status and future prospects of perturbation theory as a means of studying the equilibrium thermodynamic and near-equilibrium transport properties of deconfined QCD matter. We begin with a brief introduction to the general topic, after which we review in some detail the foundations and modern techniques of the real- and imaginary-time formalisms of thermal field theory, covering e.g. the different bases used in the real-time formalism and the resummations required to deal with soft and collinear contributions. After this, we discuss the current status of applications of these techniques, including topics such as electromagnetic rates, transport coefficients, jet quenching, heavy quarks and quarkonia, and the Equations of State of hot quark-gluon plasma as well as cold and dense quark matter. Finally, we conclude with our view of the future directions of the field, i.e. how we anticipate perturbative calculations to contribute to our collective understanding of strongly interacting matter in the coming years., Comment: 137 pages, 44 figures. v4: Version published in Physics Reports; two typographical errors corrected from v3

We study the heavy-quark momentum diffusion coefficient in far from equilibrium gluon plasma in a self-similar regime using real-time lattice techniques. We use 3 methods for the extraction: an unequal time electric field 2-point correlator integrated over the time difference, a spectral reconstruction (SR) method based on the measured equal time electric field correlator and a kinetic theory (KT) formula. The time-evolution of the momentum diffusion coefficient extracted using all methods is consistent with an approximate $t^{\frac{-1}{2}}$ power law. We also study the extracted diffusion coefficient as a function of the upper limit of the time integration and observe that including the infrared enhancement of the equal-time correlation function in the SR calculation improves the agreement with the data for transient time behavior considerably. This is a gauge invariant confirmation of the infrared enhancement previously observed only in gauge fixed correlation functions., Comment: 5 pages, 4 figures, Conference Proceedings Quark Matter 2019

We performed state-of-the-art QCD effective kinetic theory simulations of chemically equilibrating QGP in longitudinally expanding systems. We find that chemical equilibration takes place after hydrodynamization, but well before local thermalization. By relating the transport properties of QGP and the system size we estimate that hadronic collisions with final state multiplicities $dN_\text{ch}/d\eta > 10^2$ live long enough to reach approximate chemical equilibrium for all collision systems. Therefore we expect the saturation of strangeness enhancement to occur at the same multiplicity in proton-proton, proton-nucleus and nucleus-nucleus collisions., Comment: 5 pages, 1 figure, proceedings for the XVIII International Conference on Strangeness in Quark Matter (SQM 2019), based on arXiv:1811.03040, arXiv:1811.03068

Whether, how, and to what extent solutions of Bjorken-expanding systems become insensitive to aspects of their initial conditions is of importance for heavy-ion collisions. Here we study 1+1D and phenomenologically relevant boost-invariant 3+1D systems in which initial conditions approach a universal attractor solution. In Israel-Stewart theory (IS) and kinetic theory where the universal attractor extends to arbitrarily early times, we show that all initial conditions approach the attractor at early times by a power-law while their approach is exponential at late times. In these theories, the physical mechanisms of hydrodynamization operational at late times do not drive the approach to the attractor at early times, and the early-time attractor is reached prior to hydrodynamization. In marked contrast, the attractor in strongly coupled systems is realized concurrent with hydrodynamization. This qualitative difference may offer a basis for discriminating weakly and strongly coupled scenarios of heavy-ion collisions., Comment: 12 pdf figures, revised version includes analysis of AdS/CFT attractor

Motivated by the boost-invariant Glasma state in the initial stages in heavy-ion collisions, we perform classical-statistical simulations of SU(2) gauge theory in 2+1 dimensional space-time both with and without a scalar field in the adjoint representation. We show that irrespective of the details of the initial condition, the far-from-equilibrium evolution of these highly occupied systems approaches a unique universal attractor at high momenta that is the same for the gauge and scalar sectors. We extract the scaling exponents and the form of the distribution function close to this non-thermal fixed point. We find that the dynamics are governed by an energy cascade to higher momenta with scaling exponents $\alpha = 3\beta$ and $\beta = -1/5$. We argue that these values can be obtained from parametric estimates within kinetic theory indicating the dominance of small momentum transfer in the scattering processes. We also extract the Debye mass non-perturbatively from a longitudinally polarized correlator and observe an IR enhancement of the scalar correlation function for low momenta below the Debye mass., Comment: 12 pages, 7 figures; v2: accepted by PRD, references added, minor changes

We study dynamical friction in interacting relativistic systems with arbitrary mean free paths and medium constituent masses. Our novel framework recovers the known limits of ideal gas and ideal fluid when the mean free path goes to infinity or zero, respectively, and allows for a smooth interpolation between these limits. We find that in an infinite system the drag force can be expressed as a sum of ideal-gas-like and ideal-fluid-like contributions leading to a finite friction even at subsonic velocities. This simple picture receives corrections in any finite system and the corrections become especially significant for a projectile moving at a velocity $v$ close to the speed of sound $v\approx c_s$. These corrections smoothen the ideal fluid discontinuity around the speed of sound and render the drag force a continuous function of velocity. We show that these corrections can be computed to a good approximation within effective theory of viscous fluid dynamics., Comment: 33 pages, 8 figures

To study the microscopic structure of quark-gluon plasma, data from hadronic collisions must be confronted with models that go beyond fluid dynamics. Here, we study a simple kinetic theory model that encompasses fluid dynamics but contains also particle-like excitations in a boost invariant setting with no symmetries in the transverse plane and with large initial momentum asymmetries. We determine the relative weight of fluid dynamical and particle like excitations as a function of system size and energy density by comparing kinetic transport to results from the 0th, 1st and 2nd order gradient expansion of viscous fluid dynamics. We then confront this kinetic theory with data on azimuthal flow coefficients over a wide centrality range in PbPb collisions at the LHC, in AuAu collisions at RHIC, and in pPb collisions at the LHC. Evidence is presented that non-hydrodynamic excitations make the dominant contribution to collective flow signals in pPb collisions at the LHC and contribute significantly to flow in peripheral nucleus-nucleus collisions, while fluid-like excitations dominate collectivity in central nucleus-nucleus collisions at collider energies., Comment: 28 pages, 16 figures

In this conference-proceedings contribution, we review recent advances in placing model-independent constraints on the properties of cold and dense QCD matter inside neutron stars. In addition to introducing new bounds for the Equation of State, we explain how these results may be used to make robust statements about the physical phase of strongly interacting matter in the centers of neutron stars of different masses. Our findings indicate that the existence of quark-matter cores inside massive neutron stars appears to be a very common feature of the allowed Equations of State, and should not be considered an exotic or unlikely scenario., Comment: 6 pages, 5 figures. Invited plenary talk at the 12th INTEGRAL Conference, Geneva, Switzerland, 11-15 February 2019

We study the spectral properties of an overoccupied gluonic system far from equilibrium. Using classical Yang-Mills simulations and linear response theory, we determine the statistical and spectral functions. We measure dispersion relations and damping rates of transversally and longitudinally polarized excitations in the gluonic plasma, and also study further structures in the spectral function., Comment: 12 pages, 6 figures, presented by T.L. at the XXV Epiphany conference on Advances in Heavy Ion Physics, Cracow, Poland, 8-11 January 2019

The theory governing the strong nuclear force, Quantum Chromodynamics, predicts that at sufficiently high energy densities hadronic nuclear matter undergoes a deconfinement transition to a new phase of quarks and gluons. Although this has been observed in ultrarelativistic heavy-ion collisions, it is currently an open question whether quark matter exists inside neutron stars. By combining astrophysical observations and theoretical ab-initio calculations in a model-independent way, we find that the inferred properties of matter in the cores of neutron stars with mass corresponding to 1.4 solar masses are compatible with nuclear model calculations. However, the matter in the interior of maximally massive, stable neutron stars exhibits characteristics of the deconfined phase, which we interpret as evidence for the presence of quark-matter cores. For the heaviest reliably observed neutron stars with masses of about two solar masses, the presence of quark matter is found to be linked to the behaviour of the speed of sound c_s in strongly interacting matter. If the conformal bound (c_s)^2 < 1/3 is not strongly violated, massive neutron stars are predicted to have sizable quark-matter cores. This finding has important implications for the phenomenology of neutron stars, and affects the dynamics of neutron star mergers with at least one sufficiently massive participant., Comment: 34 pages, including 8 pages main text + 10 pages methods; 8 figures. 3 EoS tables included as ancillary files. v2: Substantial changes from v1. Text shortened, with major changes to analysis and structure. Analysis and conclusions made more precise; title changed. Version accepted for publication in Nature Physics

The future opportunities for high-density QCD studies with ion and proton beams at the LHC are presented. Four major scientific goals are identified: the characterisation of the macroscopic long wavelength Quark-Gluon Plasma (QGP) properties with unprecedented precision, the investigation of the microscopic parton dynamics underlying QGP properties, the development of a unified picture of particle production and QCD dynamics from small (pp) to large (nucleus--nucleus) systems, the exploration of parton densities in nuclei in a broad ($x$, $Q^2$) kinematic range and the search for the possible onset of parton saturation. In order to address these scientific goals, high-luminosity Pb-Pb and p-Pb programmes are considered as priorities for Runs 3 and 4, complemented by high-multiplicity studies in pp collisions and a short run with oxygen ions. High-luminosity runs with intermediate-mass nuclei, for example Ar or Kr, are considered as an appealing case for extending the heavy-ion programme at the LHC beyond Run 4. The potential of the High-Energy LHC to probe QCD matter with newly-available observables, at twice larger center-of-mass energies than the LHC, is investigated., Comment: 209 pages; Report from Working Group 5 of the Workshop on the Physics of the CERN HL-LHC, and Perspectives at the HE-LHC (v2: Minor updates to text and figures. Added one author.)

We present a new method to obtain spectral properties of a non-Abelian gauge theory in the region where occupation numbers are high. The method to measure the (single-particle) spectral function is based on linear response theory and classical-statistical lattice simulations. Although we apply it to a system far from equilibrium in a self-similar regime, the extracted spectral function can be understood within the hard thermal loop (HTL) formalism and can thus be connected to thermal equilibrium at high temperatures. This allows us to obtain quantities like the lifetime of quasiparticles that are beyond the leading order and difficult to compute within HTL. The approach has the potential to measure transport coefficients, to study the earliest stages of heavy-ion collisions in a controlled way and it can be employed beyond the range of validity of HTL., Comment: 5 pages, 2 figures, Proceedings of the International Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions - Hard Probes 2018

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

We study chemical equilibration in out-of-equilibrium Quark-Gluon Plasma using the first principles method of QCD effective kinetic theory, accurate at weak coupling. In longitudinally expanding systems--relevant for relativistic nuclear collisions--we find that for realistic couplings chemical equilibration takes place after hydrodynamization, but well before local thermalization. We estimate that hadronic collisions with final state multiplicities ${dN_\text{ch}}/{d\eta}\gtrsim 10^2$ live long enough to reach approximate chemical equilibrium, which is consistent with the saturation of strangeness enhancement observed in proton-proton, proton-nucleus and nucleus-nucleus collisions., Comment: 7 pages, 3 figures, see also our companion paper arXiv:1811.03068, v2 small changes, published version

We study thermalization, hydrodynamization, and chemical equilibration in out-of-equilibrium Quark-Gluon Plasma starting from various initial conditions using QCD effective kinetic theory, valid at weak coupling. In non-expanding systems gauge bosons rapidly lose information of the initial state and achieve kinetic equilibrium among themselves, while fermions approach the equilibrium distribution only at later time. In systems undergoing rapid longitudinal expansion, both gluons and quarks are kept away from equilibrium by the expansion, but the evolution is well described by fluid dynamics even before local thermal equilibrium is reached. For realistic couplings we determine the ordering between the separate hydrodynamization, chemical equilibration and thermalization time scales to be $\tau_\text{hydro}<\tau_\text{chem}<\tau_\text{therm}$., Comment: 11 pages, 9 figures, see also our companion paper arXiv:1811.03040 ; v2 minor changes, published version

At high baryon chemical potential $\mu_B$, the equation of state of QCD allows a weak-coupling expansion in the QCD coupling $\alpha_s$. The result is currently known up to and including the full next-to-next-to-leading order (NNLO) $\alpha_s^2$. Starting at this order, the computations are complicated by the modification of particle propagation in a dense medium, which necessitates non-perturbative treatment of the scale $\alpha_s^{1/2} \mu_B$. In this work, we apply a Hard-Thermal-Loop scheme for capturing the contributions of this scale to the weak-coupling expansion, and use it to determine the leading-logarithm contribution to NNNLO: $\alpha_s^3 \ln^2 \alpha_s$. This result is the first improvement to the equation of state of massless cold quark matter in 40 years. The new term is negligibly small, and thus significantly increases our confidence in the applicability of the weak-coupling expansion., Comment: 8 pages with appendix, 2 figures; v3: Added footnote to point out minus sign error in final result

Motivated by a semi-holographic approach to the dynamics of quark-gluon plasma which combines holographic and perturbative descriptions of a strongly coupled infrared and a more weakly coupled ultraviolet sector, we construct a hybrid two-fluid model where interactions between its two sectors are encoded by their effective metric backgrounds, which are determined mutually by their energy-momentum tensors. We derive the most general consistent ultralocal interactions such that the full system has a total conserved energy-momentum tensor in flat Minkowski space and study its consequences in and near thermal equilibrium by working out its phase structure and its hydrodynamic modes., Comment: 1+54 pages, 15 figures; v3: minor additions, synced with JHEP version

The observation of large azimuthal anisotropies $v_n$ in the particle spectra of proton-proton (pp) and proton-nucleus (pA) collisions challenges fluid dynamic interpretations of $v_n$, as it remains unclear how small collision systems can hydrodynamize and to what extent hydrodynamization is needed to build up $v_n$. Here, we study in a simple kinetic theory how the same physics that leads to hydrodynamization in large systems represents itself in small systems. We observe that one third to one half of the elliptic flow signal seen in fully hydrodynamized systems can be built up in collisions that extend over only one mean free path $l_{\rm mfp}$ and that do not hydrodynamize. This is qualitatively in line with observing a sizeable $v_2$ in $pp$ collisions for which other characteristics of soft multi-particle production seem well-described in a free-streaming picture. We further expose a significant system size dependence in the accuracy of hybrid approaches that match kinetic theory to viscous fluid dynamics. The implications of these findings for a reliable extraction of shear viscosity are discussed., Comment: 5 pages, 5 figures, version accepted for publication in EPJC

High-energy nuclear collisions produce a nonequilibrium plasma of quarks and gluons which thermalizes and exhibits hydrodynamic flow. There are currently no practical frameworks to connect the early particle production in classical field simulations to the subsequent hydrodynamic evolution. We build such a framework using nonequilibrium Green's functions, calculated in QCD kinetic theory, to propagate the initial energy-momentum tensor to the hydrodynamic phase. We demonstrate that this approach can be easily incorporated into existing hydrodynamic simulations, leading to stronger constraints on the energy density at early times and the transport properties of the QCD medium. Based on (conformal) scaling properties of the Green's functions, we further obtain pragmatic bounds for the applicability of hydrodynamics in nuclear collisions., Comment: 7 pages, 4 figures, v2: typos corrected and minor changes, version accepted for publication in Phys. Rev. Lett., see also our companion paper arXiv:1805.00961 for the extensive details, for the code of linear kinetic theory propagator KoMPoST used for this study see https://github.com/KMPST/KoMPoST ; v3 updated references, published version

We develop a macroscopic description of the space-time evolution of the energy-momentum tensor during the pre-equilibrium stage of a high-energy heavy-ion collision. Based on a weak coupling effective kinetic description of the microscopic equilibration process (\`a la "bottom-up"), we calculate the non-equilibrium evolution of the local background energy-momentum tensor as well as the non-equilibrium linear response to transverse energy and momentum perturbations for realistic boost-invariant initial conditions for heavy ion collisions. We demonstrate how this framework can be used on an event-by-event basis to propagate the energy momentum tensor from far-from-equilibrium initial state models, e.g. IP-Glasma, to the time $\tau_\text{hydro}$ when the system is well described by relativistic viscous hydrodynamics. The subsequent hydrodynamic evolution becomes essentially independent of the hydrodynamic initialization time $\tau_\text{hydro}$ as long as $\tau_\text{hydro}$ is chosen in an appropriate range where both kinetic and hydrodynamic descriptions overlap. We find that for $\sqrt{s_{NN}}=2.76\,\text{TeV}$ central Pb-Pb collisions, the typical time scale when viscous hydrodynamics with shear viscosity over entropy ratio $\eta/s=0.16$ becomes applicable is $\tau_\text{hydro}\sim 1\,\text{fm/c}$ after the collision., Comment: 45 pages, 27 figures, for the code of linear kinetic theory propagator KoMPoST used for this study see https://github.com/KMPST/KoMPoST

We study the spectral properties of a highly occupied non-Abelian non-equilibrium plasma appearing ubiquitously in weak coupling descriptions of QCD matter. The spectral function of this far-from-equilibrium plasma is measured by employing linear response theory in classical-statistical real-time lattice Yang-Mills simulations. We establish the existence of transversely and longitudinally polarized quasiparticles and obtain their dispersion relations, effective mass, plasmon frequency, damping rate and further structures in the spectral and statistical functions. Our new method can be interpreted as a non-perturbative generalization of hard thermal loop (HTL) effective theory. We see indications that our results approach leading order HTL in the appropriate limit. The method can also be employed beyond the range of validity of HTL., Comment: 20 pages, 15 figures; v2: published version - references and footnotes added, minor revisions in the Abstract and Secs. III.E and IV

Non-linearities in the harmonic spectra of hadron-nucleus and nucleus-nucleus collisions provide evidence for the dynamical response to azimuthal spatial eccentricities. Here, we demonstrate within the framework of transport theory that even the mildest interaction correction to a picture of free-streaming particle distributions, namely the inclusion of one perturbatively weak interaction ("one-hit dynamics"), will generically give rise to all observed linear and non-linear structures. We further argue that transport theory naturally accounts within the range of its validity for realistic signal sizes of the linear and non-linear response coefficients observed in azimuthal momentum anisotropies with a large mean free path of the order of the system size in peripheral ($\sim 50 \%$ centrality) PbPb or central pPb collisions. The shear viscosity to entropy density ratio $\eta/s$ of such a transport theory is approximately an order of magnitude larger than that of an almost perfect fluid. The phenomenological success of transport simulations thus challenges the perfect fluid paradigm of ultra-relativistic nucleus-nucleus and hadron-nucleus collisions., Comment: 8 pages, 5 figures, minor changes, typos corrected, published version

How physical systems approach hydrodynamic behavior is governed by the decay of nonhydrodynamic modes. Here, we start from a relativistic kinetic theory that encodes relaxation mechanisms governed by different timescales thus sharing essential features of generic weakly coupled nonequilib- rium systems. By analytically solving for the retarded correlation functions, we clarify how branch cuts arise generically from noncollective particle excitations, how they interface with poles arising from collective hydrodynamic excitations, and to what extent the appearance of poles remains at best an ambiguous signature for the onset of fluid dynamic behavior. We observe that processes that are slower than the hydrodynamic relaxation timescale can make a system that has already reached fluid dynamic behavior to fall out of hydrodynamics at late times. In addition, the analytical control over this model allows us to explicitly demonstrate how the hydrodynamic gradient expansion of the correlation functions can be Borel resummed such that the full nonperturbative information is recovered using perturbative input only., Comment: 21 pages, 15 figs

The LIGO/Virgo detection of gravitational waves originating from a neutron-star merger, GW170817, has recently provided new stringent limits on the tidal deformabilities of the stars involved in the collision. Combining this measurement with the existence of two-solar-mass stars, we generate a generic family of neutron-star-matter Equations of State (EoSs) that interpolate between state-of-the-art theoretical results at low and high baryon density. Comparing the results to ones obtained without the tidal-deformability constraint, we witness a dramatic reduction in the family of allowed EoSs. Based on our analysis, we conclude that the maximal radius of a 1.4-solar-mass neutron star is 13.6 km, and that smallest allowed tidal deformability of a similar-mass star is $\Lambda(1.4 M_\odot) = 120$., Comment: 6 pages, 4 figures; v2: results and discussion updated, version accepted to PRL

Classical real-time lattice simulations play an important role in understanding non-equilibrium phenomena in gauge theories and are used in particular to model the prethermal evolution of heavy-ion collisions. Above the Debye scale the classical Yang-Mills (CYM) theory can be matched smoothly to kinetic theory. First we study the limits of the quasiparticle picture of the CYM fields by determining the plasmon mass of the system using 3 different methods. Then we argue that one needs a numerical calculation of a system of classical gauge fields and small linearized fluctuations which correspond to quantum fluctuations, in a way that keeps the separation between the two manifest. We demonstrate and test an implementation of an algorithm with the linearized fluctuations showing that the linearization indeed works and the Gauss's law is conserved., Comment: 8 pages, 3 figures

We use effective kinetic theory to study the pre-equilibrium dynamics in heavy-ion collisions. We describe the evolution of linearized energy perturbations on top of out-of-equilibrium background to the energy-momentum tensor at a time when hydrodynamics becomes applicable. We apply this description to IP-Glasma initial conditions and find an overall smooth transition to hydrodynamics. In a phenomenologically favorable range of $\eta/s$ values, early time dynamics can be accurately described in terms of a few functions of a scaled time variable $\tau T/(\eta/s)$. Our framework can be readily applied to other initial state models to provide the pre-equilibrium dynamics of the energy momentum tensor., Comment: 5 pages, 5 figures, Proceedings for the 26th International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2017), February 5-11 2017, Chicago, Illinois, USA