Mass dependence of the mid-rapidity $v_1$ slope, $dv_1/dy$, for $\Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. The statistical and systematic uncertainties are presented by vertical lines and square brackets, respectively. The slopes of $p$, $d$, $t$, $^3$He and $^4$He from the same collisions are shown as black circles. The blue and dashed green lines are the results of a linear fit to the measured light nuclei and hypernuclei $v_1$ slopes, respectively. For comparison, calculations of transport models plus coalescence afterburner are shown as gold and red bars from JAM model, and blue bars from UrQMD model.

Mass dependence of the mid-rapidity $v_1$ slope, $dv_1/dy$, for $\Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. The statistical and systematic uncertainties are presented by vertical lines and square brackets, respectively. The slopes of $p$, $d$, $t$, $^3$He and $^4$He from the same collisions are shown as black circles. The blue and dashed green lines are the results of a linear fit to the measured light nuclei and hypernuclei $v_1$ slopes, respectively. For comparison, calculations of transport models plus coalescence afterburner are shown as gold and red bars from JAM model, and blue bars from UrQMD model.

Mass dependence of the mid-rapidity $v_1$ slope, $dv_1/dy$, for $\Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. The statistical and systematic uncertainties are presented by vertical lines and square brackets, respectively. The slopes of $p$, $d$, $t$, $^3$He and $^4$He from the same collisions are shown as black circles. The blue and dashed green lines are the results of a linear fit to the measured light nuclei and hypernuclei $v_1$ slopes, respectively. For comparison, calculations of transport models plus coalescence afterburner are shown as gold and red bars from JAM model, and blue bars from UrQMD model.

Mass dependence of the mid-rapidity $v_1$ slope, $dv_1/dy$, for $\Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. The statistical and systematic uncertainties are presented by vertical lines and square brackets, respectively. The slopes of $p$, $d$, $t$, $^3$He and $^4$He from the same collisions are shown as black circles. The blue and dashed green lines are the results of a linear fit to the measured light nuclei and hypernuclei $v_1$ slopes, respectively. For comparison, calculations of transport models plus coalescence afterburner are shown as gold and red bars from JAM model, and blue bars from UrQMD model.

Mass dependence of the mid-rapidity $v_1$ slope, $dv_1/dy$, for $\Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. The statistical and systematic uncertainties are presented by vertical lines and square brackets, respectively. The slopes of $p$, $d$, $t$, $^3$He and $^4$He from the same collisions are shown as black circles. The blue and dashed green lines are the results of a linear fit to the measured light nuclei and hypernuclei $v_1$ slopes, respectively. For comparison, calculations of transport models plus coalescence afterburner are shown as gold and red bars from JAM model, and blue bars from UrQMD model.

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.

The $p_{T}$ dependencies of $v_{1}$ within $-0.4<y<-0.1$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{1}$ for $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{1}$ within $-0.1<y<0$ for $^4$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{1}$ within $-0.4<y<-0.1$ for $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $p$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $d$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.1<y<0$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.4<y<-0.1$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.1<y<0$ for $^4$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.4<y<-0.1$ for $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ for $p$, $d$, and $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ for $t$ and $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ from JAM plus coalescence calculations in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ for $p$, $d$, and $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ for $t$ and $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ from JAM plus coalescence calculations in 10-40% mid-central Au+Au collisions at 3 GeV.

Light nucleus scaled $v_{1}$ slopes as a function os collision energy in 10-40 mid-cantral Au+Au collisions.

Light nucleus scaled $v_{1}$ slopes from JAM plus coalescence in 10-40 mid-cantral Au+Au collisions at 3 GeV.

Centrality dependence of the proton cumulant ratios for Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV. Protons are from $-0.5 < y < 0$ and $0.4 < p_{T} < 2.0$ GeV/$c$. Systematic uncertainties are represented by gray bars. Statistical uncertainties are smaller than marker size. CBWC is applied to all cumulant ratios. While open squares represent the data without the VFC correction, blue triangles and red circles are the results with VFC using the $\langle N_{\rm{part}} \rangle$ distributions from the UrQMD and Glauber models, respectively. UrQMD model results are represented as gold dashed line.

