We report the first measurements of cumulants, up to $4^{th}$ order, of deuteron number distributions and proton-deuteron correlations in Au+Au collisions recorded by the STAR experiment in phase-I of Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider. Deuteron cumulants, their ratios, and proton-deuteron mixed cumulants are presented for different collision centralities covering a range of center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$~=~7.7 to 200~GeV. It is found that the cumulant ratios at lower collision energies favor a canonical ensemble over a grand canonical ensemble in thermal models. An anti-correlation between proton and deuteron multiplicity is observed across all collision energies and centralities, consistent with the expectation from global baryon number conservation. The UrQMD model coupled with a phase-space coalescence mechanism qualitatively reproduces the collision-energy dependence of cumulant ratios and proton-deuteron correlations.
Cumulants of the deuteron number distributions as a function of collision energy (from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV) for most central (0 - 5$\%$), mid-central (30 - 40$\%$), and peripheral (70 - 80$\%$) Au+Au collisions. Cumulants are corrected for detector efficiency and centrality bin-width effect.
Cumulant ratios C2/C1, C3/C2, C4/C2 of deuteron number distributions and the Pearson correlation co-efficient between deuteron and proton distributions in 0-5$\%$ and 70-80$\%$ centralities in Au+Au collisions at nine energies from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV.
Cumulant ratio C4/C2 of deuteron number distributions in 0-5$\%$ centrality in Au+Au collisions at nine energies from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV.
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.