QCD Critical Point and Net-Proton Number Fluctuations at RHIC-STAR

In the search of QCD phase boundary and critical point, higher-order cumulants of conserved quantities are proposed as promising observables and have been studied extensively both experimentally and theoretically. In this paper we present cumulant ratios up to $6^{th}$-order of net-proton number distributions in Au+Au collisions at $\sqrt{\mathrm{s_{NN}}}$ = 7.7 - 200 GeV from STAR Beam Energy Scan program phase I and $\sqrt{\mathrm{s}}$ = 200 GeV $p+p$ collisions. The results are compared with various models and Lattice QCD calculations.


Introduction
In heavy-ion collision physics it is predicted that under a very high temperature and baryon density a deconfined quark-gluon plasma (QGP) phase can be created and studying the QCD phase structure is one of the main goals. A phase diagram in terms of temperature and baryon chemical potential (µ B ) is usually used to explore the QCD phase structure. Regarding the phase transition between the QGP phase and hadronic phase, first principle Lattice QCD calculation [1] at µ B = 0 MeV suggests a smooth crossover transition. While at large µ B , various QCD-based models predict first order phase transition [2]. Thermodynamically there should be an end point of the first order phase boundary which is called QCD critical point. The possible QCD critical point and the first order phase boundary have been investigated both experimentally [3][4][5] and theoretically.
Higher-order cumulants of conserved quantities like net-baryon (B), net-charge (Q) and net-strangeness (S ) are proposed as promising observables to search for the QCD critical point and the first order phase boundary. Higher-order cumulants are sensitive to the correlation length (ξ) [6] and are directly related to susceptibility (χ) of the system [7]. It is predicted that the fourth-order fluctuations will exhibit a non-monotonic energy dependence [8][9][10] when passing through the critical region. For 5 th -and 6 th -order cumulants recent calculations from Lattice QCD [11] and the functional renormalisation group approach (FRG) [12] show that they will be negative due to the crossover transition between QGP and hadronic phase. At high baryon density region, on the other hand, they are also sensitive to the first order phase boundary [13].
The data of Au+Au collisions at √ s NN = 7.7 -200 GeV and p + p collisions at √ s = 200 GeV are collected in RHIC Beam Energy Scan program phase I. Protons and antiprotons are identified by the Time Projection Chamber (TPC) and Time of flight (TOF) detectors at rapidity window −0.5 < y < 0.5 and transverse momentum window 0.4 < p T < 2.0 GeV/c.
The centrality is determined using charged particle multiplicity within |η| < 1.0 excluding protons and antiprotons to avoid auto-correlation effect [14]. The centrality bin width correction [14] is applied to suppress initial volume fluctuation effect. Cumulants are calculated at each multiplicity bin and then their weighted averages are taken for each centrality bin. The weight is the number of events at the corresponding multiplicity bin. Detector efficiency correction [15] in cumulant calculations are done by assuming binomial detector efficiency. Statistical uncertainties of cumulants are estimated by Bootstrap [16] and Delta methods [17].  Figure 1 shows energy dependence of S σ and κσ 2 of net-proton distributions from 0-5% and 70-80% centrality bins within |y| < 0.5 and 0.4 < p T < 2.0 GeV/c in Au+Au collisions at √ s NN = 7.7 -200 GeV [20,21]. The S σ (left panel) shows a decreasing trend with the increase of collision energy both in central and peripheral collisions. The decreasing trend can be qualitatively described by HRG [18] and UrQMD [19] models. The κσ 2 (right panel) shows a non-monotonic energy dependence in central collisions while the results for peripheral collisions show monotonic energy dependence. The non-monotonic trend in central collisions can not be described by different conditions (GCE, EV, and CE) of HRG and UrQMD models. Figure 2 shows energy dependence of C 5 /C 1 and C 6 /C 2 of net-proton distributions from 0-40% and 70-80% centrality bins within |y| < 0.5 and 0.4 < p T < 2.0 GeV/c in Au+Au collisions at √ s NN = 7.7 -200 GeV. It is suggested from Lattice QCD and FRG calculations that 5 th -and 6 th -order cumulants show negative signs while calculations from UrQMD and HRG models are consistent with either zero or unity. In UrQMD and HRG models, no phase transition physics is implemented. The measurements of BES-I data are shown as blue circles for 0-40% and red diamonds for 70-80%. The cumulant ratio C 5 /C 1 (left panel) deviates from zero with less than 2σ significance. In peripheral collisions C 5 /C 1 shows positive sign for all energies. The ratio C 6 /C 2 (right panel) for 0-40% decreases with negative sign with less than 2σ significance when decreasing energy while it shows positive sign in peripheral collisions (70-80%) for all energies.  Figure 3 shows multiplicity dependence of net-proton results of C 5 /C 1 and C 6 /C 2 within |y| < 0.5 and 0.4 < p T < 2.0 GeV/c in p + p collisions at √ s = 200 GeV. We see that the cumulant ratios (C 4 /C 2 , C 5 /C 1 , and C 6 /C 2 ) from p + p collisions fit into the multiplicity dependence of results from Au+Au collisions which are shown with triangles. The cumulant ratios C 5 /C 1 and C 6 /C 2 are negative for 0-40% in Au+Au collisions and positive at peripheral Au+Au collisions and p + p collisions. Pythia [22]

Summary
In this proceedings, we report the measurements of net-proton cumulant ratios up to 6 th -order in Au+Au collisions at √ s NN = 7.7 -200 GeV and p + p collisions at √ s = 200 GeV from STAR. With results from 200 GeV p + p collisions and the energy dependence of C 4 /C 2 , C 5 /C 1 , and C 6 /C 2 from the BES-I data sets, and the comparison with LQCD calculations, we conclude: 1) QCD matter is indeed created in the 200 GeV central (0-5%) Au+Au collisions at RHIC; 2) non-monotonic energy dependence of C 4 /C 2 is observed from the most central (0-5%) Au+Au collisions [20,21]. Future results from BES-II and STAR fixed-target experiment √ s NN = 3 GeV data sets will allow to answer if QCD critical point exists in the covered energy region. In 2023 to 2025, STAR plans to collect 15 to 20 billion events of Au+Au collisions at √ s NN = 200 GeV. This will allow us to perform more precise measurements of higher-order cumulants, maybe even up to 8 th -order.