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.

New measurements of directed flow for charged hadrons, characterized by the Fourier coefficient \vone, are presented for transverse momenta $\mathrm{p_T}$, and centrality intervals in Au+Au collisions recorded by the STAR experiment for the center-of-mass energy range $\mathrm{\sqrt{s_{_{NN}}}} = 7.7 - 200$ GeV. The measurements underscore the importance of momentum conservation and the characteristic dependencies on $\mathrm{\sqrt{s_{_{NN}}}}$, centrality and $\mathrm{p_T}$ are consistent with the expectations of geometric fluctuations generated in the initial stages of the collision, acting in concert with a hydrodynamic-like expansion. The centrality and $\mathrm{p_T}$ dependencies of $\mathrm{v^{even}_{1}}$, as well as an observed similarity between its excitation function and that for $\mathrm{v_3}$, could serve as constraints for initial-state models. The $\mathrm{v^{even}_{1}}$ excitation function could also provide an important supplement to the flow measurements employed for precision extraction of the temperature dependence of the specific shear viscosity.

The extreme temperatures and energy densities generated by ultra-relativistic collisions between heavy nuclei produce a state of matter with surprising fluid properties. Non-central collisions have angular momentum on the order of 1000$\hbar$, and the resulting fluid may have a strong vortical structure that must be understood to properly describe the fluid. It is also of particular interest because the restoration of fundamental symmetries of quantum chromodynamics is expected to produce novel physical effects in the presence of strong vorticity. However, no experimental indications of fluid vorticity in heavy ion collisions have so far been found. Here we present the first measurement of an alignment between the angular momentum of a non-central collision and the spin of emitted particles, revealing that the fluid produced in heavy ion collisions is by far the most vortical system ever observed. We find that $\Lambda$ and $\overline{\Lambda}$ hyperons show a positive polarization of the order of a few percent, consistent with some hydrodynamic predictions. A previous measurement that reported a null result at higher collision energies is seen to be consistent with the trend of our new observations, though with larger statistical uncertainties. These data provide the first experimental access to the vortical structure of the "perfect fluid" created in a heavy ion collision. They should prove valuable in the development of hydrodynamic models that quantitatively connect observations to the theory of the Strong Force. Our results extend the recent discovery of hydrodynamic spin alignment to the subatomic realm.

We report the energy dependence of mid-rapidity (anti-)deuteron production in Au+Au collisions at $\sqrt{s_\text{NN}} =\ $7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV, measured by the STAR experiment at RHIC. The yield of deuterons is found to be well described by the thermal model. The collision energy, centrality, and transverse momentum dependence of the coalescence parameter $B_2$ are discussed. We find that the values of $B_2$ for anti-deuterons are systematically lower than those for deuterons, indicating that the correlation volume of anti-baryons is larger than that of baryons at $\sqrt{s_\text{NN}}$ from 19.6 to 39 GeV. In addition, values of $B_2$ are found to vary with collision energy and show a broad minimum around $\sqrt{s_\text{NN}}=\ $20 to 40 GeV, which might imply a change of the equation of state of the medium in these collisions.

We present a measurement of the first-order azimuthal anisotropy, $v_1(y)$, of deuterons from Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV recorded with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The energy dependence of the $v_1(y)$ slope, $dv_{1}/dy|_{y=0}$, for deuterons, where $y$ is the rapidity, is extracted for semi-central collisions (10-40\% centrality) and compared to that of protons. While the $v_1(y)$ slopes of protons are generally negative for $\sqrt{s_{NN}} >$ 10 GeV, those for deuterons are consistent with zero, a strong enhancement of the $v_1(y)$ slope of deuterons is seen at the lowest collision energy (the largest baryon density) at $\sqrt{s_{NN}} =$ 7.7 GeV. In addition, we report the transverse momentum dependence of $v_1$ for protons and deuterons. The experimental results are compared with transport and coalescence models.

Data from the first physics run at the Relativistic Heavy-Ion Collider at Brookhaven National Laboratory, Au+Au collisions at $\sqrt{s_{NN}}=130$ GeV, have been analyzed by the STAR Collaboration using three-pion correlations with charged pions to study whether pions are emitted independently at freezeout. We have made a high-statistics measurement of the three-pion correlation function and calculated the normalized three-particle correlator to obtain a quantitative measurement of the degree of chaoticity of the pion source. It is found that the degree of chaoticity seems to increase with increasing particle multiplicity.

We present a measurement of the inclusive production of Upsilon mesons in U+U collisions at 193 GeV at mid-rapidity (|y| < 1). Previous studies in central Au+Au collisions at 200 GeV show a suppression of Upsilon(1S+2S+3S) production relative to expectations from the Upsilon yield in p+p collisions scaled by the number of binary nucleon-nucleon collisions (Ncoll), with an indication that the Upsilon(1S) state is also suppressed. The present measurement extends the number of participant nucleons in the collision (Npart) by 20% compared to Au+Au collisions, and allows us to study a system with higher energy density. We observe a suppression in both the Upsilon(1S+2S+3S) and Upsilon(1S) yields in central U+U data, which consolidates and extends the previously observed suppression trend in Au+Au collisions.

Balance functions have been measured in terms of relative pseudorapidity ($\Delta \eta$) for charged particle pairs at the Relativistic Heavy-Ion Collider (RHIC) from Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV to 200 GeV using the STAR detector. These results are compared with balance functions measured at the Large Hadron Collider (LHC) from Pb+Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV by the ALICE Collaboration. The width of the balance function decreases as the collisions become more central and as the beam energy is increased. In contrast, the widths of the balance functions calculated using shuffled events show little dependence on centrality or beam energy and are larger than the observed widths. Balance function widths calculated using events generated by UrQMD are wider than the measured widths in central collisions and show little centrality dependence. The measured widths of the balance functions in central collisions are consistent with the delayed hadronization of a deconfined quark gluon plasma (QGP). The narrowing of the balance function in central collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV implies that a QGP is still being created at this relatively low energy.

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.

We present the results of charged particle fluctuations measurements in Au + Au collisions at $\sqrt{s_{NN}}=130$ GeV using the STAR detector. Dynamical fluctuations measurements are presented for inclusive charged particle multiplicities as well as for identified charged pions, kaons, and protons. The net charge dynamical fluctuations are found to be large and negative providing clear evidence that positive and negative charged particle production is correlated within the pseudorapidity range investigated. Correlations are smaller than expected based on model-dependent predictions for a resonance gas or a quark gluon gas which undergoes fast hadronization and freeze-out. Qualitative agreement is found with comparable scaled p+p measurements and a HIJING model calculation based on independent particle collisions, although a small deviation from the 1/N scaling dependence expected from this model is observed.