Search

We propose a new method to investigate the existence and location of the conjectured high-temperature critical point of strongly interacting matter via contours of constant entropy density. By approximating these lines as a power series in the baryon chemical potential $\mu_B$, one can extrapolate them from first-principle results at zero net-baryon density, and use them to locate the QCD critical point, including the associated first-order and spinodal lines. As a proof of principle, we employ currently available continuum-extrapolated first-principle results from the Wuppertal--Budapest collaboration to find a critical point at a temperature and a baryon chemical potential of $T_c = 114.3 \pm 6.9$ MeV and $\mu_{B,c} = 602.1 \pm 62.1$ MeV, respectively. We advocate for a more precise determination of the required expansion coefficients via lattice QCD simulations as a means of pinpointing the location of the critical endpoint in the phase diagram of strongly interacting matter., Comment: 7 pages, 3 figures + supplemental material

There have long been questions about the limits to the validity of relativistic fluid dynamics, and whether it is being used outside its regime of validity in modern simulations of relativistic heavy-ion collisions. An important new tool for answering this question is a causality analysis in the nonlinear regime -- if the solutions of the evolution equations do not respect relativistic causality, they are not a faithful representation of the underlying relativistic theory (in this case, quantum chromodynamics). Using this non-linear criterion, it has recently been shown that hydrodynamics is indeed being used outside its regime of validity in simulations, at least sometimes. Here we explore the phenomenological implications, particularly the quantitative effects of demanding limits on acausality in modern Bayesian parameter estimation. We find that, while typically only a small fraction of the system's energy is initially in an acausal regime, placing strict limits on the allowed energy fraction significantly changes the preferred properties of the initial condition, which in turn alters the extracted medium properties such as bulk viscosity, where large values are no longer favored. These findings highlight the importance of developing better theoretical descriptions of the early-time, out-of-equilibrium dynamics of relativistic heavy-ion collisions., Comment: 13 pages, 9 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

The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. We confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this new framework, it is possible to separate the effects of early-time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as $p$-$p$, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable., Comment: 10 pages, 5 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 point out that there is an upper bound on the speed of sound squared given by $c_s^2 \leq 0.781$ valid for all known systems described by relativistic transient hydrodynamics where calculations of certain ratios of hydrodynamic transport coefficients can be performed from first principles. Assuming this bound is valid for ultradense matter implies that the maximum mass of isolated (non-rotating) neutron stars cannot be larger than 2.7 solar masses., Comment: 5 pages, 4 figures + appendix + references (11 pages, 6 figures)

When two nuclei collide close to the speed of light, a fluid state known as the quark-gluon plasma is formed. Attempts to understand the dynamics of this fluid have generated significant research into dissipative relativistic fluid dynamics. The fluctuation-dissipation theorem implies that any dissipative dynamical system will also experience thermal fluctuations; however, such fluctuations are not typically included in the modeling of the quark-gluon plasma. This work discusses a new method of determining whether a hydrodynamic framework is consistent with thermal fluctuations. We develop a new method for calculating the noise correlator of relativistic hydrodynamic systems and apply it to Israel-Stewart theory in a general hydrodynamic frame., Comment: 4 pages. Contribution to the proceedings of the 9th Workshop for Young Scientists on the Physics of Ultra-relativistic Nucleus-Nucleus Collisions (Hot Quarks 2022)

In neutron stars, flavor-changing weak interactions determine the equilibrium fraction of protons over neutrons. In binary neutron-star mergers, violent changes in density modify this equilibrium value at timescales of milliseconds, comparable to those required for weak interactions to take place. As a result, the fraction of protons evolves out of phase with the density oscillations, giving rise to irreversible processes. The corresponding shift in pressure leads to dissipative work that can be modeled as an effective bulk-viscous correction. In this work, we derive the relevant equations of motion of Israel-Stewart hydrodynamics within this context. Using a toy model, we compute the second-order transport coefficients. Finally, we comment on the use of a realistic equation of state. Our results are expected to be useful for the study of viscous effects in numerical simulations of binary mergers., Comment: 4 pages, 1 figure. Contribution to the Hot Quarks 2022 conference proceedings

