Source author record

Jacquelyn Noronha-Hostler

Jacquelyn Noronha-Hostler 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

43works
8topics
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

43 published item(s)

preprint2022arXiv

Initializing BSQ with Open-Source ICCING

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

preprint2022arXiv

Mirror Neutron Stars: How QCD can be used to study dark matter through gravitational waves

Given the lack of empirical evidence of weakly interacting dark matter, it is reasonable to look to other candidates such as a confining dark sector with a similar number of particles as the standard model. Twin Higgs mirror matter is one such model that is a twin of the standard model with particles masses 3--6 times heavier than the standard model that solves the hierarchy problem. This generically predicts mirror neutron stars, degenerate objects made entirely of mirror nuclear matter. We find their structure using a realistic equation of state from crust (nuclei) to core (relativistic mean-field model) and scale the particle masses using lattice QCD results. We find that mirror neutron stars have unique signatures that are detectable via gravitational waves and binary pulsars, that provides an intriguing possibility for probing dark matter.

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

Thermal-model-based characterization of heavy-ion-collision systems at chemical freeze-out

We investigate the chemical freeze-out in heavy-ion collisions (HICs) and the impact of the hadronic spectrum on thermal model analyses. Detailed knowledge of the hadronic spectrum is still an open question, which has phenomenological consequences on the study of HICs. By varying the number of resonances included in Hadron Resonance Gas (HRG) Model calculations, we can shed light on which particles may be produced. Furthermore, we study the influence of the number of states on the so-called two flavor freeze-out scenario, in which strange and light particles can freeze-out separately. We consider results for the chemical freeze-out parameters obtained from thermal model fits and from calculating net-particle fluctuations. We will show the effect of using one global temperature to fit all particles and alternatively, allowing particles with and without strange quarks to freeze-out separately.

preprint2021arXiv

Building a testable shear viscosity across the QCD phase diagram

Current experiments at the Relativistic Heavy Ion Collider (RHIC) are probing finite baryon densities where the shear viscosity to enthalpy ratio $ηT/w$ of the Quark Gluon Plasma remains unknown. We use the Hadron Resonance Gas (HRG) model with the most up-to-date hadron list to calculate $ηT/w$ at low temperatures and at finite baryon densities $ρ_B$. We then match $ηT/w$ to a QCD-based shear viscosity calculation within the deconfined phase to create a table across $\left\{T,μ_B\right\}$ for different cross-over and critical point scenarios at a specified location. We find that these new $ηT/w(T,μ_B)$ values would require initial conditions at significantly larger $ρ_B$, compared to ideal hydrodynamic trajectories, in order to reach the same freeze-out point.

preprint2021arXiv

Correlation between mean transverse momentum and anisotropic flow in heavy-ion collisions

The correlation between the mean transverse momentum of outgoing particles, $\langle p_t \rangle$, and the magnitude of anisotropic flow, $v_n$, has recently been measured in Pb+Pb collisions at the CERN Large Hadron Collider, as a function of the collision centrality. We confirm the previous observation that event-by-event hydrodynamics predicts a correlation between $v_n$ and $\langle p_t \rangle$ that is similar to that measured in data. We show that the magnitude of this correlation can be directly predicted from the initial condition of the hydrodynamic calculation, for $n=2,3$, if one replaces $v_n$ by the corresponding initial-state anisotropy, $\varepsilon_n$, and $\langle p_t\rangle$ by the total energy per unit rapidity of the fluid at the beginning of the hydrodynamic expansion.

preprint2021arXiv

Extreme Matter meets Extreme Gravity: Ultra-heavy neutron stars with crossovers and first-order phase transitions

