Search

In the paper, we discuss the development of the multi-gap resistive plate chamber Time-of-Flight (TOF) technology and the production of the STAR TOF detector in China at the beginning of the 21st century. Then we review recent experimental results from the first beam energy scan program (BES-I) at the Relativistic Heavy Ion Collider (RHIC). Topics cover measurements of collectivity, chirality, criticality, global polarization, strangeness, heavy-flavor, di-lepton and light nuclei productions., Comment: 31 pages, 33 figures. This review is dedicated to Professor Wenqing Shen on the occasion to celebrate his leadership of the Chinese STAR Collaboration, the development and production of the STAR MRPC TOF detector in China and many physics analyses

We investigate the impact of a critical end point (CEP) on the experimentally accessible baryon number fluctuations of different orders. By now, its potential location has been constrained fairly accurately within first principles functional QCD, together with the location of the chiral crossover line and further thermodynamic observables. This information is incorporated in an advanced QCD-assisted low energy effective theory which is used for the computation of baryon number fluctuations at the chemical freeze-out. This computation also takes care of global baryon number conservation at larger density, where the system changes from grand-canonical to canonical statistics. We observe a prominent peak structure, whose amplitude depends on the location of the CEP, while its position is more sensitive to the location of the freeze-out curve. Our results provide guidance for future low energy heavy-ion experiments., Comment: v2: 7+11 pages, 15 figures, 2 tables. Discussions on Polyakov loop potential (S.1.1) and the state-of-the-art estimates and predictions of CEP (S.2) with three further figures are added

We celebrate the legacies of our friend and mentor Professor Lianshou Liu who was one of the pioneers for the phenomenology of multi-particle interactions and initiated the physics of relativistic heavy-ion collisions in China. In this article, we discuss some of the recent exciting experimental observations on the collective phenomena including collectivity, chirality, criticality, strangeness production, and thermal equilibrium in high-energy nuclear collisions. Future directions, especially the physics at high baryon density, will be discussed with a focus on the first-order phase boundary and hyperon-nucleon interactions., Comment: 20 pages, 10 figures

We propose a novel method to construct the Landau thermodynamic potential directly from the fluctuations measured in heavy-ion collisions. The potential is capable of revealing the signal of the critical end-point (CEP) and the first order phase transition (FOPT) of QCD in the system even away from the phase transition region. With the available experimental data, we show that the criterion of the FOPT is negative for most of the collision energies which indicates no signal of FOPT. The data at $\sqrt{s_\mathrm{NN}}=7.7$ GeV with 0-5% centrality shows a different behavior and the mean value of the data satisfies the criterion. However, the uncertainty is still too large to make a certain conclusion. The higher order fluctuations are also required for confirming the signal. We emphasize therefore that new measurements with higher precision for the $C_{1,...,6}$ within 0-5% centrality in the vicinity of $\sqrt{s_\mathrm{NN}}=7.7$ GeV are in demand which may finally reveal the signal of QCD phase transition., Comment: 7 pages, 4 figures

Within the framework of intermittency analysis, a search for critical fluctuations is ongoing to locate the possible critical point in the quantum chromodynamics phase diagram. In this study, self-similar critical fluctuations from a critical Monte Carlo (CMC) model have been incorporated into the cascade ultrarelativistic quantum molecular dynamics (UrQMD) model. This hybrid UrQMD+CMC model exhibits a clear power-law behavior of scaled factorial moment for charged particles in Au+Au collisions at $\sqrt{s_\mathrm{NN}}$ = 7.7-200 GeV. By comparing the UrQMD+CMC model results with those from the STAR experiment, it is found that the value of a calculated scaling exponent falls in the range of the experimental measurement when 1-2 \% signal of intermittency fluctuations is added into the UrQMD sample., Comment: 8 pages, 6 figures, published version

We compute high-order baryon number fluctuations at finite temperature and density within a QCD-assisted low energy effective field theory. Quantum, thermal and density fluctuations are incorporated with the functional renormalization group approach. Quantum and in-medium fluctuations are encoded via the evolution of renormalization group flow equations. The resulting fourth- and sixth-order baryon number fluctuations meet the lattice benchmark results at vanishing density. They are consistent with experimental measurements, and in particular, the non-monotonic dependence of the kurtosis of net-baryon number distributions on the collision energy is observed in our calculations. This non-monotonicity arises from the increasingly sharpened chiral crossover with the decrease of collision energy., Comment: 6 pages, 3 figures. Proceedings for the International Conference on Critical Point and Onset of Deconfinement - CPOD2021, 15 - 19 March 2021, Online - zoom

Higher-order cumulants of conserved charge distributions are sensitive observables to probe the critical fluctuations near QCD critical point in heavy-ion collisions. Due to high interaction rate, pileup event can be one of the major sources of background in the measurements of higher-order cumulants. In this paper, we studied the effects of pileup events on higher-order cumulants of proton multiplicity distributions using UrQMD model. It is found that the proposed pileup correction fails if the correction parameters are determined by the Glauber fitting of charged particle multiplicities, which is usually done in the real heavy-ion experiment. To address this, we propose a model independent unfolding approach to determine the parameters in the pileup correction. This approach can be applied in the pileup correction for the future measurement of higher-order cumulants in heavy-ion collision experiment., Comment: 14 pages, 6 figures