We show that it is possible to define a timelike future-directed information current within relativistic first-order hydrodynamics. This constitutes the first step toward a covariantly stable and causal formulation of first-order fluctuating hydrodynamics based on thermodynamic principles. We provide several explicit examples of first-order theories with an information current, covering many physical phenomena, ranging from electric conduction to viscosity and elasticity. We use these information currents to compute the corresponding equal-time correlation functions, and we find that the physically relevant (equal-time) correlators do not depend on the choice of the hydrodynamic frame as long as the frame leads to causal and stable dynamics. In the example of chiral hydrodynamics, we find that circularly polarized shear waves have different probabilities of being excited depending on their handedness, generating net helicity in chiral fluids., Comment: 23 pages, Published on Physical Review D, see https://journals.aps.org/prd/abstract/10.1103/PhysRevD.109.085013

We reformulate the Boltzmann equation in the relaxation time approximation undergoing Bjorken flow in terms of a novel partial differential equation for the generating function of the moments of the distribution function. This is used to obtain an approximate analytic description of this system's far-from-equilibrium attractor via a series expansion at early times. This expansion possesses a finite radius of convergence and can be analytically continued to late times. We find that this procedure reproduces the known values of shear viscosity and other transport coefficients to high accuracy. We also provide a simple approximate analytic expression that describes the attractor in the entire domain of interest for studies of quark-gluon plasma dynamics.

We construct the general theory of first-order relativistic hydrodynamics for a fluid exhibiting a chiral anomaly, including all possible viscous terms allowed by symmetry. Using standard techniques, we compute the necessary and sufficient conditions for this theory to be relativistically causal in the nonlinear regime and for thermal equilibria to be linearly stable. This is the first theory of first-order chiral hydrodynamics suitable for numerical simulations., Comment: 4 pages. Contribution to the Quark Matter 2023 conference proceedings

We present results for a Bayesian analysis of the location of the QCD critical point constrained by first-principles lattice QCD results at zero baryon density. We employ a holographic Einstein-Maxwell-dilaton model of the QCD equation of state, capable of reproducing the latest lattice QCD results at zero and finite baryon chemical potential. Our analysis is carried out for two different parametrizations of this model, resulting in confidence intervals for the critical point location that overlap at one sigma. While samples of the prior distribution may not even predict a critical point, or produce critical points spread around a large region of the phase diagram, posterior samples nearly always present a critical point at chemical potentials of $\mu_{Bc} \sim 550 - 630$ MeV., Comment: Contribution to the proceedings of Quark Matter 2023. 4 pages, 1 figure

We derive necessary and sufficient conditions under which a large class of relativistic generalizations of Braginskii's magnetohydrodynamics with shear, bulk, and heat diffusion effects is causal and strongly hyperbolic in the fully nonlinear regime in curved spacetime. We find that causality severely constrains the size of nonideal effects and the onset of kinetic instabilities. Our results are crucial for assessing the regime of validity of fluid dynamical simulations of plasmas near supermassive black holes., Comment: 16 pages, 1 figure, includes Supplemental Material. Accepted for publication in Physical Review Letters

We compute first and second-order bulk-viscous transport properties due to weak-interaction processes in $npe$ matter in the neutrino transparent regime. The transport coefficients characterize the out-of-beta-equilibrium pressure corrections, which depend on the weak-interaction rates and the equation of state. We calculate these coefficients for realistic equations of state and show they are sensitive to changes in the nuclear symmetry energy $J$ and its slope $L$., Comment: 4 pages, 3 figures. Contribution to the Quark Matter 2023 conference proceedings

We present a new general formalism for introducing thermal fluctuations in relativistic hydrodynamics, which incorporates recent developments on the causality and stability of relativistic hydrodynamic theories. Our approach is based on the information current, which measures the net amount of information carried by perturbations around equilibrium in a relativistic many-body system. The resulting noise correlators are guaranteed to be observer-independent for thermodynamically stable models. We obtain an effective action within our formalism and discuss its properties., Comment: 4 pages. Contribution to the Quark Matter 2023 conference proceedings

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"

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

