We present a measurement of the first-order azimuthal anisotropy, $v_1(y)$, of deuterons from Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV recorded with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The energy dependence of the $v_1(y)$ slope, $dv_{1}/dy|_{y=0}$, for deuterons, where $y$ is the rapidity, is extracted for semi-central collisions (10-40\% centrality) and compared to that of protons. While the $v_1(y)$ slopes of protons are generally negative for $\sqrt{s_{NN}} >$ 10 GeV, those for deuterons are consistent with zero, a strong enhancement of the $v_1(y)$ slope of deuterons is seen at the lowest collision energy (the largest baryon density) at $\sqrt{s_{NN}} =$ 7.7 GeV. In addition, we report the transverse momentum dependence of $v_1$ for protons and deuterons. The experimental results are compared with transport and coalescence models.

Global hyperon polarization, $\overline{P}_\mathrm{H}$, in Au+Au collisions over a large range of collision energy, $\sqrt{s_\mathrm{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While naïve extrapolation of data trends suggests a large $\overline{P}_\mathrm{H}$ as the collision energy is reduced, the behavior of $\overline{P}_\mathrm{H}$ at small $\sqrt{s_\mathrm{NN}}<7.7$ GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of $\Lambda$ hyperons along the direction of global angular momentum in Au+Au collisions at $\sqrt{s_\mathrm{NN}}=3$ GeV. The observation of substantial polarization of $4.91\pm0.81(\rm stat.)\pm0.15(\rm syst.)$% in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

The Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider (RHIC) was extended to energies below $\sqrt{\textit{s}_{NN}}$ = 7.7 GeV in 2015 by successful implementation of the fixed-target mode of operation in the STAR (Solenoidal Track At RHIC) experiment. In the fixed-target mode, ions circulate in one ring of the collider and interact with a stationary target at the entrance of the STAR Time Projection Chamber. The first results for Au+Au collisions at $\sqrt{\textit{s}_{NN}}$ = 4.5 GeV are presented, including directed and elliptic flow of identified hadrons, and radii from pion femtoscopy. The proton flow and pion femtoscopy results agree quantitatively with earlier measurements by Alternating Gradient Synchrotron experiments at similar energies. This validates running the STAR experiment in the fixed-target configuration. Pion directed and elliptic flow are presented for the first time at this beam energy. Pion and proton elliptic flow show behavior which hints at constituent quark scaling, but large error bars preclude reliable conclusions. The ongoing second phase of BES (BES-II) will provide fixed-target data sets with 100 times more events at each of several energies down to $\sqrt{\textit{s}_{NN}}$ = 3.0 GeV.

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, $v_1$ and $v_2$, of light nuclei ($d$, $t$, $^{3}$He, $^{4}$He) produced in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured $v_1$ slopes of light nuclei at mid-rapidity. For the measured $v_2$ magnitude, a strong rapidity dependence is observed. Unlike $v_2$ at higher collision energies, the $v_2$ values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

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 report on the measurements of directed flow $v_1$ and elliptic flow $v_2$ for hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{S}^0$, $p$, $\phi$, $\Lambda$ and $\Xi^{-}$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3 GeV and $v_{2}$ for ($\pi^{\pm}$, $K^{\pm}$, $p$ and $\overline{p}$) at 27 and 54.4 GeV with the STAR experiment. While at the two higher energy midcentral collisions the number-of-constituent-quark (NCQ) scaling holds, at 3 GeV the $v_{2}$ at midrapidity is negative for all hadrons and the NCQ scaling is absent. In addition, the $v_1$ slopes at midrapidity for almost all observed hadrons are found to be positive, implying dominant repulsive baryonic interactions. The features of negative $v_2$ and positive $v_1$ slope at 3 GeV can be reproduced with a baryonic mean-field in transport model calculations. These results imply that the medium in such collisions is likely characterized by baryonic interactions.

We report on the first measurements of J/$\psi$ production at very low transverse momentum ($p_{T} <$ 0.2 GeV/c) in hadronic Au+Au collisions at $\sqrt{s_{\rm{NN}}} =$ 200 GeV and U+U collisions at $\sqrt{s_{\rm{NN}}} =$ 193 GeV. Remarkably, the inferred nuclear modification factor of J/$\psi$ at mid-rapidity in Au+Au (U+U) collisions reaches about 24 (52) for $p_{T} <$ 0.05 GeV/c in the 60-80$\%$ collision centrality class. This noteworthy enhancement cannot be explained by hadronic production accompanied by cold and hot medium effects. In addition, the $dN/dt$ distribution of J/$\psi$ for the very low $p_{T}$ range is presented for the first time. The distribution is consistent with that expected from the Au nucleus and shows a hint of interference. Comparison of the measurements to theoretical calculations of coherent production shows that the excess yield can be described reasonably well and reveals a partial disruption of coherent production in semi-central collisions, perhaps due to the violent hadronic interactions. Incorporating theoretical calculations, the results strongly suggest that the dramatic enhancement of J/$\psi$ yield observed at extremely low $p_{T}$ originates from coherent photon-nucleus interactions. In particular, coherently produced J/$\psi$'s in violent hadronic collisions may provide a novel probe of the quark-gluon-plasma.

Partons traversing the strongly interacting medium produced in heavy-ion collisions are expected to lose energy depending on their color charge and mass. We measure the nuclear modification factors for charm- and bottom-decay electrons, defined as the ratio of yields, scaled by the number of binary nucleon-nucleon collisions, in $\sqrt{s_{\rm NN}}$ = 200 GeV Au+Au collisions to $p$+$p$ collisions ($R_{\rm AA}$), or in central to peripheral Au+Au collisions ($R_{\rm CP}$). We find the bottom-decay electron $R_{\rm AA}$ and $R_{\rm CP}$ to be significantly higher than that of charm-decay electrons. Model calculations including mass-dependent parton energy loss in a strongly coupled medium are consistent with the measured data. These observations provide clear evidence of mass ordering of charm and bottom quark energy loss when traversing through the strongly coupled medium created in heavy-ion collisions.

Flow harmonics ($\textit{v}_{n}$) of the Fourier expansion for the azimuthal distributions of hadrons are commonly employed to quantify the azimuthal anisotropy of particle production relative to the collision symmetry planes. While lower order Fourier coefficients ($\textit{v}_{2}$ and $\textit{v}_{3}$) are more directly related to the corresponding eccentricities of the initial state, the higher-order flow harmonics ($\textit{v}_{n>3}$) can be induced by a mode-coupled response to the lower-order anisotropies, in addition to a linear response to the same-order anisotropies. These higher-order flow harmonics and their linear and mode-coupled contributions can be used to more precisely constrain the initial conditions and the transport properties of the medium in theoretical models. The multiparticle azimuthal cumulant method is used to measure the linear and mode-coupled contributions in the higher-order anisotropic flow, the mode-coupled response coefficients, and the correlations of the event plane angles for charged particles as functions of centrality and transverse momentum in Au+Au collisions at nucleon-nucleon center-of-mass energy \roots = 200 GeV. The results are compared to similar LHC measurements as well as to several viscous hydrodynamic calculations with varying initial conditions.

Global polarization of $\Xi$ and $\Omega$ hyperons has been measured for the first time in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The measurements of the $\Xi^-$ and $\bar{\Xi}^+$ hyperon polarization have been performed by two independent methods, via analysis of the angular distribution of the daughter particles in the parity violating weak decay $\Xi\rightarrow\Lambda+\pi$, as well as by measuring the polarization of the daughter $\Lambda$-hyperon, polarized via polarization transfer from its parent. The polarization, obtained by combining the results from the two methods and averaged over $\Xi^-$ and $\bar{\Xi}^+$, is measured to be $\langle P_\Xi \rangle = 0.47\pm0.10~({\rm stat.})\pm0.23~({\rm syst.})\,\%$ for the collision centrality 20%-80%. The $\langle P_\Xi \rangle$ is found to be slightly larger than the inclusive $\Lambda$ polarization and in reasonable agreement with a multi-phase transport model (AMPT). The $\langle P_\Xi \rangle$ is found to follow the centrality dependence of the vorticity predicted in the model, increasing toward more peripheral collisions. The global polarization of $\Omega$, $\langle P_\Omega \rangle = 1.11\pm0.87~({\rm stat.})\pm1.97~({\rm syst.})\,\%$ was obtained by measuring the polarization of daughter $\Lambda$ in the decay $\Omega \rightarrow \Lambda + K$, assuming the polarization transfer factor $C_{\Omega\Lambda}=1$.