Source author record

Jorge Noronha

Jorge Noronha 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

68works
11topics
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

68 published item(s)

preprint2026arXiv

Lattice-based equation of state with a critical point from constant entropy contours and its comparison to effective QCD approaches

In this work, we systematically assess the performance of a new method from [H. Shah et al., Phys. Rev. C 113, L012201] for locating the QCD critical point using constant-entropy contours by testing it against various effective QCD approaches. We demonstrate that, while the method yields spurious critical points in purely hadronic models (HRG) due to non-parabolic contour behavior at low temperatures ($T \lesssim 120$ MeV), it accurately reproduces the CP location in frameworks that feature a genuine phase transition and benchmarked against lattice QCD, such as Holographic Einstein-Maxwell-Dilaton, and Functional QCD approaches. Building on our previous determination of constant entropy contours using lattice data, we extend that analysis to construct a complete Lattice-based Equation of State (EoS) at finite density, which features a critical point at $(T, μ_B) \approx (114, 602)$ MeV. By integrating the extrapolated entropy density with respect to temperature, we reconstruct the pressure, baryon density, susceptibility, and speed of sound in the critical region, and analyze the focusing behavior of isentropic trajectories in the vicinity of the critical point.

preprint2026arXiv

Merging multidimensional equations of state of strongly interacting matter via a statistical mixture

We introduce a general method to merge multidimensional equations of state (EoSs) by combining them in a two-fluid equilibrium statistical mixture in the grand canonical ensemble. The merged grand potential density $ω$ is built directly from the input EoSs and the fluid fractions are fixed by minimizing $ω$ at fixed temperature $T$ and baryon chemical potential $μ_B$. Thermodynamic consistency and stability are guaranteed as all thermodynamic quantities are consistently derived from a single merged grand potential $ω(T,μ_B)$ with the correct convexity properties. Our method can accommodate a first-order phase transition and a critical endpoint with mean-field critical exponents. We use this method to merge a van der Waals Hadron-Resonance-Gas EoS with a holographic Einstein-Maxwell-Dilaton EoS that has a critical point and a first-order line. The result is a single EoS, spanning hadronic and deconfined matter over a broad range in $(T,μ_B)$, which can be readily used in heavy-ion hydrodynamic simulations. Our merging method can be generalized to consider a higher dimensional phase diagram (e.g., by considering more chemical potentials) and more than two input EoSs.

preprint2026arXiv

Universality of scaled particle spectra in ultrarelativistic heavy-ion collisions

We study the transverse momentum spectra of identified particles in ultrarelativistic collisions of large and small collision systems. In order to isolate information contained in the momentum dependence, we propose to scale the spectra by the total particle number and mean transverse momentum -- global quantities which are already well studied. We observe an interesting, nearly universal, centrality-independent shape in the scaled spectra, similar to scalings that have been studied previously. This scaling behavior breaks down at large transverse momentum and for very small systems, such as those produced in p-p collisions. We perform hybrid hydrodynamic simulations and show that, in these simulations, a centrality-independent shape is a consequence of an event-by-event independence. Our results motivate further theoretical and experimental investigations of the regime of validity of this scaling phenomenon and their physical interpretation at different collision energies and systems.

preprint2022arXiv

First-Order General-Relativistic Viscous Fluid Dynamics

We present the first generalization of Navier-Stokes theory to relativity that satisfies all of the following properties: (a) the system coupled to Einstein's equations is causal and strongly hyperbolic; (b) equilibrium states are stable; (c) all leading dissipative contributions are present, i.e., shear viscosity, bulk viscosity, and thermal conductivity; (d) non-zero baryon number is included; (e) entropy production is non-negative in the regime of validity of the theory; (f) all of the above holds in the nonlinear regime without any simplifying symmetry assumptions. These properties are accomplished using a generalization of Eckart's theory containing only the hydrodynamic variables, so that no new extended degrees of freedom are needed as in Müller-Israel-Stewart theories. Property (b), in particular, follows from a more general result that we also establish, namely, sufficient conditions that when added to stability in the fluid's rest frame imply stability in any reference frame obtained via a Lorentz transformation. All our results are mathematically rigorously established. The framework presented here provides the starting point for systematic investigations of general-relativistic viscous phenomena in neutron star mergers.

preprint2022arXiv

How Do Spherical Black Holes Grow Monopole Hair?

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.

preprint2022arXiv

Hydrodynamization times of a holographic fluid far from equilibrium

We investigate several hydrodynamization times for an ensemble of different far-from-equilibrium solutions of the strongly coupled $\mathcal{N}=4$ Supersymmetric Yang-Mills plasma undergoing Bjorken flow. For the ensemble of initial data analyzed in the present work, we find that, with typical tolerances between $3\%$ to $5\%$, the average hydrodynamization time associated with the late time convergence of the pressure anisotropy to the corresponding Borel resummed hydrodynamic attractor is approximately equal to the average hydrodynamization time associated with the Navier-Stokes result, while both are shorter than the average hydrodynamization time associated with second-order hydrodynamics. On the other hand, we find that the entropy density of the different solutions coalesces to second-order hydrodynamics long before entering in the Navier-Stokes regime. A clear hierarchy between the different average hydrodynamization times of the Bjorken expanding fluid is established for the set of analyzed initial data, comprising also some solutions which, whilst satisfying the dominant and the weak energy conditions at the initial time, evolve such as to transiently violate one or both conditions when the fluid is still far from equilibrium. In particular, solutions violating the weak energy condition are generally found to take a longer time to enter in the hydrodynamic regime than the other solutions.

preprint2022arXiv

Modeling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modeling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron-ion or electron-positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann-Vlasov equation. Crucially, this new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. In particular, we keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion disks. This new formulation will be instrumental in the construction of a large class of next-generation simulations of relativistic transient phenomena produced around black holes and neutron stars.

preprint2022arXiv

Nonlinear causality of general first-order relativistic viscous hydrodynamics

Effective theory arguments are used to derive the most general energy-momentum tensor of a relativistic viscous fluid with an arbitrary equation of state (in the absence of other conserved currents) that is first-order in the derivatives of the energy density and flow velocity and does not include extended variables such as in Mueller-Israel-Stewart-like theories. This energy-momentum tensor leads to a causal theory, provided one abandons the usual conventions for the out-of-equilibrium hydrodynamic variables put forward by Landau-Lifshitz and Eckart. In particular, causality requires nonzero out-of-equilibrium energy density corrections and heat flow. Conditions are found to ensure linear stability around equilibrium in flat space-time. We also prove local existence and uniqueness of solutions to the equations of motion. Our causality, existence, and uniqueness results hold in the full nonlinear regime, without symmetry assumptions, in four space-time dimensions, with or without coupling to Einstein's equations, and are mathematically rigorously established. Furthermore, a kinetic theory realization of this energy-momentum tensor is also provided.

