Distributions In Heavy Ion Collisions Research Articles

Impact of limited statistics on the measured hyperorder cumulants of net-proton distributions in heavy-ion collisions

Impact of limited statistics on the measured hyperorder cumulants of net-proton distributions in heavy-ion collisions

Generalized angularities measurements from STAR at √SNN = 200 GeV

Jets are produced in early stages of heavy-ion collisions and undergo modified showering in the quark-gluon plasma (QGP) medium relative to a vacuum case. These modifications can be measured using observables like jet momentum profile and generalized angularities to study the details of jetmedium interactions. Jet momentum profile (ρ(r)) encodes radially differential information about jet broadening and has shown migration of charged energy towards the jet periphery in Pb+Pb collisions at the LHC. Measurements of generalized angularities (girth g and momentum dispersion pTD) and LeSub (difference between leading and subleading constituents) from Pb+Pb collisions at the LHC show harder, or more quark-like jet fragmentation, in the presence of the medium. Measuring these distributions in heavy-ion collisions at RHIC will help us further characterize the jet-medium interactions in a phase-space region complimentary to that of the LHC. In this contribution, we present the first measurements of fully corrected g, pTD and LeSub observables using hard-core jets in Au+Au collisions at √SNN = 200 GeV, collected by the STAR experiment at RHIC.

Seventh and eighth-order cumulants of net-proton multiplicity distributions in heavy-ion collisions at RHIC-STAR

We report the first measurements of seventh and eighth-order cumulants of net-proton distributions in Au+Au collisions at √SNN=27,54.4, and 200 GeV. The measurements are performed at mid-rapidity |y| < 0.5 within 0.4 < pT < 2.0 GeV/c using the Time Projection Chamber and Time-of-Flight detector. Motivation for the measurements comes from lattice-QCD and QCD based model calculations that predict their negative signs for a crossover quark-hadron transition. While 0-40% centrality measurements at √SNN=54.4 and 200 GeV are consistent with zero within large uncertainties, at √SNN=27 GeV, they are negative with ≤ 1.4σ significance. The peripheral 70-80% measurements are either positive or consistent with zero for the three energies.

Fluctuations in heavy ion collisions and global conservation effects

Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.

In-medium hadronization of heavy quarks and its effect on charmed meson and baryon distributions in heavy-ion collisions

We present a new model for the description of heavy-quark hadronization in relativistic heavy-ion collisions in the presence of a reservoir of lighter thermal particles with which recombination can occur leading to the formation of color-singlet clusters. Color neutralization is assumed to occur locally, within the same fluid cell occupied by the heavy quark and it proceeds via the recombination of the latter with light antiquarks (quarks) or diquarks (anti-diquarks), which are assumed to be present in the medium with thermal abundance as effective degrees of freedom around the QCD crossover temperature T_c. Typically the resulting color-singlet clusters have quite low invariant mass, in most of the cases below 4 GeV, and in this case they are assumed to undergo an isotropic 2-body decay in their local rest-frame. Heavier clusters are instead fragmented as Lund strings. The possibility of recombination with light diquarks enhances the yields of charmed baryons, in qualitative agreement with recent measurements. The assumption of local color neutralization leads to a strong space-momentum correlation, which provides a substantial enhancement of the collective flow of the final-state charmed hadrons, affecting both their momentum and their angular distributions.

Deblurring for nuclei: 3D characteristics of heavy-ion collisions

Observables from nuclear and high-energy experiments can be degraded by detector performance and/or methodology in extracting the observables, such as of the final-state characteristics of heavy-ion collisions in relation to a coarsely estimated reaction-plane direction. We propose the use of deblurring methods, such as in optics, to correct for observable degradation. Our main focus is the restoration of triple-differential particle distributions in heavy-ion collisions. We demonstrate that these could be extracted from collision measurements following the Richardson-Lucy deblurring method from optics. We illustrate basic features of the restoration methodology in a schematic model assuming either ideal or more realistic particle detection. The inferred three-dimensional (3D) distributions for collisions may easier to interpret in terms of collision dynamics and sought properties of bulk matter than the currently employed Fourier coefficients, that combine information from different azimuthal angles relative to the reaction plane.

Flow-plane decorrelations in heavy-ion collisions with multiple-plane cumulants