The speed of sound of the matter within neutron stars may contain non-smooth structure related to first-order phase transitions or or crossovers. Here we investigate what are the observable consequences of structure, such as bumps, spikes, step functions, plateaus, and kinks. One of the main consequences is the possibility of ultra-heavy neutron stars, i.e.~stars with masses significantly heavier than two solar masses. These stars pass all observational and theoretical constraints, including those imposed by recent LIGO/Virgo gravitational-wave observations and NICER X-ray observations. We thoroughly investigate other consequences of this structure in the speed of sound to develop an understanding of how non-smooth features affect astrophysical observables, such as stellar radii, tidal deformability, moment of inertia, and Love number. Our results have important implications to future gravitational wave and X-ray observations of neutron stars and their impact in nuclear astrophysics.

preprint2021arXiv

Impact of Multiplicity Fluctuations on Entropy Scaling Across System Size

The initial state is one of the greatest uncertainties in heavy-ion collisions. A model-agnostic approach is taken in the phenomenological Trento framework which constrains parameters using Bayesian analysis. However, the color-glass condensate (CGC) effective theory predicts initial energy densities that lie outside the recent Bayesian analyses due, in part, to the assumption in Trento of event-by-event multiplicity fluctuations following a $Γ$ distribution. We compare the Trento-preferred $\sqrt{T_{A}T_{B}}$ scaling to CGC-like $T_{A}T_{B}$ scaling coupled with log-normal fluctuations in $AuAu$ and $dAu$ collisions and find there is a significant impact on the multiplicity distributions and on the eccentricities, which may affect the extraction of viscosity in small systems.

preprint2021arXiv

Off-of-equilibrium effects on Kurtosis Along Strangeness-Neutral Trajectories

The Beam Energy Scan program at the Relativistic Heavy Ion Collider (RHIC) is searching for the QCD critical point. The main signal for the critical point is the kurtosis of the distribution of proton yields obtained on an event-by-event basis where one expects a peak at the critical point. However, its exact behavior is still an open question due to out-of-equilibrium effects and uncertainty in the equation of state. Here we use a simplistic hydrodynamic model that enforces strangeness-neutrality, selecting trajectories that pass close to the critical point. We vary the initial conditions to estimate the effect of out-of-equilibrium hydrodynamics on the kurtosis signal.

preprint2021arXiv

Shear viscosity at finite baryon densities

We use the excluded volume Hadron Resonance Gas (HRG) model with the most up-to-date hadron list to calculate $ηT/w$ at low temperatures and at finite baryon densities $ρ_B$. This $ηT/w$ is then matched to a QCD-based shear viscosity calculation of the QGP for different profiles of $ηT/w$ across {T,$μ_{B}$} including cross-over and critical point transitions. When compared to ideal hydrodynamic trajectories across {T,$μ_{B}$}, we find that the $ηT/w(T,μ_B)$ profiles would require initial conditions at much larger baryon density to reach the same freeze-out point.

preprint2021arXiv

Skewness of mean transverse momentum fluctuations in heavy-ion collisions

We propose the skewness of mean transverse momentum, $\langle p_t \rangle$, fluctuations as a fine probe of hydrodynamic behavior in relativistic nuclear collisions. We describe how the skewness of the $\langle p_t \rangle$ distribution can be analyzed experimentally, and we use hydrodynamic simulations to predict its value. We predict in particular that $\langle p_t \rangle$ fluctuations have positive skew, which is significantly larger than if particles were emitted independently. We elucidate the origin of this result by deriving generic formulas relating the fluctuations of $\langle p_t \rangle$ to the fluctuations of the early-time thermodynamic quantities. We postulate that the large positive skewness of $\langle p_t \rangle$ fluctuations is a generic prediction of hydrodynamic models.

preprint2021arXiv

Smooth matching of $\hat{q}$ from hadronic to quark and gluon degrees of freedom