The $R$-dependence of jet observables provides a new tool in understanding the interplay between the jet energy-loss mechanism and medium response in heavy-ion collisions. This work applies the Monte Carlo events generator JEWEL and PYTHIA, coupled with $\rm T_{R}ENTo$ initial conditions and the state-of-the-art (2+1)D v-USPhydro, for the simulation of jet distributions and substructure observables for lead-lead collisions at LHC energy scales. We present the jet nuclear modification $R_{AA}$ and anisotropic flow coefficients $v_{n=2,3}$ varying the jet cone radius $R$, in the context of anti-$k_T$ jets, in addition to leading subjet fragmentation. The calculations indicate the impacts of the hydrodynamic evolution and weakly-coupled medium response, given by recoils, on the distributions. Results are compared to experimental data in a wide range of jet $p_T$ and collision centrality, and displayed along large jets ($R \ge 0.6$) predictions., Comment: 6 pages, 4 figures, contribution to the proceedings of the 11th International Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions (Hard Probes 2023), Aschaffenburg (Germany), March 26-31, 2023. Submitted to PoS

We investigate the weak-interaction-driven bulk-viscous transport properties of $npe$ matter in the neutrino transparent regime. Previous works assumed that the induced bulk viscosity correction to pressure, near beta equilibrium, is linear in deviations from the equilibrium charge fraction. We show that this is not always true for (some) realistic equations of state at densities between one and three times saturation density. This nonlinear nature of the perturbation around equilibrium motivates a far-from-beta-equilibrium description of bulk-viscous transport in neutron star mergers, which can be precisely achieved using a new Israel-Stewart formulation with resummed bulk and relaxation time transport coefficients. The computation of these transport coefficients depends on out-of-beta-equilibrium pressure corrections, which can be computed for a given equation of state. We calculate these coefficients for equations of state that satisfy the latest constraints from multi-messenger observations from LIGO/VIRGO and NICER. We show that varying the nuclear symmetry energy $J$ and its slope $L$ can significantly affect the transport coefficients and the nonlinear behavior of the out-of-equilibrium pressure corrections. Therefore, having better constraints on $J$ and $L$ will directly impact our understanding of bulk-viscous processes in neutron star mergers., Comment: 31 pages, 17 figures, 5 appendices, matches version accepted for publication in Physical Review C

A fundamental question in QCD is the existence of a phase transition at large doping of quarks over antiquarks. We present the first prediction of a QCD critical point (CP) from a Bayesian analysis constrained by first principle results at zero doping. We employ the gauge/gravity duality to map QCD onto a theory of dual black holes. Predictions for the CP location in different realizations of the model overlap at one sigma. Even if many prior samples do not include a CP, one is found in nearly 100\% of posterior samples, indicating a strong preference for a CP., Comment: 8 pages, 4 figures + appendix + references (15 pages, 9 figures), matches published version

Causality is necessary for retarded Green's functions to remain retarded in all inertial frames in relativity, which ensures that dissipation of fluctuations is a Lorentz invariant concept. For first-order BDNK theories with stochastic fluctuations, introduced via the Schwinger-Keldysh formalism, we show that imposing causality and stability leads to correlation functions of hydrodynamic fluctuations that only display the expected physical properties at small frequencies and wavenumber, i.e., within the expected regime of validity of the first-order approach. For second-order theories of Israel and Stewart type, constructed using the information current such that entropy production is always non-negative, a stochastic formulation is presented using the Martin-Siggia-Rose approach where imposing causality and stability leads to correlators with the desired properties. We also show how Green's functions can be determined from such an action. We identify a $\mathbb{Z}_2$ symmetry, analogous to the Kubo-Martin-Schwinger symmetry, under which this Martin-Siggia-Rose action is invariant. This modified Kubo-Martin-Schwinger symmetry provides a new guide for the effective action formulation of hydrodynamic systems with dynamics not solely governed by conservation laws. Furthermore, this symmetry ensures that the principle of detailed balance is valid in a covariant manner. We employ the new symmetry to further clarify the connection between the Schwinger-Keldysh and Martin-Siggia-Rose approaches, establishing a precise link between these descriptions in second-order theories of relativistic hydrodynamics. Finally, the modified Kubo-Martin-Schwinger symmetry is used to determine the corresponding action describing diffusion in Israel-Stewart theories in a general hydrodynamic frame., Comment: 28 pages