Centrality dependence of the proton cumulant ratios for Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV. Protons are from $-0.5 < y < 0$ and $0.4 < p_{T} < 2.0$ GeV/$c$. Systematic uncertainties are represented by gray bars. Statistical uncertainties are smaller than marker size. CBWC is applied to all cumulant ratios. While open squares represent the data without the VFC correction, blue triangles and red circles are the results with VFC using the $\langle N_{\rm{part}} \rangle$ distributions from the UrQMD and Glauber models, respectively. UrQMD model results are represented as gold dashed line.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Collision energy dependence of the ratios of cumulants, $C_4/C_2$, for proton (squares) and net-proton (red circles) from top 0-5% Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for proton from $\sqrt{s_{NN}}$ = 3.0 GeV collisions is shown as a filled square. HADES data of $\sqrt{s_{NN}}$ = 2.4 GeV 0-10% collisions is also shown. The vertical black and gray bars are the statistical and systematic uncertainties, respectively. In addition, results from the HRG model, based on both Canonical Ensemble (CE) and Grand-Canonical Ensemble (GCE), and transport model UrQMD are presented.

B.R. times dN/dy of $^{3}_{\Lambda}$H vs y in 3 GeV 10-50% Au+Au collisions

B.R. times dN/dy of $^{4}_{\Lambda}$H vs y in 3 GeV 10-50% Au+Au collisions

B.R. times dN/dy at |y|<0.5 of $^{3}_{\Lambda}$H vs B.R in 3 GeV 0-10% Au+Au collisions

B.R. times dN/dy at |y|<0.5 of $^{4}_{\Lambda}$H vs B.R in 3 GeV 0-10% Au+Au collisions

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 0-10%

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 0-10%

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 10-50%

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 10-50%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 0-10%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 0-10%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.75<y<-0.5, Au+Au 3 GeV, 0-10%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 10-50%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 10-50%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.75<y<-0.5, Au+Au 3 GeV, 10-50%

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\Xi^-$ (c) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

Rapidity density distributions of $K^-$ (squares), $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\phi$ meson (circles) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\Xi^-$ (diamonds) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

$\phi/K^-$ (a) and $\phi/\Xi^-$ (b) ratio as a function of collision energy, $\sqrt{s_{\mathrm{NN}}}$. The solid black circles show the measurements presented here in 0-10\% centrality bin, while empty markers in black or grey are used for data from various other energies and/or collision systems. The grey solid line represents a THERMUS calculation based on the Grand Canonical Ensemble (GCE) while the dotted lines depict calculations based on the Canonical Ensemble (CE) with different parameters of strangeness correlation radius ($r_c$). The green dashed line, green shaded band and the solid red line show transport model calculations from the public versions $\mathrm{UrQMD}^{1}$, modified $\mathrm{UrQMD}^{2}$ and SMASH, respectively.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

The centrality dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

The $p_{\rm T}$ dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

The rapidity dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

Rapidity dependence of directed flow, v1(y), for positive pions with transverse momentum pT > 0.2 GeV/c and total momentum magnitude |p| < 1.6 GeV/c from events within 10-30% centrality. Here, the BBC-based Event Plane method is used. Plotted error bars are statistical only, and systematic errors are of comparable size.

The rapidity dependence of the directed flow for the Lambda using the TPC event plane for 10-30% centrality. Plotted error bars are statistical only, while systematic errors are +/- 0.007.

The rapidity dependence of the directed flow for the K0s using the TPC event plane for 10-30% centrality. Plotted error bars are statistical only, while systematic errors are +/- 0.017.

Beam energy dependence of the directed flow slope dv1=dy at midrapidity for baryons and mesons measured by STAR.

v2 scaled by the number of constituent quarks (nq) for charged pions and protons for 0-30% central collisions. The values of v2 scaled with nq for pions and protons are consistent with each other within errors. For comparison, points from E895 are also shown.

The excitation function v2 for all charged particles or separately for protons and pions, measured by several experiments. The STAR FXT points for protons and for pions are near the region where a change in slope occurs.

Excitation function of Rout, Rside, and Rlong.

