We report the beam energy and collision centrality dependence of fifth and sixth order cumulants ($C_{5}$, $C_{6}$) and factorial cumulants ($\kappa_{5}$, $\kappa_{6}$) of net-proton and proton distributions, from $\sqrt{s_{NN}} = 3 - 200$ GeV Au+Au collisions at RHIC. The net-proton cumulant ratios generally follow the hierarchy expected from QCD thermodynamics, except for the case of collisions at $\sqrt{s_{NN}}$ = 3 GeV. $C_{6}/C_{2}$ for 0-40% centrality collisions is increasingly negative with decreasing $\sqrt{s_{NN}}$, while it is positive for the lowest $\sqrt{s_{NN}}$ studied. These observed negative signs are consistent with QCD calculations (at baryon chemical potential, $\mu_{B} \leq$ 110 MeV) that include a crossover quark-hadron transition. In addition, for $\sqrt{s_{NN}} \geq$ 11.5 GeV, the measured proton $\kappa_{n}$, within uncertainties, does not support the two-component shape of proton distributions that would be expected from a first-order phase transition. Taken in combination, the hyper-order proton number fluctuations suggest that the structure of QCD matter at high baryon density, $\mu_{B}\sim 750$ MeV ($\sqrt{s_{NN}}$ = 3 GeV) is starkly different from those at vanishing $\mu_{B}\sim 20$MeV ($\sqrt{s_{NN}}$ = 200 GeV and higher).
Event-by-event proton multiplicity distributions for 0-40$\%$, 0-5$\%$ and 50-60$\%$ Au+Au collisions at $\sqrt{s_{NN}} = 3 GeV. The distributions are not corrected for proton and antiproton detection efficiency.
Event-by-event proton multiplicity distributions for 0-40$\%$, 0-5$\%$ and 50-60$\%$ Au+Au collisions at $\sqrt{s_{NN}} = 3 GeV. The distributions are not corrected for proton and antiproton detection efficiency.
Event-by-event proton multiplicity distributions for 0-40$\%$, 0-5$\%$ and 50-60$\%$ Au+Au collisions at $\sqrt{s_{NN}} = 3 GeV. The distributions are not corrected for proton and antiproton detection efficiency.
Non-monotonic variation with collision energy ($\sqrt{s_{\rm 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 (QCD) critical point. We report the first evidence of a non-monotonic variation in kurtosis times variance of the net-proton number (proxy for net-baryon number) distribution as a function of \rootsnn with 3.1$\sigma$ significance, for head-on (central) gold-on-gold (Au+Au) collisions measured using the STAR detector at RHIC. Data in non-central Au+Au collisions and models of heavy-ion collisions without a critical point show a monotonic variation as a function of $\sqrt{s_{\rm NN}}$.
Event-by-event net-proton multiplicity distributions for central (0-5$\%$) Au+Au collisions from $\sqrt{s_{NN}} = 7.7 - 200 GeV. The distributions are normalised to total number of events. The distributions are not corrected for proton and antiproton detection efficiency.
Cumulants of net-proton distributions in Au+Au collisions for nine energies from $\sqrt{s_{NN}} = 7.7 - 200 GeV for 0-5$\%$ and 70-80$\%$ centrality.
Cumulant ratios C3/C2 and C4/C2 of net-proton distributions in Au+Au collisions for eight energies from $\sqrt{s_{NN}} = 7.7 - 62.4 GeV for 0-5$\%$ centrality. Also given are the derivative of the polynomial fits to the C3/C2 and C4/C2 vs energy at each energy and the Skellam baselines for the ratios.