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.

The polarization of $\Lambda$ and $\bar{\Lambda}$ hyperons along the beam direction has been measured relative to the second and third harmonic event planes in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. This is the first experimental evidence of the hyperon polarization by the triangular flow originating from the initial density fluctuations. The amplitudes of the sine modulation for the second and third harmonic results are comparable in magnitude, increase from central to peripheral collisions, and show a mild $p_T$ dependence. The azimuthal angle dependence of the polarization follows the vorticity pattern expected due to elliptic and triangular anisotropic flow, and qualitatively disagree with most hydrodynamic model calculations based on thermal vorticity and shear induced contributions. The model results based on one of existing implementations of the shear contribution lead to a correct azimuthal angle dependence, but predict centrality and $p_T$ dependence that still disagree with experimental measurements. Thus, our results provide stringent constraints on the thermal vorticity and shear-induced contributions to hyperon polarization. Comparison to previous measurements at RHIC and the LHC for the second-order harmonic results shows little dependence on the collision system size and collision energy.

We report a measurement of cumulants and correlation functions of event-by-event proton multiplicity distributions from fixed-target Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV measured by the STAR experiment. Protons are identified within the rapidity ($y$) and transverse momentum ($p_{\rm T}$) region $-0.9 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$ in the center-of-mass frame. A systematic analysis of the proton cumulants and correlation functions up to sixth-order as well as the corresponding ratios as a function of the collision centrality, $p_{\rm T}$, and $y$ are presented. The effect of pileup and initial volume fluctuations on these observables and the respective corrections are discussed in detail. The results are compared to calculations from the hadronic transport UrQMD model as well as a hydrodynamic model. In the most central 5% collisions, the value of proton cumulant ratio $C_4/C_2$ is negative, drastically different from the values observed in Au+Au collisions at higher energies. Compared to model calculations including Lattice QCD, a hadronic transport model, and a hydrodynamic model, the strong suppression in the ratio of $C_4/C_2$ at 3 GeV Au+Au collisions indicates an energy regime dominated by hadronic interactions.

Global polarizations ($P$) of $\Lambda$ ($\bar{\Lambda}$) hyperons have been observed in non-central heavy-ion collisions. The strong magnetic field primarily created by the spectator protons in such collisions would split the $\Lambda$ and $\bar{\Lambda}$ global polarizations ($\Delta P = P_{\Lambda} - P_{\bar{\Lambda}} < 0$). Additionally, quantum chromodynamics (QCD) predicts topological charge fluctuations in vacuum, resulting in a chirality imbalance or parity violation in a local domain. This would give rise to an imbalance ($\Delta n = \frac{N_{\text{L}} - N_{\text{R}}}{\langle N_{\text{L}} + N_{\text{R}} \rangle} \neq 0$) between left- and right-handed $\Lambda$ ($\bar{\Lambda}$) as well as a charge separation along the magnetic field, referred to as the chiral magnetic effect (CME). This charge separation can be characterized by the parity-even azimuthal correlator ($\Delta\gamma$) and parity-odd azimuthal harmonic observable ($\Delta a_{1}$). Measurements of $\Delta P$, $\Delta\gamma$, and $\Delta a_{1}$ have not led to definitive conclusions concerning the CME or the magnetic field, and $\Delta n$ has not been measured previously. Correlations among these observables may reveal new insights. This paper reports measurements of correlation between $\Delta n$ and $\Delta a_{1}$, which is sensitive to chirality fluctuations, and correlation between $\Delta P$ and $\Delta\gamma$ sensitive to magnetic field in Au+Au collisions at 27 GeV. For both measurements, no correlations have been observed beyond statistical fluctuations.