Systems with different interactions could develop the same critical behaviour due to the underlying symmetry and universality. Using this principle of universality, we can embed critical correlations modeled on the 3D Ising model into the simulated data of heavy-ion collisions, hiding weak signals of a few inter-particle correlations within a large particle cloud. Employing a point cloud network with dynamical edge convolution, we are able to identify events with critical fluctuations through supervised learning, and pick out a large fraction of signal particles used for decision-making in each single event., Comment: 10 pages, 5 figures, version accepted by Physics Letters B

Using the nucleon coalescence model based on kinetic freeze-out nucleons from the 3D MUSIC+UrQMD and the 2D VISHNU hybrid model with a crossover equation of state, we study the multiplicity dependence of deuteron ($d$) and triton ($t$) production from central to peripheral Au+Au collisions at $\sqrt{s_\mathrm{NN}}=$ 7.7, 14.5, 19.6, 27, 39, 62.4 and 200 GeV and Pb+Pb at $\sqrt{s_\mathrm{NN}}=2.76$ TeV, respectively. It is found that the ratio $N_t N_p/N_d^2$ of the proton yield $N_p$, deuteron yield $N_d$ and triton yield $N_t$ exhibits a scaling behavior in its multiplicity dependence, i.e., decreasing monotonically with increasing charged-particle multiplicity. A similar multiplicity scaling of this ratio is also found in the nucleon coalescence calculation based on kinetic freeze-out nucleons from a multiphase transport (AMPT) model. The scaling behavior of $N_t N_p/N_d^2$ can be naturally explained by the interplay between the sizes of light nuclei and the nucleon emission source. We further argue that the multiplicity scaling of $N_t N_p/N_d^2$ can be used to validate the production mechanism of light nuclei, and the collision energy dependence of this yield ratio can further serve as a baseline in the search for the QCD critical point in relativistic heavy-ion collisions., Comment: 9 pages, 5 figures

We perform a systematic study on the effect of non-uniform track efficiency correction in higher-order cumulant analysis in heavy-ion collisions. Through analytical derivation, we find that the true values of cumulants can be successfully reproduced by the efficiency correction with an average of the realistic detector efficiency for particles with the same charges within each single phase space. The theoretical conclusions are supported by a toy model simulation by tuning the non-uniformity of the efficiency employed in the track-by-track efficiency correction method. The valid averaged efficiency is found to suppress the statistical uncertainties of the reproduced cumulants dramatically. Thus, usage of the averaged efficiency requires a careful study of phase space dependence. This study is important for carrying out precision measurements of higher-order cumulants in heavy-ion collision experiments at present and in future., Comment: 21 pages and 17 figures

Local density fluctuations near the QCD critical point has been suggested to exhibit a power-law behavior which can be probed by an intermittency analysis on scaled factorial moment (SFM) in relativistic heavy-ion collisions. The collision energy and centrality dependence of the second-order SFMs are systematically investigated in Au $+$ Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27, 39, 62.4, and 200 GeV within the UrQMD model. We estimate the noncritical background in the measurement of intermittency and propose a cumulative variable method to effectively remove the background contributions. We further study the effect of particle detection efficiency by implementing the RHIC (STAR) experimental tracking efficiencies in the UrQMD events. A cell-by-cell method is proposed for experimental application of efficiency corrections on SFM. This work can provide a guidance of background subtraction and efficiency correction for the experimental measurement of intermittency in the search of the QCD critical point in heavy-ion collisions., Comment: Replaced with the published version

Mix-cumulants of conserved charge distributions are sensitive observables for probing properties of the QCD medium and phase transition in heavy-ion collisions. To perform precise measurements, efficiency correction is one of the most important step. In this study, using the binomial efficiency model, we derive efficiency correction formulas for mutually exclusive and inclusive variables. The UrQMD model is applied to verify the validity of these formulas for different types of correlations. Furthermore, we investigate the effect of the multiplicity loss and contamination emerging from the particle identifications. This study provides an important step toward future measurements of mixed-cumulants in relativistic heavy-ion collisions., Comment: 11 pages, 7 figures

Lepton scattering is an established ideal tool for studying inner structure of small particles such as nucleons as well as nuclei. As a future high energy nuclear physics project, an Electron-ion collider in China (EicC) has been proposed. It will be constructed based on an upgraded heavy-ion accelerator, High Intensity heavy-ion Accelerator Facility (HIAF) which is currently under construction, together with a new electron ring. The proposed collider will provide highly polarized electrons (with a polarization of $\sim$80%) and protons (with a polarization of $\sim$70%) with variable center of mass energies from 15 to 20 GeV and the luminosity of (2-3) $\times$ 10$^{33}$ cm$^{-2}$ s$^{-1}$. Polarized deuterons and Helium-3, as well as unpolarized ion beams from Carbon to Uranium, will be also available at the EicC. The main foci of the EicC will be precision measurements of the structure of the nucleon in the sea quark region, including 3D tomography of nucleon; the partonic structure of nuclei and the parton interaction with the nuclear environment; the exotic states, especially those with heavy flavor quark contents. In addition, issues fundamental to understanding the origin of mass could be addressed by measurements of heavy quarkonia near-threshold production at the EicC. In order to achieve the above-mentioned physics goals, a hermetical detector system will be constructed with cutting-edge technologies. This document is the result of collective contributions and valuable inputs from experts across the globe. The EicC physics program complements the ongoing scientific programs at the Jefferson Laboratory and the future EIC project in the United States. The success of this project will also advance both nuclear and particle physics as well as accelerator and detector technology in China., Comment: EicC white paper, written by the whole EicC working group

