Source author record

Paul Romatschke

Paul Romatschke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2020arXiv

Quantum Mechanical Out-Of-Time-Ordered-Correlators for the Anharmonic (Quartic) Oscillator

Out-of-time-ordered correlators (OTOCs) have been suggested as a means to study quantum chaotic behavior in various systems. In this work, I calculate OTOCs for the quantum mechanical anharmonic oscillator with quartic potential, which is classically integrable and has a Poisson-like energy-level distribution. For low temperature, OTOCs are periodic in time, similar to results for the harmonic oscillator and the particle in a box. For high temperature, OTOCs exhibit a rapid (but power-like) rise at early times, followed by saturation consistent with $2\langle x^2\rangle_T \langle p^2\rangle_T$ at late times. At high temperature, the spectral form factor decreases at early times, bounces back and then reaches a plateau with strong fluctuations.

preprint2020arXiv

Real Time Quantum Gravity Dynamics from Classical Statistical Yang-Mills Simulations

We perform microcanonical classical statistical lattice simulations of SU(N) Yang-Mills theory with eight scalars on a circle. Measuring the eigenvalue distribution of the spatial Wilson loop we find two distinct phases depending on the total energy and circle radius, which we tentatively interpret as corresponding to black hole and black string phases in a dual gravity picture. We proceed to study quenches by first preparing the system in one phase, rapidly changing the total energy, and monitoring the real-time system response. We observe that the system relaxes to the equilibrium phase corresponding to the new energy, in the process exhibiting characteristic damped oscillations. We interpret this as the topology change from black hole to black string configurations, with damped oscillations corresponding to quasi-normal mode ringing of the black hole/black string final state. This would suggest that alpha' corrections alone can resolve the singularity associated with the topology change. We extract the real and imaginary part of the lowest-lying presumptive quasinormal mode as a function of energy and N.

preprint2019arXiv

Finite temperature CFT results for all couplings: O(N) model in 2+1 dimensions

A famous example of gauge/gravity duality is the result that the entropy density of strongly coupled ${\cal N}=4$ SYM in four dimensions for large N is exactly 3/4 of the Stefan-Boltzmann limit. In this work, I revisit the massless O(N) model in 2+1 dimensions, which is analytically solvable at finite temperature $T$ for all couplings $λ$ in the large N limit. I find that the entropy density monotonically decreases from the Stefan-Boltzmann limit at $λ=0$ to exactly 4/5 of the Stefan-Boltzmann limit at $λ=\infty$. Calculating the retarded energy-momentum tensor correlator in the scalar channel at $λ=\infty$, I find that it has two logarithmic branch cuts originating at $ω=\pm 4 T \ln \frac{1+\sqrt{5}}{2}$, but no singularities in the whole complex frequency plane. I show that the ratio 4/5 and the location of the branch points both are universal within a large class of bosonic CFTs in 2+1 dimensions.

preprint2019arXiv

Simple non-perturbative resummation schemes beyond mean-field II: thermodynamics of scalar $ϕ^4$ theory in 1+1 dimensions at arbitrary coupling

Recently, non-perturbative approximate solutions were presented that go beyond the well-known mean-field resummation. In this work, these non-perturbative approximations are used to calculate finite temperature equilibrium properties for scalar $ϕ^4$ theory in two dimensions such as the pressure, entropy density and speed of sound. Unlike traditional approaches, it is found that results are well-behaved for arbitrary temperature/coupling strength, are independent of the choice of the renormalization scale $\barμ^2$, and are apparently converging as the resummation level is increased. Results also suggest the presence of a possible analytic cross-over from the high-temperature to the low-temperature regime based on the change in the thermal entropy density.

preprint2019arXiv

Strong Coupling Universality at Large N for Pure CFT Thermodynamics in 2+1 dimensions