The longitudinal and transverse spin transfers to $\Lambda$ ($\overline{\Lambda}$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $\Lambda$ and $\overline{\Lambda}$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|\eta_{\Lambda(\overline{\Lambda})}|$$<$ 1.2 and transverse momentum $p_{T,{\Lambda(\overline{\Lambda})}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions.

We report on measurements of sequential $\Upsilon$ suppression in Au+Au collisions at $\sqrt{s_{_\mathrm{NN}}}$ = 200 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC) through both the dielectron and dimuon decay channels. In the 0-60% centrality class, the nuclear modification factors ($R_{\mathrm{AA}}$), which quantify the level of yield suppression in heavy-ion collisions compared to $p$+$p$ collisions, for $\Upsilon$(1S) and $\Upsilon$(2S) are $0.40 \pm 0.03~\textrm{(stat.)} \pm 0.03~\textrm{(sys.)} \pm 0.09~\textrm{(norm.)}$ and $0.26 \pm 0.08~\textrm{(stat.)} \pm 0.02~\textrm{(sys.)} \pm 0.06~\textrm{(norm.)}$, respectively, while the upper limit of the $\Upsilon$(3S) $R_{\mathrm{AA}}$ is 0.17 at a 95% confidence level. This provides experimental evidence that the $\Upsilon$(3S) is significantly more suppressed than the $\Upsilon$(1S) at RHIC. The level of suppression for $\Upsilon$(1S) is comparable to that observed at the much higher collision energy at the Large Hadron Collider. These results point to the creation of a medium at RHIC whose temperature is sufficiently high to strongly suppress excited $\Upsilon$ states.

For the search of the chiral magnetic effect (CME), STAR previously presented the results from isobar collisions (${^{96}_{44}\text{Ru}}+{^{96}_{44}\text{Ru}}$, ${^{96}_{40}\text{Zr}}+{^{96}_{40}\text{Zr}}$) obtained through a blind analysis. The ratio of results in Ru+Ru to Zr+Zr collisions for the CME-sensitive charge-dependent azimuthal correlator ($\Delta\gamma$), normalized by elliptic anisotropy ($v_{2}$), was observed to be close to but systematically larger than the inverse multiplicity ratio. The background baseline for the isobar ratio, $Y = \frac{(\Delta\gamma/v_{2})^{\text{Ru}}}{(\Delta\gamma/v_{2})^{\text{Zr}}}$, is naively expected to be $\frac{(1/N)^{\text{Ru}}}{(1/N)^{\text{Zr}}}$; however, genuine two- and three-particle correlations are expected to alter it. We estimate the contributions to $Y$ from those correlations, utilizing both the isobar data and HIJING simulations. After including those contributions, we arrive at a final background baseline for $Y$, which is consistent with the isobar data. We extract an upper limit for the CME fraction in the $\Delta\gamma$ measurement of approximately $10\%$ at a $95\%$ confidence level on in isobar collisions at $\sqrt{s_{\text{NN}}} = 200$ GeV, with an expected $15\%$ difference in their squared magnetic fields.

We report a new measurement of the production of electrons from open heavy-flavor hadron decays (HFEs) at mid-rapidity ($|y|<$ 0.7) in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Invariant yields of HFEs are measured for the transverse momentum range of $3.5 < p_{\rm T} < 9$ GeV/$c$ in various configurations of the collision geometry. The HFE yields in head-on Au+Au collisions are suppressed by approximately a factor of 2 compared to that in $p$+$p$ collisions scaled by the average number of binary collisions, indicating strong interactions between heavy quarks and the hot and dense medium created in heavy-ion collisions. Comparison of these results with models provides additional tests of theoretical calculations of heavy quark energy loss in the quark-gluon plasma.

Azimuthal anisotropy of produced particles is one of the most important observables used to access the collective properties of the expanding medium created in relativistic heavy-ion collisions. In this paper, we present second ($v_{2}$) and third ($v_{3}$) order azimuthal anisotropies of $K_{S}^{0}$, $\phi$, $\Lambda$, $\Xi$ and $\Omega$ at mid-rapidity ($|y|<$1) in Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 54.4 GeV measured by the STAR detector. The $v_{2}$ and $v_{3}$ are measured as a function of transverse momentum and centrality. Their energy dependence is also studied. $v_{3}$ is found to be more sensitive to the change in the center-of-mass energy than $v_{2}$. Scaling by constituent quark number is found to hold for $v_{2}$ within 10%. This observation could be evidence for the development of partonic collectivity in 54.4 GeV Au+Au collisions. Differences in $v_{2}$ and $v_{3}$ between baryons and anti-baryons are presented, and ratios of $v_{3}$/$v_{2}^{3/2}$ are studied and motivated by hydrodynamical calculations. The ratio of $v_{2}$ of $\phi$ mesons to that of anti-protons ($v_{2}(\phi)/v_{2}(\bar{p})$) shows centrality dependence at low transverse momentum, presumably resulting from the larger effects from hadronic interactions on anti-proton $v_{2}$.

We report the beam energy and collision centrality dependence of fifth and sixth order cumulants ($C_{5}$, $C_{6}$) and factorial cumulants ($\kappa_{5}$, $\kappa_{6}$) of net-proton and proton distributions, from $\sqrt{s_{NN}} = 3 - 200$ GeV Au+Au collisions at RHIC. The net-proton cumulant ratios generally follow the hierarchy expected from QCD thermodynamics, except for the case of collisions at $\sqrt{s_{NN}}$ = 3 GeV. $C_{6}/C_{2}$ for 0-40% centrality collisions is increasingly negative with decreasing $\sqrt{s_{NN}}$, while it is positive for the lowest $\sqrt{s_{NN}}$ studied. These observed negative signs are consistent with QCD calculations (at baryon chemical potential, $\mu_{B} \leq$ 110 MeV) that include a crossover quark-hadron transition. In addition, for $\sqrt{s_{NN}} \geq$ 11.5 GeV, the measured proton $\kappa_{n}$, within uncertainties, does not support the two-component shape of proton distributions that would be expected from a first-order phase transition. Taken in combination, the hyper-order proton number fluctuations suggest that the structure of QCD matter at high baryon density, $\mu_{B}\sim 750$ MeV ($\sqrt{s_{NN}}$ = 3 GeV) is starkly different from those at vanishing $\mu_{B}\sim 20$MeV ($\sqrt{s_{NN}}$ = 200 GeV and higher).