Search

QTRAJ is a computer code that simulates the propagation of quarkonium in the quark-gluon plasma (QGP) based on the quantum-trajectory algorithm. This algorithm solves a master equation in which the quarkonium is treated as an open quantum system (OQS). A major advantage of this approach is that it turns a 3D spatial evolution for a density matrix into a 1D Schr\"odinger equation for a wavefunction with a non-hermitian Hamiltonian, drastically reducing the computational cost. So far, the interaction implemented in the master equation was obtained within the framework of potential non-relativistic QCD (pNRQCD), and restricted to the regime $rT \ll 1$, where $r$ is the size of the color dipole and $T$ is the temperature. In the environment produced in heavy-ion collisions (HIC's) this limit is accurate for $\Upsilon(1S)$, but the applicability to other quarkonium states is dubious. In the present study we generalize the above approach, extending it to the regime $rT\!\sim\! 1$ in the one-gluon exchange approximation, with proper Hard Thermal Loop (HTL) resummation of medium effects. This is done by implementing new jump operators connecting different color states of the $Q\bar Q$ pair and expanding them in plane waves, giving rise to a variation of the algorithm present in QTRAJ 1.0. Here we provide an overview of this approach comparing the $rT \ll 1$ and $rT\sim 1$ cases, and we discuss prospects for phenomenological application to excited states of bottomonium., Comment: Proceedings of the Conference QNP 2024

The QCD cross-over line in the temperature ($T$) -- baryo-chemical potential ($\mu_B$) plane has been computed by several lattice groups by calculating the chiral order parameter and its susceptibility at finite values of $\mu_B$. In this work we focus on the deconfinement aspect of the transition between hadronic and Quark Gluon Plasma (QGP) phases. We define the deconfinement temperature as the peak position of the static quark entropy ($S_Q(T,\mu_B)$) in $T$, which is based on the renormalized Polyakov loop. We extrapolate $S_Q(T,\mu_B)$ based on high statistics finite temperature ensembles on a $16^3\times 8$ lattice to finite density by means of a Taylor expansion to eighth order in $\mu_B$ (NNNLO) along the strangeness neutral line. For the simulations the 4HEX staggered action was used with 2+1 flavors at physical quark masses. In this setup the phase diagram can be drawn up to unprecedentedly high chemical potentials. Our results for the deconfinement temperature are in rough agreement with phenomenological estimates of the freeze-out curve in relativistic heavy ion collisions. In addition, we study the width of the deconfinement crossover. We show that up to $\mu_B \approx 400$ MeV, the deconfinement transition gets broader at higher densities, disfavoring the existence of a deconfinement critical endpoint in this range. Finally, we examine the transition line without the strangeness neutrality condition and observe a hint for the narrowing of the crossover towards large $\mu_B$., Comment: 14 pages, 8 figures

The crossover from hadronic to quark matter is understood to be both a deconfinement as well as a chiral symmetry restoring transition. Here, we study observables related to both aspects using lattice simulations: the Polyakov loop and its derivatives and the chiral condensate and its derivatives. At zero baryochemical potential, and infinite volume, the chiral and deconfinement crossover temperatures almost agree. However, chiral and deconfinement related observables have a qualitatively different chemical potential and volume dependence. In general, deconfinement related observables have a milder volume dependence. Furthermore, while the deconfinement transition appears to get broader with increasing $\mu_B$, the width as well as the strength of the chiral transition is approximately constant. Our results are based on simulations at zero and imaginary chemical potentials using 4stout-improved staggered fermions with $N_\tau=12$ time-slices and physical quark masses., Comment: 10 pages, 8 figures

An efficient way to study the QCD phase diagram at small finite density is to extrapolate thermodynamical observables from imaginary chemical potential. The phase diagram features a crossover line starting from the transition temperature already determined at zero chemical potential. In this work we focus on the Taylor expansion of this line up to $\mu^4$ contributions. We present the continuum extrapolation of the crossover temperature based on different observables at several lattice spacings., Comment: Proceedings to Quark Matter Conference 2019