We derive the set of inequalities that is necessary and sufficient for nonlinear causality and linear stability of first-order relativistic hydrodynamics with either a $U(1)_V$ conserved current or a $U(1)_A$ current with a chiral anomaly or both. Our results apply to generic hydrodynamic frames in which no relations among the transport parameters are imposed. Furthermore, our analysis yields, to the best of our knowledge, the first theory of viscous chiral hydrodynamics proven to be causal and stable. We find that causality demands the absence of vorticity-induced heat flux, forcing a departure from the thermodynamic frame in the chiral case. The inequalities for causality and stability define a hypervolume in the space of transport parameters, wherein each point corresponds to a consistent formulation. Notably, causality is determined by just three combinations of transport parameters. We present our results in a form amenable to numerical hydrodynamic simulations., Comment: 21 pages, 1 figure. Version accepted for publication in Physical Review D

We investigate the initial value problem of a very general class of $3+1$ non-Newtonian compressible fluids in which the viscous stress tensor with shear and bulk viscosity relaxes to its Navier-Stokes values. These fluids correspond to the non-relativistic limit of well-known Israel-Stewart-like theories used in the relativistic fluid dynamic simulations of high-energy nuclear and astrophysical systems. After establishing the local well-posedness of the Cauchy problem, we show for the first time in the literature that there exists a large class of initial data for which the corresponding evolution breaks down in finite time due to the formation of singularities. This implies that a large class of non-Newtonian fluids do not have finite solutions defined at all times., Comment: 15 pages, revised version accepted for publication

We show that linear superpositions of plane waves involving a single-valued, covariantly stable dispersion relation $\omega(k)$ always propagate outside the lightcone, unless $\omega(k) =a+b k$. This implies that there is no notion of causality for individual dispersion relations, since no mathematical condition on the function $\omega(k)$ (such as the front velocity or the asymptotic group velocity conditions) can serve as a sufficient condition for subluminal propagation in dispersive media. Instead, causality can only emerge from a careful cancellation that occurs when one superimposes all the excitation branches of a physical model. This is shown to happen automatically in local theories of matter that are covariantly stable. Hence, we find that the need for nonhydrodynamic modes in relativistic fluid mechanics is analogous to the need for antiparticles in relativistic quantum mechanics., Comment: 6 pages (Main Text) + 2 pages (Supplementary Material), published on Physical Review Letters, see https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.132.162301

In this review, we provide an up-to-date account of quantitative holographic descriptions of the strongly coupled quark-gluon plasma (QGP) produced in heavy-ion collisions, based on the class of gauge-gravity Einstein-Maxwell-Dilaton (EMD) models. Holography is employed to tentatively map the QCD phase diagram at finite temperature onto a dual theory of charged, asymptotically AdS black holes in 5D. With a quantitative focus on the hot QCD phase diagram, the EMD models reviewed are adjusted to describe lattice results for the finite-temperature QCD equation of state, with 2+1 flavors and physical quark masses, at zero chemical potential and vanishing electromagnetic fields. The predictive power of EMD models is tested by quantitatively comparing their predictions for the hot QCD equation of state at nonzero baryon density and the corresponding state-of-the-art lattice QCD results. The shear and bulk viscosities predicted by these EMD models are also compared to the corresponding profiles favored by the latest phenomenological multistage models describing different heavy-ion data. We report preliminary results from a Bayesian analysis which provide systematic evidence that lattice results at finite temperature and zero baryon density strongly constrains the free parameters of EMD models. Remarkably, the set of parameters constrained by lattice results at zero chemical potential produces EMD models in quantitative agreement with lattice QCD results also at finite baryon density. We also review results for equilibrium and transport properties from magnetic EMD models, describing the QGP at finite temperatures and magnetic fields. Finally, we provide a critical assessment of the main limitations and drawbacks of the holographic models reviewed in the present work, and point out some perspectives we believe are of fundamental importance for future developments., Comment: 76 pages, 14 figures, invited review accepted for publication in Progress in Particle and Nuclear Physics

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments on the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion, and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. The latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor., Comment: 45 pages, 2 figures