The azimuthal correlations between local flow planes at different (pseudo)rapidities ($\ensuremath{\eta}$) may reveal important details of the initial nuclear matter density distributions in heavy-ion collisions. Extensive experimental measurements of a factorization ratio (${r}_{2}$) and its derivative (${F}_{2}$) have shown evidence of the longitudinal flow-plane decorrelation. However, nonflow effects also affect this observable and prevent a quantitative understanding of the phenomenon. In this paper, to distinguish decorrelation and nonflow effects, we propose a new cumulant observable, ${T}_{2}$, which largely suppresses nonflow. The technique sensitivity to different initial-state scenarios and nonflow effects are tested with a simple Monte Carlo model, and in the end, the method is applied to events simulated by a multiphase transport model (AMPT) for $\mathrm{Au}+\mathrm{Au}$ collisions at $\sqrt{{s}_{NN}}=200$ GeV. We also emphasize that a distinct decorrelation signal requires not only the right sign of an observable but also its proper dependence on the $\ensuremath{\eta}$ window of the reference flow plane to be consistent with the pertinent decorrelation picture.

Nonmonotonic Energy Dependence of Net-Proton Number Fluctuations.

Nonmonotonic variation with collision energy (sqrt[s_{NN}]) of the moments of the net-baryon number distribution in heavy-ion collisions, related to the correlation length and the susceptibilities of the system, is suggested as a signature for the quantum chromodynamics critical point. We report the first evidence of a nonmonotonic variation in the kurtosis times variance of the net-proton number (proxy for net-baryon number) distribution as a function of sqrt[s_{NN}] with 3.1 σ significance for head-on (central) gold-on-gold (Au+Au) collisions measured solenoidal tracker at Relativistic Heavy Ion Collider. Data in noncentral Au+Au collisions and models of heavy-ion collisions without a critical point show a monotonic variation as a function of sqrt[s_{NN}].

QCD Critical Point and High Baryon Density Matter

We report the latest results on the search for the QCD critical point in the QCD phase diagram through high energy heavy-ion collisions. The measurements discussed are based on the higher moments of the net-proton multiplicity distributions in heavy-ion collisions. A non-monotonic variation in the product of kurtosis times the variance of the net-proton distribution is observed as a function of the collision energy with 3$\sigma$ significance. We also discuss the results of the thermal model in explaining the measured particle yield ratios in heavy-ion collisions and comparison of the different variants of hardon resonance gas model calculation to the data on higher moments of net-proton distributions. We end with a note that the upcoming programs in high baryon density regime at various experimental facilities will complete the search for the QCD critical point through heavy-ion collisions.

The dijet mass distribution in heavy ion collisions

In these proceedings, we review the production of both light and heavy flavor dijets in heavy ion collisions and highlight a promising observable to expose their distinct signatures. We propose the modification of dijet invariant mass distributions in heavy ion collisions as a new observable that exhibits striking sensitivity to the heavy quark mass dependence of in-medium parton showers. This observable has the advantage of amplifying the effects of jet quenching in contrast to conventional observables, such as the dijet momentum imbalance shift, which involve cancellations of such effects and, hence, result in less pronounced signals. Predictions are presented for Au+Au collisions at √SNN = 200 GeV to guide the future sPHENIX program at the Relativistic Heavy Ion Collider.

Bayesian Extraction of Jet Energy Loss Distributions in Heavy-Ion Collisions.

Based on the factorization in perturbative QCD, a jet cross section in heavy-ion collisions can be expressed as a convolution of the jet cross section in p+p collisions and a jet energy loss distribution. Using this simple expression and the Markov Chain MonteCarlo method, we carry out Bayesian analyses of experimental data on jet spectra to extract energy loss distributions for both single inclusive and γ-triggered jets in Pb+Pb collisions with different centralities at two colliding energies at the Large Hadron Collider. The average jet energy loss has a dependence on the initial jet energy that is slightly stronger than a logarithmic form and decreases from central to peripheral collisions. The extracted jet energy loss distributions with a scaling behavior in x=Δp_{T}/⟨Δp_{T}⟩ have a large width. These are consistent with the linear Boltzmann transport model simulations, in which the observed jet quenching is caused on the average by only a few out-of-cone scatterings.

Dynamical initialization and hydrodynamic modeling of relativistic heavy-ion collisions

