We measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au+Au collisions at the BNL Relativistic Heavy Ion Collider. A significant $v_3$ signal for protons is observed, which increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a slope at mid-rapidity, $dv_3/dy|_{(y=0)}$, opposite in sign compared to the slope for directed flow. No significant $v_3$ signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

Angular distributions of charged particles relative to jet axes are studied in $\sqrt{s_{\mathrm{NN}}}$ = 200 GeV Au+Au collisions as a function of the jet orientation with respect to the event plane. This differential study tests the expected path-length dependence of energy loss experienced by a hard-scattered parton as it traverses the hot and dense medium formed in heavy-ion collisions. A second-order event plane is used in the analysis as an experimental estimate of the reaction plane formed by the collision impact parameter and the beam direction. Charged-particle jets with $15 < p_{\rm T, jet} <$ 20 and $20 < p_{\rm T, jet} <$ 40 GeV/$c$ were reconstructed with the anti-$k_{\rm T}$ algorithm with radius parameter setting of (R=0.4) in the 20-50% centrality bin to maximize the initial-state eccentricity of the interaction region. The reaction plane fit method is implemented to remove the flow-modulated background with better precision than prior methods. Yields and widths of jet-associated charged-hadron distributions are extracted in three angular bins between the jet axis and the event plane. The event-plane (EP) dependence is further quantified by ratios of the associated yields in different EP bins. No dependence on orientation of the jet axis with respect to the event plane is seen within the uncertainties in the kinematic regime studied. This finding is consistent with a similar experimental observation by ALICE in $\sqrt{s_{\mathrm{NN}}}$ = 2.76 TeV Pb+Pb collision data.

Quark interactions with topological gluon configurations can induce local chirality imbalance and parity violation in quantum chromodynamics, which can lead to the chiral magnetic effect (CME) -- an electric charge separation along the strong magnetic field in relativistic heavy-ion collisions. The CME-sensitive azimuthal correlator observable ($\Delta\gamma$) is contaminated by background arising, in part, from resonance decays coupled with elliptic anisotropy ($v_{2}$). We report here differential measurements of the correlator as a function of the pair invariant mass ($m_{\rm inv}$) in 20-50% centrality Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$= 200 GeV by the STAR experiment at RHIC. Strong resonance background contributions to $\Delta\gamma$ are observed. At large $m_{\rm inv}$ where this background is significantly reduced, the $\Delta\gamma$ value is found to be significantly smaller. An event-shape-engineering technique is deployed to determine the $v_{2}$ background shape as a function of $m_{\rm inv}$. We extract a $v_2$-independent and $m_{\rm inv}$-averaged signal $\Delta\gamma_{\rm sig}$ = (0.03 $\pm$ 0.06 $\pm$ 0.08) $\times10^{-4}$, or $(2\pm4\pm5)\%$ of the inclusive $\Delta\gamma(m_{\rm inv}>0.4$ GeV/$c^2$)$ =(1.58 \pm 0.02 \pm 0.02) \times10^{-4}$, within pion $p_{T}$ = 0.2 - 0.8~\gevc and averaged over pseudorapidity ranges of $-1 < \eta < -0.05$ and $0.05 < \eta < 1$. This represents an upper limit of $0.23\times10^{-4}$, or $15\%$ of the inclusive result, at $95\%$ confidence level for the $m_{\rm inv}$-integrated CME contribution.

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).

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.

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.

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.

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 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.