Fluctuations of conserved charges are sensitive to the QCD phase transition and a possible critical endpoint in the phase diagram at finite density. In this work, we compute the baryon number fluctuations up to tenth order at finite temperature and density. This is done in a QCD-assisted effective theory that accurately captures the quantum- and in-medium effects of QCD at low energies. A direct computation at finite density allows us to assess the applicability of expansions around vanishing density. By using different freeze-out scenarios in heavy-ion collisions, we translate these results into baryon number fluctuations as a function of collision energy. We show that a non-monotonic energy dependence of baryon number fluctuations can arise in the non-critical crossover region of the phase diagram. Our results compare well with recent experimental measurements of the kurtosis and the sixth-order cumulant of the net-proton distribution from the STAR collaboration. They indicate that the experimentally observed non-monotonic energy dependence of fourth-order net-proton fluctuations is highly non-trivial. It could be an experimental signature of an increasingly sharp chiral crossover and may indicate a QCD critical point. The physics implications and necessary upgrades of our analysis are discussed in detail., Comment: 23 pages, 14 figures, 2 tables, published version

We study the chiral phase transition inside a rotating cylinder within the framework of the Namb--Jona-Lasinio model. A spectral boundary condition is imposed to avoid faster than light. We investigate how the geometry of the cylinder and rotation influence the chiral phase transition at finite temperature and chemical potential. The inhomogeneous effects caused by the finite size and rotation are also taken into account. It is found that finite size will reduce the chiral transition temperature and raises the chiral transition chemical potential, while the rotation reduces both the chiral transition temperature and chemical potential. In addition, we discuss the implications of our results in heavy-ion collisions and equation of states of neutron star.

We studied the effects of centrality fluctuation and deuteron formation on the cumulants ($C_n$) and correlation functions ($\kappa_n$) of protons up to sixth order in most central ($b<3$ fm) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$ = 3 GeV from a microscopic transport model (JAM). The results are presented as a function of rapidity acceptance within transverse momentum $0.4

With the aim of understanding the phase structure of nuclear matter created in high-energy nuclear collisions at finite baryon density, a beam energy scan program has been carried out at Relativistic Heavy Ion Collider (RHIC). In this mini-review, most recent experimental results on collectivity, criticality and heavy flavor productions will be discussed. The goal here is to establish the connection between current available data and future heavy-ion collision experiments in a high baryon density region., Comment: Invited review, 30 pages, 18 figures

We study the chiral phase transition of the two-flavor Nambu-Jona-Lasinio (NJL) model in a rotating sphere, which includes both rotation and finite size effects. We find that rotation leads to a suppression of the chiral condensate at a finite temperature, while its effects are smaller than the finite size effects. Our work can be helpful to study the effects relevant to rotation in heavy-ion collisions in a more realistic way.

We employ a new self-consistent mean field approximation of NJL model, which introduces a free parameter $\alpha$ ($\alpha$ reflects the weight of different interaction channels), to study the effects of the chiral chemical potential $\mu_5$ on QCD phase structure, especially the location of the QCD critical end point (CEP). We find that, at a high temperature, the critical temperature of QCD phase transition smoothly increases with $\mu_5$ at the beginning, and then decreases rapidly. At low temperature and high baryon density region, the increase of the chiral chemical potential will reduce the critical chemical potential of phase transition. The temperature of the CEP shows a non-monotonic dependence on the chiral chemical potential with a long plateau around the maximum. At $\mu_5=0$, we found that the CEP will disappear when the $\alpha$ value is larger than 0.71 and will reappear when the $\mu_5$ increases. This study is important for exploring the QCD phase structure and the location of CEP in chiral imbalanced systems., Comment: 6 pages, 6 figures

We studied the centrality selection effect on cumulants (up to fourth order) and the cumulants ratios of net-proton multiplicity distributions in Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 7.7, 19.6 and 200 GeV from UrQMD model. The net-proton cumulants are calculated with collision centralities by using charged particle multiplicity from different pesudorapidity ($\eta$) region. By comparing the results from various collision centralities, we found that the autocorrelation effects are not significant in the results with collision centralities "refmult-3" and "refmult-2", which are using mid-rapidity charged particles but excluding (anti-)protons and analysis region, respectively. Furthermore, due to the contributions of spectator protons, we observed poor centrality resolution when using charged particles at forward $\eta$ region at low energies. This work can serve as a baseline for centrality selection of future fluctuations analysis in relativistic heavy-ion collisions., Comment: 11 pages and 12 figures, Minor update on Fig. 5, 6, 12 and references