Pure CFTs have vanishing $β$-function at any value of the coupling. One example of a pure CFT is the O(N) Wess-Zumino model in 2+1 dimensions in the large N limit. This model can be analytically solved at finite temperature for any value of the coupling, and we find that its entropy density at strong coupling is exactly equal to 31/35 of the non-interacting Stefan-Boltzmann result. We show that a large class of theories with equal numbers of N-component fermions and bosons, supersymmetric or not, for a large class of interactions, exhibit the same universal ratio. For unequal numbers of fermions and bosons we find that the strong-weak thermodynamic ratio is bounded to lie in between 4/5 and 1.

preprint2016arXiv

On high-order perturbative calculations at finite density

We discuss the prospects of performing high-order perturbative calculations in systems characterized by a vanishing temperature but finite density. In particular, we show that the determination of generic Feynman integrals containing fermionic chemical potentials can be reduced to the evaluation of three-dimensional phase space integrals over vacuum on-shell amplitudes - a result reminiscent of a previously proposed "naive real-time formalism" for vacuum diagrams. Applications of these rules are discussed in the context of the thermodynamics of cold and dense QCD, where it is argued that they facilitate an extension of the Equation of State of cold quark matter to higher perturbative orders.

preprint2016arXiv

Retarded Correlators in Kinetic Theory: Branch Cuts, Poles and Hydrodynamic Onset Transitions

In this work the collective modes of an effective kinetic theory description based on the Boltzmann equation in relaxation time approximation applicable to gauge theories at weak but finite coupling and low frequencies are studied. Real time retarded two-point correlators of the energy-momentum tensor and the R-charge current are calculated at finite temperature in flat space-times for large N gauge theories. It is found that the real time correlators possess logarithmic branch cuts which in the limit of large coupling disappear and give rise to non-hydrodynamic poles that are reminiscent of quasi-normal modes in black holes. In addition to branch cuts, correlators can have simple hydrodynamic poles, generalizing the concept of hydrodynamic modes to intermediate wavelength. Surprisingly, the hydrodynamic poles cease to exist for some critical value of the wavelength and coupling reminiscent of the properties of onset transitions.

preprint2015arXiv

Collective flow without hydrodynamics: simulation results for relativistic ion collisions

Flow signatures in experimental data from relativistic ion collisions are usually interpreted as a fingerprint of the presence of a hydrodynamic phase during the evolution of these systems. In this work, flow signatures arising from event-by-event viscous hydrodynamics are compared to those arising from event-by-event non-interacting particle dynamics (free-streaming), both followed by a late-stage hadronic cascade, in d+Au, 3He+Au at sqrt(s)=200 GeV and p+Pb collisions at sqrt(s)=5 TeV, respectively. For comparison, also Pb+Pb collisions at sqrt(s)=2.76 TeV are simulated. It is found that non-hydrodynamic evolution can give rise to equal or larger radial flow than hydrodynamics with eta/s=0.08 in all simulated collision systems. In light-on-heavy-ion collisions, free-streaming gives rise to triangular and quadrupolar flow comparable to or larger than that from hydrodynamics, but it generally leads to considerably smaller elliptic flow. As expected, free-streaming leads to considerably less elliptic, triangular and quadrupolar flow than hydrodynamics in nucleus-nucleus collisions, such as event-by-event Pb+Pb collisions at sqrt(s)=2.76 TeV.

preprint2015arXiv

Lattice Boltzmann simulations of a two-dimensional Fermi gas at unitarity

We present fully nonlinear dissipative fluid dynamics simulations of a trapped two-dimensional Fermi gas at unitarity using a Lattice Boltzmann algorithm. We are able to simulate non-harmonic trapping potentials, temperature-dependent viscosities as well as a discretized version of the ballistic (non-interacting) behavior. Our approach lends itself to direct comparison with experimental data, opening up the possibility of a precision determination of transport coefficients in the unitary Fermi gas. Furthermore, we predict the presence of a non-hydrodynamic component in the quadrupole mode, which should be observable experimentally.

preprint2015arXiv

Non-hydrodynamic transport in trapped unitary Fermi gases

Many strongly coupled fluids are known to share similar hydrodynamic transport properties. In this work we argue that this similarity could extend beyond hydrodynamics to transient dynamics through the presence of non-hydrodynamic modes. We review non-hydrodynamic modes in kinetic theory and gauge/gravity duality and discuss their signatures in trapped Fermi gases close to unitarity. Reanalyzing previously published experimental data, we find hints of non-hydrodynamic modes in cold Fermi gases in two and three dimensions.