Transverse mass dependence of Rout, Rside, and Rlong for three experiments: E895, STAR, and E866. The corresponding center-of-mass collision energies are 4.3 GeV, 4.5 GeV, and 4.9 GeV, respectively. Pairs for the STAR points are created from negative pion tracks in the momentum range 0.15 < pT < 0.8 GeV/c from events in the 0-15% centrality range.

The centrality dependence of Rout, Rside, and Rlong. Errors are statistical only. Here positive and negative pion pairs in the momentum range 0.15 < pT < 0.8 GeV/c are used.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_4\{\Psi_{4}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_5\{\Psi_{5}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_7\{\Psi_{7}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_4\{\Psi_{4}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_5\{\Psi_{5}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_7\{\Psi_{7}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of centrality.

The cumulant $c_2\{8\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in pPb collisions.

The elliptic flow $v_2\{6\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in PbPb collisions.

The elliptic flow $v_2\{8\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in PbPb collisions.

The elliptic flow $v_2\{\text{LYZ}\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in PbPb collisions.

The elliptic flow $v_2\{6\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in pPb collisions.

The elliptic flow $v_2\{8\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in pPb collisions.

The elliptic flow $v_2\{\text{LYZ}\}$ extracted for all charged particles with $0.3 < p_T < 3.0$ GeV/c as a function of $N_{trk}^{offline}$ in pPb collisions.

The cumulant ratio $c_2\{6\}/v_2\{4\}$ as a function of $c_2\{4\}/v_2\{2\}$ in pPb collisions.

The cumulant ratio $c_2\{8\}/v_2\{6\}$ as a function of $c_2\{4\}/v_2\{2\}$ in pPb collisions.

The cumulant ratio $c_2\{6\}/v_2\{4\}$ as a function of $c_2\{4\}/v_2\{2\}$ in PbPb collisions.

The cumulant ratio $c_2\{8\}/v_2\{6\}$ as a function of $c_2\{4\}/v_2\{2\}$ in PbPb collisions.

Results for G1(DEUT)/F1(DEUT) for the deuteron in bins of (W;Q**2), along with average kinematic values and correction factors for each bin. All values are averaged over the event distribution.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV $p$+$p$ collisions.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of centrality for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The scaled variance as a funciton of $N_{part}$ for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

Cross sections for final states consisting of two charm mesons and an extra pion.

Total charm cross sections for exclusive and inclusive final states.

The measured total hadronic cross section and R value with correction for initial state radiation.

The cross section for E+ E- --> ETA K+ K- excluding the contribution from PHI --> K+ K- with statistical errors only.

The cross section for E+ E- --> PHI ETA with statistical errors only.

Isoscalar (I=0) component of the cross section for E+ E- --> K*(892) K with statistical errors only.

Isovector (I=1) component of the cross section for E+ E- --> K*(892) K with statistical errors only.

Isoscalar (I=0) component of the cross section for E+ E- --> K2*(1430) K with statistical errors only.

Isovector (I=1) component of the cross section for E+ E- --> K2*(1430) K with statistical errors only.

Measured cross section for E+ E- --> OMEGA F0(975) with statistical errors only.

Measured cross section for E+ E- --> 2(PI+ PI-) ETA with statistical errorsonly.

Measured cross section for E+ E- --> ETAPRIME PI+ PI- with statistical errors only.

Measured cross section for E+ E- --> F1(1285) PI+ PI- with statistical errors only.

Measured cross section for E+ E- --> K+ K- PI+ PI- PI0 with statistical errors only.

Measured cross section for E+ E- --> PHI ETA with statistical errors only.

Measured cross section for E+ E- --> OMEGA K+ K- with statistical errors only.

Measured cross section for E+ E- --> K+ K- PI+ PI- ETA with statistical errors only.

Measurement of the E+ E- --> PHI F(980)0 cross section.

Measurement of the E+ E- --> K+ K- PI0 PI0 cross section. Statistical errors only.

Measurement of the E+ E- --> PHI F0 cross section for the F0 --> PI0 PI0 final state. Statistical errors only.

Measurement of the E+ E- --> K+ K- K+ K- cross section. Statistical errors only.

(Color Online) Variation of $N_{ch}$ normalized to the number of collisions in the FTPC coverage $(2.9 \leq \eta \leq 3.9)$ and $N_{\gamma}$ normalized to number of collisions, in the PMD coverage $(2.3 \leq \eta \leq 3.7)$ as a function of $N_{coll}$. The lower band shows the uncertainty in the ratio due to uncertainties in $N_{coll}$ calculations.

(Color Online) $dN/d\eta$ for charged particles and photons for Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for various event centrality classes.

(Color Online) $dN/d\eta$ for charged particles and photons for Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for various event centrality classes.

(Color Online) Half width at half maximum of the pseudorapidity distributions ($\eta_{h}$) of charged particles as a function of total charged particle multiplicity ($N_{T}$) normalized to the center of mass energy. The Au + Au collision data are from the PHOBOS [8] experiment and p + p collision data are from the ISR [31] experiments.

(Color Online) Half width at half maximum of the pseudorapidity distributions ($\eta_{h}$) of charged particles as a function of total charged particle multiplicity ($N_{T}$) normalized to the center of mass energy. The Au + Au collision data are from the PHOBOS [8] experiment and p + p collision data are from the ISR [31] experiments.

(Color Online) Variation of $dN_{ch}/d\eta$ normalized to $N_{part}$ with $\eta – y_{beam}$ for central and peripheral collisions for positively charged hadrons ($h^{+}$) and negatively charged hadrons ($h^{−}$).

(Color Online) The top panel shows the variation of pion rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies. Also shown is the estimated $dN_{\pi^{0}}/dy$ obtained from $dN\_{\gamma}/dy$ normalized to $N_{part}$. The bottom panel shows the variation of net proton rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies.

(Color Online) The top panel shows the variation of pion rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies. Also shown is the estimated $dN_{\pi^{0}}/dy$ obtained from $dN\_{\gamma}/dy$ normalized to $N_{part}$. The bottom panel shows the variation of net proton rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies.

Parameter $\alpha$ from the fit to the data.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 200 GeV.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 130 GeV.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 19.6 GeV.

Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of $N_P$ is the average between two energies.

Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of $N_P$ is the average between two energies.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.975 and azimuthal angle PHI = 105 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.975 and azimuthal angle PHI = 135 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.975 and azimuthal angle PHI = 165 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.875 and azimuthal angle PHI = 15 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.875 and azimuthal angle PHI = 45 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.875 and azimuthal angle PHI = 75 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.875 and azimuthal angle PHI = 105 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.875 and azimuthal angle PHI = 135 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.875 and azimuthal angle PHI = 165 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.650 and azimuthal angle PHI = 15 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.650 and azimuthal angle PHI = 45 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.650 and azimuthal angle PHI = 75 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.650 and azimuthal angle PHI = 105 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.650 and azimuthal angle PHI = 135 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Cross section for the reaction E P --> E P GAMMA at a polar angle given by COS(THETA) = -0.650 and azimuthal angle PHI = 165 degrees both in the centre-of-mass frame of the GAMMA* P --> GAMMA* P reaction.

Measured xi =-ln(2p/sqrt(s)) spectra for centre of mass energy 3.2 GeV.. Errors are statistical and systematic added in quadrature.

Measured xi =-ln(2p/sqrt(s)) spectra for centre of mass energy 4.6 GeV.. Errors are statistical and systematic added in quadrature.

Measured xi =-ln(2p/sqrt(s)) spectra for centre of mass energy 4.8 GeV.. Errors are statistical and systematic added in quadrature.

Charged particle momentum sprecrum for center of mass energy 2.2 GeV.. Errors are statistical and systematic added in quadrature.

Charged particle momentum sprecrum for center of mass energy 2.6 GeV.. Errors are statistical and systematic added in quadrature.

Charged particle momentum sprecrum for center of mass energy 3.0 GeV.. Errors are statistical and systematic added in quadrature.

Charged particle momentum sprecrum for center of mass energy 3.2 GeV.. Errors are statistical and systematic added in quadrature.

Charged particle momentum sprecrum for center of mass energy 4.6 GeV.. Errors are statistical and systematic added in quadrature.

Charged particle momentum sprecrum for center of mass energy 4.8 GeV.. Errors are statistical and systematic added in quadrature.

Two-jet rapidity distribution.