Balance functions have been measured in terms of relative pseudorapidity ($\Delta \eta$) for charged particle pairs at the Relativistic Heavy-Ion Collider (RHIC) from Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV to 200 GeV using the STAR detector. These results are compared with balance functions measured at the Large Hadron Collider (LHC) from Pb+Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV by the ALICE Collaboration. The width of the balance function decreases as the collisions become more central and as the beam energy is increased. In contrast, the widths of the balance functions calculated using shuffled events show little dependence on centrality or beam energy and are larger than the observed widths. Balance function widths calculated using events generated by UrQMD are wider than the measured widths in central collisions and show little centrality dependence. The measured widths of the balance functions in central collisions are consistent with the delayed hadronization of a deconfined quark gluon plasma (QGP). The narrowing of the balance function in central collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV implies that a QGP is still being created at this relatively low energy.

The extreme temperatures and energy densities generated by ultra-relativistic collisions between heavy nuclei produce a state of matter with surprising fluid properties. Non-central collisions have angular momentum on the order of 1000$\hbar$, and the resulting fluid may have a strong vortical structure that must be understood to properly describe the fluid. It is also of particular interest because the restoration of fundamental symmetries of quantum chromodynamics is expected to produce novel physical effects in the presence of strong vorticity. However, no experimental indications of fluid vorticity in heavy ion collisions have so far been found. Here we present the first measurement of an alignment between the angular momentum of a non-central collision and the spin of emitted particles, revealing that the fluid produced in heavy ion collisions is by far the most vortical system ever observed. We find that $\Lambda$ and $\overline{\Lambda}$ hyperons show a positive polarization of the order of a few percent, consistent with some hydrodynamic predictions. A previous measurement that reported a null result at higher collision energies is seen to be consistent with the trend of our new observations, though with larger statistical uncertainties. These data provide the first experimental access to the vortical structure of the "perfect fluid" created in a heavy ion collision. They should prove valuable in the development of hydrodynamic models that quantitatively connect observations to the theory of the Strong Force. Our results extend the recent discovery of hydrodynamic spin alignment to the subatomic realm.

Notwithstanding decades of progress since Yukawa first developed a description of the force between nucleons in terms of meson exchange, a full understanding of the strong interaction remains a major challenge in modern science. One remaining difficulty arises from the non-perturbative nature of the strong force, which leads to the phenomenon of quark confinement at distances on the order of the size of the proton. Here we show that in relativistic heavy-ion collisions, where quarks and gluons are set free over an extended volume, two species of produced vector (spin-1) mesons, namely $\phi$ and $K^{*0}$, emerge with a surprising pattern of global spin alignment. In particular, the global spin alignment for $\phi$ is unexpectedly large, while that for $K^{*0}$ is consistent with zero. The observed spin-alignment pattern and magnitude for the $\phi$ cannot be explained by conventional mechanisms, while a model with a connection to strong force fields, i.e. an effective proxy description within the Standard Model and Quantum Chromodynamics, accommodates the current data. This connection, if fully established, will open a potential new avenue for studying the behaviour of strong force fields.

We present results from a harmonic decomposition of two-particle azimuthal correlations measured with the STAR detector in Au+Au collisions for energies ranging from $\sqrt{s_{NN}}=7.7$ GeV to 200 GeV. The third harmonic $v_3^2\{2\}=\langle \cos3(\phi_1-\phi_2)\rangle$, where $\phi_1-\phi_2$ is the angular difference in azimuth, is studied as a function of the pseudorapidity difference between particle pairs $\Delta\eta = \eta_1-\eta_2$. Non-zero {\vthree} is directly related to the previously observed large-$\Delta\eta$ narrow-$\Delta\phi$ ridge correlations and has been shown in models to be sensitive to the existence of a low viscosity Quark Gluon Plasma (QGP) phase. For sufficiently central collisions, $v_3^2\{2\}$ persist down to an energy of 7.7 GeV suggesting that QGP may be created even in these low energy collisions. In peripheral collisions at these low energies however, $v_3^2\{2\}$ is consistent with zero. When scaled by pseudorapidity density of charged particle multiplicity per participating nucleon pair, $v_3^2\{2\}$ for central collisions shows a minimum near {\snn}$=20$ GeV.

Local parity-odd domains are theorized to form inside a Quark-Gluon-Plasma (QGP) which has been produced in high-energy heavy-ion collisions. The local parity-odd domains manifest themselves as charge separation along the magnetic field axis via the chiral magnetic effect (CME). The experimental observation of charge separation has previously been reported for heavy-ion collisions at the top RHIC energies. In this paper, we present the results of the beam-energy dependence of the charge correlations in Au+Au collisions at midrapidity for center-of-mass energies of 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV from the STAR experiment. After background subtraction, the signal gradually reduces with decreased beam energy, and tends to vanish by 7.7 GeV. The implications of these results for the CME will be discussed.

