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.

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.

Density fluctuations near the QCD critical point can be probed via an intermittency analysis in relativistic heavy-ion collisions. We report the first measurement of intermittency in Au$+$Au collisions at $\sqrt{s_\mathrm{_{NN}}}$ = 7.7-200 GeV measured by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The scaled factorial moments of identified charged hadrons are analyzed at mid-rapidity and within the transverse momentum phase space. We observe a power-law behavior of scaled factorial moments in Au$+$Au collisions and a decrease in the extracted scaling exponent ($\nu$) from peripheral to central collisions. The $\nu$ is consistent with a constant for different collisions energies in the mid-central (10-40%) collisions. Moreover, the $\nu$ in the 0-5% most central Au$+$Au collisions exhibits a non-monotonic energy dependence that reaches a possible minimum around $\sqrt{s_\mathrm{_{NN}}}$ = 27 GeV. The physics implications on the QCD phase structure are discussed.

We report the energy dependence of mid-rapidity (anti-)deuteron production in Au+Au collisions at $\sqrt{s_\text{NN}} =\ $7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV, measured by the STAR experiment at RHIC. The yield of deuterons is found to be well described by the thermal model. The collision energy, centrality, and transverse momentum dependence of the coalescence parameter $B_2$ are discussed. We find that the values of $B_2$ for anti-deuterons are systematically lower than those for deuterons, indicating that the correlation volume of anti-baryons is larger than that of baryons at $\sqrt{s_\text{NN}}$ from 19.6 to 39 GeV. In addition, values of $B_2$ are found to vary with collision energy and show a broad minimum around $\sqrt{s_\text{NN}}=\ $20 to 40 GeV, which might imply a change of the equation of state of the medium in these collisions.

We present two-particle $p_{\rm t}$ correlations as a function of event centrality for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV at the Relativistic Heavy Ion Collider using the STAR detector. These results are compared to previous measurements from CERES at the Super Proton Synchrotron and from ALICE at the Large Hadron Collider. The data are compared with UrQMD model calculations and with a model based on a Boltzmann-Langevin approach incorporating effects from thermalization. The relative dynamical correlations for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV show a power law dependence on the number of participant nucleons and agree with the results for Pb+Pb collisions at $\sqrt{s_{\rm NN}} = 2.76~ {\rm TeV}$ from ALICE. As the collision energy is lowered from $\sqrt{s_{\rm NN}}$ = 200 GeV to 7.7 GeV, the centrality dependence of the relative dynamical correlations departs from the power law behavior observed at the higher collision energies. In central collisions, the relative dynamical correlations increase with collision energy up to $\sqrt{s_{\rm NN}}$ = 200 GeV in contrast to previous measurements that showed little dependence on the collision energy.

The two-particle angular correlation functions, $R_2$, of pions, kaons, and protons in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV were measured by the STAR experiment at RHIC. These correlations were measured for both like-sign and unlike-sign charge combinations and versus the centrality. The correlations of pions and kaons show the expected near-side ({\it i.e.}, at small relative angles) peak resulting from short-range mechanisms. The amplitudes of these short-range correlations decrease with increasing beam energy. However, the proton correlation functions exhibit strong anticorrelations in the near-side region. This behavior is observed for the first time in an A+A collision system. The observed anticorrelation is $p_{T}$-independent and decreases with increasing beam energy and centrality. The experimental results are also compared to the Monte Carlo models UrQMD, Hijing, and AMPT.

The linear and mode-coupled contributions to higher-order anisotropic flow are presented for Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 27, 39, 54.4, and 200 GeV and compared to similar measurements for Pb+Pb collisions at the Large Hadron Collider (LHC). The coefficients and the flow harmonics' correlations, which characterize the linear and mode-coupled response to the lower-order anisotropies, indicate a beam energy dependence consistent with an influence from the specific shear viscosity ($\eta/s$). In contrast, the dimensionless coefficients, mode-coupled response coefficients, and normalized symmetric cumulants are approximately beam-energy independent, consistent with a significant role from initial-state effects. These measurements could provide unique supplemental constraints to (i) distinguish between different initial-state models and (ii) delineate the temperature ($T$) and baryon chemical potential ($\mu_{B}$) dependence of the specific shear viscosity $\frac{\eta}{s} (T, \mu_B)$.

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.