We present a fully three-dimensional model providing initial conditions for energy and conserved charge density distributions in heavy ion collisions at RHIC Beam Energy Scan (BES) collision energies [C. Shen, B. Schenke, Dynamical initial state model for relativistic heavy-ion collisions, Phys. Rev. C97 (2) (2018) 024907. arXiv:1710.00881, doi:10.1103/PhysRevC.97.024907; C. Shen, B. Schenke, Initial state and hydrodynamic modeling of heavy-ion collisions at RHIC BES energies, PoS CPOD2017 (2018) 006. arXiv:1711.10544]. The model includes the dynamical deceleration of participating nucleons or valence quarks. It provides a realistic estimation of the initial baryon stopping during the early stage of collisions. We also present the implementation of the model with 3+1 dimensional hydrodynamics, which involves the addition of source terms that deposit energy and net-baryon densities produced by the initial state model at proper times greater than the initial time for the hydrodynamic simulation. The importance of this dynamical initialization stage on hadronic flow observables at the RHIC BES is quantified.

Strong directed flow of heavy flavor as a probe of matter distribution in heavy-ion collisions

The breaking of longitudinal boost invariance in non-central relativistic heavy ion collisions due to asymmetric local participant densities gives rise to a tilt in the reaction plane of the thermalized fireball. A direct consequence of this is the observed rapidity-odd directed flow v1 of charged particles. We study the directed flow of D and D‾ mesons by evolving the charm quark phase space distribution within Langevin dynamics dragged by matter in the tilted bulk, that itself is evolved according to relativistic hydrodynamics. This gives rise to a directed flow of D and D‾ mesons that is of the same sign but several times larger than that of the observed charged particle v1. Hence, rapidity-odd heavy flavor v1 is an ideal candidate to unravel the longitudinal profile of initial entropy deposition. Further, the electromagnetic fields of the initial stage can also give rise to rapidity-odd v1, however unlike the tilt mechanism the electromagnetic field driven v1 is of opposite signs for D and D‾. We study the interplay of the tilt and the electromagnetic fields on heavy flavor directed flow and compare our model predictions with the preliminary STAR data presented in this conference.

Heavy flavor azimuthal correlations in cold nuclear matter

Background: It has been proposed that the azimuthal distributions of heavy flavor quark-antiquark pairs may be modified in the medium of a heavy-ion collision. Purpose: This work tests this proposition through next-to-leading order (NLO) calculations of the azimuthal distribution, $d\sigma/d\phi$, including transverse momentum broadening, employing $<k_T^2>$ and fragmentation in exclusive $Q \bar Q$ pair production. While these studies were done for $p+p$, $p + \bar p$ and $p+$Pb collisions, understanding azimuthal angle correlations between heavy quarks in these smaller, colder systems is important for their interpretation in heavy-ion collisions. Methods: First, single inclusive $p_T$ distributions calculated with the exclusive HVQMNR code are compared to those calculated in the fixed-order next-to-leading logarithm approach. Next the azimuthal distributions are calculated and sensitivities to $<k_T^2>$, $p_T$ cut, and rapidity are studied at $\sqrt{s} = 7$ TeV. Finally, calculations are compared to $Q \bar Q$ data in elementary $p+p$ and $p + \bar p$ collisions at $\sqrt{s} = 7$ TeV and 1.96 TeV as well as to the nuclear modification factor $R_{p {\rm Pb}}(p_T)$ in $p+$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV measured by ALICE. Results: The low $p_T$ ($p_T < 10$ GeV) azimuthal distributions are very sensitive to the $k_T$ broadening and rather insensitive to the fragmentation function. The NLO contributions can result in an enhancement at $\phi \sim 0$ absent any other effects. Agreement with the data was found to be good. Conclusions: The NLO calculations, assuming collinear factorization and introducing $k_T$ broadening, result in significant modifications of the azimuthal distribution at low $p_T$ which must be taken into account in calculations of these distributions in heavy-ion collisions.

Dynamical initial-state model for relativistic heavy-ion collisions

We present a fully three-dimensional model providing initial conditions for energy and net-baryon density distributions in heavy ion collisions at arbitrary collision energy. The model includes the dynamical deceleration of participating nucleons or valence quarks, depending on the implementation. The duration of the deceleration continues until the string spanned between colliding participants is assumed to thermalize, which is either after a fixed proper time, or a fluctuating time depending on sampled final rapidities. Energy is deposited in space-time along the string, which in general will span a range of space-time rapidities and proper times. We study various observables obtained directly from the initial state model, including net-baryon rapidity distributions, 2-particle rapidity correlations, as well as the rapidity decorrelation of the transverse geometry. Their dependence on the model implementation and parameter values is investigated. We also present the implementation of the model with 3+1 dimensional hydrodynamics, which involves the addition of source terms that deposit energy and net-baryon densities produced by the initial state model at proper times greater than the initial time for the hydrodynamic simulation.