Measurements of the elliptic flow, $v_{2}$, of identified hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{s}^{0}$, $p$, $\bar{p}$, $\phi$, $\Lambda$, $\bar{\Lambda}$, $\Xi^{-}$, $\bar{\Xi}^{+}$, $\Omega^{-}$, $\bar{\Omega}^{+}$) in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV are presented. The measurements were done at mid-rapidity using the Time Projection Chamber and the Time-of-Flight detectors of the STAR experiment during the Beam Energy Scan program at RHIC. A significant difference in the $v_{2}$ values for particles and the corresponding anti-particles was observed at all transverse momenta for the first time. The difference increases with decreasing center-of-mass energy, $\sqrt{s_{NN}}$ (or increasing baryon chemical potential, $\mu_{B}$) and is larger for the baryons as compared to the mesons. This implies that particles and anti-particles are no longer consistent with the universal number-of-constituent quark (NCQ) scaling of $v_{2}$ that was observed at $\sqrt{s_{NN}}=$ 200 GeV. However, for the group of particles NCQ scaling at $(m_{T}-m_{0})/n_{q}>$ 0.4 GeV/$c^{2}$ is not violated within $\pm$10%. The $v_{2}$ values for $\phi$ mesons at 7.7 and 11.5 GeV are approximately two standard deviations from the trend defined by the other hadrons at the highest measured $p_{T}$ values.

Elliptic flow (v_2) values for identified particles at midrapidity in Au + Au collisions measured by the STAR experiment in the Beam Energy Scan at the Relativistic Heavy Ion Collider at sqrt{s_{NN}}= 7.7--62.4 GeV are presented for three centrality classes. The centrality dependence and the data at sqrt{s_{NN}}= 14.5 GeV are new. Except at the lowest beam energies we observe a similar relative v_2 baryon-meson splitting for all centrality classes which is in agreement within 15% with the number-of-constituent quark scaling. The larger v_2 for most particles relative to antiparticles, already observed for minimum bias collisions, shows a clear centrality dependence, with the largest difference for the most central collisions. Also, the results are compared with A Multiphase Transport Model and fit with a Blast Wave model.

We present measurements of $\pi^-$ and $\pi^+$ elliptic flow, $v_2$, at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}} =$ 200, 62.4, 39, 27, 19.6, 11.5 and 7.7 GeV, as a function of event-by-event charge asymmetry, $A_{ch}$, based on data from the STAR experiment at RHIC. We find that $\pi^-$ ($\pi^+$) elliptic flow linearly increases (decreases) with charge asymmetry for most centrality bins at $\sqrt{s_{_{\rm NN}}} = \text{27 GeV}$ and higher. At $\sqrt{s_{_{\rm NN}}} = \text{200 GeV}$, the slope of the difference of $v_2$ between $\pi^-$ and $\pi^+$ as a function of $A_{ch}$ exhibits a centrality dependence, which is qualitatively similar to calculations that incorporate a chiral magnetic wave effect. Similar centrality dependence is also observed at lower energies.

Elliptic flow ($v_{2}$) values for identified particles at mid-rapidity in Au+Au collisions, measured by the STAR experiment in the Beam Energy Scan at RHIC at $\sqrt{s_{NN}}=$ 7.7--62.4 GeV, are presented. A beam-energy dependent difference of the values of $v_{2}$ between particles and corresponding anti-particles was observed. The difference increases with decreasing beam energy and is larger for baryons compared to mesons. This implies that, at lower energies, particles and anti-particles are not consistent with the universal number-of-constituent-quark (NCQ) scaling of $v_{2}$ that was observed at $\sqrt{s_{NN}}=$ 200 GeV.

The inclusive $J/\psi$ transverse momentum ($p_{T}$) spectra and nuclear modification factors are reported at midrapidity ($|y|<1.0$) in Au+Au collisions at $\sqrt{s_{NN}}=$ 39, 62.4 and 200 GeV taken by the STAR experiment. A suppression of $J/\psi$ production, with respect to {\color{black}the production in $p+p$ scaled by the number of binary nucleon-nucleon collisions}, is observed in central Au+Au collisions at these three energies. No significant energy dependence of nuclear modification factors is found within uncertainties. The measured nuclear modification factors can be described by model calculations that take into account both suppression of direct $J/\psi$ production due to the color screening effect and $J/\psi$ regeneration from recombination of uncorrelated charm-anticharm quark pairs.