One of the key signatures of the Quark Gluon Plasma (QGP) is the energy loss of high momentum particles as they traverse the strongly interacting medium. The energy loss of these particles is governed by the jet transport coefficient $\hat{q}/T^3$, wherein it has been thought that there is a large jump as the system transitions between the hadron gas and Quark Gluon Plasma phases. Here we calculate $\hat{q}/T^3$ within the Hadron Resonance Gas (HRG) model with the particle list PDG16+ and find that, if one incorporates the experimental error in the hadronic calculation of $\hat{q}/T^3$ and assumes a higher pseudo-critical temperature, then a smooth transition from the hadron gas phase into the Quark Gluon Plasma phase is possible. We also find a significant enhancement in $\hat{q}/T^3$ at finite baryon chemical potential and find issues with the relationship between the shear viscosity and the jet transport coefficient within a hadron gas phase.

preprint2021arXiv

The slope, the hill, the drop, and the swoosh: Learning about the nuclear matter equation of state from the binary Love relations

Analyses that connect astrophysical observations of neutron stars with nuclear matter properties sometimes rely on equation-of-state insensitive relations. We show that the slope of the binary Love relations (i.e.~between the tidal deformabilities of binary neutron stars) encodes the rate of change of the nuclear matter speed of sound below three times nuclear saturation density. Twin stars lead to relations that present a signature ''hill'', ''drop'', and ''swoosh'' due to the second (mass-radius) stable branch, requiring a new description of the binary love relations.

preprint2020arXiv

Chemical freeze-out parameters of net-kaons in heavy-ion collisions

We study chemical freeze-out parameters for heavy-ion collisions by performing two different thermal analyses. We analyze results from thermal fits for particle yields, as well as, net-charge fluctuations in order to characterize the chemical freeze-out. The Hadron Resonance Gas (HRG) model is employed for both methods. By separating the light hadrons from the strange hadrons in thermal fits, we study the proposed flavor hierarchy. For the net-charge fluctuations, we calculate the mean-over-variance ratio of the net-kaon fluctuations in the HRG model at the five highest energies of the RHIC Beam Energy Scan (BES) for different particle data lists. We compare these results with recent experimental data from the STAR collaboration in order to extract sets of chemical freeze-out parameters for each list. We focused on particle lists which differ largely in the number of resonant states. By doing so, our analysis determines the effect of the amount of resonances included in the HRG model on the freeze-out conditions. Our findings have potential impact on various other models in the field of relativistic heavy-ion collisions.

preprint2020arXiv

Importance of Multiplicity Fluctuations in Entropy Scaling

One of the greatest uncertainties in heavy-ion collisions is the description of the initial state. Different models predict a wide range of initial energy density distributions based on their underlying assumptions. Final flow harmonics are sensitive to these differences in the initial state due to the nearly linear mapping between eccentricities and anisotropic flow harmonics. The Trento code uses a model-agnostic approach by phenomenologically parameterizing the initial state and constraining those parameters from a Bayesian analysis. There the multiplicity fluctuations were determined by a one parameter $Γ$ distribution. However, initial-state models arising from the Color-Glass Condensate (CGC) framework lead to an initial energy density which is outside the functional form considered in Trento and its later Bayesian analyses because they rely on log-normal multiplicity fluctuations. We compare $T_{A}T_{B}$ scaling (CGC-like) to $\sqrt{T_{A}T_{B}}$ scaling (preferred from a Trento Bayesian analysis) and find that the $T_A T_B$ form together with log-normal fluctuations is a reasonable candidate to describe the multiplicity fluctuations but leads to larger eccentricities, which would affect the extraction of viscosity in small systems.

preprint2020arXiv

Neutron Star Equation of State in light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing non-trivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

preprint2020arXiv

New opportunities at the photon energy frontier

Ultra-peripheral collisions (UPCs) involving heavy ions and protons are the energy frontier for photon-mediated interactions. UPC photons can be used for many purposes, including probing low-$x$ gluons via photoproduction of dijets and vector mesons, probes of beyond-standard-model processes, such as those enabled by light-by-light scattering, and studies of two-photon production of the Higgs.

preprint2020arXiv

Off-diagonal correlators of conserved charges from lattice QCD and how to relate them to experiment

