Nonmonotonic Energy Dependence of Net-Proton Number Fluctuations

Simulations of Stochastic Fluid Dynamics near a Critical Point in the Phase Diagram

We present simulations of stochastic fluid dynamics in the vicinity of a critical endpoint belonging to the universality class of the Ising model. This study is motivated by the challenge of modeling the dynamics of critical fluctuations near a conjectured critical endpoint in the phase diagram of quantum chromodynamics (QCD). We focus on the interaction of shear modes with a conserved scalar density, which is known as model H. We show that the observed dynamical scaling behavior depends on the correlation length and the shear viscosity of the fluid. As the correlation length is increased or the viscosity is decreased we observe a crossover from the dynamical exponent of critical diffusion, z≃4, to the expected scaling exponent of model H, z≃3. We use our method to investigate the time-dependent correlation function of non-Gaussian moments M^{n}(t) of the order parameter. We find that the relaxation time depends in a nontrivial manner on the power n.

Beam Energy Dependence of Triton Production and Yield Ratio (N_{t}×N_{p}/N_{d}^{2}) in Au+Au Collisions at RHIC

We report the triton (t) production in midrapidity (|y|<0.5) Au+Au collisions at sqrt[s_{NN}]=7.7-200  GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider. The nuclear compound yield ratio (N_{t}×N_{p}/N_{d}^{2}), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity (dN_{ch}/dη) and follows a scaling behavior. The dN_{ch}/dη dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3σ and 3.4σ, respectively, giving a combined significance of 4.1σ. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller p_{T} acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

Beam Energy Dependence of Fifth- and Sixth-Order Net-Proton Number Fluctuations in Au+Au Collisions at RHIC

We report the beam energy and collision centrality dependence of fifth and sixth order cumulants (C_{5}, C_{6}) and factorial cumulants (κ_{5}, κ_{6}) of net-proton and proton number distributions, from center-of-mass energy (sqrt[s_{NN}]) 3 GeV to 200 GeV Au+Au collisions at RHIC. Cumulant ratios of net-proton (taken as proxy for net-baryon) distributions generally follow the hierarchy expected from QCD thermodynamics, except for the case of collisions at 3 GeV. The measured values of C_{6}/C_{2} for 0%-40% centrality collisions show progressively negative trend with decreasing energy, while it is positive for the lowest energy studied. These observed negative signs are consistent with QCD calculations (for baryon chemical potential, μ_{B}≤110  MeV) which contains the crossover transition range. In addition, for energies above 7.7 GeV, the measured proton κ_{n}, within uncertainties, does not support the two-component (Poisson+binomial) shape of proton number distributions that would be expected from a first-order phase transition. Taken in combination, the hyperorder proton number fluctuations suggest that the structure of QCD matter at high baryon density, μ_{B}∼750  MeV at sqrt[s_{NN}]=3  GeV is starkly different from those at vanishing μ_{B}∼24  MeV at sqrt[s_{NN}]=200  GeV and higher collision energies.

Deciphering the phases of QCD matter with fluctuations and correlations of conserved charges

A review is given on recent experimental and theoretical/phenomenological developments regarding the phase structure of the strongly interacting matter. Specifically, evolution with collision energy of netproton number fluctuations as measured by several experiments are presented and their implications for the QCD phase diagram are outlined. In addition, theoretical calculations on correlations between conserved charges are presented and prospects for their experimental explorations are addressed.

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 of conserved charges in hydrodynamics and molecular dynamics

We present an overview of recent theoretical results on fluctuations of conserved charges in heavy-ion collisions obtained in relativistic hydrodynamics and molecular dynamics frameworks. In particular, we discuss the constraints on the location of the QCD critical point based on comparisons of experimental data on proton number cumulants with precision calculations of noncritical contributions. We also cover recent developments on critical fluctuations in molecular dynamics simulations.

Particle Multiplicity Fluctuations and Spatiotemporal Properties of Particle-Emitting Source of Strongly Interacting Matter for NICA and RHIC Energies

The results of the model analysis of hadron femtoscopic correlations and factorial moments of particle multiplicity in heavy ion collisions for the energy range of the Beam Energy Scan (BES) program at RHIC and future NICA collider are presented. For this purpose, the simulation of Au+Au collisions at center-of-mass energies 7.7 and 11.5 GeV per nucleon pair using the UrQMD, vHLLE+UrQMD (with the crossover and first-order equation of states), and HYDJET++ event generators was performed. The sensitivity of pion and kaon correlation radii and the dependence of the factorial moments on heavy ion beam energy to quark–hadron phase transition details was studied. In addition, the possible influence of some relevant detector effects on the corresponding experimental observables is discussed.

Measurements of Proton High-Order Cumulants in sqrt[s_{NN}]=3 GeV Au+Au Collisions and Implications for the QCD Critical Point

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au+Au collisions at sqrt[s_{NN}]=3.0  GeV, measured by the STAR experiment in the kinematic acceptance of rapidity (y) and transverse momentum (p_{T}) within -0.5<y<0 and 0.4<p_{T}<2.0  GeV/c. In the most central 0%-5% collisions, a proton cumulant ratio is measured to be C_{4}/C_{2}=-0.85±0.09 (stat)±0.82 (syst), which is 2σ below the Poisson baseline with respect to both the statistical and systematic uncertainties. The hadronic transport UrQMD model reproduces our C_{4}/C_{2} in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in C_{4}/C_{2} is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV.

Heavy-Ion Collisions toward High-Density Nuclear Matter

In the present paper, the current efforts in heavy-ion collisions toward high-density nuclear matter will be discussed. First, the essential points learned from RHIC and LHC will be reviewed. Then, the present data from the STAR Beam Energy Scan are discussed. Finally, the current efforts, NICA, FAIR, HIAF, and J-PARC-HI (heavy ion) are described. In particular, the efforts of the J-PARC-HI project are described in detail.

Probing Compressed Baryonic Matter

The objective of the compressed baryonic matter (CBM) experiment at the future Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany, is the investigation of the fundamental properties of strongly interacting matter. Of particular interest for our understanding of compact stellar objects is the determination of the equation-of-state (EOS) at high baryon densities and the exploration of the microscopic degrees-of-freedom under these conditions. The results of these laboratory experiments will complement astronomical observations, which also constrain the high-density EOS. Recent results of QCD-based calculations suggest that a possible first-order chiral phase transition should be observable in heavy-ion collisions at FAIR energies. This article reviews relevant observables from heavy-ion collisions and describes the detector configuration and the physics performance of the CBM experiment.

Measurement of the Sixth-Order Cumulant of Net-Proton Multiplicity Distributions in Au+Au Collisions at sqrt[s_{NN}]=27, 54.4, and 200 GeV at RHIC

According to first-principle lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region μ_{B}≤T_{c}. In this range the ratio, C_{6}/C_{2}, of net-baryon distributions are predicted to be negative. In this Letter, we report the first measurement of the midrapidity net-proton C_{6}/C_{2} from 27, 54.4, and 200 GeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). The dependence on collision centrality and kinematic acceptance in (p_{T}, y) are analyzed. While for 27 and 54.4 GeV collisions the C_{6}/C_{2} values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C_{6}/C_{2} ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict mostly positive values except for the most central collisions within uncertainties. These observations seem to favor a smooth crossover in the high-energy nuclear collisions at top RHIC energy.