We present a novel construction of the QCD equation of state (EoS) at finite baryon density. Our work combines a recently proposed resummation scheme for lattice QCD results with the universal critical behavior at the QCD critical point. This allows us to obtain a family of equations of state in the range $0 \leq \mu_B \leq 700$ MeV and 25 MeV $\leq T \leq 800$ MeV, which match lattice QCD results near $\mu_B=0$ while featuring a critical point in the 3D Ising model universality class. The position of the critical point can be chosen within the range accessible to beam-energy scan heavy-ion collision experiments. The strength of the singularity and the shape of the critical region are parameterized using a standard parameter set. We impose stability and causality constraints and discuss the available ranges of critical point parameter choices, finding that they extend beyond earlier parametric QCD EoS proposals. We present thermodynamic observables, including baryon density, pressure, entropy density, energy density, baryon susceptibility and speed of sound, that cover a wide range in the QCD phase diagram relevant for experimental exploration., Comment: 16 pages, 18 figures

The effect of a finite volume presents itself both in heavy ion experiments as well as in recent model calculations. The magnitude is sensitive to the proximity of a nearby critical point. We calculate the finite volume effects at finite temperature in continuum QCD using lattice simulations and set the focus on the vicinity of the chiral crossover. We investigate the impact of finite volumes at zero and small chemical potentials on the QCD transition through the chiral observables.

The BEST Collaboration equation of state combining lattice data with the 3D Ising critical point encounters limitations due to the truncated Taylor expansion up to $\frac{\mu_B}{T} \sim 2.5$. This truncation consequently restricts its applicability at high densities. Through a resummation scheme, the lattice results have been extended to $\frac{\mu_B}{T} = 3.5$. In this article, we amalgamate these ideas with the 3D-Ising model, yielding a family of equations of state valid up to $\mu_B=700 \text{MeV}$ with the correct critical behavior. Our equations of state feature tunable parameters, providing a stable and causal framework-a crucial tool for hydrodynamics simulations., Comment: 4 pages, 2 figures, Contribution to Quark Matter2023 (Houston, TX, 3-9 Sep. 2023)

Fluctuations play a key role in the study of QCD phases. Lattice QCD is a valuable tool to calculate them, but going to high orders is challenging. Up to the fourth order, continuum results are available since 2015. We present the first continuum results for sixth order baryon fluctuations for temperatures between $T=130 - 200$ MeV, and eighth order at $T=145$ MeV in a fixed volume. We show that for $T \leq 145$ MeV, relevant for criticality search, finite volume effects are under control. Our results are in sharp contrast with well known results in the literature obtained at finite lattice spacing., Comment: 5 pages, 2 figures (main text) + 5 pages, 7 figures (supplemental material)

Hadron lists based on experimental studies summarized by the Particle Data Group (PDG) are a crucial input for the equation of state and thermal models used in the study of strongly-interacting matter produced in heavy-ion collisions. Modeling of these strongly-interacting systems is carried out via hydrodynamical simulations, which are followed by hadronic transport codes that also require a hadronic list as input. To remain consistent throughout the different stages of modeling of a heavy-ion collision, the same hadron list with its corresponding decays must be used at each step. It has been shown that even the most uncertain states listed in the PDG from 2016 are required to reproduce partial pressures and susceptibilities from Lattice Quantum Chromodynamics with the hadronic list known as the PDG2016+. Here, we update the hadronic list for use in heavy-ion collision modeling by including the latest experimental information for all states listed in the Particle Data Booklet in 2021. We then compare our new list, called PDG2021+, to Lattice Quantum Chromodynamics results and find that it achieves even better agreement with the first principles calculations than the PDG2016+ list. Furthermore, we develop a novel scheme based on intermediate decay channels that allows for only binary decays, such that PDG2021+ will be compatible with the hadronic transport framework SMASH. Finally, we use these results to make comparisons to experimental data and discuss the impact on particle yields and spectra., Comment: 17 pages, 16 figures, 2 tables

Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document., Comment: 190 pages, 69 figures

QCD with infinite heavy quark masses exhibits a first-order thermal transition which is driven by the spontaneous breaking of the global $\mathcal{Z}_3$ center symmetry. We analyze the corresponding order parameter, namely the Polyakov loop and its moments, and show, with a rigorous finite size scaling, that in the continuum limit the transition is of first order. We show that the use of a parallel tempering algorithm can significantly reduce the large auto-correlation times which are mainly caused by the supercritical slowing down. As a result, we calculate the transition temperature $w_0 T_c$ with per-mill precision, and the latent heat, carrying out controlled continuum and infinite volume extrapolations., Comment: 8 pages, presentation at the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, University of Bonn, Germany

The first order transition between the confining and the center symmetry breaking phases of the SU(3) Yang-Mills theory is marked by discontinuities in various thermodynamics functions, such as the energy density or the value of the Polyakov loop. We investigate the non-analytical behaviour of the topological susceptibility and its higher cumulant around the transition temperature and make the connection to the curvature of the phase diagram in the $T-\theta$ plane and to the latent heat., Comment: 13 pages, 10 figures