We present measurements of elliptic flow ($v_2$) of electrons from the decays of heavy-flavor hadrons ($e_{HF}$) by the STAR experiment. For Au+Au collisions at $\sqrt{s_{\rm NN}} = $ 200 GeV we report $v_2$, for transverse momentum ($p_T$) between 0.2 and 7 GeV/c using three methods: the event plane method ($v_{2}${EP}), two-particle correlations ($v_2${2}), and four-particle correlations ($v_2${4}). For Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 62.4 and 39 GeV we report $v_2${2} for $p_T< 2$ GeV/c. $v_2${2} and $v_2${4} are non-zero at low and intermediate $p_T$ at 200 GeV, and $v_2${2} is consistent with zero at low $p_T$ at other energies. The $v_2${2} at the two lower beam energies is systematically lower than at $\sqrt{s_{\rm NN}} = $ 200 GeV for $p_T < 1$ GeV/c. This difference may suggest that charm quarks interact less strongly with the surrounding nuclear matter at those two lower energies compared to $\sqrt{s_{\rm NN}} = 200$ GeV.

The acceptance-corrected dielectron excess mass spectra, where the known hadronic sources have been subtracted from the inclusive dielectron mass spectra, are reported for the first time at mid-rapidity $|y_{ee}|<1$ in minimum-bias Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 200 GeV. The excess mass spectra are consistently described by a model calculation with a broadened $\rho$ spectral function for $M_{ee}<1.1$ GeV/$c^{2}$. The integrated dielectron excess yield at $\sqrt{s_{NN}}$ = 19.6 GeV for $0.4<M_{ee}<0.75$ GeV/$c^2$, normalized to the charged particle multiplicity at mid-rapidity, has a value similar to that in In+In collisions at $\sqrt{s_{NN}}$ = 17.3 GeV. For $\sqrt{s_{NN}}$ = 200 GeV, the normalized excess yield in central collisions is higher than that at $\sqrt{s_{NN}}$ = 17.3 GeV and increases from peripheral to central collisions. These measurements indicate that the lifetime of the hot, dense medium created in central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV is longer than those in peripheral collisions and at lower energies.

A search for the quantum chromodynamics (QCD) critical point was performed by the STAR experiment at the Relativistic Heavy Ion Collider, using dynamical fluctuations of unlike particle pairs. Heavy-ion collisions were studied over a large range of collision energies with homogeneous acceptance and excellent particle identification, covering a significant range in the QCD phase diagram where a critical point may be located. Dynamical $K\pi$, $p\pi$, and $Kp$ fluctuations as measured by the STAR experiment in central 0-5\% Au+Au collisions from center-of-mass collision energies $\rm \sqrt{s_{NN}}$ = 7.7 to 200 GeV are presented. The observable $\rm \nu_{dyn}$ was used to quantify the magnitude of the dynamical fluctuations in event-by-event measurements of the $K\pi$, $p\pi$, and $Kp$ pairs. The energy dependences of these fluctuations from central 0-5\% Au+Au collisions all demonstrate a smooth evolution with collision energy.