Like fluctuations, non-diagonal correlators of conserved charges provide a tool for the study of chemical freeze-out in heavy ion collisions. They can be calculated in thermal equilibrium using lattice simulations, and be connected to moments of event-by-event net-particle multiplicity distributions. We calculate them from continuum extrapolated lattice simulations at $μ_B=0$, and present a finite-$μ_B$ extrapolation, comparing two different methods. In order to relate the grand canonical observables to the experimentally available net-particle fluctuations and correlations, we perform a Hadron Resonance Gas (HRG) model analysis, which allows us to completely break down the contributions from different hadrons. We then construct suitable hadronic proxies for fluctuations ratios, and study their behavior at finite chemical potentials. We also study the effect of introducing acceptance cuts, and argue that the small dependence of certain ratios on the latter allows for a direct comparison with lattice QCD results, provided that the same cuts are applied to all hadronic species. Finally, we perform a comparison for the constructed quantities for experimentally available measurements from the STAR Collaboration. Thus, we estimate the chemical freeze-out temperature to 165 MeV using a strangeness-related proxy. This is a rather high temperature for the use of the Hadron Resonance Gas, thus, further lattice studies are necessary to provide first principle results at intermediate $μ_B$.

preprint2020arXiv

QCD equation of state matched to lattice data and exhibiting a critical point singularity

We construct a family of equations of state for QCD in the temperature range 30 MeV $\leq T\leq$ 800 MeV and in the chemical potential range $0\leq μ_B \leq$ 450 MeV. These equations of state match available lattice QCD results up to $\mathcal{O}(μ_B^4)$ and in each of them we place a critical point in the 3D Ising model universality class. The position of this critical point can be chosen in the range of chemical potentials covered by the second Beam Energy Scan at RHIC. We discuss possible choices for the free parameters, which arise from mapping the Ising model onto QCD. Our results for the pressure, entropy density, baryon density, energy density and speed of sound can be used as inputs in the hydrodynamical simulations of the fireball created in heavy ion collisions. We also show our result for the second cumulant of the baryon number in thermal equilibrium, displaying its divergence at the critical point. In the future, comparisons between RHIC data and the output of the hydrodynamic simulations, including calculations of fluctuation observables, built upon the model equations of state that we have constructed may be used to locate the critical point in the QCD phase diagram, if there is one to be found.

preprint2019arXiv

D meson sensitivity to a system size scan at LHC

Experimental measurements in pA collisions indicate no D meson suppression ($R_{\rm pPb} \sim 1$) but a surprisingly large $v_2$. To better understand these results we propose a system size scan at the LHC involving $^{16}$OO, $^{40}$ArAr, $^{129}$XeXe and $^{208}$PbPb collisions. Using Trento+ v-USPhydro+DAB-MOD to make predictions, we find that the $R_{\rm AA}$ tends towards unity when the system size is decreased, but nonetheless, in the most central collisions $v_2\{2\}$ is almost independent of the colliding system. These results are analyzed in light of path length and initial eccentricity variations.

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.

preprint2019arXiv

Determination of Chemical Freeze-out Parameters from Net-kaon Fluctuations at RHIC

We calculate the mean-over-variance ratio of the net-kaon fluctuations in the Hadron Resonance Gas (HRG) Model for the five highest energies of the RHIC Beam Energy Scan (BES) for different particle data lists. We compare these results with the latest experimental data from the STAR collaboration in order to extract sets of chemical freeze-out parameters for each list. We focused on the PDG2012 and PDG2016+ particle lists, which differ largely in the number of resonant states. Our analysis determines the effect of the amount of resonances included in the HRG on the freeze-out conditions.

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

Hydrodynamic predictions for 5.02 TeV Pb-Pb collisions

We make predictions for momentum-integrated elliptic and triangular flow as well as mean transverse momentum for 5.02 TeV Pb-Pb collisions, as planned at the Large Hadron Collider. We use hydrodynamic calculations to predict the change of these observables as the center-of-mass collision energy evolves from 2.76 TeV to 5.02 TeV per nucleon pair. By using previously measured values as a baseline, we are able to make a robust prediction without relying on a particular model for initial conditions and without precise knowledge of medium properties such as viscosity. Thus, though the predicted changes are small, they can provide a significant test of the current hydrodynamic picture of heavy-ion collisions.