The equation of state (EoS) of QCD is a crucial input for the modeling of heavy-ion-collision (HIC) and neutron-star-merger systems. Calculations of the fundamental theory of QCD, which could yield the true EoS, are hindered by the infamous Fermi sign problem which only allows direct simulations at zero or imaginary baryonic chemical potential. As a direct consequence, the current coverage of the QCD phase diagram by lattice simulations is limited. In these proceedings, two different equations of state based on first-principle lattice QCD (LQCD) calculations are discussed. The first is solely informed by the fundamental theory by utilizing all available diagonal and non-diagonal susceptibilities up to $\mathcal{O}(\mu_B^4)$ in order to reconstruct a full EoS at finite baryon number, electric charge and strangeness chemical potentials. For the second, we go beyond information from the lattice in order to explore the conjectured phase structure, not yet determined by LQCD methods, to assist the experimental HIC community in their search for the critical point. We incorporate critical behavior into this EoS by relying on the principle of universality classes, of which QCD belongs to the 3D Ising Model. This allows one to study the effects of a singularity on the thermodynamical quantities that make up the equation of state used for hydrodynamical simulations of HICs. Additionally, we ensure that these EoSs are valid for applications to HICs by enforcing conditions of strangeness neutrality and fixed charge-to-baryon-number ratio., Comment: Contribution to the 37th Winter Workshop on Nuclear Dynamics. arXiv admin note: text overlap with arXiv:2103.08146

In this work, we study the lensing effect of the QCD critical point on hydrodynamic trajectories, and its consequences on the net-proton kurtosis $\kappa_4$. Including critical behavior by means of the BEST Collaboration equation of state (EoS), we first consider a scenario in equilibrium, then compare with hydrodynamic 0+1D simulations with Bjorken expansion, including both shear and bulk viscous terms. We find that, both in and out-of-equilibrium, the size and shape of the critical region directly affect if the signal will survive through the dynamical evolution.

The Beam Energy Scan Theory (BEST) collaboration's equation of state (EoS) incorporates a 3D Ising model critical point into the Quantum Chromodynamics (QCD) equation of state from lattice simulations. However, it contains 4 free parameters related to the size and location of the critical region in the QCD phase diagram. Certain combinations of the free parameters lead to acausal or unstable realizations of the EoS that should not be considered. In this work, we use an active learning framework to rule out pathological EoS efficiently. We find that checking stability and causality for a small portion of the parameters' range is sufficient to construct algorithms that perform with $>$96% accuracy across the entire parameter space. Though in this work we focus on a specific case, our approach can be generalized to any EoS containing a parameter space-class correspondence., Comment: 14 pages, 9 figures

We calculate a resummed equation of state with lattice QCD simulations at imaginary chemical potentials. This work presents a generalization of the scheme introduced in 2102.06660 to the case of non-zero $\mu_S$, focusing on the line of strangeness neutrality. We present results up to $\mu_B/T \leq 3.5$ on the strangeness neutral line $\left\langle S \right\rangle = 0$ in the temperature range $130 \rm{MeV} \leq T \leq 280 \rm{MeV}$. We also extrapolate the finite baryon density equation of state to small non-zero values of the strangeness-to-baryon ratio $R=\left\langle S \right\rangle / \left\langle B \right\rangle$. We perform a continuum extrapolation using lattice simulations of the 4stout-improved staggered action with 8, 10, 12 and 16 timeslices., Comment: 15 pages, 12 figures; v2: contains ancillary files with tabulated data

We perform large scale simulations to characterize the transition in quenched QCD. It is shown by a rigorous finite size scaling that the transition is of first order. After this qualitative feature quantitative results are obtained with unprecedented precision: we calculate the transition temperature $w_0T_c$=0.25384(23), -- which is the first per-mill accurate result in QCD thermodynamics -- and the latent heat $\Delta E/T_c^4$=1.025(21)(27) in both cases carrying out controlled continuum and infinite volume extrapolations. As it is well known the cost of lattice simulations explodes in the vicinity of phase transitions, a phenomenon called critical slowing down for second order phase transitions and supercritical slowing down for first order phase transitions. We show that a generalization of the parallel tempering algorithm of Marinari and Parisi [Europhys. Lett. 19, 451 (1992)] originally for spin systems can efficiently overcome these difficulties even if the transition is of first order, like in the case of QCD without quarks, or with very heavy quarks. We also report on our investigations on the autocorrelation times and other details., Comment: 15 pages, 12 figures, PRD version, minor changes

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., Comment: 5 pages, 2 figures, 1 table, Proceedings of SQM 2021 - Strangeness in Quark Matter