preprint2015arXiv

Weak and strong coupling equilibration in nonabelian gauge theories

We present a direct comparison studying equilibration through kinetic theory at weak coupling and through holography at strong coupling in the same set-up. The set-up starts with a homogeneous thermal state, which then smoothly transitions through an out-of-equilibrium phase to an expanding system undergoing boost-invariant flow. This first apples-to-apples comparison of equilibration provides a benchmark for similar equilibration processes in heavy-ion collisions, where the equilibration mechanism is still under debate. We find that results at weak and strong coupling can be smoothly connected by simple, empirical power-laws for the viscosity, equilibration time and entropy production of the system.

preprint2014arXiv

Absence of a local rest frame in far from equilibrium quantum matter

In a collision of strongly coupled quantum matter we find that the dynamics of the collision produces regions where a local rest frame cannot be defined because the energy-momentum tensor does not have a real time-like eigenvector. This effect is purely quantum mechanical, since for classical systems, a local rest frame can always be defined. We study the relation with the null and weak energy condition, which are violated in even larger regions, and compare with previously known examples. While no pathologies or instabilities arise, it is interesting that regions without a rest frame are possibly present in heavy ion collisions.

preprint2014arXiv

Extracting the shear viscosity of a high temperature hadron gas

Quark-Gluon plasmas produced in relativistic heavy-ion collisions quickly expand and cool, entering a phase consisting of multiple interacting hadronic resonances just below the QCD deconfinement temperature. The transport properties of this hot hadron gas are poorly understood, yet they play an important role in our ability to infer transport properties of the quark-gluon plasma, because experimental measurements integrate over the whole system evolution. Assuming that the hot hadron gas can be modeled by a hadron cascade code based on kinetic theory assuming binary interactions, the shear viscosity over entropy ratio of a hot hadron gas for temperatures in between 120 MeV and 170 MeV is extracted. Furthermore, we present estimates for a second order transport coefficient, the shear viscous relaxation time at a temperature of 165 MeV.

preprint2014arXiv

Onset of cavitation in the quark-gluon plasma

We study the onset of bubble formation (cavitation) in the quark-gluon plasma as a result of the reduction of the effective pressure from bulk-viscous corrections. By calculating velocity gradients in typical models for quark-gluon plasma evolution in heavy-ion collisions, we obtain results for the critical bulk viscosity above which cavitation occurs. Since present experimental data for heavy-ion collisions seems inconsistent with the presence of bubbles above the phase transition temperature of QCD, our results may be interpreted as an upper limit of the bulk viscosity in nature. Our results indicate that bubble formation is consistent with the expectation of hadronisation in low-temperature QCD.

preprint2014arXiv

Precision studies of v_n fluctuations

The power spectrum of heavy ion collisions is investigated by studying initial state fluctuations on top of a smooth hydrodynamic flow. In particular, the stability of the location of the first minimum of the power spectrum and the dependence of the hydrodynamic response v_n / epsilon_n on n and on eta / s are discussed. In our study we develop a new Green's function method for the analytic hydrodynamic flow by S. Gubser and make use of a fully non-linear hydrodynamics code. We find that there will be no well-defined first minimum of the response for n < 10, due to the fact that all minima in that region are found to be sensitive to the location of the initial perturbations. Also, we find that the often proposed form of the hydrodynamical response, that ln ( v_n / epsilon_n ) depend quadratically on n and linearly on eta / s, should not hold once many events have been averaged over.

preprint2014arXiv

Simulation of Black Hole Collisions in Asymptotically AdS Spacetimes