We develop the first causal and stable theory of a bulk-viscous relativistic pseudoplastic (or dilatant) fluid. This new formalism brings to light the rheological properties of several relativistic physical systems. Neutron star collisions can behave as a relativistic pseudoplastic material with viscous properties dictated by the non-conservation of lepton currents due to weak decay. Two-temperature relativistic plasmas, such as those surrounding supermassive galactic black holes, are predominantly pseudoplastic. Our framework can also be employed to construct novel viscous models for the evolution of the Universe with pseudoplastic or dilatant features., Comment: 9 pages (Main Text) + 5 pages (Supplementary Material); 3 figures; Published on Physical Review D, see https://journals.aps.org/prd/abstract/10.1103/PhysRevD.109.096040

This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided.

A general organizing principle is proposed that can be used to derive the equations of motion describing the near-equilibrium dynamics of causal and thermodynamically stable relativistic systems. The latter are found to display some new type of universal behavior near equilibrium that allows them to be grouped into universality classes defined by their degrees of freedom, information content, and conservation laws. The universality classes expose a number of surprising equivalences between different theories, shedding new light on the near-equilibrium behavior of relativistic systems., Comment: 8 pages (Main Text) + 9 pages (Supplementary Material), published on PRL, see https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.132.222302

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS., Comment: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

We study the dynamical formation of scalar monopole and dipole hair in scalar Gauss-Bonnet theory and dynamical Chern-Simons theory. We prove that the spherically-symmetric mode of the dipole hair is completely determined by the product of the mass of the spacetime and the value of the monopole hair. We then show that the dynamics of the $\ell=1$ mode of the dipole hair is intimately tied to the appearance of the event horizon during axisymmetric collapse, which results in the radiation of certain modes that could have been divergent in the future of the collapse. We confirm these analytical predictions by simulating the gravitational collapse of a rapidly rotating neutron star in the decoupling limit, both in scalar Gauss-Bonnet and dynamical Chern-Simons theory. Our results, combined with those of Ref.~\cite{R:2022cwe}, provide a clear physical picture of the dynamics of scalar monopole and dipole radiation in axisymmetric and spherical gravitational collapse in these theories., Comment: v2-matches published version in PRD

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

Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science, Comment: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

We investigate the out-of-equilibrium dynamics of viscous fluids in a spatially flat Friedmann-Lema\^itre-Robertson-Walker cosmology using the most general causal and stable viscous energy-momentum tensor defined at first order in spacetime derivatives. In this new framework a pressureless viscous fluid having density $\rho$ can evolve to an asymptotic future solution in which the Hubble parameter approaches a constant while $\rho \rightarrow 0$, even in the absence of a cosmological constant (i.e., $\Lambda = 0$). Thus, while viscous effects in this model drive an accelerated expansion of the universe, the density of the viscous component itself vanishes, leaving behind only the acceleration. This behavior emerges as a consequence of causality in first-order theories of relativistic fluid dynamics and it is fully consistent with Einstein's equations., Comment: 14 pages, no figures, references updated, discussion added

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

By using gravity/gauge correspondence, we construct a holographic model, constrained to mimic the lattice QCD equation of state at zero density, to investigate the temperature and baryon chemical potential dependence of the equation of state. We also obtained the energy loss of light and heavy partons within the hot and dense plasma represented by the heavy quark drag force, Langevin diffusion coefficients and jet quenching parameter at the critical point and across the first-order transition line predicted by the model., Comment: 4 pages, proceedings for the 20th International Conference on Strangeness in Quark Matter (SQM2022)

We combine, for the first time, event-by-event $\rm T_RENTo$ initial conditions with the relativistic viscous hydrodynamic model v-USPhydro and the Monte Carlo event generator JEWEL to make predictions for the nuclear modification factor $R_{AA}$ and jet azimuthal anisotropies $v_n\left\{2\right\}$ in $\sqrt{s_{NN}}=5.02 \, \rm TeV$ PbPb collisions for multiple centralities and values of the jet cone radius $R$. The $R$-dependence of $R_{AA}$ and $v_2\left\{2\right\}$ strongly depends on the presence of recoiling scattering centers. We find a small jet $v_3\left\{2\right\}$ in central collisions, which is consistent with ATLAS results., Comment: 7 pages, 4 figures, higher number of generated events, refs. added

