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.

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.

We present three-particle mixed-harmonic correlations $\la \cos (m\phi_a + n\phi_b - (m+n) \phi_c)\ra$ for harmonics $m,n=1-3$ for charged particles in $\sqrt{s_{NN}}=$200 GeV Au+Au collisions at RHIC. These measurements provide information on the three-dimensional structure of the initial collision zone and are important for constraining models of a subsequent low-viscosity quark-gluon plasma expansion phase. We investigate correlations between the first, second and third harmonics predicted as a consequence of fluctuations in the initial state. The dependence of the correlations on the pseudorapidity separation between particles show hints of a breaking of longitudinal invariance. We compare our results to a number of state-of-the art hydrodynamic calculations with different initial states and temperature dependent viscosities. These measurements provide important steps towards constraining the temperature dependent transport and the longitudinal structure of the initial state at RHIC.

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.

Measurements of two-particle correlations on transverse momentum $p_t$ for Au-Au collisions at $\sqrt{s_{NN}} = 130$ GeV are presented. Significant large-momentum-scale correlations are observed for charged primary hadrons with $0.15 \leq p_t \leq 2$ GeV/$c$ and pseudorapidity $|\eta| \leq 1.3$. Such correlations were not observed in a similar study at lower energy and are not predicted by theoretical collision models. Their direct relation to mean-$p_t$ fluctuations measured in the same angular acceptance is demonstrated. Positive correlations are observed for pairs of particles which have large $p_t$ values while negative correlations occur for pairs in which one particle has large $p_t$ and the other has much lower $p_t$. The correlation amplitudes per final state particle increase with collision centrality. The observed correlations are consistent with a scenario in which the transverse momentum of hadrons associated with initial-stage semi-hard parton scattering is dissipated by the medium to lower $p_t$.

We report measurements of single-particle inclusive spectra and two-particle azimuthal distributions of charged hadrons at high transverse momentum (high $p_T$) in minimum bias and central d+Au collisions at $\sqrt{s_{NN}}$=200 GeV. The inclusive yield is enhanced in d+Au collisions relative to binary-scaled p+p collisions, while the two-particle azimuthal distributions are very similar to those observed in p+p collisions. These results demonstrate that the strong suppression of the inclusive yield and back-to-back correlations at high $p_T$ previously observed in central Au+Au collisions are due to final-state interactions with the dense medium generated in such collisions.

In relativistic heavy-ion collisions, a global spin polarization, $P_\mathrm{H}$, of $\Lambda$ and $\bar{\Lambda}$ hyperons along the direction of the system angular momentum was discovered and measured across a broad range of collision energies and demonstrated a trend of increasing $P_\mathrm{H}$ with decreasing $\sqrt{s_{NN}}$. A splitting between $\Lambda$ and $\bar{\Lambda}$ polarization may be possible due to their different magnetic moments in a late-stage magnetic field sustained by the quark-gluon plasma which is formed in the collision. The results presented in this study find no significant splitting at the collision energies of $\sqrt{s_{NN}}=19.6$ and $27$ GeV in the RHIC Beam Energy Scan Phase II using the STAR detector, with an upper limit of $P_{\bar{\Lambda}}-P_{\Lambda}<0.24$% and $P_{\bar{\Lambda}}-P_{\Lambda}<0.35$%, respectively, at a 95% confidence level. We derive an upper limit on the na\"ive extraction of the late-stage magnetic field of $B<9.4\cdot10^{12}$ T and $B<1.4\cdot10^{13}$ T at $\sqrt{s_{NN}}=19.6$ and $27$ GeV, respectively, although more thorough derivations are needed. Differential measurements of $P_\mathrm{H}$ were performed with respect to collision centrality, transverse momentum, and rapidity. With our current acceptance of $|y|<1$ and uncertainties, we observe no dependence on transverse momentum and rapidity in this analysis. These results challenge multiple existing model calculations following a variety of different assumptions which have each predicted a strong dependence on rapidity in this collision-energy range.

Results on high transverse momentum charged particle emission with respect to the reaction plane are presented for Au+Au collisions at $\sqrt{s_{_{NN}}}$= 200 GeV. Two- and four-particle correlations results are presented as well as a comparison of azimuthal correlations in Au+Au collisions to those in $p+p$ at the same energy. Elliptic anisotropy, $v_2$, is found to reach its maximum at $p_t \sim 3$ GeV/c, then decrease slowly and remain significant up to $p_t\approx 7$ -- 10 GeV/c. Stronger suppression is found in the back-to-back high-$p_t$ particle correlations for particles emitted out-of-plane compared to those emitted in-plane. The centrality dependence of $v_2$ at intermediate $p_t$ is compared to simple models based on jet quenching.

We present the first measurements of charge-dependent correlations on angular difference variables $\eta_1 - \eta_2$ (pseudorapidity) and $\phi_1 - \phi_2$ (azimuth) for primary charged hadrons with transverse momentum $0.15 \leq p_t \leq 2$ GeV/$c$ and $|\eta| \leq 1.3$ from Au-Au collisions at $\sqrt{s_{NN}} = 130$ GeV. We observe correlation structures not predicted by theory but consistent with evolution of hadron emission geometry with increasing centrality from one-dimensional fragmentation of color strings along the beam direction to an at least two-dimensional hadronization geometry along the beam and azimuth directions of a hadron-opaque bulk medium.

First measurements of charge-independent correlations on angular difference variables $\eta_1 - \eta_2$ (pseudorapidity) and $\phi_1 - \phi_2$ (azimuth) are presented for primary charged hadrons with transverse momentum $0.15 \leq p_t \leq 2$ GeV/$c$ and $|\eta| \leq 1.3$ from Au-Au collisions at $\sqrt{s_{NN}} = 130$ GeV. Strong charge-independent angular correlations are observed associated with jet-like structures and elliptic flow. The width of the jet-like peak on $\eta_1 - \eta_2$ increases by a factor 2.3 from peripheral to central collisions, suggesting strong coupling of semi-hard scattered partons to a longitudinally-expanding medium. New methods of jet analysis introduced here provide evidence for nonperturbative QCD medium effects in heavy ion collisions.