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.

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.

Measurements of the elliptic flow, $v_{2}$, of identified hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{s}^{0}$, $p$, $\bar{p}$, $\phi$, $\Lambda$, $\bar{\Lambda}$, $\Xi^{-}$, $\bar{\Xi}^{+}$, $\Omega^{-}$, $\bar{\Omega}^{+}$) in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV are presented. The measurements were done at mid-rapidity using the Time Projection Chamber and the Time-of-Flight detectors of the STAR experiment during the Beam Energy Scan program at RHIC. A significant difference in the $v_{2}$ values for particles and the corresponding anti-particles was observed at all transverse momenta for the first time. The difference increases with decreasing center-of-mass energy, $\sqrt{s_{NN}}$ (or increasing baryon chemical potential, $\mu_{B}$) and is larger for the baryons as compared to the mesons. This implies that particles and anti-particles are no longer consistent with the universal number-of-constituent quark (NCQ) scaling of $v_{2}$ that was observed at $\sqrt{s_{NN}}=$ 200 GeV. However, for the group of particles NCQ scaling at $(m_{T}-m_{0})/n_{q}>$ 0.4 GeV/$c^{2}$ is not violated within $\pm$10%. The $v_{2}$ values for $\phi$ mesons at 7.7 and 11.5 GeV are approximately two standard deviations from the trend defined by the other hadrons at the highest measured $p_{T}$ values.

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.

We report on K*0 production at mid-rapidity in Au+Au and Cu+Cu collisions at \sqrt{s_{NN}} = 62.4 and 200 GeV collected by the Solenoid Tracker at RHIC (STAR) detector. The K*0 is reconstructed via the hadronic decays K*0 \to K+ pi- and \bar{K*0} \to K-pi+. Transverse momentum, pT, spectra are measured over a range of pT extending from 0.2 GeV/c to 5 GeV/c. The center of mass energy and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, <pT>, are presented. The measured N(K*0)/N(K) and N(\phi)/N(K*0) ratios favor the dominance of re-scattering of decay daughters of K*0 over the hadronic regeneration for the K*0 production. In the intermediate pT region (2.0 < pT < 4.0 GeV/c), the elliptic flow parameter, v2, and the nuclear modification factor, RCP, agree with the expectations from the quark coalescence model of particle production.

We present STAR measurements of strange hadron ($\mathrm{K}^{0}_{\mathrm S}$, $\Lambda$, $\overline{\Lambda}$, $\Xi^-$, $\overline{\Xi}^+$, $\Omega^-$, $\overline{\Omega}^+$, and $\phi$) production at mid-rapidity ($|y| < 0.5$) in Au+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 7.7 - 39 GeV from the Beam Energy Scan Program at the Relativistic Heavy Ion Collider (RHIC). Transverse momentum spectra, averaged transverse mass, and the overall integrated yields of these strange hadrons are presented versus the centrality and collision energy. Antibaryon-to-baryon ratios ($\overline{\Lambda}$/$\Lambda$, $\overline{\Xi}^+$/$\Xi^-$, $\overline{\Omega}^+$/$\Omega^-$) are presented as well, and used to test a thermal statistical model and to extract the temperature normalized strangeness and baryon chemical potentials at hadronic freeze-out ($\mu_{B}/T_{\rm ch}$ and $\mu_{S}/T_{\rm ch}$) in central collisions. Strange baryon-to-pion ratios are compared to various model predictions in central collisions for all energies. The nuclear modification factors ($R_{\textrm{CP}}$) and antibaryon-to-meson ratios as a function of transverse momentum are presented for all collision energies. The $\mathrm{K}^{0}_{\mathrm S}$$R_{\textrm{CP}}$ shows no suppression for $p_{\rm T}$ up to 3.5 $\mathrm{GeV} / c$ at energies of 7.7 and 11.5 GeV. The $\overline{\Lambda}$/$\mathrm{K}^{0}_{\mathrm S}$ ratio also shows baryon-to-meson enhancement at intermediate $p_{\rm T}$ ($\approx$2.5 $\mathrm{GeV} / c$) in central collisions at energies above 19.6 GeV. Both observations suggest that there is likely a change of the underlying strange quark dynamics at collision energies below 19.6 GeV.

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.

We report on new measurements of elliptic flow ($v_2$) of electrons from heavy-flavor hadron decays at mid-rapidity ($|y|<0.8$) in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 27 and 54.4 GeV from the STAR experiment. Heavy-flavor decay electrons ($e^{\rm HF}$) in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 54.4 GeV exhibit a non-zero $v_2$ in the transverse momentum ($p_{\rm T}$) region of $p_{\rm T}<$ 2 GeV/$c$ with the magnitude comparable to that at $\sqrt{s_{_{\rm NN}}}=200$ GeV. The measured $e^{\rm HF}$$v_2$ at 54.4 GeV is also consistent with the expectation of their parent charm hadron $v_2$ following number-of-constituent-quark scaling as other light and strange flavor hadrons at this energy. These suggest that charm quarks gain significant collectivity through the evolution of the QCD medium and may reach local thermal equilibrium in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=54.4$ GeV. The measured $e^{\rm HF}$$v_2$ in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=$ 27 GeV is consistent with zero within large uncertainties. The energy dependence of $v_2$ for different flavor particles ($\pi,\phi,D^{0}/e^{\rm HF}$) shows an indication of quark mass hierarchy in reaching thermalization in high-energy nuclear collisions.

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.