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.

'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.'

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Efficiency corrected $\phi$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.1 in the method section.

Efficiency and acceptance corrected $K^{*0}$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.4 in the method section.

$\phi$-meson $\rho_{00}$ obtained from 1st- and 2nd-order event planes. The red stars (gray squares) show the $\phi$-meson $\rho_{00}$ as a function of beam energy, obtained with the 2nd-order (1st-order) EP.

$\phi$-meson $\rho_{00}$ with respect to different quantization axes. $\phi$-meson $\rho_{00}$ as a function of beam energy, for the out-of-plane direction (stars) and the in-plane direction (diamonds). Curves are fits based on theoretical calculations with a $\phi$-meson field. The corresponding $G_s^{(y)}$ values obtained from the fits are shown in the legend.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Results for the spin-exotic $1^{-+}1^+[\pi\pi]_{1^{-\,-}}\pi P$ wave from the free-isobar partial-wave analysis performed in the fourth $t^\prime$ bin from $0.326$ to $1.000\;(\text{GeV}/c)^2$. The plotted values represent the intensity of the coherent sum of the dynamic isobar amplitudes $\{\mathcal{T}_k^\text{fit}\}$ as a function of $m_{3\pi}$, where the coherent sums run over all $m_{\pi^-\pi^+}$ bins indexed by $k$. These intensity values are given in number of events per $40\;\text{MeV}/c^2$ $m_{3\pi}$ interval and correspond to the orange points in Fig. 8(b). In the "Resources" section of this $t^\prime$ bin, we provide the JSON file named <code>transition_amplitudes_tBin_3.json</code> for download, which contains for each $m_{3\pi}$ bin the values of the transition amplitudes $\{\mathcal{T}_k^\text{fit}\}$ for all $m_{\pi^-\pi^+}$ bins, their covariances, and further information. The data in this JSON file are organized in independent bins of $m_{3\pi}$. The information in these bins can be accessed via the key <code>m3pi_bin_<#>_t_prime_bin_3</code>. Each independent $m_{3\pi}$ bin contains <ul> <li>the kinematic ranges of the $(m_{3\pi}, t^\prime)$ cell, which are accessible via the keys <code>m3pi_lower_limit</code>, <code>m3pi_upper_limit</code>, <code>t_prime_lower_limit</code>, and <code>t_prime_upper_limit</code>.</li> <li>the $m_{\pi^-\pi^+}$ bin borders, which are accessible via the keys <code>m2pi_lower_limits</code> and <code>m2pi_upper_limits</code>.</li> <li>the real and imaginary parts of the transition amplitudes $\{\mathcal{T}_k^\text{fit}\}$ for all $m_{\pi^-\pi^+}$ bins, which are accessible via the keys <code>transition_amplitudes_real_part</code> and <code>transition_amplitudes_imag_part</code>, respectively.</li> <li>the covariance matrix of the real and imaginary parts of the $\{\mathcal{T}_k^\text{fit}\}$ for all $m_{\pi^-\pi^+}$ bins, which is accessible via the key <code>covariance_matrix</code>. Note that this matrix is real-valued and that its rows and columns are indexed such that $(\Re,\Im)$ pairs of the transition amplitudes are arranged with increasing $k$.</li> <li>the normalization factors $\mathcal{N}_a$ in Eq. (13) for all $m_{\pi^-\pi^+}$ bins, which are accessible via the key <code>normalization_factors</code>.</li> <li>the shape of the zero mode, i.e., the values of $\tilde\Delta_k$ for all $m_{\pi^-\pi^+}$ bins, which is accessible via the key <code>zero_mode_shape</code>.</li> <li>the reference wave, which is accessible via the key <code>reference_wave</code>. Note that this is always the $4^{++}1^+\rho(770)\pi G$ wave.</li> </ul>

The centrality dependencies of the $\Delta\gamma\{\psi_\mathrm{ZDC}\}$ for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the a for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A/a using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A/a using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $f_\mathrm{CME}$ using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $f_\mathrm{CME}$ using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $\Delta\gamma_\mathrm{CME}$ using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $\Delta\gamma_\mathrm{CME}$ using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $f_\mathrm{CME}$ from different methods for 20-50% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $f_\mathrm{CME}$ from different methods for 50-80% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $\Delta\gamma_\mathrm{CME}$ from different methods for 20-50% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $\Delta\gamma_\mathrm{CME}$ from different methods for 50-80% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

Charged particle mean multiplicities measured in 4x4 kinematic bins of $Q^2$ and $y$.

Charged particle mean multiplicities measured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

Variance of the charged particle multiplicity measured in 4x4 kinematic bins of $Q^2$ and $y$.

Variance of the charged particle multiplicity measured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

KNO function for charged particles neasured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

Hadron entropy for charged particles selected in a sliding pseudorapidity window of 1.4 units size with transverse momenta $P_{T lab}>150$ MeV, in the accessible kinematic bins of momentum transfer $Q^2$ and $y$. The data are presented as a function of $Q^2$ and the corresponding average $x$-Bjorken.

Hadron entropy for charged particles neasured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$. The data are presented as a function of $Q^2$ and the corresponding average $x$-Bjorken.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons with a Soeding-inspired analytic function --- statistical correlations only. Only one half of the (symmetric) matrix is stored.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic $\rho^0(770)$ photoproduction cross section off protons as a function of energy (various experiments)

Fit of elastic $\rho^0(770)$ photoproduction cross sections off protons as a function of energy (various experiments) --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single differential in the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross section off protons, single differential in the momentum transfer $\vert t\vert$ and in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

The $\gamma_{SS}$ correlators in du+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

The correlation difference variable $\Delta \gamma = \gamma_{OS}-\gamma_{SS}$ in p+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

The correlation difference variable $\Delta \gamma = \gamma_{OS}-\gamma_{SS}$ in d+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

The measured two-particle cumulant $v_2\{2\}$ in p+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

The measured two-particle cumulant $v_2\{2\}$ in d+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

The scaled observable $\Delta{\gamma}_{scaled} = \Delta{\gamma} \times dN_{ch}/d\eta/v_{2}$ in p+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

The scaled observable $\Delta{\gamma}_{scaled} = \Delta{\gamma} \times dN_{ch}/d\eta/v_{2}$ in d+Au collisions at $\sqrt{s_{NN}}=200$ GeV at RHIC as a function of multiplicity.

Raw $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 27 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 39 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 200 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collisions energy dependence of $M/\sigma^2$ for $\Delta N_k$ multiplicity distributions from 0–5% most central and 70–80% peripheral collisions in Au+Au collisions at \sqrt{s_{NN}} = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collisions energy dependence of $S\sigma$ for $\Delta N_k$ multiplicity distributions from 0–5% most central and 70–80% peripheral collisions in Au+Au collisions at \sqrt{s_{NN}} = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collisions energy dependence of $\kappa\sigma^2$ for $\Delta N_k$ multiplicity distributions from 0–5% most central and 70–80% peripheral collisions in Au+Au collisions at \sqrt{s_{NN}} = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.