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.

Parity-odd domains, corresponding to non-trivial topological solutions of the QCD vacuum, might be created during relativistic heavy-ion collisions. These domains are predicted to lead to charge separation of quarks along the orbital momentum of the system created in non-central collisions. To study this effect, we investigate a three particle mixed harmonics azimuthal correlator which is a \P-even observable, but directly sensitive to the charge separation effect. We report measurements of this observable using the STAR detector in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$=200 and 62~GeV. The results are presented as a function of collision centrality, particle separation in rapidity, and particle transverse momentum. A signal consistent with several of the theoretical expectations is detected in all four data sets. We compare our results to the predictions of existing event generators, and discuss in detail possible contributions from other effects that are not related to parity violation.

We study the beam-energy and system-size dependence of \phi meson production (using the hadronic decay mode \phi -- K+K-) by comparing the new results from Cu+Cu collisions and previously reported Au+Au collisions at \sqrt{s_NN} = 62.4 and 200 GeV measured in the STAR experiment at RHIC. Data presented are from mid-rapidity (|y|<0.5) for 0.4 < pT < 5 GeV/c. At a given beam energy, the transverse momentum distributions for \phi mesons are observed to be similar in yield and shape for Cu+Cu and Au+Au colliding systems with similar average numbers of participating nucleons. The \phi meson yields in nucleus-nucleus collisions, normalised by the average number of participating nucleons, are found to be enhanced relative to those from p+p collisions with a different trend compared to strange baryons. The enhancement for \phi mesons is observed to be higher at \sqrt{s_NN} = 200 GeV compared to 62.4 GeV. These observations for the produced \phi(s\bar{s}) mesons clearly suggest that, at these collision energies, the source of enhancement of strange hadrons is related to the formation of a dense partonic medium in high energy nucleus-nucleus collisions and cannot be alone due to canonical suppression of their production in smaller systems.

Identified charged particle spectra of $\pi^{\pm}$, $K^{\pm}$, $p$ and $\pbar$ at mid-rapidity ($|y|<0.1$) measured by the $\dedx$ method in the STAR-TPC are reported for $pp$ and d+Au collisions at $\snn = 200$ GeV and for Au+Au collisions at 62.4 GeV, 130 GeV, and 200 GeV. ... [Shortened for arXiv list. Full abstract in manuscript.]

We present results on strange and multi-strange particle production in Au+Au collisions at $\sqrt{s_{NN}}=62.4$ GeV as measured with the STAR detector at RHIC. Mid-rapidity transverse momentum spectra and integrated yields of $K^{0}_{S}$, $\Lambda$, $\Xi$, $\Omega$ and their anti-particles are presented for different centrality classes. The particle yields and ratios follow a smooth energy dependence. Chemical freeze-out parameters, temperature, baryon chemical potential and strangeness saturation factor obtained from the particle yields are presented. Intermediate transverse momentum ($p_T$) phenomena are discussed based on the ratio of the measured baryon-to-meson spectra and nuclear modification factor. The centrality dependence of various measurements presented show a similar behavior as seen in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

A systematic study is presented for centrality, transverse momentum ($p_T$) and pseudorapidity ($\eta$) dependence of the inclusive charged hadron elliptic flow ($v_2$) at midrapidity($|\eta| < 1.0$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27 and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and 4-particle cumulants ($v_2{4}$), are presented in order to investigate non-flow correlations and $v_2$ fluctuations. We observe that the difference between $v_2{2}$ and $v_2{4}$ is smaller at the lower collision energies. Values of $v_2$, scaled by the initial coordinate space eccentricity, $v_{2}/\varepsilon$, as a function of $p_T$ are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider ($\sqrt{s_{NN}}$ = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV). The $v_2(p_T)$ values for fixed $p_T$ rise with increasing collision energy within the $p_T$ range studied ($< 2 {\rm GeV}/c$). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of $v_{2}(p_{T})$. We also compare the $v_2$ results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at Beam Energy Scan (BES) energies are discussed.

Local parity-odd domains are theorized to form inside a Quark-Gluon-Plasma (QGP) which has been produced in high-energy heavy-ion collisions. The local parity-odd domains manifest themselves as charge separation along the magnetic field axis via the chiral magnetic effect (CME). The experimental observation of charge separation has previously been reported for heavy-ion collisions at the top RHIC energies. In this paper, we present the results of the beam-energy dependence of the charge correlations in Au+Au collisions at midrapidity for center-of-mass energies of 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV from the STAR experiment. After background subtraction, the signal gradually reduces with decreased beam energy, and tends to vanish by 7.7 GeV. The implications of these results for the CME will be discussed.

Elliptic flow (v_2) values for identified particles at midrapidity in Au + Au collisions measured by the STAR experiment in the Beam Energy Scan at the Relativistic Heavy Ion Collider at sqrt{s_{NN}}= 7.7--62.4 GeV are presented for three centrality classes. The centrality dependence and the data at sqrt{s_{NN}}= 14.5 GeV are new. Except at the lowest beam energies we observe a similar relative v_2 baryon-meson splitting for all centrality classes which is in agreement within 15% with the number-of-constituent quark scaling. The larger v_2 for most particles relative to antiparticles, already observed for minimum bias collisions, shows a clear centrality dependence, with the largest difference for the most central collisions. Also, the results are compared with A Multiphase Transport Model and fit with a Blast Wave model.

We present results from a harmonic decomposition of two-particle azimuthal correlations measured with the STAR detector in Au+Au collisions for energies ranging from $\sqrt{s_{NN}}=7.7$ GeV to 200 GeV. The third harmonic $v_3^2\{2\}=\langle \cos3(\phi_1-\phi_2)\rangle$, where $\phi_1-\phi_2$ is the angular difference in azimuth, is studied as a function of the pseudorapidity difference between particle pairs $\Delta\eta = \eta_1-\eta_2$. Non-zero {\vthree} is directly related to the previously observed large-$\Delta\eta$ narrow-$\Delta\phi$ ridge correlations and has been shown in models to be sensitive to the existence of a low viscosity Quark Gluon Plasma (QGP) phase. For sufficiently central collisions, $v_3^2\{2\}$ persist down to an energy of 7.7 GeV suggesting that QGP may be created even in these low energy collisions. In peripheral collisions at these low energies however, $v_3^2\{2\}$ is consistent with zero. When scaled by pseudorapidity density of charged particle multiplicity per participating nucleon pair, $v_3^2\{2\}$ for central collisions shows a minimum near {\snn}$=20$ GeV.

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.