QCD with heavy dynamical quarks exhibits a first order thermal transition which is driven by the spontaneous breaking of the global $\mathcal{Z}_3$ center symmetry. Decreasing the quark masses weakens the transition until the corresponding latent heat vanishes at the critical mass. We explore the heavy mass region with three flavors of staggered quarks and analyze the Polyakov loop and its moments in a finite volume scaling study. We calculate the heavy critical mass in the three flavor theory in the infinite volume limit for $N_t=8$., Comment: 8 pages, presentation at the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Massachusetts Institute of Technology, USA

The hadron resonance gas (HRG) model and its extensions are often used to describe the hadronic phase of strongly interacting matter. In our work we use lattice-QCD simulations with temporal extents of $N_\tau=8,10$ and $12$ to quantify corrections to the ideal HRG. Firstly, we determine a number of subleading fugacity expansion coefficients of the QCD free energy via a two-dimensional scan on the imaginary baryon number chemical potential ($\mu_B$) - strangeness chemical potential ($\mu_S$) plane. Using the aforementioned coefficients, we also extrapolate ratios of baryon number and strangeness fluctuations and correlations to finite chemical potentials via a truncated fugacity expansion. Our results extrapolated along the crossover line $T_\mathrm{c}(\mu_B)$ at strangeness neutrality are able to reproduce trends of experimental net-proton fluctuations measured by the STAR Collaboration., Comment: 10 pages, 4 figures, Contribution to the 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021

The equation of state of Quantum Chromodynamics (QCD) at finite density is currently known only in a limited range in the baryon chemical potential $\mu_B$. This is due to fundamental shortcomings of traditional methods such as Taylor expansion around $\mu_B=0$. In this contribution, we present an alternative scheme that displays substantially improved convergence over the Taylor expansion method. We calculate the alternative expansion coefficients in the continuum, and show our results for the thermodynamic observables up to $\mu_B/T\le3.5$., Comment: 8 pages, presentation at the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Massachusetts Institute of Technology, USA

We formulate a family of equations of state for Quantum Chromodynamics that exhibit critical features and obey the charge conservation conditions present in heavy-ion collisions (HICs). This construction utilizes the first-principle Lattice QCD (LQCD) results up to $\mathcal{O}(\mu_B^4)$ by matching the Taylor coefficients at each order. The criticality of the equation of state (EoS) is implemented based on the well-established principle of universal scaling behavior, where QCD belongs to the 3D Ising Model universality class. The critical point can be placed in a range of temperature and baryonic chemical potential relevant for the Beam Energy Scan II at RHIC. Furthermore, the strength of the critical features can be varied as well. This flexibility is embedded in four free parameters, that could potentially be constrained by the experimental data. We will discuss the features and versatility of our EoS, as well as present the relevant thermodynamic quantities which are important for hydrodynamical simulations of HICs., Comment: 6 pages, 4 figures, Proceedings of the International Conference on Critical Point and Onset of Deconfinement 2021

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.

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Finally, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world., Comment: 103 pages

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

We present a strangeness-neutral equation of state for QCD that exhibits critical behavior and matches lattice QCD results for the Taylor-expanded thermodynamic variables up to fourth-order in $\mu_B/T$. It is compatible with the SMASH hadronic transport approach and has a range of temperatures and baryonic chemical potentials relevant for phase II of the Beam Energy Scan at RHIC. We provide an updated version of the software BES-EoS, which produces an equation of state for QCD that includes a critical point in the 3D Ising model universality class. This new version also includes isentropic trejectories and the critical contribution to the correlation length. Since heavy-ion collisions have zero global net-strangeness density and a fixed ratio of electric charge to baryon number, the BES-EoS is more suitable to describe this system. Comparison with the previous version of the EoS is thoroughly discussed., Comment: 11 pages, 11 figures

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.

Current experiments at the Relativistic Heavy Ion Collider (RHIC) are probing finite baryon densities where the shear viscosity to enthalpy ratio $\eta 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 $\eta T/w$ at low temperatures and at finite baryon densities $\rho_B$. We then match $\eta T/w$ to a QCD-based shear viscosity calculation within the deconfined phase to create a table across $\left\{T,\mu_B\right\}$ for different cross-over and critical point scenarios at a specified location. We find that these new $\eta T/w(T,\mu_B)$ values would require initial conditions at significantly larger $\rho_B$, compared to ideal hydrodynamic trajectories, in order to reach the same freeze-out point.