Light nuclei production is sensitive to the baryon density fluctuations and can be used to probe the QCD phase transition in relativistic heavy-ion collisions. In this work, we studied the production of proton, deuteron, triton in central Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5, 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4 and 200 GeV from a transport model (JAM). Based on the coalescence production of light nuclei, we calculated the energy dependence of rapidity density $dN/dy$ and particle ratios ($d/p$, $t/p$, and $t/d$). More importantly, the yield ratio $N_{{t}} \times N_{{p}} / N_{{d}}^{2}$, which is sensitive to the neutron density fluctuations, shows a flat energy dependence and cannot describe the non-monotonic trend observed by the STAR experiment. Based on the nucleon coalescence, this work can provide constraint and reference to search for the QCD critical point and/or first order phase transition with light nuclei production in future heavy-ion collision experiments., Comment: 5 pages, 4 figures

The system of a quarkonium confined by an impenetrable spherical cavity filled with a hot quantum chromodynamics (QCD) medium is studied by solving the Schr\"{o}dinger equation. This is the first time this issue has been raised for discussion. The Schr\"{o}dinger equation with an appropriate boundary condition of a quarkonium in an impenetrable cavity filled with a hot medium is derived. The numerical results are obtained with the help of Gaussian Expansion Method. Binding energies and radii of the ground and low-excited states are obtained as a function of the medium temperature and the cavity radius. We find the behaviour of quarkonium in this cavity is different from that in infinite space. Our results show that the quarkonium dissociation temperature decreases as the cavity radius decreases and the finite volume effects on the ground state are more obvious than on the excited states. We also find that the less mass of the constituents and the bigger radius of the quarkonium lead the finite volume effects to become more obvious., Comment: 9 pages, 16 figures

Proton number fluctuation is sensitive observable to search for the QCD critical point in heavy-ion collisions. In this paper, we studied rapidity acceptance dependence of the proton cumulants and correlation functions in most central Au+Au collisions at $\sqrt{s_\mathrm{NN}} = 5$ GeV from a microscopic hadronic transport model (JAM). At mid-rapidity, we found the effects of resonance weak decays and hadronic re-scattering on the proton cumulants and correlation functions are small, but those effects get larger when further increasing the rapidity acceptance. On the other hand, we found the baryon number conservation is a dominant background effect on the rapidity acceptance dependence of proton number fluctuations. It leads to a strong suppression of cumulants and cumulant ratios, as well as the negative proton correlation functions. We also studied those two effects on the energy dependence of cumulant ratios of net-proton distributions in most central Au+Au collisions at $\sqrt{s_\mathrm{NN}} = 5-200$ GeV from JAM model. This work can serve as a non-critical baseline for future QCD critical point search in heavy-ion collisions at high baryon density region., Comment: 8 pages, 7 figures, To be published in Physical Review C

We propose the transverse velocity ($\beta_T$) dependence of the anti-deuteron to deuteron ratio as a new observable to search for the QCD critical point in heavy-ion collisions. The QCD critical point can attract the system evolution trajectory in the QCD phase diagram, which is known as focusing effect. To quantify this effect, we employ thermal model and hadronic transport model to simulate the dynamical particle emission along a hypothetical focusing trajectory near critical point. We found the focusing effect can lead to anomalous $\beta_T$ dependence of $\bar{p}/p$, $\bar{d}/d$ and $^3\overline{\text{He}}/^3\text{He}$ ratios. We examined the $\beta_T$ dependence of $\bar{p}/p$ and $\bar{d}/d$ ratios of central Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 to 200 GeV measured by the STAR experiment at RHIC. Surprisingly, we only observe a negative slope in $\beta_T$ dependence of $\bar{d}/d$ ratio at $\sqrt{s_{NN}}$ = 19.6 GeV, which indicates the trajectory evolution has passed through the critical region. In the future, we could constrain the location of the critical point and/or width of the critical region by making precise measurements on the $\beta_T$ dependence of $\bar{d}/d$ ratio at different energies and rapidity., Comment: 4 pages, 5 figures, To be published in Chinese Physics C

We present a detailed analysis of the beam energy dependence of the mechanisms for the generation of directed and elliptic flows in Au+Au collisions focusing on the role of hadronic rescattering and spectator shadowing within a microscopic transport model JAM with different equation of state. A systematic study of the beam energy dependence is performed for Au+Au collisions at $\sqrt{s_{NN}} =2.3 - 62.4$ GeV. The transition of the dynamical origin of the directed flow is observed. We find that the initial Glauber type nucleon-nucleon collisions generate negative $v_1$ for nucleons at midrapidity due to the presence of spectator matter, and this negative nucleon $v_1$ is turned to be positive by the meson-baryon interactions at the beam energy region of $\sqrt{s_{NN}} < 30$ GeV. In contrast, above 30 GeV there is no spectator shadowing at midrapidity, and initial nucleon-nucleon collisions do not generate directed flow, but subsequent rescatterings among produced particles generate negative $v_1$ for nucleons. It is demonstrated that negative pion-directed flows are mostly generated by the interaction with the spectator matter. It is also shown that the squeeze-out effect is largely suppressed in the case of softening, which leads to the enhancement of elliptic flow around $\sqrt{s_{NN}}=5-7$ GeV. The elliptic flow at midrapidity above 10 GeV is not influenced by the squeeze-out due to spectator matter, while its effect is seen at the forward rapidity range of $y/y_\mathrm{c.m.}>0.5$, which decreases as beam energy increases., Comment: 14 pages, 17 figures