preprint2022arXiv

Perturbative approaches in relativistic kinetic theory and the emergence of first-order hydrodynamics

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.

preprint2022arXiv

QCD Equation of State and Phase Diagram from Holographic Black Holes

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.

preprint2022arXiv

QCD Equilibrium and Dynamical Properties from Holographic Black Holes

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.

preprint2022arXiv

Transport coefficients of quasi-particle models within a new relaxation time approximation of the Boltzmann equation

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.

preprint2020arXiv

Connecting far-from-equilibrium hydrodynamics to resummed transport coefficients and attractors

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

preprint2020arXiv

Equivalence between first-order causal and stable hydrodynamics and Israel-Stewart theory for boost-invariant systems with a constant relaxation time

We show that the recently formulated causal and stable first-order hydrodynamics has the same dynamics as Israel-Stewart theory for boost-invariant, Bjorken expanding systems with a conformal equation of state and a regulating sector determined by a constant relaxation time. In this case, the general solution of the new first-order formulation can be determined analytically.

preprint2020arXiv

Exact hydrodynamic attractor of an ultrarelativistic gas of hard spheres

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

preprint2020arXiv

Measuring Momentum-Dependent Flow Fluctuations in Heavy-Ion Collisions

In heavy-ion collisions, momentum-dependent pair correlations can be characterized by a principal component analysis (PCA), in which subleading modes are expected to reveal new information on flow fluctuations. However, we find that, as currently measured, these modes can be dominated by multiplicity fluctuations, which serve as an unwanted background. Here, we propose new PCA observables that are robust against multiplicity fluctuations and isolate novel sources of flow fluctuations, thus being suited to provide fresh insight into the initial stages of the system at small length scales.

preprint2020arXiv

Momentum-dependent flow fluctuations as a hydrodynamic response to initial geometry

We propose a redefinition of the principal component analysis (PCA) of anisotropic flow that makes it more directly connected to fluctuations of the initial geometry of the system. Then, using state-of-the-art hydrodynamic simulations, we make an explicit connection between flow fluctuations and a cumulant expansion of the initial transverse geometry. In particular, we show that the second principal component of elliptic flow is generated by higher-order cumulants, and therefore probes smaller length scales of the initial state. With this information, it will be possible to put new constraints on properties of the early-time dynamics of a heavy-ion collision, including small-scale structure, as well as properties of the quark-gluon plasma.

preprint2020arXiv

Probing the structure of the initial state of heavy-ion collisions with $p_T$-dependent flow fluctuations

The connection between initial-state geometry and anisotropic flow can be quantified through a well-established mapping between $p_T$-integrated flow harmonics and cumulants of the initial transverse energy distribution. In this paper we successfully extend this mapping to also include $p_T$-differential flow. In doing so, we find that subleading principal components of anisotropic flow can reveal previously unobserved details of the hydrodynamic response, in both the linear and the nonlinear regimes. Most importantly, we show that they provide novel information on the small-scale structures present in the initial stage of relativistic heavy-ion collisions.

preprint2020arXiv

System response to the initial energy-momentum tensor in relativistic heavy-ion collisions

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 summarize a proposed extension of this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. Numerical tests validate this proposal.

preprint2019arXiv

DAB-MOD sensitivity study of heavy flavor $R_{AA}$ and azimuthal anisotropies based on beam energy, initial conditions, hadronization, and suppression mechanisms

Heavy flavor probes provide important information about the in-medium properties of the quark-gluon plasma produced in heavy-ion collisions. In this work, we investigate the effects of 2D+1 event-by-event fluctuating hydrodynamic backgrounds on the nuclear suppression factor and momentum anisotropies of heavy flavor mesons and non-photonic electrons. Using the state-of-the-art D and B mesons modular simulation code (called "DAB-MOD"), we perform a systematic comparison of different transport equations in the same background, including a few energy loss models --- with and without energy loss fluctuations --- and a relativistic Langevin model with different drag parametrizations. We present the resulting D and B mesons $R_{AA}$, $v_2$, $v_3$, and $v_4$ as well as multi-particle cumulants, in AuAu collisions at $\sqrt{s_\text{NN}}=200\, \text{GeV}$ and PbPb collisions at $\sqrt{s_\text{NN}}=2.76 \, \text{TeV}$ and $\sqrt{s_\text{NN}}=5.02 \, \text{TeV}$, and compare them to the available experimental data. The $v_2\{4\}/v_2\{2\}$ ratio, which is known to be a powerful probe of the initial conditions and flow fluctuations in the soft sector, is also studied in the context of heavy flavor. We also investigate the correlations between the transverse anisotropies of heavy mesons and all charged particles to better understand how heavy quarks couple to the hydrodynamically expanding quark-gluon plasma. We study the influence that different initial conditions and the implementation of heavy-light quark coalescence has on our results.

preprint2016arXiv

Divergence of the Chapman-Enskog expansion in relativistic kinetic theory

In this letter we show for the first time that the relativistic Chapman-Enskog series for a massless gas undergoing Bjorken expansion diverges. In order to fix this problem, we propose a novel type of expansion that includes non-perturbative contributions in the Knudsen number that are not considered in Chapman-Enskog theory. This approach is in good agreement with exact solutions of the Boltzmann equation for a wide range of values of Knudsen number and does not display the clear signs of divergence exhibited by the Chapman-Enskog series.

preprint2016arXiv

Energy loss, equilibration, and thermodynamics of a baryon rich strongly coupled quark-gluon plasma

