Integrated $D^0$ cross section per nucleon-nucleon collision at mid-rapidity for $p_T >0$ (a) and $p_T >4$ GeV/c (b) as a function of centrality $N_{part}$. The statistical and systematic uncertainties are shown as error bars and brackets on the data points. The green boxes on the data points depict the overall normalization uncertainties in p+p and Au+Au data respectively.

$D^0$ $R_{AA}$ in Au+Au collisions at $\sqrt{s_{NN}} =200$ GeV for $0-10\%$ (a), $10-40\%$ (b) and $40-80\%$ (c) centrality bins, respectively. The first two and last two data points are presented as empty circles, indicating that the p+p reference is extrapolated into these $p_T$ ranges. The statistical and systematic uncertainties are shown as error bars and brackets on the data points. The light and dark green boxes on the right depict the normalization uncertainties in determining the $N_{bin}$ in Au+Au collisions and the total inelastic cross section in p+p collisions, respectively.

$D^0$ and $\bar{D}^0$ invariant yields at mid-rapidity $(|y|<1)$ vs transverse momentum for different centrality classes. Error bars indicate statistical uncertainties and brackets depict systematic uncertainties. Global systematic uncertainties in B.R. and $N_{bin}$ are not plotted. Solid lines depict Levy function fits.

$\bar{D}^0/D^0$ invariant yield ratio at mid-rapidity $(|y|<1)$ vs transverse momentum for different centrality classes. Error bars indicate statistical uncertainties and brackets depict systematic uncertainties. Dashed lines depict constant function fits to the $\bar{D}^0/D^0$ ratios.

$D_0$ $R_{CP}$ with the $60-80\%$ spectrum as the reference for different centrality classes compared to the LBT model calculations shown by dashed lines [48, 49].

$D^0$ $R_{CP}$ with the $40-60\%$ spectrum as the reference for different centrality classes compared to model calculations from LBT (black dotted lines) and the Duke (blue dashedlines) groups [48-50].

BW fit parameters for $D^0$

$D^0$ invariant yield at mid-rapidity $(|y|<1)$ vs transverse kinetic energy ($m_T -m_0$) for different centrality classes. Error bars indicate statistical uncertainties and brackets depict systematic uncertainties. Global systematic uncertainties in B.R. are not plotted. Solid and dashed black lines depict exponential function fits and the dot-dashed line depicts a power-law function fit to the spectrum in the $60-80\%$ centrality bin.

mT fit parameters Teff $D^0$

mT fit parameters Teff Jpsi

mT fit parameters Teff K

mT fit parameters Teff Lambda

mT fit parameters Teff Omega

mT fit parameters Teff p

mT fit parameters Teff phi

mT fit parameters Teff pi

mT fit parameters Teff Xi

Nbin and Npart for different centralities in AuAu collisions

Vertex resolution uncertainties for different centralities in Au+Au collisions (for $R_{CP}$ plot)

Integrated $D^0$ $R_{\rm AA}$ as a function of $N_{\mathrm{part}}$ in different $p_{\rm T}$ regions, 0$-$8 ${\mathrm{GeV/}}c$ (squares), 0.7$-$2.2 ${\mathrm{GeV/}}c$ (diamonds) and 3$-$8 ${\mathrm{GeV/}}c$ (circles). Open symbols are for the 0$-$80% minimum bias events. The vertical bar around unity denotes the overall normalization uncertainty from p+p reference.

v2 ratio of 10-20%HeAu/0-10%dAu and UC pAu/0-10%dAu

v3 ratio of 10-20%HeAu/0-10%dAu and UC pAu/0-10%dAu

Inclusive Y(2S) $R_{AA}$ as a function of centrality in Au+Au collisions at 200 GeV. The bin corresponding to $N_{part}$ = 162 is for 0-60% centrality. Global uncertainty of 20.5% not shown.

Upper limit of inclusive Y(3S) $R_{AA}$ with 95% confidence level in 0-60% Au+Au collisions at 200 GeV

Upper limit of inclusive Y(3S) $R_{AA}$ with 95% confidence level in 0-60% Au+Au collisions at 200 GeV

Inclusive Y(1S) $R_{AA}$ as a function of $p_{T}$ in Au+Au collisions at 200 GeV. Global uncertainty of 6.8% not shown.