We performed systematic studies on the effects of event-by-event efficiency fluctuations on efficiency correction for cumulant analysis in relativistic heavy-ion collision experiments. Experimentally, particle efficiencies of events measured under different experimental conditions should be different. For fluctuation measurements, the final event-by-event multiplicity distributions should be the superposed distributions of various type of events measured under different conditions. We demonstrate efficiency fluctuation effects using numerical simulation, in which we construct an event ensemble consisting of events with two different efficiencies. By using the mean particle efficiencies, we find that the efficiency corrected cumulants show large deviations from the original inputs when the discrepancy between the two efficiencies is large. We further studied the effects of efficiency fluctuations for the cumulants of net-proton distributions by implementing the UrQMD events of Au+Au collisions at $\sqrt{s_{NN}}=7.7$ GeV in a realistic STAR detector acceptance. We consider the unequal efficiency in two sides of the Time Projection Chamber (TPC), multiplicity dependent efficiency, and the event-by-event variations of the collision vertex position along the longitudinal direction ($V_\mathrm{z}$). When the efficiencies fluctuate dramatically within the studied event sample, the effects of efficiency fluctuations have significant impacts on the efficiency corrections of cumulants with the mean efficiency. We find that this effect can be effectively suppressed by binning the entire event ensemble into various sub-event samples, in which the efficiency variations are relatively small. The final efficiency corrected cumulants can be calculated from the weighted average of the corrected factorial moments of the sub-event samples with the mean efficiency., Comment: 21 pages, 14 figures

The beam energy dependence of the elliptic flow,$v_2$, is studied in mid-central Au+Au collisions in the energy range of $3\leq \sqrt{s_{NN}} \leq 30$ GeV within the microscopic transport model JAM. The results of three different modes of JAM are compared; cascade-,hadronic mean field-, and a new mode with modified equations of state, with a first order phase transition (1.O.P.T.) and with a crossover transition. The standard hadronic mean field suppresses $v_2$, while the inclusion of the effects of a 1.O.P.T. (and also of a crossover transition) does enhance $v_2$ at $\sqrt{s_{NN}}<30$ GeV. The enhancement or suppression of the scaled energy flow, dubbed "elliptic flow"is understood as being due to out of plane- flow, i.e. $v_2<0$, dubbed out of plane - "squeeze-out", which occurs predominantly in the early, compression stage. Subsequently, the in-plane flow dominates, in the expansion stage, $v_2 > 0$. The directed flow, dubbed "bounce- off", is an independent measure of the pressure, which quickly builds up the transverse momentum transfer in the reaction plane. When the spectator matter leaves the participant fireball region, where the highest compression occurs, a hard expansion leads to larger $v_2$. A combined analysis of the three transverse flow coefficients, radial $v_0$-, directed $v_1$- and elliptic $v_2$- flow, in the beam energy range of $3\leq\sqrt{s_{NN}}\leq10$ GeV, distinguishes the different compression and expansion scenarios: a characteristic dependence on the early stage equation of state is observed. The enhancement of both the elliptic and the transverse radial flow and the simultaneous collapse of the directed flow of nucleons offers a clear signature if 1.O.P.T. is realized at the highest baryon densities created in high energy heavy-ion collisions., Comment: 12 pages,10 figures

We discuss properties and applications of factorial cumulants of various particle numbers and for their mixed channels measured by the event-by-event analysis in relativistic heavy-ion collisions. After defining the factorial cumulants for systems with multi-particle species, their properties are elucidated. The uses of the factorial cumulants in the study of critical fluctuations are discussed. We point out that factorial cumulants play useful roles in understanding fluctuation observables when they have underlying physics approximately described by the binomial distribution. As examples, we suggest novel utilization methods of the factorial cumulants in the study of the momentum cut and rapidity window dependences of fluctuation observables., Comment: 17 pages, 4 figures

We studied the acceptance dependence of proton cumulants (up to fourth order) and correlation functions in 0--5\% most central Au+Au collisions at $\sqrt{s_\text{NN}}$=7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV from UrQMD model. We found that high order proton cumulants show suppressions at large acceptance. By decomposing the proton cumulants into linear combination of multi-proton correlation functions, we observed the two-proton correlation functions always show negative values due to the effects of baryon number conservations. The three and four-proton correlation functions are close to zero and show negligible acceptance dependence. We further observed that the proton cumulants and correlation functions follow similar trends and show a scaling behavior when plotting the results versus mean number of protons. The comparisons between experimental data and the UrQMD calculations show that the non-monotonic energy dependence of proton correlation functions measured by STAR experiment cannot be described by the UrQMD model. The UrQMD calculations can provide us baselines for the experimental studies of the proton cumulants and correlation functions. Finally, we propose to measure the rapidity dependence of the reduced proton correlation functions to search for the QCD critical point in heavy-ion collisions., Comment: 9 pages, 8 figures, Accepted by Physics Letters B