Lattice data for the QCD equation of state and the baryon susceptibility near the crossover phase transition (at zero baryon density) are used to determine the input parameters of a 5-dimensional Einstein-Maxwell-Dilaton holographic model that provides a consistent holographic framework to study both equilibrium and out-of-equilibrium properties of a hot and {\it baryon rich} strongly coupled quark-gluon plasma (QGP). We compare our holographic equation of state computed at nonzero baryon chemical potential, $μ_B$, with recent lattice calculations and find quantitative agreement for the pressure and the speed of sound for $μ_B \leq 400$ MeV. This holographic model is used to obtain holographic predictions for the temperature and $μ_B$ dependence of the drag force and the Langevin diffusion coefficients associated with heavy quark jet propagation as well as the jet quenching parameter $\hat{q}$ and the shooting string energy loss of light quarks in the baryon dense plasma. We find that the energy loss of heavy and light quarks generally displays a nontrivial, fast-varying behavior as a function of the temperature near the crossover. Moreover, energy loss is also found to generally increase due to nonzero baryon density effects even though this strongly coupled liquid cannot be described in terms of well defined quasiparticle excitations. Furthermore, to get a glimpse of how thermalization occurs in a hot and baryon dense QGP, we study how the lowest quasinormal mode of an external massless scalar disturbance in the bulk is affected by a nonzero baryon charge. We find that the equilibration time associated with the lowest quasinormal mode decreases in a dense medium.

preprint2016arXiv

Entropy Production, Hydrodynamics, and Resurgence in the Primordial Quark-Gluon Plasma from Holography

Microseconds after the Big Bang quarks and gluons formed a strongly-coupled non-conformal liquid driven out-of-equilibrium by the expansion of the Universe. We use holography to determine the non-equilibrium behavior of this liquid in a Friedmann-Lemaitre-Robertson-Walker Universe and develop an expansion for the corresponding entropy production in terms of the derivatives of the cosmological scale factor. We show that the resulting series has zero radius of convergence and we discuss its resurgent properties. Finally, we compute the resummed entropy production rate in de Sitter Universe at late times and show that the leading order approximation given by bulk viscosity effects can strongly overestimate/underestimate the rate depending on the microscopic parameters.

preprint2016arXiv

Event-by-event hydrodynamics $+$ jet energy loss: A solution to the $R_{AA} \otimes v_2$ puzzle

High $p_T > 10$ GeV elliptic flow, which is experimentally measured via the correlation between soft and hard hadrons, receives competing contributions from event-by-event fluctuations of the low $p_T$ elliptic flow and event plane angle fluctuations in the soft sector. In this paper, a proper account of these event-by-event fluctuations in the soft sector, modeled via viscous hydrodynamics, is combined with a jet energy loss model to reveal that the positive contribution from low $p_T$ $v_2$ fluctuations overwhelms the negative contributions from event plane fluctuations. This leads to an enhancement of high $p_T > 10$ GeV elliptic flow in comparison to previous calculations and provides a natural solution to the decade long high $p_T$ $R_{AA} \otimes v_2$ puzzle. We also present the first theoretical calculation of high $p_T$ $v_3$, which is shown to be compatible with current LHC data. Furthermore, we discuss how short wavelength jet-medium physics can be deconvoluted from the physics of soft, bulk event-by-event flow observables using event shape engineering techniques.

preprint2016arXiv

Holographic calculation of the QCD crossover temperature in a magnetic field

Lattice data for the QCD equation of state and the magnetic susceptibility computed near the crossover transition at zero magnetic field are used to determine the input parameters of a five dimensional Einstein-Maxwell-Dilaton holographic model. Once the model parameters are fixed at zero magnetic field, one can use this holographic construction to study the effects of a magnetic field on the equilibrium and transport properties of the quark-gluon plasma. In this paper we use this model to study the dependence of the crossover temperature with an external magnetic field. Our results for the pressure of the plasma and the crossover temperature are in quantitative agreement with current lattice data for values of the magnetic field $0 \le eB \lesssim 0.3$ GeV$^2$, which is the relevant range for ultrarelativistic heavy ion collision applications.

preprint2016arXiv

Polyakov loop and heavy quark entropy in strong magnetic fields from holographic black hole engineering

We investigate the temperature and magnetic field dependence of the Polyakov loop and heavy quark entropy in a bottom-up Einstein-Maxwell-dilaton (EMD) holographic model for the strongly coupled quark-gluon plasma (QGP) that quantitatively matches lattice data for the $(2+1)$-flavor QCD equation of state at finite magnetic field and physical quark masses. We compare the holographic EMD model results for the Polyakov loop at zero and nonzero magnetic fields and the heavy quark entropy at vanishing magnetic field with the latest lattice data available for these observables and find good agreement for temperatures $T\gtrsim 150$ MeV and magnetic fields $eB\lesssim 1$ GeV$^2$. Predictions for the behavior of the heavy quark entropy at nonzero magnetic fields are made that could be readily tested on the lattice.

preprint2015arXiv

Heavy flavor electron $R_\text{AA}$ and $v_2$ in event-by-event relativistic hydrodynamics

In this work we investigate how event-by-event hydrodynamics fluctuations affect the nuclear suppression factor and elliptic flow of heavy flavor mesons and non-photonic electrons. We use a 2D+1 Lagrangian ideal hydrodynamic code on an event-by-event basis in order to compute local temperature and flow profiles. Using a strong coupling inspired energy loss parametrization on top of the evolving space-time energy density distributions we are able to propagate the heavy quarks inside the medium until the freeze-out temperature is reached and a Pythia modeling of hadronization takes place. The resulting D$^0$ and heavy-flavor electron yield is compared with recent experimental data for $R_\text{AA}$ and $v_2$ from the STAR and Phenix collaborations. In addition we present preditions for the higher order Fourier harmonic coefficients $v_3(p_T)$ of heavy-flavor electrons at RHIC's $\sqrt{S_\text{NN}} = 200$ GeV collisions.

preprint2015arXiv

Hydrodynamic transport coefficients for the non-conformal quark-gluon plasma from holography

In this paper we obtain holographic formulas for the transport coefficients $κ$ and $τ_π$ present in the second-order derivative expansion of relativistic hydrodynamics in curved spacetime associated with a non-conformal strongly coupled plasma described holographically by an Einstein+Scalar action in the bulk. We compute these coefficients as functions of the temperature in a bottom-up non-conformal model that is tuned to reproduce lattice QCD thermodynamics at zero baryon chemical potential. We directly compute, besides the speed of sound, 6 other transport coefficients that appear at second-order in the derivative expansion. We also give an estimate for the temperature dependence of 11 other transport coefficients taking into account the simplest contribution from non-conformal effects that appear near the QCD crossover phase transition. Using these results, we construct an Israel-Stewart-like theory in flat spacetime containing 13 of these 17 transport coefficients that should be suitable for phenomenological applications in the context of numerical hydrodynamic simulations of the strongly-coupled, non-conformal quark-gluon plasma. Using several different approximations, we give parametrizations for the temperature dependence of all the second-order transport coefficients that appear in this theory in a format that can be easily implemented in existing numerical hydrodynamic codes.