Inclusive Y(1S) $R_{AA}$ as a function of $p_{T}$ in 0-60% Au+Au collisions at 200 GeV. Global uncertainty of 6.8% not shown.

Inclusive Y(2S) R_{AA} as a function of p_{T} in Au+Au collisions at 200 GeV. Global uncertainty of 6.8% not shown.

Inclusive Y(2S) $R_{AA}$ as a function of $p_{T}$ in 0-60% Au+Au collisions at 200 GeV. Global uncertainty of 6.8% not shown.

Inclusive Y(1S) yield as a function of centrality in Au+Au collisions at 200 GeV. The bin corresponding to $N_{part}$ = 162 is for 0-60% centrality.

Inclusive Y(2S) yield as a function of centrality in Au+Au collisions at 200 GeV. The bin corresponding to $N_{part}$ = 162 is for 0-60% centrality.

Upper limit of inclusive Y(3S) yield with 95% confidence level in 0-60% Au+Au collisions at 200 GeV

Inclusive Y(1S) yield as a function of $p_{T}$ in 0-60% Au+Au collisions at 200 GeV.

Inclusive Y(2S) yield as a function of p_{T} in 0-60% Au+Au collisions at 200 GeV.

Y(2S)/Y(1S) ratio as a function of centrality in Au+Au collisions at 200 GeV. The bin corresponding to $N_{part}$ = 162 is for 0-60% centrality.

Upper limit of Y(3S)/Y(1S) ratio with 95% confidence level in 0-60% Au+Au collisions at 200 GeV

Y(2S)/Y(1S) ratio as a function of $p_{T}$ in 0-60% Au+Au collisions at 200 GeV.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Extrapolated average transverse momenta ⟨pT ⟩ as a function of dNch/dy for different particle species in Au+Au collisions at 62.4 GeV. Statistical uncertainties are represented by the error bars at the points while the systematic uncertainties are represented by the gray bars. The π, charged K and p data were extracted from Ref. [14].

Extrapolated average transverse momenta ⟨pT ⟩ as a function of dNch/dy for different particle species in Au+Au collisions at 62.4 GeV. Statistical uncertainties are represented by the error bars at the points while the systematic uncertainties are represented by the gray bars. The π, charged K and p data were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Anti-baryon to baryon yield ratios for strange baryons versus the center of mass energy √sNN. Λ/Λ is shown in the top panel while the multi-strange baryons are on the bottom panel. The data from AGS are not corrected for the weak decay feed-down from the multistrange baryons while the data from SPS and RHIC are corrected. The lines are the results of a thermal model calculation (see text section IV A). The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi- strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Anti-baryon to baryon yield ratios for strange baryons versus the center of mass energy √sNN. Λ/Λ is shown in the top panel while the multi-strange baryons are on the bottom panel. The data from AGS are not corrected for the weak decay feed-down from the multistrange baryons while the data from SPS and RHIC are corrected. The lines are the results of a thermal model calculation (see text section IV A). The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi- strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Antibaryon-to-baryon yield ratios for strange particles and protons as a function of dNch/dy at √sNN=62.4 and 200 GeV. The p data were extracted from Ref. [14]. The √sNN=200 GeV strange hadron data were extracted from Ref. [15].

Antibaryon-to-baryon yield ratios for strange particles and protons as a function of dNch/dy at √sNN=62.4 and 200 GeV. The p data were extracted from Ref. [14]. The √sNN=200 GeV strange hadron data were extracted from Ref. [15].

Particle-yield ratios as obtained by measurements (black dots) for the most central (0–5%) Au+Au collisions at 62.4 GeV and statistical model predictions (lines). The ratios indicated by the dashed lines (blue) were obtained by using only π, K, and protons, whereas the ratios indicated by the full lines (green) were obtained by also using the hyperons in the fit.

Particle-yield ratios as obtained by measurements (black dots) for the most central (0–5%) Au+Au collisions at 62.4 GeV and statistical model predictions (lines). The ratios indicated by the dashed lines (blue) were obtained by using only π, K, and protons, whereas the ratios indicated by the full lines (green) were obtained by also using the hyperons in the fit.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π− as a function of √sNN. The lines are the results of the thermal model calculation (see text Sec. 4a). The SPS values are from NA49 [20] (centrality 0–7%) and the RHIC values are from STAR [4, 15] (centrality 0–5%). For the multistrange baryons Ξ and Ω (bottom), the SPS results are from NA57 [2] (centrality 0–11%) and the RHIC values are from STAR [15, 21] (centrality 0–20%).

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π− as a function of √sNN. The lines are the results of the thermal model calculation (see text Sec. 4a). The SPS values are from NA49 [20] (centrality 0–7%) and the RHIC values are from STAR [4, 15] (centrality 0–5%). For the multistrange baryons Ξ and Ω (bottom), the SPS results are from NA57 [2] (centrality 0–11%) and the RHIC values are from STAR [15, 21] (centrality 0–20%).