Fluctuations of conserved quantities, such as baryon, electric charge and strangeness number, are sensitive observables in heavy-ion collisions to search for the QCD phase transition and critical point. In this paper, we performed a systematical analysis on the various cumulants and cumulant ratios of event-by-event net-strangeness distributions in Au+Au collisions at $\sqrt{s_{NN}}$=7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV from UrQMD model. We performed a systematical study on the contributions from various strange baryons and mesons to the net-strangeness fluctuations. The results demonstrate that the cumulants and cumulant ratios of net-strangeness distributions extracted from different strange particles show very different centrality and energy dependence behavior. By comparing with the net-kaon fluctuations, we found that the strange baryons play an important role in the fluctuations of net-strangeness. This study can provide useful baselines to study the QCD phase transition and search for the QCD critical point by using the fluctuations of net-strangeness in heavy-ion collisions experiment. It can help us to understand non-critical physics contributions to the fluctuations of net-strangeness., Comment: 8 pages, 8 figures

Fluctuations of conserved quantities, such as baryon, electric charge and strangeness number, are sensitive observables in relativistic heavy-ion collisions to probe the QCD phase transition and search for the QCD critical point. In this paper, we review the experimental measurements of the cumulants (up to fourth order) of event-by-event net-proton (proxy for net-baryon), net-charge and net-kaon (proxy for net-strangeness) multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}=7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, 200$ GeV from the first phase of beam energy scan program at the Relativistic Heavy-Ion Collider (RHIC). We also summarize the data analysis methods of suppressing the volume fluctuations, auto-correlations and the unified description of efficiency correction and error estimation. Based on theoretical and model calculations, we will discuss the characteristic signatures of critical point as well as backgrounds for the fluctuation observables in heavy-ion collisions. The physics implications and the future second phase of the beam energy scan (2019-2020) at RHIC will be also discussed., Comment: 42 pages, 40 figures, Invited Review for Nuclear Science and Techniques

We perform a systematic study of elliptic flow ($v_2$) in Au+Au collision at $\sqrt{s_{NN}}=5$ GeV by using a microscopic transport model JAM. The centrality, pseudorapidity, transverse momentum and beam energy dependence of $v_2$ for charged as well as identified hadrons are studied. We investigate the effects of both the hadronic mean-field and the softening of equation of state (EoS) on elliptic flow. The softening of EoS is realized by imposing attractive orbits in two body scattering, which can reduce the pressure of the system. We found that the softening of EoS leads to the enhancement of $v_2$, while the hadronic mean-field suppresses $v_2$ relative to the cascade mode. It indicates that elliptic flow at high baryon density regions is highly sensitive to the EoS and the enhancement of $v_2$ may probe the signature of a first-order phase transition in heavy-ion collisions at beam energies of a strong baryon stopping region., Comment: 7 pages, 9 figures, to appear in Chinese Physics C

Fluctuations and correlations of conserved charges are sensitive observables for studying the QCD phase transition and critical point in high-energy heavy-ion collisions. We have studied the centrality and energy dependence of mixed-cumulants (up to fourth order) between net-baryon and net-strangeness in Au+Au collisions at $\sqrt{s_{NN}}$= 7.7, 11.5, 19.6, 27, 39, 62.4, 200 GeV from UrQMD model. To compare with other theoretical calculations, we normalize these mixed-cumulants by various order cumulants of net-strangeness distributions. We found that the results obtained from UrQMD calculations are comparable with the results from Lattice QCD at low temperature and hadron resonance gas model. The ratios of mixed-cumulants ($R_{11}^{BS},R_{13}^{BS},R_{22}^{BS}$,~$R_{31}^{BS}$) from UrQMD calculations show weak centrality dependence. However, the mixed-cumulant ratios $R_{11}^{BS}$ and $R_{31}^{BS}$ show strong increase at low energy, while the $R_{13}^{BS}$ snd $R_{22}^{BS}$ are similar at different energies. Furthermore, we have also studied the correlations between different hadron species and their contributions to the net-baryon and net-strangeness correlations. These models studies can provide baselines for searching for the signals of QCD phase transition and critical point in heavy-ion collisions., Comment: 11 pages, 5 figures

Baryon number fluctuations are sensitive to the QCD phase transition and QCD critical point. According to the Feynman rules of finite-temperature field theory, we calculated various order moments and cumulants of the baryon number distributions in the quasi-particle model of quark gluon plasma. Furthermore, we compared our results with the experimental data measured by the STAR experiment at RHIC. It is found that the experimental data can be well described by the model for the colliding energies above 30 GeV and show large discrepancies at low energies. It can put new constraint on qQGP model and also provide a baseline for the QCD critical point search in heavy-ion collisions at low energies., Comment: 13 pages, 5 figures