preprint2015arXiv

Perfect fluidity of a dissipative system: Analytical solution for the Boltzmann equation in $\mathrm{AdS}_{2}\otimes \mathrm{S}_{2}$

In this paper we obtain an analytical solution of the relativistic Boltzmann equation under the relaxation time approximation that describes the out-of-equilibrium dynamics of a radially expanding massless gas. This solution is found by mapping this expanding system in flat spacetime to a static flow in the curved spacetime $\mathrm{AdS}_{2}\otimes \mathrm{S}_{2}$. We further derive explicit analytic expressions for the momentum dependence of the single particle distribution function as well as for the spatial dependence of its moments. We find that this dissipative system has the ability to flow as a perfect fluid even though its entropy density does not match the equilibrium form. The non-equilibrium contribution to the entropy density is shown to be due to higher order scalar moments (which possess no hydrodynamical interpretation) of the Boltzmann equation that can remain out of equilibrium but do not couple to the energy-momentum tensor of the system. Thus, in this system the slowly moving hydrodynamic degrees of freedom can exhibit true perfect fluidity while being totally decoupled from the fast moving, non-hydrodynamical microscopic degrees of freedom that lead to entropy production.

preprint2015arXiv

Sensitivity of flow harmonics to sub-nucleon scale fluctuations in heavy ion collisions

In this paper a new procedure to smooth out the initial energy densities of hydrodynamics is employed to show that the initial spatial eccentricities $\varepsilon_{m,n}$, which drive the final state flow harmonics $v_n$, are remarkably robust with respect to variations of the underlying scale of initial energy density spatial gradients, $λ$, in nucleus-nucleus collisions. For $\sqrt{s}=2.76$ TeV Pb+Pb collisions, the $\varepsilon_{m,n}$'s (across centrality classes) change by less than $10\%$ if the scale of fluctuations is varied from $0.1$ to $1$ fm. We show, using the 2+1 Lagrangian hydrodynamic code, v-USPhydro, that this robustness is transferred to the final $v_n$'s computed within event by event viscous hydrodynamics. This indicates that the flow harmonics in nucleus-nucleus collisions are not yet particularly sensitive to the underlying microscopic sub-nucleon physics below the confinement scale. On the other hand, the eccentricities of top $1\%$ high multiplicity $\sqrt{s}=5.02$ TeV p+Pb collisions are found to be very sensitive to sub-nucleonic scale fluctuations, which should be contrasted with the robustness found in peripheral Pb+Pb collisions with the same multiplicity.

preprint2015arXiv

Suppression of baryon diffusion and transport in a baryon rich strongly coupled quark-gluon plasma

Five dimensional black hole solutions that describe the QCD crossover transition seen in $(2+1)$-flavor lattice QCD calculations at zero and nonzero baryon densities are used to obtain predictions for the baryon susceptibility, baryon conductivity, baryon diffusion constant, and thermal conductivity of the strongly coupled quark-gluon plasma in the range of temperatures $130\,\textrm{MeV}\le T\le 300\,\textrm{MeV}$ and baryon chemical potentials $0\le μ_B \le 400\,\textrm{MeV}$. Diffusive transport is predicted to be suppressed in this region of the QCD phase diagram, which is consistent with the existence of a critical end point at larger baryon densities. We also calculate the fourth-order baryon susceptibility at zero baryon chemical potential and find quantitative agreement with recent lattice results. The baryon transport coefficients computed in this paper can be readily implemented in state-of-the-art hydrodynamic codes used to investigate the dense QGP currently produced at RHIC's low energy beam scan.

preprint2015arXiv

The unreasonable effectiveness of hydrodynamics in heavy ion collisions

Event-by-event hydrodynamic simulations of AA and pA collisions involve initial energy densities with large spatial gradients. This is associated with the presence of large Knudsen numbers ($K_n\approx 1$) at early times, which may lead one to question the validity of the hydrodynamic approach in these rapidly evolving, largely inhomogeneous systems. A new procedure to smooth out the initial energy densities is employed to show that the initial spatial eccentricities, $\varepsilon_n$, are remarkably robust with respect to variations in the underlying scale of initial energy density spatial gradients, $λ$. For $\sqrt{s_{NN}}=2.76$ TeV LHC initial conditions generated by the MCKLN code, $\varepsilon_n$ (across centralities) remains nearly constant if the fluctuation scale varies by an order of magnitude, i.e., when $λ$ varies from 0.1 to 1 fm. Given that the local Knudsen number $K_n\approx \frac{1}λ$, the robustness of the initial eccentricities with respect to changes in the fluctuation scale suggests that the vn's cannot be used to distinguish between events with large $K_n$ from events where $K_n$ is in the hydrodynamic regime. We use the 2+1 Lagrangian hydrodynamic code v-USPhydro to show that this is indeed the case: anisotropic flow coefficients computed within event-by-event viscous hydrodynamics are only sensitive to long wavelength scales of order $\frac{1}{Λ_{QCD}}\approx 1$ fm and are incredibly robust with respect to variations in the initial local Knudsen number. This robustness can be used to justify the somewhat unreasonable effectiveness of the nearly perfect fluid paradigm in heavy ion collisions.

preprint2015arXiv

Thermal suppression of moving heavy quark pairs in strongly coupled plasma

The imaginary part of the heavy quark-antiquark potential experienced by moving heavy quarkonia in strongly coupled plasmas dual to theories of gravity is computed by considering thermal worldsheet fluctuations of the holographic Nambu-Goto string. General results for a wide class of gravity duals are presented and an explicit formula for $\mathrm{Im} \, V_{Q\bar{Q}}$ is found in the case where the axis of the moving $\bar{Q}Q$ pair has an arbitrary orientation with respect to its velocity in the plasma. These results are applied to the study of heavy quarkonia propagating through a strongly coupled $\mathcal{N} = 4$ SYM plasma. Our results indicate that the onset of $\mathrm{Im} \, V_{Q\bar{Q}}$ decreases with increasing rapidity (though our analysis is limited to slowly moving quarkonia) and that, in general, a $Q\bar{Q}$ pair is more strongly bound if its axis is aligned with its direction of motion through the strongly coupled plasma.

preprint2015arXiv

Vanishing DC holographic conductivity from a magnetic monopole condensate