Nuclear modification factor RCP, calculated as the ratio between 0–10% central spectra and 40–80% peripheral spectra, for π, K0S, Λ, and Ξ particles in Au+Au collisions at 62.4 GeV. The π RCP values were extracted from Ref. [10]. The gray band on the right side of the plot shows the uncertainties on the estimation of the number of binary collisions and the gray band on the lower left side indicates the uncertainties on the number of participants.

Nuclear modification factor RCP, calculated as the ratio between 0–10% central spectra and 40–80% peripheral spectra, for π, K0S, Λ, and Ξ particles in Au+Au collisions at 62.4 GeV. The π RCP values were extracted from Ref. [10]. The gray band on the right side of the plot shows the uncertainties on the estimation of the number of binary collisions and the gray band on the lower left side indicates the uncertainties on the number of participants.

Nuclear modification factor RCP, calculated as the ratio between 0–5% central spectra and 40–60% peripheral spectra, for Λ and Ξ particles measured in Au + Au collisions at 62.4 GeV. The gray band corresponds to the equivalent RCP curve for the Λ particles measured in Au+Au collisions at 200 GeV [15].

Nuclear modification factor RCP, calculated as the ratio between 0–5% central spectra and 40–60% peripheral spectra, for Λ and Ξ particles measured in Au + Au collisions at 62.4 GeV. The gray band corresponds to the equivalent RCP curve for the Λ particles measured in Au+Au collisions at 200 GeV [15].

Λ/K0S ratio as a function of transverse momentum for different centrality classes. 0–5% (solid circles), 40–60% (open squares), and 60–80% (solid triangles) in Au+Au collisions at 62.4 GeV.

Λ/K0S ratio as a function of transverse momentum for different centrality classes. 0–5% (solid circles), 40–60% (open squares), and 60–80% (solid triangles) in Au+Au collisions at 62.4 GeV.

Maximum value of the Λ/K0S ratio from Au+Au collisions at 62.4 GeV (solid circles) and 200 GeV (open circles) [11] as a function of ⟨Npart⟩ for different centrality classes. The lowest ⟨Npart⟩ point corresponds to p+p collisions at 200 GeV [44]. The maximum of the Λ––/K0S from Au+Au collisions at 62.4 GeV is shown as solid triangles.

Maximum value of the Λ/K0S ratio from Au+Au collisions at 62.4 GeV (solid circles) and 200 GeV (open circles) [11] as a function of ⟨Npart⟩ for different centrality classes. The lowest ⟨Npart⟩ point corresponds to p+p collisions at 200 GeV [44]. The maximum of the Λ––/K0S from Au+Au collisions at 62.4 GeV is shown as solid triangles.

This dataset corresponds to Figure 2, the v2 ratio value estimated between epd (\Psi_1) and tpc (\Psi_2) in the paper

This dataset corresponds to Figure 2, the v2 ratio value estimated between epd (\Psi_1) and epd (\Psi_2) in the paper

This dataset corresponds to Figure 3, the \Delta\gamma_{112}multiply N_{part} value estimated by tpc (\Psi_2) in the paper

This dataset corresponds to Figure 3, the \Delta\gamma_{112}multiply N_{part} value estimated by epd (\Psi_2) in the paper

This dataset corresponds to Figure 3, the \Delta\gamma_{1111}multiply N_{part} value estimated by epd (\Psi_1) in the paper

This dataset corresponds to Figure 3, the ratio between \Delta\gamma_{1111} (by epd \Psi_1) and \Delta\gamma(112) (by tpc \Psi_2) in the paper

This dataset corresponds to Figure 3, the ratio between \Delta\gamma_{1111} (by epd \Psi_1) and \Delta\gamma(112) (by epd \Psi_2) in the paper

This dataset corresponds to Figure 4, the \Delta\gamma_{112}multiply N_{part} and scaled by v_2 estimated by tpc (\Psi_2) in the paper