We analyze the rapidity and transverse momentum dependence for the cumulants of the net-proton and net-baryon distributions in Au+Au collisions at $\sqrt{s_{\text{NN}}} = 5\,\text{GeV}$ with a microscopic hadronic transport (JAM) model. To study the effects of mean field potential and softening of equation of state (EoS) on the fluctuations of net-proton (baryon) in heavy-ion collisions, the calculations are performed with two different modes. The softening of EoS is realized in the model by implementing the attractive orbit in the two-body scattering to introduce a reduction pressure of the system. By comparing the results from the two modes with the results from default cascade, we find the mean field potential and softening of EoS have strong impacts on the rapidity distributions ($\text{d}N/\text{d}y$) and the shape of the net-proton (baryon) multiplicity distributions. The net-proton (baryon) cumulants and their ratios calculated from all of the three modes are with similar trends and show significant suppression with respect to unity, which can be explained by the presence of baryon number conservations. It indicates that the effects of mean field potential and softening of EoS might be not the ingredients that are responsible to the observed strong enhancement in the most central Au+Au collisions at 7.7 GeV measured by the STAR experiment at RHIC., Comment: 6 pages, 5 figures

Fluctuations of conserved quantities are sensitive observables to probe the signature of QCD phase transition and critical point in heavy-ion collisions. With the UrQMD model, we have studied the centrality and energy dependence of various order cumulants and cumulant ratios (up to fourth order) of net-proton,net-charge and net-kaon multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}$= 7.7, 11.5, 19.6, 27, 39, 62.4, 200 GeV. The model results show that the production mechanism of the particles and anti-particles have significant impacts on the cumulants of net-particles multiplicity distributions and show strong energy dependence. We also made comparisons between model calculations and experimental data measured in the first phase of the beam energy scan (BES) program by the STAR experiment at RHIC. The comparisons indicate that the baryon conservation effect strongly suppress the cumulants of net-proton distributions at low energies and the non-monotonic energy dependence for the net-proton {\KV} at the most central Au+Au collisions measured by the STAR experiment can not be described by the UrQMD model. Since there has no physics of QCD phase transition and QCD critical point implemented in the UrQMD, the model results provide us baselines and qualitative estimates about the non-critical background contributions to the fluctuations observables in heavy-ion collisions., Comment: 9 pages, 10 figures

We show the importance of separated efficiency corrections between positively and negatively charged particles for cumulant calculations by Monte Carlo toy models and analytical calculations. Our results indicate that S{\sigma} in published net-proton results from the STAR experiment will be suppressed about 5 to 10% in central collisions, and 10 to 20% in peripheral collisions at the beam energy of \sqrt s_{NN} = 62.4 and 200 GeV if the separated efficiencies are used to efficiency correction., Comment: 10 pages and 6 figures

Beam energy scan programs in heavy-ion collisions aim to explore the QCD phase structure at high baryon density. Sensitive observables are applied to probe the signatures of the QCD phase transition and critical point in heavy-ion collisions at RHIC and SPS. Intriguing structures, such as dip, peak and oscillation, have been observed in the energy dependence of various observables. In this paper, an overview is given and corresponding physics implications will be discussed for the experimental highlights from the beam energy scan programs at the STAR, PHENIX and NA61/SHINE experiments. Furthermore, the beam energy scan phase II at RHIC (2019-2020) and other future experimental facilities for studying the physics at low energies will be also discussed., Comment: 8 pages, 14 figures. Proceedings of the Quark Matter 2015 conference, Kobe, Japan

We present the energy dependence of moments of net-proton and net-charge multiplicity distributions in Au+Au collisions measured by the STAR experiment in the first phase of the Beam Energy Scan (BES) at the Relativistic Heavy Ion Collider (RHIC). By using the time of flight detector for particle identification, the upper transverse momentum ($p_{T}$) limit for proton and anti-proton is extended from 0.8 GeV/c up to 2 GeV/c. The $p_{T}$ and rapidity acceptance dependence study for the moments of net-proton distribution show that the larger the acceptance is, the greater the deviation from unity. The most pronounced structure is found in the energy dependence of {\KV} of net-proton distributions from the $0\sim5%$ most central collisions within $0.4

Fluctuations of conserved quantities in heavy-ion collisions are used to probe the phase transition and the QCD critical point for the strongly interacting hot and dense nuclear matter. The STAR experiment has carried out moment analysis of net-proton (proxy for net-baryon (B)), net-kaon (proxy for net-strangeness (S)), and net-charge (Q). These measurements are important for understanding the quantum chromodynamics phase diagram. We present the analysis techniques used in the moment analysis by the STAR experiment and discuss the moments of net-proton and net-charge distributions from the first phase of the Beam Energy Scan program at the Relativistic Heavy Ion Collider., Comment: 5 pages, 2 figures, Proceedings of Workshop for Young Scientists with Research Interests Focused on Physics at FAIR, Italy, Sept. 22-27, 2014 (FAIRNESS Workshop 2014)

We provide a unified description of efficiency correction and error estimation for moments of conserved quantifies in heavy-ion collisions. Moments and cumulants are expressed in terms of the factorial moments, which can be easily corrected for the efficiency effect. By deriving the covariance between factorial moments, one can obtain the general error formula for the efficiency corrected moments based on the error propagation derived from the Delta theorem. The skellam distribution based Monto Carlo simulation is used to test the Delta theorem and Bootstrap error estimation methods. The statistical errors calculated from the two methods can well reflect the statistical fluctuations of the efficiency corrected moments., Comment: 23 pages, 5 figures