We show how to obtain a vanishing DC conductivity in 3-dimensional strongly coupled QFT's using a massive 2-form field in the bulk that satisfies a special kind of boundary condition. The real and imaginary parts of the AC conductivity are evaluated in this holographic setup and we show that the DC conductivity identically vanishes even for an arbitrarily small (though nonzero) value of the 2-form mass in the bulk. We identify the bulk action of the massive 2-form with an effective theory describing a phase in which magnetic monopoles have condensed in the bulk. Our results indicate that a condensate of magnetic monopoles in a 4-dimensional bulk leads to a vanishing DC holographic conductivity in 3-dimensional strongly coupled QFT's.

preprint2014arXiv

A new exact solution of the relativistic Boltzmann equation and its hydrodynamic limit

We present an exact solution of the relativistic Boltzmann equation for a system undergoing boost-invariant longitudinal and azimuthally symmetric transverse flow ("Gubser flow"). The resulting exact non-equilibrium dynamics is compared to 1st- and 2nd-order relativistic hydrodynamic approximations for various shear viscosity to entropy density ratios. This novel solution can be used to test the validity and accuracy of different hydrodynamic approximations in conditions similar to those generated in relativistic heavy-ion collisions.

preprint2014arXiv

A systematic study of exact solutions in second-order conformal hydrodynamics

In this paper we present the details of our previous work on exact solutions in second-order conformal hydrodynamics together with a number of new solutions found by mapping Minkowski space onto various curved spacetimes such as anti-de Sitter space and hyperbolic space. We analytically show how the solutions of ideal hydrodynamics are modified by the second-order effects including vorticity. We also find novel boost-invariant exact solutions which exist only in second-order hydrodynamics and have an unusual dependence on the proper time.

preprint2014arXiv

Azimuthal asymmetries in high-energy collisions of protons with holographic shockwaves

Large azimuthal quadrupole and octupole asymmetries have recently been found in p+Pb collisions at the LHC. We argue that these might arise from a projectile dipole scattering off fluctuations in the target with a size on the order of the dipole. In a holographic scenario, parity even angular moments $v_{2n}$ are generated by the real part of the light-like Wilson loop due to the contribution from the background metric to the Nambu-Goto action. On the other hand, parity odd moments $v_{2n+1}$ must arise from the imaginary part of the light-like Wilson loop, which is naturally induced by a fluctuating Neveu-Schwarz 2-form.

preprint2014arXiv

Elliptic Flow Suppression due to Hadron Mass Spectrum

Hadron resonance gas models provide a good description of the equation of state of quantum chromodynamics determined by lattice QCD calculations at temperatures $T \sim 100-155$ MeV. In this paper we investigate the effects of an exponentially increasing hadron mass spectrum (Hagedorn spectrum) on the azimuthal anisotropy of the rapidly expanding matter formed in ultrarelativistic heavy ion collisions. If the temperature at which the conversion from fluid degrees of freedom to hadrons is sufficiently close to the Hagedorn temperature, the production of Hagedorn resonances suppresses the differential elliptic flow of all hadron species.

preprint2014arXiv

Exact analytical solutions of second-order conformal hydrodynamics

We present some exact solutions of relativistic second-order hydrodynamic equations in theories with conformal symmetry. Starting from a spherically expanding solution in ideal hydrodynamics, we take into account general conformal second-order corrections, and construct, for the first time, fully analytical axisymmetric exact solutions including the case with nonzero vorticity. These solutions are time-reversible despite having a nonvanishing shear stress tensor and provide a useful quantitative measure of the second-order effects in relativistic hydrodynamics.

preprint2014arXiv

Flow harmonics within an analytically solvable viscous hydrodynamic model

Based on a viscous hydrodynamic model with anisotropically perturbed Gubser flow and isothermal Cooper-Frye freezeout at early times, we analytically compute the flow harmonics $v_n(p_T)$ and study how they scale with the harmonic number $n$ and transverse momentum, as well as the system size, shear and bulk viscosity coefficients, and collision energy. In particular, we find that the magnitude of shear viscous corrections grows linearly with $n$. The mixing between different harmonics is also discussed. While this model is rather simple as compared to realistic heavy-ion collisions, we argue that the scaling results presented here may be meaningfully compared to experimental data collected over many energies, system sizes, and geometries.

preprint2014arXiv

Jet quenching effects on the anisotropic flow at RHIC

In this article we investigate how the energy and momentum deposited by partonic dijets in the quark-gluon plasma may perturb the geometry-induced hydrodynamic expansion of the bulk nuclear matter created in heavy ion collisions at RHIC. The coupling between the jets and the medium is done through a source term in the energy-momentum conservation equations for ideal hydrodynamics. We concentrate our attention at mid-rapidity and solve the equations event-by-event imposing boost-invariance. For $p_T \gtrsim 1$ GeV the anisotropic flow is found to be considerably enhanced, if the dijets deposit on average more than 12 GeV in the medium (or equivalently 6 GeV for each jet of the pair), which corresponds, in our model, to an average suppression greater than $65\%$ of the initial jet transverse energy.

preprint2014arXiv

Jet quenching effects on the direct, elliptic, and triangular flow at RHIC

In this paper we investigate how the energy and momentum deposited by partonic dijets in the quark-gluon plasma may affect the direct, elliptic and triangular flow of low (and intermediate) $p_T$ hadrons in central Au+Au collisions at RHIC. The dijets are modeled as external sources in the energy-momentum conservation equations for hydrodynamics, which are solved on an event-by-event basis within the ideal fluid approximation. We focus our investigation at mid-rapidity and solve the hydrodynamic equations imposing boost invariance. Differential anisotropic flow coefficients for $p_T \gtrsim 1$ GeV are found to be significantly enhanced if the dijets deposit on average more than 12 GeV in the QGP (or more than 6 GeV per jet). Because this jet-induced extra anisotropic flow is not related to the fluctuations of the initial geometry of the collision, the correlation between the $v_2$ and $v_3$ coefficients and their corresponding eccentricities is considerably weakened. In addition, we argue that the extra amount of direct flow induced by dijets may be quantified by comparing the azimuthal dependence of dihadron correlations in dijet events with the corresponding quantity obtained in events without dijets. This comparison could be used to give a rough estimate of the magnitude of the effective coupling between the jets and the medium.

preprint2014arXiv

Studying the validity of relativistic hydrodynamics with a new exact solution of the Boltzmann equation

