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.