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.

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.

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%

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) k- in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) k+ in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) pbar in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 7.7 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 11.5 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 19.6 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 27 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 39 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

The midrapidity (|y| < 0.1) dN/dy normalized by ⟨Npart⟩/2 as a function of √sNN for π±, K±, and p and p ̄ in 0–5% Au+Au collisions at BES energies. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

⟨mT⟩ − m of π±, K±, and p and p ̄ as a function of √sNN . Midrapidity (|y| < 0.1) results are shown for 0–5% central Au+Au collisions at BES energies. The errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

π−/π+, K−/K+, and p ̄/p ratios at midrapidity (|y| < 0.1) in central 0–5% Au+Au collisions at √sNN = 7.7, 11.5, 19.6, 27, and 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

K/π ratio at midrapidity (|y| < 0.1) for central 0–5% Au+Au collisions at √sNN = 7.7, 11.5, 19.6, 27, and 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

The GCE model particle yields fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The GCE model particle ratios fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The SCE model particle yields fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The SCE model particle ratios fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

Chemical freeze-out parameter γS plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μB plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μS plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter Tch plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter R plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μS between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Chemical freeze-out parameter γS plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μB plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter Tch plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter R plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter R between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

" (a) Energy dependence of kinetic and chemical freezeout temperatures for central heavy-ion collisions. The curves represent various theoretical predictions [81,82]. (b) Energy dependence of average transverse radial flow velocity for central heavy-ion collisions. The data points other than BES energies are taken from Refs. [43,53–64,66] and references therein. The BES data points are for 0–5% central collisions, AGS energies are mostly for 0–5%, SPS energies are for mostly 0–7%, and top RHIC and LHC energies are for 0–5% central collisions. Uncertainties represent systematic uncertainties."

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\phi_{c})\rangle$ as a function of reference multiplicity for different charge combinations, after corrections for acceptance effects. In the legend the signs indicate the charge of particles $\alpha$, $\beta$, and c. The results shown are for Au+Au collisions at 200 GeV obtained in the Full Field.

$\langle cos(\phi_{\alpha}-\phi_{\beta})\rangle$ as a function of centrality for different charge combinations and FF and RFF configurations. The data points corresponding to different charge and field configurations are slightly shifted in the horizontal direction with respect to each other for clarity. The error bars are statistical. Also shown are model predictions described in Section VII.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs. The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs after scaling by the flow of the third particle. The shaded areas represent the uncertainty from $v_{2,c}$ scaling (see text for details). The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations, with lower (in magnitude) limit obtained with elliptic flow from two-particle correlations and upper limit from four-particle cumulants. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of centrality. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of number of participants. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Three-particle correlator in Au+Au and Cu+Cu collisions compared to HIJING calculations shown as thick lines. All three particles are taken in the main TPC region, |$\eta$| < 1.0. The correlator has been scaled with number of participants and are plotted versus number of participants. In this representation HIJING results for Au+Au and Cu+Cu collisions coincide in the region of overlap.

$\langle cos(\phi_{\alpha} + \phi_{\beta} − 2\Psi_{RP})\rangle$ calculated for 200 GeV Au+Au events with event generators HIJING (with and without an “elliptic flow afterburner”), UrQMD, and MEVSIM. Blue symbols mark opposite-charge correlations, and red are same-charge. Solid stars represent the values from the data to facilitate comparison. Acceptance cuts of 0.15 < p$_{t}$ < 2 GeV/c and |$\eta$| < 1.0 were used in all cases. For MEVSIM, HIJING, and UrQMD points the true reaction plane from the generated event was used for $\Psi_{RP}$ . Thick solid lighter colored lines represent possible non-reactionplane dependent contribution from many-particle clusters as estimated by HIJING and discussed in Section VII A. Corresponding estimates from UrQMD are about factor of two smaller.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Event plane $\Phi_{2}$ resolution as a function of centrality in Au + Au 200 GeV collisions, where the vertical axis starts from 0.3 for clarity.

$\phi−\Phi_{2}$ distribution for $\phi$ meson at 1.5 < $p_{T}$ < 2.0 GeV/c in Au + Au collisions (0–80$\%$) at 200 GeV. The line is the fitting result. Error bars are statistical only.

<cos2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (full red points) and <sin2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (open blue points) as a function of $m_{inv}$ of $K^{+}K^{-}$ pairs at 0.5 < $p_{T}$ < 1.0 GeV/c in Au + Au 200 GeV collisions (0–80\%), where the solid curve is the result of fitting by Eq. (10). The arrow shows the position of the $\phi$ invariant mass peak. The dashed line shows the zero horizontal line.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the elliptic flow $v_{2}$ of $\phi$, $\Lambda$, and $K_{s}^{0}$ in Au + Au collisions (0–80$\%$) at 200 GeV. Data points for $\phi$ are from the reaction plane method (full up-triangles) and invariant mass method (full circles), where data points from the reaction plane method are shifted slightly along the $x$ axis for clarity. Vertical error bars represent statistical errors, while the square bands represent systematic uncertainties. The magenta curved band represents the $v_{2}$ of $\phi$ meson from the AMPT model with a string melting mechanism [116]. The dash and dot curves represent parametrizations inspired by number-of-quark scaling ideas from Ref. [117] for NQ = 2 and NQ = 3, respectively.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

Upper panel, Azimuthal distributions of associated particles for trigger particles in-plane (squares) and out-of-plane (triangles) for Au+Au collisions at centrality 20-60%. Open symbols are reflections of solid symbols around $\Delta \phi$ = 0 and $\Delta \phi$ = $\pi$. Elliptic flow contribution is shown by dashed lines. Lower panel, Distributions after substracting elliptic flow, and the corresponding measurement in p + p collisions (histogram).

$v_{2}$ at 3 ≤ $p_{t}$ ≤ 6 GeV/c versus impact parameter, b, compared to models of particle emission by a static source (see text).

Efficiency corrected two-particle azimuthal distributions for minimum bias and central d+Au collisions, and for p+p collisions. Curves are fits using Eq. 3, with parameters given in Table I. Trigger particles are selected with 4 GeV/c < $p_{T}$ (trig) < 6 GeV/c, while associated particles are selected within 2 GeV/c < $p_{T}$ < $p_{T}$ (trig), with |$\eta$| < 0.7 for both sets. Only statistical errors are listed. Normalization uncertainties are less than 5%.

Comparison of two-particle azimuthal distributions for central d+Au collisions to those seen in p+p and central Au+Au collisions. The respective pedestals have been subtracted. Trigger particles are selected with 4 GeV/c < $p_{T}$ (trig) < 6 GeV/c, while associated particles are selected within 2 GeV/c < $p_{T}$ < $p_{T}$ (trig), with |$\eta$| < 0.7 for both sets. Only statistical errors are listed. Normalization uncertainties are less than 5%.