Taylor expansion of the equation of state of QCD suffers from shortcomings at chemical potentials $\mu_B \geq (2-2.5)T$. First, one faces difficulties inherent in performing such an expansion with a limited number of coefficients; second, higher order coefficients determined from lattice calculations suffer from a poor signal-to-noise ratio. In this work, we present a novel scheme for extrapolating the equation of state of QCD to finite, real chemical potential that can extend its reach further than previous methods. We present continuum extrapolated lattice results for the new expansion coefficients and show the thermodynamic observables up to $\mu_B/T\le3.5$., Comment: 13 pages, 10 figures

The hadron resonance gas (HRG) model is often believed to correctly describe the confined phase of QCD. This assumption is the basis of many phenomenological works on QCD thermodynamics and of the analysis of hadron yields in relativistic heavy ion collisions. We use first-principle lattice simulations to calculate corrections to the ideal HRG. Namely, we determine the sub-leading fugacity expansion coefficients of the grand canonical free energy, receiving contributions from processes like kaon-kaon or baryon-baryon scattering. We achieve this goal by performing a two dimensional scan on the imaginary baryon number chemical potential ($\mu_B$) - strangeness chemical potential ($\mu_S$) plane, where the fugacity expansion coefficients become Fourier coefficients. We carry out a continuum limit estimation of these coefficients by performing lattice simulations with temporal extents of $N_\tau=8,10,12$ using the 4stout-improved staggered action. We then use the truncated fugacity expansion to extrapolate ratios of baryon number and strangeness fluctuations and correlations to finite chemical potentials. Evaluating the fugacity expansion along the crossover line, we reproduce the trend seen in the experimental data on net-proton fluctuations by the STAR collaboration., Comment: 13 pages, 6 figures

In the context of the ongoing search for the QCD critical point at the Relativistic Heavy-Ion Collider, we study the equation of state near the critical point in the temperature and baryon chemical potential plane. We use the parametric representation introduced in earlier literature, which maps the universal 3D Ising equation of state onto the QCD phase diagram using several non-universal parameters. We focus on the quartic cumulant of the baryon number, or baryon number susceptibility~$\chi_4^B$, which can be accessed experimentally via net-proton fluctuation kurtosis measurements. It was originally predicted, through universality arguments based on the {\em leading} singular contribution, that $\chi_4^B$ and net-proton kurtosis should show a specific non-monotonic behavior due to the critical point. In particular, when following the freeze-out curve on the phase diagram by decreasing beam energy, the kurtosis is expected to dip, and then peak, when the beam energy scan passes close to the critical point. We study the effects of the non-universal and thus far unknown parameters of the Ising-to-QCD mapping on the behavior of~$\chi_4^B$. We find that, while the peak remains a solid feature, the presence of the critical point does not necessarily cause a dip in $\chi_4^B$ on the freezeout line {\em below} the transition temperature. The critical point contribution to the dip appears only for a narrow set of mapping parameters, when subleading singular terms are sufficiently suppressed., Comment: 10 pages, 4 figures

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., Comment: 4 pages, 2 figures, 1 table, contribution to the Quark Matter 2019 proceedings

Detailed knowledge of the hadronic spectrum is still an open question, which has phenomenological consequences on the study of heavy-ion collisions. A previous lattice QCD study concluded that additional strange resonances are missing in the currently tabulated lists provided by the Particle Data Group (PDG). That study identified the list labeled PDG2016+ as the ideal spectrum to be used as an input in thermal-model-based analyses. In this work, we determine the effect of additional resonances on the freeze-out parameters of systems created in heavy-ion collisions. These parameters are obtained from thermal fits of particle yields and net-particle fluctuations. For a complete picture, we compare several hadron lists including both experimentally discovered and theoretically predicted states. We find that the inclusion of additional resonances mildly influences the extracted parameters -- with a general trend of progressively lowering the temperature -- but is not sufficient to close the gap in temperature between light and strange hadrons previously observed in the literature., Comment: 13 pages, 8 figures, v2: typos corrected, references updated