We present an exact solution to the Boltzmann equation which describes a system undergoing boost-invariant longitudinal and azimuthally symmetric radial expansion for arbitrary shear viscosity to entropy density ratio. This new solution is constructed by considering the conformal map between Minkowski space and the direct product of three dimensional de Sitter space with a line. The resulting solution respects SO(3)_q x SO(1,1) x Z_2 symmetry. We compare the exact kinetic solution with exact solutions of the corresponding macroscopic equations that were obtained from the kinetic theory in ideal and second-order viscous hydrodynamic approximations. The macroscopic solutions are obtained in de Sitter space and are subject to the same symmetries used to obtain the exact kinetic solution.

preprint2013arXiv

A holographic calculation of the electric conductivity of the strongly coupled quark-gluon plasma near the deconfinement transition

The frequency dependent conductivity $σ(ω)$ of the strongly coupled Quark-Gluon Plasma (QGP) is estimated using a bottom up holographic model that can adequately describe recent lattice data for QCD thermodynamics at zero chemical potential. Different choices for the coupling between the bulk gauge field and the other bulk fields that define the background (the metric and a scalar field) are used in order to fit the lattice data for the electric charge susceptibility $χ_2^Q/T^2$. The ratio $σ_{DC}/T$ is found to vary near the deconfinement transition in a way that is similar to recent lattice results. This model is used to compute the charge diffusion coefficient $D$ of the strongly coupled plasma. We find that the dimensionless combination $DT$ has the same type of temperature dependence displayed by $σ_{DC}/T$ and, thus, charge diffusion is suppressed at low temperatures. The frequency dependent conductivity $σ(ω)$ reveals some nontrivial structure for values of the temperature near the phase transition. None of these structures appear in the associated Euclidean correlator, which we also compute. Our results suggest that the conformal invariance violation near the QCD deconfinement phase transition may be seen in the Euclidean correlator through a downward shift of its value at the minimum, which gives a rough estimate of the temperature dependence of the DC conductivity in the plasma.

preprint2013arXiv

Bulk Viscosity Effects in Event-by-Event Relativistic Hydrodynamics

Bulk viscosity effects on the collective flow harmonics in heavy ion collisions are investigated, on an event by event basis, using a newly developed 2+1 Lagrangian hydrodynamic code named v-USPhydro which implements the Smoothed Particle Hydrodynamics (SPH) algorithm for viscous hydrodynamics. A new formula for the bulk viscous corrections present in the distribution function at freeze-out is derived starting from the Boltzmann equation for multi-hadron species. Bulk viscosity is shown to enhance the collective flow Fourier coefficients from $v_2(p_T)$ to $v_5(p_T)$ when $% p_{T}\sim 1-3$ GeV even when the bulk viscosity to entropy density ratio, $% ζ/s$, is significantly smaller than $1/(4π)$.

preprint2013arXiv

Estimates for the Thermal Width of Heavy Quarkonia in Strongly Coupled Plasmas from Holography

The gauge/gravity duality is used to investigate the imaginary part of the heavy quark potential (defined via the rectangular Wilson loop) in strongly coupled plasmas. This quantity can be used to estimate the width of heavy quarkonia in a plasma at strong coupling. In this paper the thermal worldsheet fluctuation method, proposed in [J.Noronha and A.Dumitru, Phys.\ Rev.\ Lett.\ {\bf 103}, 152304 (2009)], is revisited and general conditions for the existence of an imaginary part for the heavy quark potential computed within classical gravity models are obtained. We prove a general result that establishes the connection between this imaginary part of the potential determined holographically and the area law displayed by the Wilson loop in the vacuum of confining gauge theories. We also determine the imaginary part of the heavy quark potential in a strongly coupled plasma dual to Gauss-Bonnet gravity. This provides an estimate of how the thermal width of heavy quarkonia changes with the shear viscosity to entropy density ratio, $η/s$, at strong coupling.

preprint2013arXiv

Large Nc Deconfinement Transition in the Presence of a Magnetic Field

We investigate the effect of a homogeneous magnetic field on the thermal deconfinement transition of QCD in the large $N_c$ limit. First we discuss how the critical temperature decreases due to the inclusion of $N_f \ll N_c$ flavors of massless quarks. Then we study the equivalent correction in the presence of an external Abelian magnetic field. To leading order in $N_{f}/N_{c}$, the fact that the deconfinement critical temperature decreases with the magnetic field depends solely on quarks behaving paramagnetically. Finally, we discuss the effects from a finite quark mass and its competition with magnetic effects.

preprint2013arXiv

Modification of the B Meson Mass in a Magnetic Field from QCD Sum Rules

In this paper we extend the well known QCD sum rules used in the calculation of the mass of heavy mesons to estimate the modification of the charged B meson mass, m_B, in the presence of an external Abelian magnetic field, eB. Two simplifying limits were considered: the weak field limit in which the external field satisfies eB << m^2 (with m being any of the masses involved) and the strong field limit in which the field strength is small in comparison to the bottom quark mass (or the B meson mass) squared but it is large compared to the mass of the light quarks, i.e., m^2_{u,d} << eB << m^2_{b,B}. We found that m_B decreases with the magnetic field in the both of these limits.

preprint2013arXiv

On the fluid behavior of a baryon rich hadron resonance gas

We investigate the effects of finite baryon chemical potential on the transport properties of a hadron resonance gas. We find that a hadron resonance gas with large baryon number density is closer to the ideal fluid limit than the corresponding gas with zero baryon number. This suggests that the system created at the Relativistic Heavy Ion Collider (RHIC) at lower collision energies may behave as a fluid, with an effective fluidity close to the one found at RHIC's highest energy near phase transition. This might explain why the differential elliptic flow coefficient measured at lower collisional energies at RHIC is similar to the one observed at high energies.

preprint2013arXiv

Solutions of Conformal Israel-Stewart Relativistic Viscous Fluid Dynamics

We use symmetry arguments developed by Gubser to construct the first radially-expanding explicit solutions of the Israel-Stewart formulation of hydrodynamics. Along with a general semi-analytical solution, an exact analytical solution is given which is valid in the cold plasma limit where viscous effects from shear viscosity and the relaxation time coefficient are important. The radially expanding solutions presented in this paper can be used as nontrivial checks of numerical algorithms employed in hydrodynamic simulations of the quark-gluon plasma formed in ultra-relativistic heavy ion collisions. We show this explicitly by comparing such analytic and semi-analytic solutions with the corresponding numerical solutions obtained using the MUSIC viscous hydrodynamics simulation code.

preprint2013arXiv