This dataset corresponds to Figure 4, the \Delta\gamma_{112}multiply N_{part} and scaled by v_2 estimated by epd (\Psi_2) in the paper

This dataset corresponds to Figure 4, the \Delta\gamma_{1111}multiply N_{part} and scaled by v_2 estimated by epd (\Psi_1) in the paper

This dataset corresponds to Figure 4, the ratio between \Delta\gamma_{1111} scaled by v_{211} (by epd \Psi_1) and \Delta\gamma(112) scaled by v_2 (by tpc \Psi_2) in the paper

This dataset corresponds to Figure 4, the ratio between \Delta\gamma_{1111} scaled by v_{211} (by epd \Psi_1) and \Delta\gamma(112) scaled by v_2 (by epd \Psi_2) in the paper

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV (Multiplicity class 0-10%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV (Multiplicity class 20-30%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV (Multiplicity class 30-40%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV (Multiplicity class 0-10%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV (Multiplicity class 20-30%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV (Multiplicity class 30-40%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV (Multiplicity class 0-10%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV (Multiplicity class 20-30%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV (Multiplicity class 30-40%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV (Multiplicity class 0-10%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV (Multiplicity class 20-30%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV (Multiplicity class 30-40%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV (Multiplicity class 0-10%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV (Multiplicity class 20-30%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV (Multiplicity class 30-40%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV (Multiplicity class 60-80%).

$p_{\mathrm T}$- integrated yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$p_{\mathrm T}$- integrated yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$p_{\mathrm T}$- integrated yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$p_{\mathrm T}$- integrated yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$p_{\mathrm T}$- integrated yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$p_{\mathrm T}$- integrated yield of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

<$p_{\mathrm T}$> of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

<$p_{\mathrm T}$> of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

<$p_{\mathrm T}$> of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

<$p_{\mathrm T}$> of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

<$p_{\mathrm T}$> of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

<$p_{\mathrm T}$> of $\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. $<(dN_{ch}/dy)^{1/3}>$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. $<(dN_{ch}/dy)^{1/3}>$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~11.5 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. $<(dN_{ch}/dy)^{1/3}>$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. $<(dN_{ch}/dy)^{1/3}>$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. $<(dN_{ch}/dy)^{1/3}>$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{\mathrm{K^{*0}} + \bar{\mathrm{K^{*0}}}}{K^{+} + K^{-}}$ vs. $<(dN_{ch}/dy)^{1/3}>$ in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV. Total systematic error is the quadrature sum of the correlated and uncorrelated systematic errors

$\frac{2\phi}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV

$\frac{2\phi}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV

$\frac{2\phi}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV

$\frac{2\phi}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV

$\frac{2\phi}{K^{+} + K^{-}}$ vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$7.7 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$11.5 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$14.5 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$19.6 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$27 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$39 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$62.4 GeV

lower limit of hadronic phase lifetime vs. < Npart > in AuAu collisions at $\sqrt{s_{\mathrm{NN}}}~=~$200 GeV

Reference multiplicity distributions obtained from Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV data (black markers), Glauber model (red histogram), and unfolding approach to separate single and pileup contributions. Vertical lines represent statistical uncertainties. Single, pileup, and single+pileup collisions are shown in solid blue markers, dashed green, and dashed pink lines, respectively. The 0–5% central events and 5–60% mid-central to peripheral events are indicated by black arrows. The ratio of the single+pileup to the measured multiplicity distribution is shown in the lower panel.

Reference multiplicity distributions obtained from Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV data (black markers), Glauber model (red histogram), and unfolding approach to separate single and pileup contributions. Vertical lines represent statistical uncertainties. Single, pileup, and single+pileup collisions are shown in solid blue markers, dashed green, and dashed pink lines, respectively. The 0–5% central events and 5–60% mid-central to peripheral events are indicated by black arrows. The ratio of the single+pileup to the measured multiplicity distribution is shown in the lower panel.

Proton cumulants as a function of reference multiplicity (black circles) from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Centrality-binned results with and without centrality bin width corrections are represented by red circles and blue squares, respectively. Vertical dashed lines indicate the centrality classes, from right to left: 0–5%, 5–10%, 10–20%. Data points in this figure are only corrected for detector efficiency but not for the pileup effect, which will be discussed in a later section.

Proton cumulants as a function of reference multiplicity (black circles) from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Centrality-binned results with and without centrality bin width corrections are represented by red circles and blue squares, respectively. Vertical dashed lines indicate the centrality classes, from right to left: 0–5%, 5–10%, 10–20%. Data points in this figure are only corrected for detector efficiency but not for the pileup effect, which will be discussed in a later section.

Proton cumulants as a function of reference multiplicity (black circles) from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Centrality-binned results with and without centrality bin width corrections are represented by red circles and blue squares, respectively. Vertical dashed lines indicate the centrality classes, from right to left: 0–5%, 5–10%, 10–20%. Data points in this figure are only corrected for detector efficiency but not for the pileup effect, which will be discussed in a later section.

Proton cumulants as a function of reference multiplicity from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Pileup corrected and uncorrected cumulants as a function of reference multiplicity are represented by black circles and blue open squares, respectively. Red circles and blue-filled squares represent the results of centrality binned data.

Proton cumulants as a function of reference multiplicity from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Pileup corrected and uncorrected cumulants as a function of reference multiplicity are represented by black circles and blue open squares, respectively. Red circles and blue-filled squares represent the results of centrality binned data.

Ratios of proton cumulants as a function of reference multiplicity from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Pileup corrected and uncorrected cumulants are represented by black circles and blue open squares, respectively. Red circles and blue-filled squares represent the results of centrality binned data.

Ratios of proton cumulants as a function of reference multiplicity from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Pileup corrected and uncorrected cumulants are represented by black circles and blue open squares, respectively. Red circles and blue-filled squares represent the results of centrality binned data.

UrQMD results of the proton cumulant ratios up to sixth order in Au+Au collisions at $\sqrt{s_{\rm NN}}$= 3 GeV. The black circles are without VF correction while blue squares and red triangles are results with VFC which used $N_{\rm part}$ distributions from UrQMD and Glauber models, respectively. The blue crosses are calculations using UrQMD events with b $\leq$ 3 fm. The above results are applied CBWC except for the one (blue crosses) using b $\leq$3 fm events.

UrQMD results of the proton cumulant ratios up to sixth order in Au+Au collisions at $\sqrt{s_{\rm NN}}$= 3 GeV. The black circles are without VF correction while blue squares and red triangles are results with VFC which used $N_{\rm part}$ distributions from UrQMD and Glauber models, respectively. The blue crosses are calculations using UrQMD events with b $\leq$ 3 fm. The above results are applied CBWC except for the one (blue crosses) using b $\leq$3 fm events.

Proton cumulants up to sixth order in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Data without volume fluctuation correction is shown as grey open squares while data with volume fluctuation correction using $N_{\rm part}$ distributions from Glauber and UrQMD models are shown as black circles and black open triangles, respectively. The corresponding centrality binned cumulants are shown in blue squares, red circles, and orange triangles, respectively. Similarly to Fig. 6, the vertical dashed lines indicate the centrality classes.

Proton cumulants up to sixth order in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Data without volume fluctuation correction is shown as grey open squares while data with volume fluctuation correction using $N_{\rm part}$ distributions from Glauber and UrQMD models are shown as black circles and black open triangles, respectively. The corresponding centrality binned cumulants are shown in blue squares, red circles, and orange triangles, respectively. Similarly to Fig. 6, the vertical dashed lines indicate the centrality classes.

Proton cumulant ratios up to sixth order in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Data without volume fluctuation correction are shown as grey open squares while data with volume fluctuation correction using $N_{\rm part}$ distributions from Glauber and UrQMD models are shown as black circles and black open triangles, respectively. The corresponding centrality binned cumulants are shown in blue squares, red circles, and orange triangles, respectively. Similarly to Fig. 6, the vertical dashed lines indicate the centrality classes.

Proton cumulant ratios up to sixth order in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. Data without volume fluctuation correction are shown as grey open squares while data with volume fluctuation correction using $N_{\rm part}$ distributions from Glauber and UrQMD models are shown as black circles and black open triangles, respectively. The corresponding centrality binned cumulants are shown in blue squares, red circles, and orange triangles, respectively. Similarly to Fig. 6, the vertical dashed lines indicate the centrality classes.

UrQMD results of proton cumulant ratios up to sixth order in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. The vertical dashed lines indicate the centrality classes.

Experimental results on centrality dependence of cumulants (left panels) and cumulant ratios (right panels) up to sixth order of the proton multiplicity distributions in Au+Au collisions at $N_{\rm part}$ = 3 GeV. The open squares are data without VF correction while red circles and blue triangles are results with VF correction with $N_{\rm part}$ distributions from Glauber and UrQMD models, respectively.

Experimental results on centrality dependence of cumulants (left panels) and cumulant ratios (right panels) up to sixth order of the proton multiplicity distributions in Au+Au collisions at $N_{\rm part}$ = 3 GeV. The open squares are data without VF correction while red circles and blue triangles are results with VF correction with $N_{\rm part}$ distributions from Glauber and UrQMD models, respectively.

Same as Fig. 14 but for correlation function (left panels) and their normalized ratios (right panels).

Same as Fig. 14 but for correlation function (left panels) and their normalized ratios (right panels).

Cumulants and cumulant ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. The transverse momentum window is $p_{\rm T}$ from $0.4<p_{\rm T}<2.0$ GeV/$c$ and the rapidity window is $−0.5<y<0$. Statistical and systematic uncertainties are represented by black and gray bars, respectively. UrQMD predictions are depicted by gold bands.

Cumulants and cumulant ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. The transverse momentum window is $p_{\rm T}$ from $0.4<p_{\rm T}<2.0$ GeV/$c$ and the rapidity window is $−0.5<y<0$. Statistical and systematic uncertainties are represented by black and gray bars, respectively. UrQMD predictions are depicted by gold bands.

Same as Fig. 16 but for correlation functions and correlation function ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV.

Same as Fig. 16 but for correlation functions and correlation function ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV.

The transverse-momentum and rapidity dependence of cumulant ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. In the left column, the $p_{\rm T}$ analysis window is $0.4<p_{\rm T}<2.0$ GeV/$c$ while the rapidity window is varied in the range $y_{\rm min}<y<0$. In the right column, the rapidit$y$ analysis window is $−0.5<y<0$ while the $p_{\rm T}$ is varied in the range $0.4<p_{\rm T}<p_{\rm T}^{\rm max}$ GeV/$c$. The most central (0–5%) and peripheral (50–60%) events are depicted by black squares and blue triangles, respectively. Statistical and systematic uncertainties are represented by black and gray bars, respectively. UrQMD simulations for the top 0–5% and 50–60% are shown by gold and blue bands, respectively.

The transverse-momentum and rapidity dependence of cumulant ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. In the left column, the $p_{\rm T}$ analysis window is $0.4<p_{\rm T}<2.0$ GeV/$c$ while the rapidity window is varied in the range $y_{\rm min}<y<0$. In the right column, the rapidit$y$ analysis window is $−0.5<y<0$ while the $p_{\rm T}$ is varied in the range $0.4<p_{\rm T}<p_{\rm T}^{\rm max}$ GeV/$c$. The most central (0–5%) and peripheral (50–60%) events are depicted by black squares and blue triangles, respectively. Statistical and systematic uncertainties are represented by black and gray bars, respectively. UrQMD simulations for the top 0–5% and 50–60% are shown by gold and blue bands, respectively.

The transverse-momentum and rapidity dependence of cumulant ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. In the left column, the $p_{\rm T}$ analysis window is $0.4<p_{\rm T}<2.0$ GeV/$c$ while the rapidity window is varied in the range $y_{\rm min}<y<0$. In the right column, the rapidit$y$ analysis window is $−0.5<y<0$ while the $p_{\rm T}$ is varied in the range $0.4<p_{\rm T}<p_{\rm T}^{\rm max}$ GeV/$c$. The most central (0–5%) and peripheral (50–60%) events are depicted by black squares and blue triangles, respectively. Statistical and systematic uncertainties are represented by black and gray bars, respectively. UrQMD simulations for the top 0–5% and 50–60% are shown by gold and blue bands, respectively.

The transverse-momentum and rapidity dependence of cumulant ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV. In the left column, the $p_{\rm T}$ analysis window is $0.4<p_{\rm T}<2.0$ GeV/$c$ while the rapidity window is varied in the range $y_{\rm min}<y<0$. In the right column, the rapidit$y$ analysis window is $−0.5<y<0$ while the $p_{\rm T}$ is varied in the range $0.4<p_{\rm T}<p_{\rm T}^{\rm max}$ GeV/$c$. The most central (0–5%) and peripheral (50–60%) events are depicted by black squares and blue triangles, respectively. Statistical and systematic uncertainties are represented by black and gray bars, respectively. UrQMD simulations for the top 0–5% and 50–60% are shown by gold and blue bands, respectively.

As in Fig. 18 but for transverse-momentum and rapidity dependence of correlation function ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV.

As in Fig. 18 but for transverse-momentum and rapidity dependence of correlation function ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV.

As in Fig. 18 but for transverse-momentum and rapidity dependence of correlation function ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV.

As in Fig. 18 but for transverse-momentum and rapidity dependence of correlation function ratios of proton multiplicity distributions for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV.

Collision energy dependence of the cumulant ratios: $C_2/C_1=\sigma/M$, $C_3/C_2=S\sigma$, and $C_4/C_2=\kappa\sigma^2$, for protons (open squares) and net protons (red circles) from top 0–5% (top panels) and 50–60% (bottom panels) Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for protons from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions is shown as a filled square. UrQMD results with $|y|<0.5$ for protons are shown as gold bands while those for net protons are shown as green dashed lines or green bands. At 3GeV, the model results for protons (−0.5) are shown as blue crosses. UrQMD results of proton and net-proton $C_4/C_2$, see right panels, are almost totally overlapped. The open cross is the result of the model with a fixed impact parameter $b < 3$ fm. The hydrodynamic calculations, for 5% central Au+Au collisions, for protons from $|y|<0.5$ are shown as dashed red linea and the result of the 3 GeV protons from $−0.5<y<0$ is shown as an open red star.

Collision energy dependence of the cumulant ratios: $C_2/C_1=\sigma/M$, $C_3/C_2=S\sigma$, and $C_4/C_2=\kappa\sigma^2$, for protons (open squares) and net protons (red circles) from top 0–5% (top panels) and 50–60% (bottom panels) Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for protons from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions is shown as a filled square. UrQMD results with $|y|<0.5$ for protons are shown as gold bands while those for net protons are shown as green dashed lines or green bands. At 3GeV, the model results for protons (−0.5) are shown as blue crosses. UrQMD results of proton and net-proton $C_4/C_2$, see right panels, are almost totally overlapped. The open cross is the result of the model with a fixed impact parameter $b < 3$ fm. The hydrodynamic calculations, for 5% central Au+Au collisions, for protons from $|y|<0.5$ are shown as dashed red linea and the result of the 3 GeV protons from $−0.5<y<0$ is shown as an open red star.

Collision energy dependence of the cumulant ratios: $C_2/C_1=\sigma/M$, $C_3/C_2=S\sigma$, and $C_4/C_2=\kappa\sigma^2$, for protons (open squares) and net protons (red circles) from top 0–5% (top panels) and 50–60% (bottom panels) Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for protons from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions is shown as a filled square. UrQMD results with $|y|<0.5$ for protons are shown as gold bands while those for net protons are shown as green dashed lines or green bands. At 3GeV, the model results for protons (−0.5) are shown as blue crosses. UrQMD results of proton and net-proton $C_4/C_2$, see right panels, are almost totally overlapped. The open cross is the result of the model with a fixed impact parameter $b < 3$ fm. The hydrodynamic calculations, for 5% central Au+Au collisions, for protons from $|y|<0.5$ are shown as dashed red linea and the result of the 3 GeV protons from $−0.5<y<0$ is shown as an open red star.

Collision energy dependence of the cumulant ratios: $C_2/C_1=\sigma/M$, $C_3/C_2=S\sigma$, and $C_4/C_2=\kappa\sigma^2$, for protons (open squares) and net protons (red circles) from top 0–5% (top panels) and 50–60% (bottom panels) Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for protons from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions is shown as a filled square. UrQMD results with $|y|<0.5$ for protons are shown as gold bands while those for net protons are shown as green dashed lines or green bands. At 3GeV, the model results for protons (−0.5) are shown as blue crosses. UrQMD results of proton and net-proton $C_4/C_2$, see right panels, are almost totally overlapped. The open cross is the result of the model with a fixed impact parameter $b < 3$ fm. The hydrodynamic calculations, for 5% central Au+Au collisions, for protons from $|y|<0.5$ are shown as dashed red linea and the result of the 3 GeV protons from $−0.5<y<0$ is shown as an open red star.

Collision energy dependence of the cumulant ratios: $C_2/C_1=\sigma/M$, $C_3/C_2=S\sigma$, and $C_4/C_2=\kappa\sigma^2$, for protons (open squares) and net protons (red circles) from top 0–5% (top panels) and 50–60% (bottom panels) Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for protons from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions is shown as a filled square. UrQMD results with $|y|<0.5$ for protons are shown as gold bands while those for net protons are shown as green dashed lines or green bands. At 3GeV, the model results for protons (−0.5) are shown as blue crosses. UrQMD results of proton and net-proton $C_4/C_2$, see right panels, are almost totally overlapped. The open cross is the result of the model with a fixed impact parameter $b < 3$ fm. The hydrodynamic calculations, for 5% central Au+Au collisions, for protons from $|y|<0.5$ are shown as dashed red linea and the result of the 3 GeV protons from $−0.5<y<0$ is shown as an open red star.

Collision energy dependence of the cumulant ratios: $C_2/C_1=\sigma/M$, $C_3/C_2=S\sigma$, and $C_4/C_2=\kappa\sigma^2$, for protons (open squares) and net protons (red circles) from top 0–5% (top panels) and 50–60% (bottom panels) Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for protons from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions is shown as a filled square. UrQMD results with $|y|<0.5$ for protons are shown as gold bands while those for net protons are shown as green dashed lines or green bands. At 3GeV, the model results for protons (−0.5) are shown as blue crosses. UrQMD results of proton and net-proton $C_4/C_2$, see right panels, are almost totally overlapped. The open cross is the result of the model with a fixed impact parameter $b < 3$ fm. The hydrodynamic calculations, for 5% central Au+Au collisions, for protons from $|y|<0.5$ are shown as dashed red linea and the result of the 3 GeV protons from $−0.5<y<0$ is shown as an open red star.

$v_{2}(p_{T})$ for $K_{S}^{0}$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $K_{S}^{0}$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $K_{S}^{0}$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $K_{S}^{0}$ (Centrality:40-80%)

$v_{3}(p_{T})$ for $K_{S}^{0}$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\Lambda$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\Lambda$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\Lambda$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\Lambda$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\Lambda$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\Lambda$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\Lambda$ (Centrality:40-80%)

$v_{3}(p_{T})$ for $\Lambda$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\bar{\Lambda}$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\bar{\Lambda}$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\bar{\Lambda}$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\bar{\Lambda}$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\bar{\Lambda}$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\bar{\Lambda}$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\bar{\Lambda}$ (Centrality:40-80%)

$v_{3}(p_{T})$ for $\bar{\Lambda}$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\phi$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\phi$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\phi$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\phi$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\phi$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\phi$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\phi$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\Xi^{-}$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\Xi^{-}$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\Xi^{-}$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\Xi^{-}$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\Xi^{-}$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\Xi^{-}$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\Xi^{-}$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\bar{\Xi^{+}}$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\Omega^{-}$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\Omega^{-}$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\Omega^{-}$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\Omega^{-}$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\Omega^{-}$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\Omega^{-}$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\Omega^{-}$ (Centrality:0-80%)

$v_{2}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:0-10%)

$v_{2}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:10-40%)

$v_{2}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:40-80%)

$v_{2}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:0-80%)

$v_{3}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:0-10%)

$v_{3}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:10-40%)

$v_{3}(p_{T})$ for $\bar{\Omega^{+}}$ (Centrality:0-80%)

$v_{2}$ of $\phi$ to $\bar{p}$ ratio

Difference of $v_{2}$ between particle and anti-particle

Difference of $v_{3}$ between particle and anti-particle

'Identified transverse momentum spectra of $\pi^{-}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'

'Identified transverse momentum spectra of $K^{-}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'

'Identified transverse momentum spectra of $\bar{p}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'

'Average transverse momentum as a function of $\langle$N$_{part}$$\rangle$ for $\pi^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'Average transverse momentum as a function of $\langle$N$_{part}$$\rangle$ for $K^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'Average transverse momentum as a function of $\langle$N$_{part}$$\rangle$ for p at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for $\pi^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for $K^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for p at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for $\bar{p}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$\pi^{−}$/$\pi^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$K^{−}$/$K^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$\bar{p}$/p ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$K^{+}$/$\pi^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$K^{-}$/$\pi^{-}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'p/$\pi^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$\bar{p}$/$\pi^{-}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'