We present results from the evolution of spacetimes that describe the merger of asymptotically global AdS black holes in 5D with an SO(3) symmetry. Prompt scalar field collapse provides us with a mechanism for producing distinct trapped regions on the initial slice, associated with black holes initially at rest. We evolve these black holes towards a merger, and follow the subsequent ring-down. The boundary stress tensor of the dual CFT is conformally related to a stress tensor in Minkowski space which inherits an axial symmetry from the bulk SO(3). We compare this boundary stress tensor to its hydrodynamic counterpart with viscous corrections of up to second order, and compare the conformally related stress tensor to ideal hydrodynamic simulations in Minkowski space, initialized at various time slices of the boundary data. Our findings reveal far-from-hydrodynamic behavior at early times, with a transition to ideal hydrodynamics at late times.

preprint2014arXiv

The equation of state in two-, three-, and four-color QCD at non-zero temperature and density

We calculate the equation of state at non-zero temperature and density from first principles in two-, three- and four-color QCD with two fermion flavors in the fundamental and two-index, antisymmetric representation. By matching low-energy results (from a `hadron resonance gas') to high-energy results from (resummed) perturbative QCD, we obtain results for the pressure and trace anomaly that are in quantitative agreement with full lattice-QCD studies for three colors at zero chemical potential. Our results for non-zero chemical potential at zero temperature constitute predictions for the equation of state in QCD-like theories that can be tested by traditional lattice studies for two-color QCD with two fundamental fermions and four-color QCD with two two-index, antisymmetric fermions. We find that the speed of sound squared at zero temperature can exceed one third, which may be relevant for the phenomenology of high-mass neutron stars.

preprint2013arXiv

A fully dynamical simulation of central nuclear collisions

We present a fully dynamical simulation of central nuclear collisions around mid-rapidity at LHC energies. Unlike previous treatments, we simulate all phases of the collision, including the equilibration of the system. For the simulation, we use numerical relativity solutions to AdS/CFT for the pre-equilibrium stage, viscous hydrodynamics for the plasma equilibrium stage and kinetic theory for the low density hadronic stage. Our pre-equilibrium stage provides initial conditions for hydrodynamics, resulting in sizable radial flow. The resulting light particle spectra reproduce the measurements from the ALICE experiment at all transverse momenta.

preprint2013arXiv

Implications of hydrodynamic fluctuations on the minimum shear viscosity of the dilute Fermi gas at unitarity

We confirm and expand on work by Chafin and Schaefer on hydrodynamic fluctuations in the unitary Fermi gas. Using the result for the equation of state from a recent MIT experiment, we derive lower bounds for η/n and η/s as a function of temperature. Re-analyzing recent quantum Monte Carlo data for the shear-viscosity spectral function we point out a possible resolution for the tension between the viscosity bound η/n> 0.3 from Chafin and Schaefer and the quantum Monte Carlo results η/n<0.2$ from Wlazlowski et al. near the critical temperature

preprint2013arXiv

Pre-Equilibrium Radial Flow from Central Shock-Wave Collisions in AdS5

Using gauge/gravity duality, central ultrarelativistic nucleus-nucleus collisions are modelled as collisions of shock waves in five-dimensional asymptotic AdS space. For early times after the collision, it is possible to analytically match the metric from the past to the future light-cone. This allows extraction of the pre-equilibrium energy-momentum tensor of the strongly coupled, large N gauge theory. For central collisions, this allows qualitative statements concerning the build-up of radial flow at mid-rapidity in AA and pA collisions. We find that the early-time radial flow buildup is identical to that expected from ideal hydrodynamics with an entropy density proportional to the square root of the product of the matter densities in the individual "nuclei".

preprint2013arXiv

Tests of the Quark-Gluon Plasma Coupling Strength at Early Times with Heavy Quarks

The redistribution in momentum space of heavy quarks via their interactions in the quark-gluon plasma is an excellent probe of the heavy quark coupling strength to the medium. We utilize a Monte Carlo Langevin calculation for tracking heavy quark - antiquark pairs embedded in a viscous hydrodynamic space-time evolution. We find that the nuclear modification factor ($R_{AA}$) for charm quarks is relatively insensitive to the coupling to the quark-gluon plasma at early times where the highest temperatures are acheived. In contrast the azimuthal angular correlation of charm and anticharm quarks is extremely sensitive to the early time evolution. For beauty quarks the situation is reversed in terms of sensitivity. Future measurements of heavy quarks have the potential to map out the temperature dependence of the shear viscosity to entropy density ratio (η/s).

preprint2012arXiv

Relativistic (Lattice) Boltzmann Equation with Non-Ideal Equation of State

The relativistic Boltzmann equation for a single particle species generally implies a fixed, unchangeable equation of state that corresponds to that of an ideal gas. Real-world systems typically have more complicated equation of state which cannot be described by the Boltzmann equation. The present work derives a 'Boltzmann-like' equation that gives rise to a conserved energy-momentum tensor with an arbitrary (but thermodynamically consistent) equation of state. Using this, a Lattice Boltzmann scheme for diagonal metric tensors and arbitrary equations of state is constructed. The scheme is verified for QCD in the Milne metric by comparing to viscous fluid dynamics.

preprint2011arXiv

Shock wave collisions in AdS5: approximate numerical solutions

We numerically study the evolution of a boost-invariant N=4 SYM medium using AdS/CFT. We consider a toy model for the collision of gravitational shock waves, finding that the energy density first increases, reaches a maximum and then starts to decrease, matching hydrodynamics for late times. For the initial conditions we consider, the hydrodynamic scale governing the late time behaviour is to very good approximation determined by the area of the black hole horizon at initial times. Our results provide a toy model for the early time evolution of the bulk system in heavy-ion collisions at RHIC and the LHC.

preprint2011arXiv

The stickiness of sound: An absolute lower limit on viscosity and the breakdown of second order relativistic hydrodynamics

Hydrodynamics predicts long-lived sound and shear waves. Thermal fluctuations in these waves can lead to the diffusion of momentum density, contributing to the shear viscosity and other transport coefficients. Within viscous hydrodynamics in 3+1 dimensions, this leads to a positive contribution to the shear viscosity, which is finite but inversely proportional to the microscopic shear viscosity. Therefore the effective infrared viscosity is bounded from below. The contribution to the second-order transport coefficient $τ_π$ is divergent, which means that second-order relativistic viscous hydrodynamics is inconsistent below some frequency scale. We estimate the importance of each effect for the Quark-Gluon Plasma, finding them to be minor if $η/s = 0.16$ but important if $η/s = 0.08$.

preprint2010arXiv

Cold Quark Matter

We perform an O(alpha_s^2) perturbative calculation of the equation of state of cold but dense QCD matter with two massless and one massive quark flavor, finding that perturbation theory converges reasonably well for quark chemical potentials above 1 GeV. Using a running coupling constant and strange quark mass, and allowing for further non-perturbative effects, our results point to a narrow range where absolutely stable strange quark matter may exist. Absent stable strange quark matter, our findings suggest that quark matter in compact star cores becomes confined to hadrons only slightly above the density of atomic nuclei. Finally, we show that equations of state including quark matter lead to hybrid star masses up to M~2M_solar, in agreement with current observations. For strange stars, we find maximal masses of M~2.75M_solar and conclude that confirmed observations of compact stars with M>2M_solar would strongly favor the existence of stable strange quark matter.

preprint2010arXiv

Looking inside neutron stars: Microscopic calculations confront observations

While QCD appears not to be accurately solvable in the regime of interest for neutron star physics, microscopic calculations are feasible at both low and very high densities. In this work, we propose using the most realistic calculations in these two regimes of nuclear physics and perturbative QCD, and construct equations of state by matching the results requiring thermodynamic consistency. We find that the resulting equations of state --- in contrast to several hadronic ones --- are able to reproduce current observational data on neutron stars without any fine tuning, and allow stable hybrid stars with masses up to 2.1M_{sun}. Using recent observations of star radii, we perform a maximum likelihood analysis to further constrain the equation of state, and in addition show that the effects of rotation on radii and masses should not be neglected in future precision studies.

preprint2010arXiv

Relativistic Viscous Fluid Dynamics and Non-Equilibrium Entropy

Fluid dynamics corresponds to the dynamics of a substance in the long wavelength limit. Writing down all terms in a gradient (long wavelength) expansion up to second order for a relativistic system at vanishing charge density, one obtains the most general (causal) equations of motion for a fluid in the presence of shear and bulk viscosity, as well as the structure of the non-equilibrium entropy current. Requiring positivity of the divergence of the non-equilibrium entropy current relates some of its coefficients to those entering the equations of motion. I comment on possible applications of these results for conformal and non-conformal fluids.

preprint2009arXiv

Conformal Relativistic Viscous Hydrodynamics: Applications to RHIC results at sqrt(s_NN) = 200 GeV

A new set of equations for relativistic viscous hydrodynamics that captures both weak-coupling and strong-coupling physics to second order in gradients has been developed recently. We apply this framework to bulk physics at RHIC, both for standard (Glauber-type) as well as for Color-Glass-Condensate initial conditions and show that the results do not depend strongly on the values for the second-order transport coefficients. Results for multiplicity, radial flow and elliptic flow are presented and we quote the ratio of viscosity over entropy density for which our hydrodynamic model is consistent with experimental data. For Color-Glass-Condensate initial conditions, early thermalization does not seem to be required in order for hydrodynamics to describe charged hadron elliptic flow.

preprint2009arXiv

New Developments in Relativistic Viscous Hydrodynamics

Starting with a brief introduction into the basics of relativistic fluid dynamics, I discuss our current knowledge of a relativistic theory of fluid dynamics in the presence of (mostly shear) viscosity. Derivations based on the generalized second law of thermodynamics, kinetic theory, and a complete second-order gradient expansion are reviewed. The resulting fluid dynamic equations are shown to be consistent for all these derivations, when properly accounting for the respective region of applicability, and can be applied to both weakly and strongly coupled systems. In its modern formulation, relativistic viscous hydrodynamics can directly be solved numerically. This has been useful for the problem of ultrarelativistic heavy-ion collisions, and I will review the setup and results of a hydrodynamic description of experimental data for this case.

preprint2008arXiv

On the collision of two shock waves in AdS5

We consider two ultrarelativistic shock waves propagating and colliding in five-dimensional Anti-de-Sitter spacetime. By transforming to Rosen coordinates, we are able to find the form of the metric shortly after the collision. Using holographic renormalization, we calculate the energy-momentum tensor on the boundary of AdS space for early times after the collision. Via the gauge-gravity duality, this gives some insights on bulk dynamics of systems created by high energy scattering in strongly coupled gauge theories. We find that Bjorken boost-invariance is explicitely violated at early times and we obtain an estimate for the thermalization time in this simple system.

preprint2006arXiv

Collective non-Abelian instabilities in a melting Color Glass Condensate

We present first results for 3+1-D simulations of SU(2) Yang-Mills equations for matter expanding into the vacuum after a heavy ion collision. Violations of boost invariance cause a Weibel instability leading soft modes to grow with proper time $τ$ as $\exp(Γ\sqrt{g^2μτ})$, where $g^2μ$ is a scale arising from the saturation of gluons in the nuclear wavefunction. The scale for the growth rate $Γ$ is set by a plasmon mass, defined as $ω_{\rm pl}= κ_0 \sqrt{\frac{g^2μ}τ}$, generated dynamically in the collision. We compare the numerical ratio $Γ/κ_0$ to the corresponding value predicted by the Hard Thermal Loop formalism for anisotropic plasmas.

preprint2006arXiv

The Unstable Glasma

We discuss results from 3+1-D numerical simulations of SU(2) Yang--Mills equations for an unstable Glasma expanding into the vacuum after a high energy heavy ion collision. We expand on our earlier work on a non-Abelian Weibel instability in such a system and study the behavior of the instability in greater detail on significantly larger lattices than previously. We establish the time scale for the onset of the instability and demonstrate that the growth rate is robust as one approaches the continuum limit. For large violations of boost invariance, non-Abelian effects cause the growth of soft modes to saturate. At late times, we observe significant creation of longitudinal pressure and a systematic trend towards isotropy. These time scales however are significantly larger than those required for early thermalization in heavy ion collisions. We discuss additional effects in the produced Glasma that may speed up thermalization.