Transient Fluid Dynamics of the Quark-Gluon Plasma According to AdS/CFT

We argue, using the AdS/CFT correspondence, that the transient dynamics of the shear stress tensor in a strongly coupled $\mathcal{N}=4$ SYM plasma is not described by relaxation-type, fluid dynamical equations: at long times the equations of motion should contain a \textit{second-order} comoving derivative of the shear stress tensor. This occurs because in this strongly-coupled system the lowest "non-hydrodynamical" quasinormal modes associated with shear stress possess a nonzero real part at zero wavenumber. We use Weyl invariance to obtain the most general equations of motion containing 2 comoving derivatives of the shear stress tensor in the transient regime that are compatible with the symmetries. We show that the asymptotic solution of this theory valid at times much larger than the timescale associated with the "non-hydrodynamical" modes reproduces the well-known results previously obtained directly from the AdS/CFT correspondence. If the QGP formed in heavy ion collisions can be at least qualitatively understood in terms of strongly-coupled $\mathcal{N}=4$ SYM theory, the second time derivative present in the equations of motion of the fluid may lead to an unexpected dependence on the initial conditions for the shear stress tensor needed in numerical hydrodynamic simulations.

preprint2012arXiv

Falling Strings and Light Quark Jet Quenching at LHC

We explore phenomenological signatures of light quark jet quenching within the AdS/CFT correspondence. Numerical studies of the instantaneous energy loss of light quarks, modeled as falling strings, suggest a linear path dependence. We propose a phenomenological model for the energy loss and use it to compute the nuclear modification factor $R_{AA}$ for light quarks in an expanding plasma with Glauber initial conditions. The results are compared to the light hadron $R_{AA}$ data at the LHC and, although they show qualitative agreement, the quantitative disagreement we found motivated the exploration of effects from higher order derivative corrections to $AdS_5$ on this observable.

preprint2012arXiv

Further results on peripheral-tube model for ridge correlation

Peripheral one-tube model has shown to be a nice tool for dynamically understanding several aspects of ridge structures in long-range two-particle correlations, observed experimentally and obtained also in our model calculations using NexSPheRIO code. Here, we study an extension of the model, to initial configurations with several peripheral tubes distributed randomly in azimuth. We show that the two-particle correlation is almost independent of the number of tubes, although the flow distribution becomes indeed strongly event dependent. In our picture, the ridge structures are causally connected not only in the longitudinal direction but also in azimuth.

preprint2012arXiv

Hadron Mass Spectrum and the Shear Viscosity to Entropy Density Ratio of Hot Hadronic Matter

Lattice calculations of the QCD trace anomaly at temperatures $T<160$ MeV have been shown to match hadron resonance gas model calculations, which include an exponentially rising hadron mass spectrum. In this paper we perform a more detailed comparison of the model calculations to lattice data that confirms the need for an exponentially increasing density of hadronic states. Also, we find that the lattice data is compatible with a hadron density of states that goes as $ρ(m) \sim m^{-a}\exp(m/T_H) $ at large $m$ with $a> 5/2$ (where $T_H \sim 167$ MeV). With this specific subleading contribution to the density of states, heavy resonances are most likely undergo 2-body decay (instead of multi-particle decay), which facilitates their inclusion into hadron transport codes. Moreover, estimates for the shear viscosity and the shear relaxation time coefficient of the hadron resonance model computed within the excluded volume approximation suggest that these transport coefficients are sensitive to the parameters that define the hadron mass spectrum.

preprint2011arXiv

Hagedorn States and Thermalization in Heavy Ion Collisions

In recent years, Hagedorn states have been used to explain the equilibrium and transport properties of a hadron gas close to the QCD critical temperature. These massive resonances are shown to lower $η/s$ to near the AdS/CFT limit close to the phase transition. A comparison of the Hagedorn model to recent lattice results is made and it is found that the hadrons can reach chemical equilibrium almost immediately, well before the chemical freeze-out temperatures found in thermal fits for a hadron gas without Hagedorn states.

preprint2011arXiv

Jet Quenching in Non-Conformal Holography

We use our non-conformal holographic bottom-up model for QCD described in 1012.0116 to further study the effect of the QCD trace anomaly on the energy loss of both light and heavy quarks in a strongly coupled plasma. We compute the nuclear modification factor $R_{AA}$ for bottom and charm quarks in an expanding plasma with Glauber initial conditions. We find that the maximum stopping distance of light quarks in a non-conformal plasma scales with the energy with a temperature (and energy) dependent effective power.

preprint2011arXiv

Non-conformal Holography of Heavy Quark Quenching

We develop a holographic (bottom-up) gravity model for QCD to understand the connection between the peak in the trace anomaly and the magnitude of heavy quark energy loss in a strongly-coupled plasma. The potential of the scalar field on the gravity side is constructed to reproduce some properties of QCD at finite temperature and its parameters are constrained by fitting lattice gauge theory results. The energy loss of heavy quarks (as predicted by the holographic model) is found to be strongly sensitive to the medium properties.

preprint2010arXiv

Connecting Polyakov Loops to the Thermodynamics of $SU(N_c)$ Gauge Theories Using the Gauge-String Duality

We show that in four-dimensional gauge theories dual to five-dimensional Einstein gravity coupled to a single scalar field in the bulk the derivative of the single heavy quark free energy in the deconfined phase is $dF_{Q}(T)/dT \sim -1/c_s^2(T)$, where $c_s(T)$ is the speed of sound. This general result provides a direct link between the softest point in the equation of state of strongly-coupled plasmas and the deconfinement phase transition described by the expectation value of the Polyakov loop. We give an explicit example of a gravity dual with black hole solutions that can reproduce the lattice results for the expectation value of the Polyakov loop and the thermodynamics of SU(3) Yang-Mills theory in the (non-perturbative) temperature range between $T_c$ and $3T_c$.

preprint2010arXiv

Flavoring the Quark-Gluon Plasma with Charm

We demonstrate that a Quark-Gluon Plasma (QGP) with a dilute admixture of heavy quarks has, in general, a lower speed of sound than a "pure" QGP without effects from heavy flavors. The change in the speed of sound is sensitive to the details of the theory, making the hydrodynamic response to "flavoring" a sensitive probe of the underlying microscopic dynamics. We suggest that this effect may be measured in ultrarelativistic heavy ion collisions by relating the event-by-event number of charm quarks to flow observables such as the average transverse momentum.

preprint2010arXiv

Heavy quarks and the collective properties of hot QCD

