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.

The system created in non-central relativistic nucleus-nucleus collisions possesses large orbital angular momentum. Due to spin-orbit coupling, particles produced in such a system could become globally polarized along the direction of the system angular momentum. We present the results of Lambda and anti-Lambda hyperon global polarization measurements in Au+Au collisions at sqrt{s_NN}=62.4 GeV and 200 GeV performed with the STAR detector at RHIC. The observed global polarization of Lambda and anti-Lambda hyperons in the STAR acceptance is consistent with zero within the precision of the measurements. The obtained upper limit, |P_{Lambda,anti-Lambda}| <= 0.02, is compared to the theoretical values discussed recently in the literature.

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 a systematic measurement of cumulants, $C_{n}$, for net-proton, proton and antiproton multiplicity distributions, and correlation functions, $\kappa_n$, for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{\mathrm {NN}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4 and 200 GeV. The $C_{n}$ and $\kappa_n$ are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, $y$, and transverse momentum, $p_{T}$. The data were taken during the first phase of the Beam Energy Scan (BES) program (2010 -- 2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity ($|y| <$ 0.5) and transverse momentum 0.4 $<$$p_{\rm T}$$<$ 2.0 GeV/$c$, using the STAR detector at RHIC. We observe a non-monotonic energy dependence ($\sqrt{s_{\mathrm {NN}}}$ = 7.7 -- 62.4 GeV) of the net-proton $C_{4}$/$C_{2}$ with the significance of 3.1$\sigma$ for the 0-5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{\mathrm {NN}}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, $\kappa_2$, of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, $\kappa_4$, of protons plays a role in determining the energy dependence of proton $C_4/C_1$ below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.

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.

According to first-principle lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region μB ≤ T c. In this range the ratio, C6=C2, of net-baryon distributions are predicted to be negative. In this Letter, we report the first measurement of the midrapidity net-proton C6=C2 from 27, 54.4, and 200 GeV Au þ Au collisions at the Relativistic Heavy Ion Collider (RHIC). The dependence on collision centrality and kinematic acceptance in (p T , y) are analyzed. While for 27 and 54.4 GeV collisions the C6=C2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C6=C2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict mostly positive values except for the most central collisions within uncertainties. These observations seem to favor a smooth crossover in the high-energy nuclear collisions at top RHIC energy.

Identified charged pion, kaon, and proton spectra are used to explore the system size dependence of bulk freeze-out properties in Cu+Cu collisions at $\sqrt{s_{NN}}$=200 and 62.4 GeV. The data are studied with hydrodynamically-motivated Blast-wave and statistical model frameworks in order to characterize the freeze-out properties of the system. The dependence of freeze-out parameters on beam energy and collision centrality is discussed. Using the existing results from Au+Au and $pp$ collisions, the dependence of freeze-out parameters on the system size is also explored. This multi-dimensional systematic study furthers our understanding of the QCD phase diagram revealing the importance of the initial geometrical overlap of the colliding ions. The analysis of Cu+Cu collisions, which expands the system size dependence studies from Au+Au data with detailed measurements in the smaller system, shows that the bulk freeze-out properties of charged particles studied here scale with the total charged particle multiplicity at mid-rapidity, suggesting the relevance of initial state effects.

We present a systematic analysis of two-pion interferometry in Au+Au collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 GeV and Cu+Cu collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 and 200 GeV using the STAR detector at RHIC. The multiplicity and transverse momentum dependences of the extracted femtoscopic radii are studied. The scaling of the apparent freeze-out volume with charged particle multiplicity is studied for the RHIC energy domain. The multiplicity scaling of the measured radii is found to be independent of colliding system and collision energy.

We present results for the measurement of $\phi$ meson production via its charged kaon decay channel $\phi \to K^+K^-$ in Au+Au collisions at $\sqrt{s_{_{NN}}}=62.4$, 130, 200 GeV, and in p+p and d+Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV from the STAR experiment at RHIC. The mid-rapidity ($|y|<0.5$) $\phi$ meson spectra in central Au+Au collisions are found to be well described by a single exponential distribution. On the other hand, the spectra from p+p, d+Au and peripheral Au+Au collisions show power-law tails at intermediate and high transverse momenta ($p_{T}$) and are described better by Levy distributions. The constant $\phi/K^-$ yield ratio vs. beam species, collision centrality and colliding energy is in contradiction with expectations from models having kaon coalescence as the dominant mechanism for $\phi$ production at RHIC. The $\Omega/\phi$ yield ratio as a function of $p_{T}$ is consistent with a model based on the recombination of thermal $s$ quarks up to $p_{T}\sim 4$ GeV/c, but disagrees at higher transverse momenta. The measured nuclear modification factor, $R_{dAu}$, for the $\phi$ meson increases above unity at intermediate $p_{T}$, similar to that for pions and protons, while $R_{AA}$ is suppressed due to jet quenching in central Au+Au collisions. Number of constituent quark scaling of both $R_{cp}$ and $v_{2}$ for the $\phi$ meson with respect to other hadrons in Au+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV at intermediate $p_{T}$ is observed. These observations support quark coalescence as being the dominant mechanism of hadronization in the intermediate $p_{T}$ region at RHIC.

We present the directed flow ($v_1$) measured in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV in the mid-pseudorapidity region $|\eta|<1.3$ and in the forward pseudorapidity region $2.5 < |\eta| < 4.0$. The results are obtained using the three-particle cumulant method, the event plane method with mixed harmonics, and for the first time at the Relativistic Heavy Ion Collider (RHIC), the standard method with the event plane reconstructed from spectator neutrons. Results from all three methods are in good agreement. Over the pseudorapidity range studied, charged particle directed flow is in the direction opposite to that of fragmentation neutrons.