We present a systematic comparison between the recently measured cumulants of the net-proton distributions by STAR for 0-5% central Au+Au collisions at $\sqrt{s_{NN}}$=7.7-200 GeV and two kinds of possible baseline measures, the Poisson and Binomial baselines. These baseline measures are assuming that the proton and anti-proton distributions independently follow Poisson statistics or Binomial statistics. The higher order cumulant net-proton data are observed to deviate from all the baseline measures studied at 19.6 and 27 GeV. We also compare the net-proton with net-baryon fluctuations in UrQMD and AMPT model, and convert the net-proton fluctuations to net-baryon fluctuations in AMPT model by using a set of formula., Comment: 4 pages, 5 figures, Proceedings of XXIV International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2014), May 19-24, 2014, Darmstadt, Germany

We present various order moments (Variance ($\sigma^2$), Skewness($S$), Kurtosis($\kappa$)) of net-proton multiplicity distributions in Au+Au collisions from the first phase of RHIC beam energy scan ({\sNN}=7.7$-$200 GeV). The measurements are carried out by the STAR detector at mid-rapidity ($|y|<0.5$) and within transverse momentum range $0.4

Moments (Variance ($\sigma^2$), Skewness($S$), Kurtosis($\kappa$)) of multiplicity distributions of conserved quantities, such as net-baryon,net-charge and net-strangeness, are predicted to be sensitive to the correlation length of the system and connected to the thermodynamic susceptibilities computed in Lattice QCD and Hadron Resonance Gas (HRG) model. In this paper, we present several measurement artifacts that could lead to volume fluctuation and auto-correlation effects in the moment analysis of net-proton multiplicity distributions in heavy-ion collisions using the UrQMD model. We discuss methods to overcome these artifacts so that the extracted moments could be used to obtain physical conclusions. In addition we present methods to properly estimate the statistical errors in moment analysis., Comment: 8 pages, 12 figures

We present the measurements of the higher moments of the net-proton multiplicity distributions for Au+Au collisions from the first phase of the beam energy scan program at RHIC. The measurements are carried out at midrapidity ($|y|<0.5$) and within the transverse momentum range $0.4

Higher moments of net-proton multiplicity distributions are applied to search for the QCD critical point. In this paper, we will present measurements for kurtosis ($\kappa$), skewness ($S$) and variance $\sigma^2$) of net-proton multiplicity distributions at the mid-rapidity ($|y|<0.5$) and transverse momentum range $0.4

Fluctuations of conserved quantities are predicted to be sensitive to the correlation length and connected to the thermodynamic susceptibility. Thus, moments of net-baryon, net-charge and net-strangeness have been extensively studied theoretically and experimentally to explore phase structure and bulk properties of QCD matters created in heavy ion collision experiment. As the moments analysis is statistics hungry study, the error estimation is crucial to extract physics information from the limited experimental data. In this paper, we will derive the limit distributions and error formula based on Delta theorem in statistics for various order moments used in the experimental data analysis. The Monte Carlo simulation is also applied to test the error formula., Comment: 14 pages, 10 figures

Higher moments of event-by-event net-proton multiplicity distributions are applied to search for the QCD critical point in the heavy ion collisions. It has been demonstrated that higher moments as well as moment products are sensitive to the correlation length and directly connected to the thermodynamic susceptibilities computed in the Lattice QCD and Hadron Resonance Gas (HRG) model. In this paper, we will present measurements for kurtosis ($\kappa$), skewness ($S$) and variance ($\sigma^{2}$) of net-proton multiplicity distributions at the mid-rapidity ($|y|<0.5$) and $0.4

Higher moments of event-by-event net-proton multiplicity distributions have been applied to search for the QCD critical point. Model results are used to provide a baseline for this search. The measured moment products, $\kappa\sigma^2$ and $S\sigma$ of net-proton distributions, which are directly connected to the thermodynamical baryon number susceptibility ratio in Lattice QCD and Hadron Resonance Gas (HRG) model, are compared to the transport and thermal model results. We argue that a non-monotonic dependence of $\kappa \sigma^2$ and $S \sigma$ as a function of beam energy can be used to search for the QCD critical point., Comment: 7 pages, 3 figures. CPOD 2010 Proceedings

Matter described by Quantum Chromodynamics (QCD), the theory of strong interactions, may undergo phase transitions when its temperature and the chemical potentials are varied. QCD at finite temperature is studied in the laboratory by colliding heavy-ions at varying beam energies. We present a test of QCD in the non-perturbative domain through a comparison of thermodynamic fluctuations predicted in lattice computations with the experimental data of baryon number distributions in high-energy heavy-ion collisions. This study provides evidence for thermalization in these collisions, and allows us to find the crossover temperature between normal nuclear matter and a deconfined phase called the quark gluon plasma. This value allows us to set a scale for the phase diagram of QCD., Comment: 16 pages and 4 figures. Accepted for publication in Science