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

We report on the first measurement of charm-strange meson $D_s^{\pm}$ production at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV from the STAR experiment. The yield ratio between strange ($D_{s}^{\pm}$) and non-strange ($D^{0}$) open-charm mesons is presented and compared to model calculations. A significant enhancement, relative to a PYTHIA simulation of $p$+$p$ collisions, is observed in the $D_{s}^{\pm}/D^0$ yield ratio in Au+Au collisions over a large range of collision centralities. Model calculations incorporating abundant strange-quark production in the quark-gluon plasma (QGP) and coalescence hadronization qualitatively reproduce the data. The transverse-momentum integrated yield ratio of $D_{s}^{\pm}/D^0$ at midrapidity is consistent with a prediction from a statistical hadronization model with the parameters constrained by the yields of light and strange hadrons measured at the same collision energy. These results suggest that the coalescence of charm quarks with strange quarks in the QGP plays an important role in $D_{s}^{\pm}$ meson production in heavy-ion collisions.

Quark interactions with topological gluon configurations can induce chirality imbalance and local parity violation in quantum chromodynamics. This can lead to electric charge separation along the strong magnetic field in relativistic heavy-ion collisions -- the chiral magnetic effect (CME). We report measurements by the STAR collaboration of a CME-sensitive observable in $p$+Au and $d$+Au collisions at 200 GeV, where the CME is not expected, using charge-dependent pair correlations relative to a third particle. We observe strong charge-dependent correlations similar to those measured in heavy-ion collisions. This bears important implications for the interpretation of the heavy-ion data.

We report the first measurement of the inclusive jet and the dijet longitudinal double-spin asymmetries, $A_{LL}$, at midrapidity in polarized $pp$ collisions at a center-of-mass energy $\sqrt{s} = 510$ GeV. The inclusive jet $A_{LL}$ measurement is sensitive to the gluon helicity distribution down to a gluon momentum fraction of $x\approx 0.015$, while the dijet measurements, separated into four jet-pair topologies, provide constraints on the $x$ dependence of the gluon polarization. Both results are consistent with previous measurements made at $\sqrt{s}= 200$ GeV in the overlapping kinematic region, $x > 0.05$, and show good agreement with predictions from recent next-to-leading order global analyses.

We present the first measurements of transverse momentum spectra of $\pi^{\pm}$, $K^{\pm}$, $p(\bar{p})$ at midrapidity ($|y| < 0.1$) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The centrality dependence of particle yields, average transverse momenta, particle ratios and kinetic freeze-out parameters are discussed. The results are compared with the published results from Au+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV in STAR. The results are also compared to those from A Multi Phase Transport (AMPT) model.

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in \pp and \AuAu collisions at \sqrtsn $= 200$ GeV. The substructure observables include SoftDrop groomed momentum fraction (\zg), groomed jet radius (\rg), and subjet momentum fraction (\zsj) and opening angle (\tsj). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for \pp collisions and are compared to leading order Monte Carlo models. The subjet \tsj~distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in \AuAu collisions. We compare data from \AuAu collisions to those from \pp which are embedded in minimum-bias \AuAu events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than \zg~and \rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuum-like. We also report measurements of the differential di-jet momentum imbalance ($A_{\rm{J}}$) for jets of varying \tsj. We find no qualitative differences in energy loss signatures for varying angular scales in the range $0.1 < $\tsj $ < 0.3$, leading to the possible interpretation that energy loss in this population of high momentum di-jet pairs, is due to soft medium-induced gluon radiation from a single color-charge as it traverses the medium.