We provide the most accurate results for the QCD transition line so far. We optimize the definition of the crossover temperature $T_c$, allowing for its very precise determination, and extrapolate from imaginary chemical potential up to real $\mu_B \approx 300$ MeV. The definition of $T_c$ adopted in this work is based on the observation that the chiral susceptibility as a function of the condensate is an almost universal curve at zero and imaganiary $\mu_B$. We obtain the parameters $\kappa_2=0.0153(18)$ and $\kappa_4=0.00032(67)$ as a continuum extrapolation based on $N_t=10,12$ and $16$ lattices with physical quark masses. We also extrapolate the peak value of the chiral susceptibility and the width of the chiral transition along the crossover line. In fact, both of these are consistent with a constant function of $\mu_B$. We see no sign of criticality in the explored range., Comment: 5 pages main text + 7 pages supplementary material

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., Comment: 5 pages, 2 figures, contribution to the proceedings from the 18th International Conference on Strangeness in Quark Matter (SQM 2019)

We present cross-correlators of QCD conserved charges at $\mu_B=0$ from lattice simulations and perform a Hadron Resonance Gas (HRG) model analysis to break down the hadronic contributions to these correlators. We construct a suitable hadronic proxy for the ratio $-\chi_{11}^{BS}/\chi_2^S$ and discuss the dependence on the chemical potential and experimental cuts. We then perform a comparison to preliminary STAR results and comment on a possible direct comparison of lattice and experiment., Comment: 6 pages, 3 figures, contribution to the proceedings from the 18th International Conference on Strangeness in Quark Matter (SQM 2019)

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 $\mu_B=0$, and present a finite-$\mu_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 $\mu_B$., Comment: 20 pages, 17 figures, 1 table, version published in Phys.Rev.D

Using the moments of the net-kaon distribution calculated within a state of-the-art hadron resonance gas model compared to experimental data from STAR's Beam Energy Scan, we find that the extracted strange freeze-out temperature is incompatible with the light one extracted from net-proton and net-charge fluctuations. Additionally predictions for net-$Lambda$ fluctuations are made that also appear to be consistent with a higher freeze-out temperature for strange particles. This strangeness freeze-out temperature is roughly $10-15$ MeV higher than the corresponding light freeze-out temperature. We also discuss cross-susceptibilities using different identified particles, which may be a further test of this two freeze-out temperature picture. Finally, we lay out the necessary updates needed in relativistic hydrodynamic models to take into account for this two freeze-out temperature scenario and present preliminary results of $\Lambda$ spectra at RHIC for AuAu $\sqrt{s_{NN}}=200$ GeV collisions that indicate a higher freeze-out temperature is preferred., Comment: 11 pages, 6 figure, Proceedings for CPOD2018

We construct an equation of state for Quantum Chromodynamics (QCD) at finite temperature and chemical potentials for baryon number $B$, electric charge $Q$ and strangeness $S$. We use the Taylor expansion method, up to the fourth power for the chemical potentials. This requires the knowledge of all diagonal and non-diagonal $BQS$ correlators up to fourth order: these results recently became available from lattice QCD simulations, albeit only at a finite lattice spacing $N_t=12$. We smoothly merge these results to the Hadron Resonance Gas (HRG) model, to be able to reach temperatures as low as 30 MeV; in the high temperature regime, we impose a smooth approach to the Stefan-Boltzmann limit. We provide a parameterization for each one of these $BQS$ correlators as functions of the temperature. We then calculate pressure, energy density, entropy density, baryonic, strangeness, electric charge densities and compare the two cases of strangeness neutrality and $\mu_S=\mu_Q=0$. We also calculate the isentropic trajectories and compare them in the two cases. Our equation of state can be readily used as an input of hydrodynamical simulations of matter created at the Relativistic Heavy Ion Collider (RHIC)., Comment: 14 pages, 6 figures; version published in Phys.Rev.C; clarified discussion of range of validity, added treatment of the speed of sound

Current knowledge of the finite-density QCD equation of state from first principles is limited to a Taylor expansion in the baryonic chemical potential around $\mu_B=0$. By means of a scaling form for the equation of state of the 3D Ising model and a non-universal, parametrized map to QCD coordinates, we construct a family of equations of state matching state of the art first principle Lattice QCD calculations and including the correct critical behavior, which can be readily employed in hydrodynamical simulations of heavy ion collisions at finite density, covering most of the BES range at RHIC. This contribution reports on work done within the Fluctuations/Equation of State working group of the BEST Collaboration., Comment: 3 pages, 1 figure, Quark Matter 2018 conference proceeding

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 \mu_B \leq$ 450 MeV. These equations of state match available lattice QCD results up to $\mathcal{O}(\mu_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., Comment: 39 pages, 18 figures, 2 tables, version published in Phys.Rev.C