A systematic study is presented for centrality, transverse momentum ($p_T$) and pseudorapidity ($\eta$) dependence of the inclusive charged hadron elliptic flow ($v_2$) at midrapidity($|\eta| < 1.0$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27 and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and 4-particle cumulants ($v_2{4}$), are presented in order to investigate non-flow correlations and $v_2$ fluctuations. We observe that the difference between $v_2{2}$ and $v_2{4}$ is smaller at the lower collision energies. Values of $v_2$, scaled by the initial coordinate space eccentricity, $v_{2}/\varepsilon$, as a function of $p_T$ are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider ($\sqrt{s_{NN}}$ = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV). The $v_2(p_T)$ values for fixed $p_T$ rise with increasing collision energy within the $p_T$ range studied ($< 2 {\rm GeV}/c$). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of $v_{2}(p_{T})$. We also compare the $v_2$ results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at Beam Energy Scan (BES) energies are 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 bulk properties of the matter produced in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, and 39 GeV using identified hadrons ($\pi^\pm$, $K^\pm$, $p$ and $\bar{p}$) from the STAR experiment in the Beam Energy Scan (BES) Program at the Relativistic Heavy Ion Collider (RHIC). Midrapidity ($|y|<$0.1) results for multiplicity densities $dN/dy$, average transverse momenta $\langle p_T \rangle$ and particle ratios are presented. The chemical and kinetic freeze-out dynamics at these energies are discussed and presented as a function of collision centrality and energy. These results constitute the systematic measurements of bulk properties of matter formed in heavy-ion collisions over a broad range of energy (or baryon chemical potential) at RHIC.

We report systematic measurements of bulk properties of the system created in Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 14.5 GeV recorded by the STAR detector at the Relativistic Heavy Ion Collider (RHIC).The transverse momentum spectra of $\pi^{\pm}$, $K^{\pm}$ and $p(\bar{p})$ are studied at mid-rapidity ($|y| < 0.1$) for nine centrality intervals. The centrality, transverse momentum ($p_T$),and pseudorapidity ($\eta$) dependence of inclusive charged particle elliptic flow ($v_2$), and rapidity-odd charged particles directed flow ($v_{1}$) results near mid-rapidity are also presented. These measurements are compared with the published results from Au+Au collisions at other energies, and from Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 2.76 TeV. The results at $\sqrt{s_{\mathrm{NN}}}$ = 14.5 GeV show similar behavior as established at other energies and fit well in the energy dependence trend. These results are important as the 14.5 GeV energy fills the gap in $\mu_B$, which is of the order of 100 MeV,between $\sqrt{s_{\mathrm{NN}}}$ =11.5 and 19.6 GeV. Comparisons of the data with UrQMD and AMPT models show poor agreement in general.

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.

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 STAR measurements of strange hadron ($\mathrm{K}^{0}_{\mathrm S}$, $\Lambda$, $\overline{\Lambda}$, $\Xi^-$, $\overline{\Xi}^+$, $\Omega^-$, $\overline{\Omega}^+$, and $\phi$) production at mid-rapidity ($|y| < 0.5$) in Au+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 7.7 - 39 GeV from the Beam Energy Scan Program at the Relativistic Heavy Ion Collider (RHIC). Transverse momentum spectra, averaged transverse mass, and the overall integrated yields of these strange hadrons are presented versus the centrality and collision energy. Antibaryon-to-baryon ratios ($\overline{\Lambda}$/$\Lambda$, $\overline{\Xi}^+$/$\Xi^-$, $\overline{\Omega}^+$/$\Omega^-$) are presented as well, and used to test a thermal statistical model and to extract the temperature normalized strangeness and baryon chemical potentials at hadronic freeze-out ($\mu_{B}/T_{\rm ch}$ and $\mu_{S}/T_{\rm ch}$) in central collisions. Strange baryon-to-pion ratios are compared to various model predictions in central collisions for all energies. The nuclear modification factors ($R_{\textrm{CP}}$) and antibaryon-to-meson ratios as a function of transverse momentum are presented for all collision energies. The $\mathrm{K}^{0}_{\mathrm S}$$R_{\textrm{CP}}$ shows no suppression for $p_{\rm T}$ up to 3.5 $\mathrm{GeV} / c$ at energies of 7.7 and 11.5 GeV. The $\overline{\Lambda}$/$\mathrm{K}^{0}_{\mathrm S}$ ratio also shows baryon-to-meson enhancement at intermediate $p_{\rm T}$ ($\approx$2.5 $\mathrm{GeV} / c$) in central collisions at energies above 19.6 GeV. Both observations suggest that there is likely a change of the underlying strange quark dynamics at collision energies below 19.6 GeV.

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.

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 on new measurements of elliptic flow ($v_2$) of electrons from heavy-flavor hadron decays at mid-rapidity ($|y|<0.8$) in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 27 and 54.4 GeV from the STAR experiment. Heavy-flavor decay electrons ($e^{\rm HF}$) in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 54.4 GeV exhibit a non-zero $v_2$ in the transverse momentum ($p_{\rm T}$) region of $p_{\rm T}<$ 2 GeV/$c$ with the magnitude comparable to that at $\sqrt{s_{_{\rm NN}}}=200$ GeV. The measured $e^{\rm HF}$$v_2$ at 54.4 GeV is also consistent with the expectation of their parent charm hadron $v_2$ following number-of-constituent-quark scaling as other light and strange flavor hadrons at this energy. These suggest that charm quarks gain significant collectivity through the evolution of the QCD medium and may reach local thermal equilibrium in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=54.4$ GeV. The measured $e^{\rm HF}$$v_2$ in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=$ 27 GeV is consistent with zero within large uncertainties. The energy dependence of $v_2$ for different flavor particles ($\pi,\phi,D^{0}/e^{\rm HF}$) shows an indication of quark mass hierarchy in reaching thermalization in high-energy nuclear collisions.

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.