Effects of ρ-meson width on pion distributions in heavy-ion collisions

The influence of the finite width of ρ meson on the pion momentum distribution is studied quantitatively in the framework of the S-matrix approach combined with a blast-wave model to describe particle emissions from an expanding fireball. We find that the proper treatment of resonances which accounts for their production dynamics encoded in data for partial wave scattering amplitudes can substantially modify spectra of daughter particles originating in their two body decays. In particular, it results in an enhancement of the low-pT pions from the decays of ρ mesons which improves the quantitative description of the pion spectra in heavy ion collisions obtained by the ALICE collaboration at the LHC energy.

A description of the pseudorapidity distributions in heavy ion collisions at RHIC and LHC energies

The charged particles produced in nucleus–nucleus collisions are classified into two parts: One is from the hot and dense matter created in collisions. The other is from leading particles. The hot and dense matter is assumed to expand and generate particles according to BJP hydrodynamics, a theory put forward by A. Bialas, R.A. Janik and R. Peschanski. The leading particles are argued to possess a Gaussian rapidity distribution with the normalization constant equaling the number of participants. A comparison is made between the theoretical results and the experimental measurements performed by BRAHMS and PHOBOS Collaborations at BNL–RHIC in Au–Au and Cu–Cu collisions at sNN=200 GeV and by ALICE Collaboration at CERN–LHC in Pb–Pb collisions at sNN=2.76 TeV. The theoretical results are well consistent with experimental data.

Higher order fluctuations of strangeness and flavour hierarchy

We present a preliminary analysis on the sensitivity to the chemical freeze-out temperature of higher order moments of strange particle multiplicity distributions in heavy ion collisions. Within the Hadron Resonance Gas (HRG) model we evaluate ratios of cumulants for kaons (K±) and hyperons (Λ, Σ±, Ξ−,0, Ω−) as a function of the temperature and compare them to the sensitivity profiles obtained from ratios of particle yields. We show that ratios of higher order fluctuations of strange baryons could provide a useful tool to extract the range of freeze-out temperature, once experimental data are available. Finally, a connection to lattice data through the fourth to second cumulant ratio is made. The deconfinement transition on the lattice seems to indicate the possibility of a flavour hierarchy, namely strange quarks seem to deconfine at a higher temperature. We would like to test the possibility for the same scenario to occur at the chemical freeze-out and we show how the inclusion of multi-strange baryons in the evaluation of higher order cumulants might provide a sensitive observable to extract the freeze-out temperature.

The Unified Hydrodynamics and the Pseudorapidity Distributions in Heavy Ion Collisions at BNL-RHIC and CERN-LHC Energies

The charged particles produced in nucleus-nucleus collisions are divided into two parts. One is from the hot and dense matter created in collisions. The other is from leading particles. The hot and dense matter is assumed to expand according to unified hydrodynamics and freezes out into charged particles from a space-like hypersurface with a fixed proper time ofτFO. The leading particles are conventionally taken as the particles which inherit the quantum numbers of colliding nucleons and carry off most of incident energy. The rapidity distributions of the charged particles from these two parts are formulated analytically, and a comparison is made between the theoretical results and the experimental measurements performed in Au-Au and Pb-Pb collisions at the respective BNL-RHIC and CERN-LHC energies. The theoretical results are well consistent with experimental data.

Momentum imbalance ofDmesons in ultrarelativistic heavy-ion collisions at the CERN Large Hadron Collider

As a new observable for heavy flavor correlations the momentum imbalance $A_D$ of D mesons is proposed. It is defined analogously to the jet momentum imbalance $A_J$ of fully reconstructed jets. However, since D mesons are flavor tagged particles, no jet reconstruction is necessary. $A_D$ quantifies the influence of the medium created in heavy-ion collisions on correlated charm pairs. We present results with the partonic transport model Boltzmann Approach to MultiParton Scatterings (BAMPS), which describes well the nuclear modification factor and elliptic flow of all heavy flavor particles at RHIC and LHC. The $A_D$ distribution in heavy-ion collisions at LHC is shifted to larger values of $A_D$ compared to proton-proton collisions. We argue that this shift is due to medium effects and can be explained partially by a path length imbalance of charm pairs and partially by momentum fluctuations in the initial charm pair distribution.