Data from Figure 2, panel b, U+U, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel a, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel b, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel c, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel d, 0-5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel e, 0-5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel f, 0-5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel a, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel a, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel b, Au+Au, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel b, U+U, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel c, Au+Au, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel c, U+U, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel d, Au+Au, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel d, U+U, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel a, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel a, Au+Au, standard method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel a, Au+Au, difference, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel b, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel b, Au+Au, standard, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel b, Au+Au, difference, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel c, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel c, Au+Au, standard, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel c, Au+Au, difference, 0.2<p_{T}<3 GeV/c

Upper limits at 90% CL of $B(K^+\to\pi^+a)\times B(a\to\gamma\gamma)$ in the prompt ALP decay assumption.

See caption of Figure 9.

See caption of Figure 9.

See caption of Figure 9.

See caption of Figure 9.

See caption of Figure 9.

See caption of Figure 9.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The frequentist $3\sigma$ confidence region in $\sin^2(2\theta)$ $\Delta m^2$ for a 2-neutrino muon-to-electron oscillation fit.

The test-statistic $-2\ln\mathcal{L}$ as a function of $\Delta m^2$ and $\sin^2(2\theta)$ for a 2-neutrino muon-to-electron oscillation fit.

Observed NuE data and background prediction for arXiv:2006.16883

Observed NuMu data and prediction for arXiv:2006.16883

NuMu to NuE oscillation simulation for arXiv:2006.16883. Used used for the combined fit and the neutrino-mode analysis. ntuple has 17,204 predicted muon-to-electron neutrino and antineutrino full transmutation events simulated for the neutrino-mode beam, and contains information on reconstructed neutrino energy ($E_\text{vis}$), true neutrino energy, neutrino baseline, and event weight for each event.

Observed NuEBar data and background prediction for arXiv:2006.16883

Observed NuMuBar data and prediction for arXiv:2006.16883

NuMuBar to NuEBar oscillation simulation for arXiv:2006.16883. Used only for the combined fit. ntuple has 117,949 predicted muon-to-electron neutrino and antineutrino full transmutation events simulated for the antineutrino-mode beam, and contains information on reconstructed neutrino energy ($E_\text{vis}$), true neutrino energy, neutrino baseline, and event weight for each event. This simulation is only used for the combined fit.

NuMuBar to NuEBar oscillation simulation for arXiv:1207.4809. Used only for the antineutrino-mode analysis. ntuple has 86,403 predicted muon-to-electron neutrino and antineutrino full transmutation events simulated for the antineutrino-mode beam, and contains information on reconstructed neutrino energy ($E_\text{vis}$), true neutrino energy, neutrino baseline, and event weight for each event. This simulation is only used for the antineutrino-mode analysis.

The combined fractional covariance matrix for arXiv:2006:16883 includes bins for six classes of events in the following order: 1) neutrino-mode muon-to-electron full transmutation events, 2) neutrino-mode electron neutrino background events, 3) neutrino-mode muon neutrino CCQE events, 4) antineutrino-mode muon-to-electron full transmutation events, 5) antineutrino-mode electron antineutrino background events, 6) antineutrino-mode muon antineutrino CCQE events. Of these categories, the fractional covariance matrix includes the data statistical and systematic errors for four categories: (anti)neutrino-mode electron (anti)neutrino background events, and (anti)neutrino-mode muon (anti)neutrino CCQE events. The (anti)neutrino-mode muon-to-electron full transmutation event portions of the fractional covariance matrix only include systematic errors. The six event classes map to four observable event classifications: 1) NuE-like which includes classes 1 and 2, 2) NuMu-like which includes class 3, 3) NuEBar-like which includes classes 4 and 5, and 4) NuMuBar-like which includes class 6.

The neutrino-mode fractional covariance matrix for arXiv:2006:16883 includes bins for three classes of events in the following order: 1) neutrino-mode muon-to-electron full transmutation events, 2) neutrino-mode electron neutrino background events, 3) neutrino-mode muon neutrino CCQE events. Of these categories, the fractional covariance matrix includes the data statistical and systematic errors for two categories: neutrino-mode electron neutrino background events, and neutrino-mode muon neutrino CCQE events. The neutrino-mode muon-to-electron full transmutation event portion of the fractional covariance matrix only includes systematic errors. The three event classes map to two observable event classifications: 1) NuE-like which includes classes 1 and 2, 2) NuMu-like which includes class 3.

The antineutrino fractional covariance matrix for arXiv:1207.4809 includes bins for three classes of events in the following order: 1) antineutrino-mode muon-to-electron full transmutation events, 2) antineutrino-mode electron antineutrino background events, 3) antineutrino-mode muon antineutrino CCQE events. Of these categories, the fractional covariance matrix includes the data statistical and systematic errors for two categories: antineutrino-mode electron antineutrino background events, and antineutrino-mode muon antineutrino CCQE events. The antineutrino-mode muon-to-electron full transmutation event portions of the fractional covariance matrix only include systematic errors. The three event classes map to two observable event classifications: 1) NuEBar-like which includes classes 1 and 2, and 2) NuMuBar-like which includes class 3.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Dielectron invariant mass spectra in 2.5-3.0 GeV/c.

Dielectron invariant mass spectra in 5.0-6.0 GeV/c.

Dielectron invariant mass spectra in 6.0-7.0 GeV/c.

Dielectron invariant mass spectra in 7.0-8.0 GeV/c.

Dielectron invariant mass spectra in 8.0-9.0 GeV/c.

Dielectron invariant mass spectra in 9.0-10.0 GeV/c.

Dielectron invariant mass spectra with combined uncertainties, 2.0-2.5 GeV/C.

Data versus fit function, 2.0-2.5 GeV/C.

Data versus cocktail, 2.0-2.5 GeV/C.

Fraction of direct photon versus inclusive photon.

Direct virtual photon invariant yield in 60-80% centrality.

Direct virtual photon invariant yield in 40-60% centrality.

Direct virtual photon invariant yield in 0-80% centrality.

Direct virtual photon invariant yield in 20-40% centrality.

Direct virtual photon invariant yield in 0-20% centrality.

Direct virtual photon excess yield in 1.0-3.0GeV.

Direct virtual photon excess yield in 1.5-3.0GeV.

Direct virtual photon total yield in 1.0-3.0GeV.

Direct virtual photon total yield in 1.5-3.0GeV.