preprint2016arXiv

Implications of Missing Resonances in Heavy Ions Collisions

Between 2004 until today many new resonances in the QCD spectrum have been added by the Particle Data Group. However, it is still possible that there are further states that have yet to be measured, as predicted by Lattice QCD calculations and Quark Models. In this paper I review how these missing states would influence the field of heavy-ion collisions. Additionally, the quantum numbers that characterize these missing states are also discussed.

preprint2016arXiv

Kaon fluctuations from lattice QCD

We show that it is possible to isolate a set of kaon fluctuations in lattice QCD. By means of the Hadron Resonance Gas (HRG) model, we calculate the actual kaon second-to-first fluctuation ratio, which receives contribution from primordial kaons and resonance decays, and show that it is very close to the one obtained for primordial kaons in the Boltzmann approximation. The latter only involves the strangeness and electric charge chemical potentials, which are functions of $T$ and $μ_B$ due to the experimental constraint on strangeness and electric charge, and can therefore be calculated on the lattice. This provides an unambiguous method to extract the kaon freeze-out temperature, by comparing the lattice results to the experimental values for the corresponding fluctuations.

preprint2016arXiv

Linear and cubic response to the initial eccentricity in heavy-ion collisions

We study the relation between elliptic flow, $v_2$ and the initial eccentricity, $\varepsilon_2$, in heavy-ion collisions, using hydrodynamic simulations. Significant deviations from linear eccentricity scaling are seen in more peripheral collisions. We identify the mechanism responsible for these deviations as a cubic response, which we argue is a generic property of the hydrodynamic response to the initial density profile. The cubic response increases elliptic flow fluctuations, thereby improving agreement of initial condition models with experimental data.

preprint2016arXiv

Resolving the $R_{AA}$ to $v_n$ puzzle

After 10 years of struggling to simultaneously describe the nuclear modification factor $R_{AA}$ and flow harmonics $v_n$'s at high $p_T$, now theoretical models are able to reproduce experimental data well. The necessary theoretical developments such as event-by-event fluctuations, choice of initial conditions, and the scalar product method to calculate flow harmonics at high $p_T$ are reviewed. Additionally, a discussion of new proposed experimental observables known as Soft Hard Event Engineering (SHEE) that are sensitive to the path length dependence of the energy loss is included.

preprint2016arXiv

Solving the $R_{AA}\otimes v_2$ puzzle

For the past ten years $R_{AA}(p_T)$, the nuclear modification factor that encodes the suppression of high $p_T$ particles due to energy loss within the medium was fairly well described by many theoretical models. However, the same models systematically under-predicted the high $p_T$ elliptic flow, $v_2$, which is experimentally measured as the correlation between soft and hard hadrons. All previous calculations neglected the effect of event-by-event fluctuations of an expanding viscous hydrodynamical background as well as the soft-hard flow harmonic correlations in the experimentally measured $v_2$. In this talk I show how event-by-event viscous hydrodynamics (computed using the v-USPhydro code) coupled to an energy loss model (BBMG) is able to simultaneously describe soft physics observables as well as the high-$p_T$ $R_{AA}$ and $v_2$. Suggestions for future more differential calculations at the LHC run2 are made to explore soft-hard flow correlations.

preprint2016arXiv

Symmetric cumulants and event-plane correlations in Pb+Pb collisions

The ALICE Collaboration has recently measured the correlations between amplitudes of anisotropic flow in different Fourier harmonics, referred to as symmetric cumulants. We derive approximate relations between symmetric cumulants involving $v_4$ and $v_5$ and the event-plane correlations measured by ATLAS. The validity of these relations is tested using event-by-event hydrodynamic calculations. The corresponding results are in better agreement with ALICE data than existing hydrodynamic predictions. We make quantitative predictions for three symmetric cumulants which are not yet measured.