We report measurements of the nuclear modification factor, $R_{ \mathrm{CP}}$, for charged hadrons as well as identified $\pi^{+(-)}$, $K^{+(-)}$, and $p(\overline{p})$ for Au+Au collision energies of $\sqrt{s_{_{ \mathrm{NN}}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, and 62.4 GeV. We observe a clear high-$p_{\mathrm{T}}$ net suppression in central collisions at 62.4 GeV for charged hadrons which evolves smoothly to a large net enhancement at lower energies. This trend is driven by the evolution of the pion spectra, but is also very similar for the kaon spectra. While the magnitude of the proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ does depend on collision energy, neither the proton nor the anti-proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ exhibit net suppression at any energy. A study of how the binary collision scaled high-$p_{\mathrm{T}}$ yield evolves with centrality reveals a non-monotonic shape that is consistent with the idea that jet-quenching is increasing faster than the combined phenomena that lead to enhancement.

New measurements of directed flow for charged hadrons, characterized by the Fourier coefficient \vone, are presented for transverse momenta $\mathrm{p_T}$, and centrality intervals in Au+Au collisions recorded by the STAR experiment for the center-of-mass energy range $\mathrm{\sqrt{s_{_{NN}}}} = 7.7 - 200$ GeV. The measurements underscore the importance of momentum conservation and the characteristic dependencies on $\mathrm{\sqrt{s_{_{NN}}}}$, centrality and $\mathrm{p_T}$ are consistent with the expectations of geometric fluctuations generated in the initial stages of the collision, acting in concert with a hydrodynamic-like expansion. The centrality and $\mathrm{p_T}$ dependencies of $\mathrm{v^{even}_{1}}$, as well as an observed similarity between its excitation function and that for $\mathrm{v_3}$, could serve as constraints for initial-state models. The $\mathrm{v^{even}_{1}}$ excitation function could also provide an important supplement to the flow measurements employed for precision extraction of the temperature dependence of the specific shear viscosity.

We report a systematic measurement of cumulants, $C_{n}$, for net-proton, proton and antiproton multiplicity distributions, and correlation functions, $\kappa_n$, for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{\mathrm {NN}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4 and 200 GeV. The $C_{n}$ and $\kappa_n$ are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, $y$, and transverse momentum, $p_{T}$. The data were taken during the first phase of the Beam Energy Scan (BES) program (2010 -- 2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity ($|y| <$ 0.5) and transverse momentum 0.4 $<$$p_{\rm T}$$<$ 2.0 GeV/$c$, using the STAR detector at RHIC. We observe a non-monotonic energy dependence ($\sqrt{s_{\mathrm {NN}}}$ = 7.7 -- 62.4 GeV) of the net-proton $C_{4}$/$C_{2}$ with the significance of 3.1$\sigma$ for the 0-5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{\mathrm {NN}}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, $\kappa_2$, of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, $\kappa_4$, of protons plays a role in determining the energy dependence of proton $C_4/C_1$ below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.

Rapidity-odd directed flow measurements at midrapidity are presented for $\Lambda$, $\bar{\Lambda}$, $K^\pm$, $K^0_s$ and $\phi$ at $\sqrt{s_{NN}} =$ 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV in Au+Au collisions recorded by the STAR detector at the Relativistic Heavy Ion Collider. These measurements greatly expand the scope of data available to constrain models with differing prescriptions for the equation of state of quantum chromodynamics. Results show good sensitivity for testing a picture where flow is assumed to be imposed before hadron formation and the observed particles are assumed to form via coalescence of constituent quarks. The pattern of departure from a coalescence-inspired sum-rule can be a valuable new tool for probing the collision dynamics.

We report the first measurements of a complete second-order cumulant matrix of net-charge, net-proton, and net-kaon multiplicity distributions for the first phase of the beam energy scan program at RHIC. This includes the centrality and, for the first time, the pseudorapidity window dependence of both diagonal and off-diagonal cumulants in Au+Au collisions at \sNN~= 7.7-200 GeV. Within the available acceptance of $|\eta|<0.5$, the cumulants grow linearly with the pseudorapidity window. Relative to the corresponding measurements in peripheral collisions, the ratio of off-diagonal over diagonal cumulants in central collisions indicates an excess correlation between net-charge and net-kaon, as well as between net-charge and net-proton. The strength of such excess correlation increases with the collision energy. The correlation between net-proton and net-kaon multiplicity distributions is observed to be negative at \sNN~= 200 GeV and change to positive at the lowest collision energy. Model calculations based on non-thermal (UrQMD) and thermal (HRG) production of hadrons cannot explain the data. These measurements will help map the QCD phase diagram, constrain hadron resonance gas model calculations, and provide new insights on the energy dependence of baryon-strangeness correlations. An erratum has been added to address the issue of self-correlation in the previously considered efficiency correction for off-diagonal cumulant measurement. Previously considered unidentified (net-)charge correlation results ($\sigma^{11}_{Q,p}$ and $\sigma^{11}_{Q,k})$ are now replaced with identified (net-)charge correlation ($\sigma^{11}_{Q^{PID},p}$ and $\sigma^{11}_{Q^{PID},k}$)

We report the measurement of $K^{*0}$ meson at midrapidity ($|y|<$ 1.0) in Au+Au collisions at $\sqrt{s_{\rm NN}}$~=~7.7, 11.5, 14.5, 19.6, 27 and 39 GeV collected by the STAR experiment during the RHIC beam energy scan (BES) program. The transverse momentum spectra, yield, and average transverse momentum of $K^{*0}$ are presented as functions of collision centrality and beam energy. The $K^{*0}/K$ yield ratios are presented for different collision centrality intervals and beam energies. The $K^{*0}/K$ ratio in heavy-ion collisions are observed to be smaller than that in small system collisions (e+e and p+p). The $K^{*0}/K$ ratio follows a similar centrality dependence to that observed in previous RHIC and LHC measurements. The data favor the scenario of the dominance of hadronic re-scattering over regeneration for $K^{*0}$ production in the hadronic phase of the medium.

Invariant cross-sections are presented for the inclusive reaction p + p → π o + anything, Measurements of large transverse momentum π o 's (2.5 GeV/ c < p ⊥ <9 GeV/ c ) were made near 90° at the CERN ISR at five centre-of-mass energies (√ s = 23.5, 30.6, 44.8, 52.7 and 62.4 GeV. At large p ⊥ , the invariant cross-sections are seem to vary with s and p ⊥ , in good agreement with a fit of the form Ap ⊥ − n F ( p ⊥ /√ s ), with n ≈8 and F ( p ⊥ /√ s )≈ exp (−26 p ⊥ /√ s ).