After reviewing the evidence that QCD matter at ultrarelativistic energies behaves as a very good fluid, we describe the connection of QCD fluidity to heavy quark observables. We review the way in which heavy quark spectra can place tighter limits on the viscosity of QCD matter. Finally, we show that correlations between flow observables and the event-by-event charm quark abundance ("flavoring") can shed light on the system's equation of state.

preprint2010arXiv

The Heavy Quark Free Energy in QCD and in Gauge Theories with Gravity Duals

We show that the Polyakov loop approaches its asymptotic value at high temperatures {\it from below} in any strongly-coupled gauge theory at large $N_c$ with a dual supergravity description. In these theories the Polyakov loop can be used to define a quantity that behaves as a free energy of a heavy quark in the medium. This should be compared to the recent lattice results in pure glue SU(3) theory at high temperatures where the Polyakov loop approaches its non-interacting limit {\it from above}. We show that the "heavy quark free energy" obtained from the renormalized loop computed on the lattice does not behave as a true thermodynamic free energy. We conjecture that this should always be the case in asymptotically free gauge theories.

preprint2010arXiv

Transport Properties of QCD at Large $N_c$ and the Gauge/String Duality

Below the deconfinement phase transition large $N_c$ QCD is expected to be a very viscous hadronic fluid because both the shear and bulk viscosity to entropy density ratio, $η/s,ζ/s \sim N_c^2$. In this letter I show that $η/s \sim N_c^2$ in the confined phase of holographic models of QCD at large $N_c$ defined in the supergravity approximation. Our results show that the gauge/string duality can be used to describe not only nearly perfect fluids but also extremely viscous systems such as a cold gas of glueballs.

preprint2010arXiv

Universal Flow-Driven Conical Emission in Ultrarelativistic Heavy-Ion Collisions

The double-peak structure observed in soft-hard hadron correlations is commonly interpreted as a signature for a Mach cone generated by a supersonic jet interacting with the hot and dense medium created in ultrarelativistic heavy-ion collisions. We show that it can also arise due to averaging over many jet events in a transversally expanding background. We find that the jet-induced away-side yield does not depend on the details of the energy-momentum deposition in the plasma, the jet velocity, or the system size. Our claim can be experimentally tested by comparing soft-hard correlations induced by heavy-flavor jets with those generated by light-flavor jets.

preprint2009arXiv

Anomalous Conical Di-jet Correlations in pQCD vs AdS/CFT

We propose an identified heavy quark jet observable to discriminate between weakly coupled pQCD and strongly coupled AdS/CFT models of quark gluon plasma dynamics in ultra-relativistic nuclear collisions at RHIC and LHC energies. These models are shown to predict qualitatively different associated hadron correlations with respect to tagged heavy quark jets. While both models feature similar far zone Mach and diffusion wakes, the far zone stress features are shown to be too weak to survive thermal broadening at hadron freeze-out. However, these models differ significantly in a near zone "Neck" region where strong chromo-fields sourced by the heavy quark jet couple to the polarizable plasma. Conical associated correlations, if any, are shown to be dominated by the jet induced transverse flow in the Neck zone and unrelated to the weak far zone wakes. Unlike in AdS/CFT, we show that the induced transverse flow in the Neck zone is too weak in pQCD to produce conical correlations after Cooper-Frye freeze-out. The observation of conical correlations violating Mach's law would favor the strongly-coupled AdS/CFT string drag dynamics, while their absence would favor weakly-coupled pQCD-based chromo-hydrodynamics.

preprint2009arXiv

Constraints on AdS/CFT Gravity Dual Models of Heavy Ion Collisions

We show that the five-fold constraints due to (1) the observed nuclear modification of heavy quark jets measured via non-photonic electrons $R_{AA}^e(p_T \sim 6 {\rm GeV})$ in central Au+Au collisions at 200 AGeV, (2) the "perfect fluid' elliptic transverse flow of low transverse momenta pions, $v_2(p_T\sim 1 {\rm GeV})$ reported for noncentral collisions, (3) the global pion rapidity density $dN_π/dy$, (4) the lattice QCD entropy density deficiency, $S/S_{SB}$, of strongly coupled Quark-Gluon Plasmas (sQGP), and (5) a causal requirement are analytically correlated in a class of gauge/string gravity dual models of sQGP dynamics. Current RHIC/BNL and lattice QCD data are found to be remarkably compatible with these models if the t'Hooft and Gauss-Bonnet coupling parameters lie in the range $λ\approx 10-25$ and $0< λ_{GB}< 0.09$. In addition, the observed five-fold correlation appears to favor color glass condensate over Glauber initial conditions within current systematic errors.

preprint2009arXiv

Thermal Width of the $Υ$ at Large t' Hooft Coupling

We use the AdS/CFT correspondence to show that the heavy quark (static) potential in a strongly-coupled plasma develops an imaginary part at finite temperature. Thus, deeply bound heavy quarkonia states acquire a small nonzero thermal width when the t'Hooft coupling $λ=g^2 N_c \gg 1$ and the number of colors $N_c \to \infty$. In the dual gravity description, this imaginary contribution comes from thermal fluctuations around the bottom of the classical sagging string in the bulk that connects the heavy quarks located at the boundary. We predict a strong suppression of $Υ$'s in heavy-ion collisions and discuss how this may be used to estimate the initial temperature.

preprint2009arXiv

Universality of the Diffusion Wake from Stopped and Punch-Through Jets in Heavy-Ion Collisions

We solve (3+1)-dimensional ideal hydrodynamical equations with source terms that describe punch-through and fully stopped jets in order to compare their final away-side angular correlations in a static medium. For fully stopped jets, the backreaction of the medium is described by a simple Bethe-Bloch-like model which leads to an explosive burst of energy and momentum (Bragg peak) close to the end of the jet's evolution through the medium. Surprisingly enough, we find that the medium's response and the corresponding away-side angular correlations are largely insensitive to whether the jet punches through or stops inside the medium. This result is also independent of whether momentum deposition is longitudinal (as generally occurs in pQCD energy loss models) or transverse (as the Bethe-Bloch formula implies). The existence of the diffusion wake is therefore shown to be universal to all scenarios where momentum as well as energy is deposited into the medium, which can readily be understood in ideal hydrodynamics through vorticity conservation. The particle yield coming from the strong forward moving diffusion wake that is formed in the wake of both punch-through and stopped jets largely overwhelms their weak Mach cone signal after freeze-out.