By using the AdS/CFT correspondence, we construct an Einstein-Maxwell-Dilaton model to map the thermodynamics of strongly interacting matter. The holographic model, constrained to reproduce the lattice QCD equation of state at zero baryon chemical potential, predicts a critical end point and a first order phase transition line. We also obtain the equation of state of the model for a large region of the phase diagram. We characterize the crossover transition region by two lines, one associated to the inflection of the second order baryon susceptibility, and another one associated to the minimum of the square of the speed of sound along trajectories of constant entropy per baryon number. We observe that both lines merge at the critical point., Comment: 6 pages, 3 figures, contribution to the proceedings of The International conference on Critical Point and Onset of Deconfinement - CPOD2021

By using gravity/gauge correspondence, we employ an Einstein-Maxwell-Dilaton model to compute the equilibrium and out-of-equilibrium properties of a hot and baryon rich strongly coupled quark-gluon plasma. The family of 5-dimensional holographic black holes, which are constrained to mimic the lattice QCD equation of state at zero density, is used to investigate the temperature and baryon chemical potential dependence of the equation of state. We also obtained the baryon charge conductivity, and the bulk and shear viscosities with a particular focus on the behavior of these observables on top of the critical end point and the line of first order phase transition predicted by the model., Comment: Proceedings for the 37th Winter Workshop on Nuclear Dynamics

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

In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ($\beta$-)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in post-merger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact on the post-merger gravitational wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the post-merger gravitational wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation., Comment: 15 pages, 7 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

A bottom-up Einstein-Maxwell-Dilaton holographic model is used to compute, for the first time, the behavior of several transport coefficients of the hot and baryon-rich strongly coupled quark-gluon plasma at the critical point and also across the first-order phase transition line in the phase diagram. The observables under study are of the shear and bulk viscosities, baryon diffusion, thermal conductivity, the jet quenching parameter $\hat{q}$, as well as the heavy-quark drag force and Langevin diffusion coefficients. These calculations provide a phenomenologically promising estimate for these coefficients, given that our model quantitatively reproduces lattice QCD thermodynamics results, both at zero and finite baryon density, besides naturally incorporating the nearly-perfect fluidity of the quark-gluon plasma. We find that the diffusion of baryon charge, and also the shear and bulk viscosities, are suppressed with increasing baryon density, indicating that the medium becomes even closer to perfect fluidity at large densities. On the other hand, the jet quenching parameter and the heavy-quark momentum diffusion are enhanced with increasing density. The observables display a discontinuity gap when crossing the first-order phase transition line, while developing an infinite slope at the critical point. The transition temperatures associated with different transport coefficients differ in the crossover region but are found to converge at the critical point., Comment: 18 pages, 13 figures

Black holes in certain modified gravity theories that contain a scalar field coupled to curvature invariants are known to possess (monopole) scalar hair while non-black-hole spacetimes (like neutron stars) do not. Therefore, as a neutron star collapses to a black hole, scalar hair must grow until it settles to the stationary black hole solution with (monopole) hair. In this paper, we study this process in detail and show that the growth of scalar hair is tied to the appearance and growth of the event horizon (before an apparent horizon forms), which forces scalar modes that would otherwise (in the future) become divergent to be radiated away. We prove this result rigorously in general first for a large class of modified theories, and then we exemplify the results by studying the temporal evolution of the scalar field in scalar Gauss-Bonnet gravity in two backgrounds: (i) a collapsing Oppenheimer-Snyder background, and (ii) a collapsing neutron star background. In case (i), we find an exact scalar field solution analytically, while in case (ii) we solve for the temporal evolution of the scalar field numerically, with both cases supporting the conclusion presented above. Our results suggest that the emission of a burst of scalar field radiation is a necessary condition for black hole formation in a large class of modified theories of gravity., Comment: v2 : 23 pages, minor corrections in text, matches published version

We compute the radius of convergence of the linearized relativistic hydrodynamic expansion around a non-trivially rotating strongly coupled N=4 Super-Yang-Mills plasma. Our results show that the validity of hydrodynamics is sustained and can even get enhanced in a highly vortical quark-gluon plasma, such as the one produced in heavy-ion collisions. The hydrodynamic dispersion relations are computed using a rotating background that is an analytic solution of the ideal hydrodynamic equations of motion with non-vanishing angular momentum and large vorticity gradients, giving rise to a particular boost symmetry. Analytic equations for the transport coeffcients of the rotating plasma as a function of their values in a plasma at rest are given., Comment: 8 pages, 2 figures, comments welcome!