preprint2015arXiv

Extracting shear viscosity of the Quark Gluon Plasma in the presence of bulk viscosity

One of the most remarkable features of the Quark Gluon Plasma is its nearly perfect fluidity behavior indicated by the small shear viscosity to entropy density ratio obtained from fitting relativistic viscous hydrodynamics flow harmonics to experimental data. In recent years, bulk viscosity has also been considered in the context of event-by-event relativistic hydrodynamics and it has been found to have a non-trivial interplay with shear viscosity. In this paper some of the issues are discussed that require further work when extracting the shear viscosity to entropy density ratio in the presence of a non-zero bulk viscosity.

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.

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

Suppression of the LHC $p/π$ ratio due to the QCD mass spectrum

Recent measurements of the proton to pion ratio at $\sqrt{s}_{NN}=2.7$ TeV by the ALICE collaboration at the LHC have found it to be lower than predictions from thermal fits. In this paper we investigate the role that the extended mass spectrum via Hagedorn states- massive resonances that follow an exponential mass spectrum and possess large decay widths- play in the determination of particle ratios at LHC through a scenario of multi-particle reactions and dynamical chemical equilibrium within the hadron gas phase. We show that it is possible to describe the lower $p/π$ ratio at LHC while still obtaining the experimental ratio of $K/π$ and $Λ/π^+$ in the Hagedorn state scenario if the protons are underpopulated at the switching temperature from hydrodynamics to the hadron gas phase.

preprint2014arXiv

Understanding the $p/π$ ratio at LHC due to QCD mass spectrum

Thermal fits have consistently reproduced the experimental particles yields of heavy ion collisions, however, the proton to pion ratio from ALICE Pb+Pb $\sqrt{s_{NN}}=2.76$ TeV is over-predicted by thermal models- known at the $p/π$ puzzle. Here we test the relevance of the extended mass spectrum, i.e., include Hagedorn states (resonances that follow an exponential mass spectrum and have very short life times) on the $p/π$ puzzle. We find that the extended mass spectrum is able to reproduce particle ratios at both RHIC and the LHC as well as being able to match the lower $p/π$ ratio at the LHC through dynamical chemical equilibration.

preprint2014arXiv

Viscous Effects on the Mapping of the Initial to Final State in Heavy Ion Collisions

We investigate the correlation between various aspects of the initial geometry of heavy ion collisions at the Relativistic Heavy Ion Collider energies and the final anisotropic flow, using v-USPhydro, a 2+1 event-by-event viscous relativistic hydrodynamical model. We test the extent of which shear and bulk viscosity affect the prediction of the final flow harmonics, $v_n$, from the initial eccentricities, $\varepsilon_n$. We investigate in detail the flow harmonics $v_1$ through $v_5$ where we find that $v_1$, $v_4$, and $v_5$ are dependent on more complicated aspects of the initial geometry that are especially important for the description of peripheral collisions, including a non-linear dependence on eccentricities as well as a dependence on shorter-scale features of the initial density. Furthermore, we compare our results to previous results from NeXSPheRIO, a 3+1 relativistic ideal hydrodynamical model that has a non-zero initial flow contribution, and find that the combined contribution from 3+1 dynamics and non-zero, fluctuating initial flow decreases the predictive ability of the initial eccentricities, in particular for very peripheral collisions, but also disproportionately in central collisions.

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π)$.

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.

preprint2010arXiv

Hagedorn States and Thermalization

In recent years Hagedorn states have been used to explain the physics close to the critical temperature within a hadron gas. Because of their large decay widths these massive resonances lower $η/s$ to near the AdS/CFT limit within the hadron gas phase. A comparison of the Hagedorn model to recent lattice results is made and it is found that for both Tc =176 MeV and Tc=196 MeV, 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. In this paper we also observe the effects of Hagedorn States on the $K^+/π^+$ horn seen at AGS, SPS, and RHIC.