We present high precision measurements of elliptic flow near midrapidity ($|y|<1.0$) for multi-strange hadrons and $\phi$ meson as a function of centrality and transverse momentum in Au+Au collisions at center of mass energy $\sqrt{s_{NN}}=$ 200 GeV. We observe that the transverse momentum dependence of $\phi$ and $\Omega$ $v_{2}$ is similar to that of $\pi$ and $p$, respectively, which may indicate that the heavier strange quark flows as strongly as the lighter up and down quarks. This observation constitutes a clear piece of evidence for the development of partonic collectivity in heavy-ion collisions at the top RHIC energy. Number of constituent quark scaling is found to hold within statistical uncertainty for both 0-30$\%$ and 30-80$\%$ collision centrality. There is an indication of the breakdown of previously observed mass ordering between $\phi$ and proton $v_{2}$ at low transverse momentum in the 0-30$\%$ centrality range, possibly indicating late hadronic interactions affecting the proton $v_{2}$.

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.

Collisions between prolate uranium nuclei are used to study how particle production and azimuthal anisotropies depend on initial geometry in heavy-ion collisions. We report the two- and four-particle cumulants, $v_2\{2\}$ and $v_2\{4\}$, for charged hadrons from U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV and Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. Nearly fully overlapping collisions are selected based on the amount of energy deposited by spectators in the STAR Zero Degree Calorimeters (ZDCs). Within this sample, the observed dependence of $v_2\{2\}$ on multiplicity demonstrates that ZDC information combined with multiplicity can preferentially select different overlap configurations in U+U collisions. An initial-state model with gluon saturation describes the slope of $v_2\{2\}$ as a function of multiplicity in central collisions better than one based on Glauber with a two-component multiplicity model.

The STAR collaboration presents for the first time two-dimensional di-hadron correlations with identified leading hadrons in 200 GeV central Au+Au and minimum-bias d+Au collisions to explore hadronization mechanisms in the quark gluon plasma. The enhancement of the jet-like yield for leading pions in Au+Au data with respect to the d+Au reference and the absence of such an enhancement for leading non-pions (protons and kaons) are discussed within the context of a quark recombination scenario. The correlated yield at large angles, specifically in the \emph{ridge region}, is found to be significantly higher for leading non-pions than pions. The consistencies of the constituent quark scaling, azimuthal harmonic model and a mini-jet modification model description of the data are tested, providing further constraints on hadronization.

The acceptance-corrected dielectron excess mass spectra, where the known hadronic sources have been subtracted from the inclusive dielectron mass spectra, are reported for the first time at mid-rapidity $|y_{ee}|<1$ in minimum-bias Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 200 GeV. The excess mass spectra are consistently described by a model calculation with a broadened $\rho$ spectral function for $M_{ee}<1.1$ GeV/$c^{2}$. The integrated dielectron excess yield at $\sqrt{s_{NN}}$ = 19.6 GeV for $0.4<M_{ee}<0.75$ GeV/$c^2$, normalized to the charged particle multiplicity at mid-rapidity, has a value similar to that in In+In collisions at $\sqrt{s_{NN}}$ = 17.3 GeV. For $\sqrt{s_{NN}}$ = 200 GeV, the normalized excess yield in central collisions is higher than that at $\sqrt{s_{NN}}$ = 17.3 GeV and increases from peripheral to central collisions. These measurements indicate that the lifetime of the hot, dense medium created in central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV is longer than those in peripheral collisions and at lower energies.

A search for the quantum chromodynamics (QCD) critical point was performed by the STAR experiment at the Relativistic Heavy Ion Collider, using dynamical fluctuations of unlike particle pairs. Heavy-ion collisions were studied over a large range of collision energies with homogeneous acceptance and excellent particle identification, covering a significant range in the QCD phase diagram where a critical point may be located. Dynamical $K\pi$, $p\pi$, and $Kp$ fluctuations as measured by the STAR experiment in central 0-5\% Au+Au collisions from center-of-mass collision energies $\rm \sqrt{s_{NN}}$ = 7.7 to 200 GeV are presented. The observable $\rm \nu_{dyn}$ was used to quantify the magnitude of the dynamical fluctuations in event-by-event measurements of the $K\pi$, $p\pi$, and $Kp$ pairs. The energy dependences of these fluctuations from central 0-5\% Au+Au collisions all demonstrate a smooth evolution with collision energy.

We present measurements of $\Omega$ and $\phi$ production at mid-rapidity from Au+Au collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27 and 39 GeV by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). Motivated by the coalescence formation mechanism for these strange hadrons, we study the ratios of $N(\Omega^{-}+\Omega^{+})/(2N(\phi))$. These ratios as a function of transverse momentum ($p_T$) fall on a consistent trend at high collision energies, but start to show deviations in peripheral collisions at $\sqrt{s_{NN}}$ = 19.6, 27 and 39 GeV, and in central collisions at 11.5 GeV in the intermediate $p_T$ region of 2.4-3.6 GeV/c. We further evaluate empirically the strange quark $p_T$ distributions at hadronization by studying the $\Omega/\phi$ ratios scaled by the number of constituent quarks. The NCQ-scaled $\Omega/\phi$ ratios show a suppression of strange quark production in central collisions at 11.5 GeV compared to $\sqrt{s_{NN}} >= 19.6$ GeV. The shapes of the presumably thermal strange quark distributions in 0-60% most central collisions at 7.7 GeV show significant deviations from those in 0-10% most central collisions at higher energies. These features suggest that there is likely a change of the underlying strange quark dynamics in the transition from quark-matter to hadronic matter at collision energies below 19.6 GeV.

We present measurements of $\pi^-$ and $\pi^+$ elliptic flow, $v_2$, at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}} =$ 200, 62.4, 39, 27, 19.6, 11.5 and 7.7 GeV, as a function of event-by-event charge asymmetry, $A_{ch}$, based on data from the STAR experiment at RHIC. We find that $\pi^-$ ($\pi^+$) elliptic flow linearly increases (decreases) with charge asymmetry for most centrality bins at $\sqrt{s_{_{\rm NN}}} = \text{27 GeV}$ and higher. At $\sqrt{s_{_{\rm NN}}} = \text{200 GeV}$, the slope of the difference of $v_2$ between $\pi^-$ and $\pi^+$ as a function of $A_{ch}$ exhibits a centrality dependence, which is qualitatively similar to calculations that incorporate a chiral magnetic wave effect. Similar centrality dependence is also observed at lower energies.

We report on measurements of dielectron ($e^+e^-$) production in Au$+$Au collisions at a center-of-mass energy of 200 GeV per nucleon-nucleon pair using the STAR detector at RHIC. Systematic measurements of the dielectron yield as a function of transverse momentum ($p_{\rm T}$) and collision centrality show an enhancement compared to a cocktail simulation of hadronic sources in the low invariant-mass region ($M_{ee}<$ 1 GeV/$c^2$). This enhancement cannot be reproduced by the $\rho$-meson vacuum spectral function. In minimum-bias collisions, in the invariant-mass range of 0.30 $-$ 0.76 GeV/$c^2$, integrated over the full $p_{\rm T}$ acceptance, the enhancement factor is 1.76 $\pm$ 0.06 (stat.) $\pm$ 0.26 (sys.) $\pm$ 0.29 (cocktail). The enhancement factor exhibits weak centrality and $p_{\rm T}$ dependence in STAR's accessible kinematic regions, while the excess yield in this invariant-mass region as a function of the number of participating nucleons follows a power-law shape with a power of 1.44 $\pm$ 0.10. Models that assume an in-medium broadening of the $\rho$ meson spectral function consistently describe the observed excess in these measurements. Additionally, we report on measurements of $\omega$ and $\phi$-meson production through their $e^+e^-$ decay channel. These measurements show good agreement with Tsallis Blast-Wave model predictions as well as, in the case of the $\phi$-meson, results through its $K^+K^-$ decay channel. In the intermediate invariant-mass region (1.1$<M_{ee}<$ 3 GeV/$c^2$), we investigate the spectral shapes from different collision centralities. Physics implications for possible in-medium modification of charmed hadron production and other physics sources are discussed.

We report the observation of transverse polarization-dependent azimuthal correlations in charged pion pair production with the STAR experiment in $p^\uparrow+p$ collisions at RHIC. These correlations directly probe quark transversity distributions. We measure signals in excess of five standard deviations at high transverse momenta, at high pseudorapidities eta>0.5, and for pair masses around the mass of the rho-meson. This is the first direct transversity measurement in p+p collisions. Comparing the results to data from lepton-nucleon scattering will test the universality of these spin-dependent quantities.

Dihadron correlations are analyzed in $\sqrt{s_{_{\rm NN}}} = 200$ GeV $d$+Au collisions classified by forward charged particle multiplicity and zero-degree neutral energy in the Au-beam direction. It is found that the jetlike correlated yield increases with the event multiplicity. After taking into account this dependence, the non-jet contribution on the away side is minimal, leaving little room for a back-to-back ridge in these collisions.

Dihadron angular correlations in $d$+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV are reported as a function of the measured zero-degree calorimeter neutral energy and the forward charged hadron multiplicity in the Au-beam direction. A finite correlated yield is observed at large relative pseudorapidity ($\Delta\eta$) on the near side (i.e. relative azimuth $\Delta\phi\sim0$). This correlated yield as a function of $\Delta\eta$ appears to scale with the dominant, primarily jet-related, away-side ($\Delta\phi\sim\pi$) yield. The Fourier coefficients of the $\Delta\phi$ correlation, $V_{n}=\langle\cos n\Delta\phi\rangle$, have a strong $\Delta\eta$ dependence. In addition, it is found that $V_{1}$ is approximately inversely proportional to the mid-rapidity event multiplicity, while $V_{2}$ is independent of it with similar magnitude in the forward ($d$-going) and backward (Au-going) directions.

We report measurements of Upsilon meson production in p+p, d+Au, and Au+Au collisions using the STAR detector at RHIC. We compare the Upsilon yield to the measured cross section in p+p collisions in order to quantify any modifications of the yield in cold nuclear matter using d+Au data and in hot nuclear matter using Au+Au data separated into three centrality classes. Our p+p measurement is based on three times the statistics of our previous result. We obtain a nuclear modification factor for Upsilon(1S+2S+3S) in the rapidity range |y|<1 in d+Au collisions of R_dAu = 0.79 +/- 0.24 (stat.) +/- 0.03 (sys.) +/- 0.10 (pp sys.). A comparison with models including shadowing and initial state parton energy loss indicates the presence of additional cold-nuclear matter suppression. Similarly, in the top 10% most-central Au+Au collisions, we measure a nuclear modification factor of R_AA=0.49 +/- 0.1 (stat.) +/- 0.02 (sys.) +/- 0.06 (pp sys.), which is a larger suppression factor than that seen in cold nuclear matter. Our results are consistent with complete suppression of excited-state Upsilon mesons in Au+Au collisions. The additional suppression in Au+Au is consistent with the level expected in model calculations that include the presence of a hot, deconfined Quark-Gluon Plasma. However, understanding the suppression seen in d+Au is still needed before any definitive statements about the nature of the suppression in Au+Au can be made.

A simple C++ class structure for construction of a Monte Carlo event generators which can produce unweighted events within relativistic phase space is presented. The generator is self-adapting to the provided matrix element and acceptance cuts. The program is designed specially for exclusive processes and includes, as an example of such an application, implementation of the model for exclusive production of meson pairs $pp \rightarrow p M^+M^- p $ in high energy proton-proton collisions.

We report on a polarization measurement of inclusive $J/\psi$ mesons in the di-electron decay channel at mid-rapidity at 2 $<p_{T}<$ 6 GeV/$c$ in $p+p$ collisions at $\sqrt{s}$ = 200 GeV. Data were taken with the STAR detector at RHIC. The $J/\psi$ polarization measurement should help to distinguish between different models of the $J/\psi$ production mechanism since they predict different $p_{T}$ dependences of the $J/\psi$ polarization. In this analysis, $J/\psi$ polarization is studied in the helicity frame. The polarization parameter $\lambda_{\theta}$ measured at RHIC becomes smaller towards high $p_{T}$, indicating more longitudinal $J/\psi$ polarization as $p_{T}$ increases. The result is compared with predictions of presently available models.

The $\jpsi$ $\pt$ spectrum and nuclear modification factor ($\raa$) are reported for $\pt < 5 \ \gevc$ and $|y|<1$ from 0\% to 60\% central Au+Au and Cu+Cu collisions at $\snn = 200 \ \gev$ at STAR. A significant suppression of $\pt$-integrated $\jpsi$ production is observed in central Au+Au events. The Cu+Cu data are consistent with no suppression, although the precision is limited by the available statistics. $\raa$ in Au+Au collisions exhibits a strong suppression at low transverse momentum and gradually increases with $\pt$. The data are compared to high-$\pt$ STAR results and previously published BNL Relativistic Heavy Ion Collider results. Comparing with model calculations, it is found that the invariant yields at low $\pt$ are significantly above hydrodynamic flow predictions but are consistent with models that include color screening and regeneration.

We report the first measurement of charmed-hadron ($D^0$) production via the hadronic decay channel ($D^0\rightarrow K^- + \pi^+$) in Au+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 200\,GeV with the STAR experiment. The charm production cross-section per nucleon-nucleon collision at mid-rapidity scales with the number of binary collisions, $N_{bin}$, from $p$+$p$ to central Au+Au collisions. The $D^0$ meson yields in central Au+Au collisions are strongly suppressed compared to those in $p$+$p$ scaled by $N_{bin}$, for transverse momenta $p_{T}>3$ GeV/$c$, demonstrating significant energy loss of charm quarks in the hot and dense medium. An enhancement at intermediate $p_{T}$ is also observed. Model calculations including strong charm-medium interactions and coalescence hadronization describe our measurements.

A data-driven method was applied to measurements of Au+Au collisions at $\sqrt{s_{_{\rm NN}}} =$ 200 GeV made with the STAR detector at RHIC to isolate pseudorapidity distance $\Delta\eta$-dependent and $\Delta\eta$-independent correlations by using two- and four-particle azimuthal cumulant measurements. We identified a component of the correlation that is $\Delta\eta$-independent, which is likely dominated by anisotropic flow and flow fluctuations. It was also found to be independent of $\eta$ within the measured range of pseudorapidity $|\eta|<1$. The relative flow fluctuation was found to be $34\% \pm 2\% (stat.) \pm 3\% (sys.)$ for particles of transverse momentum $p_{T}$ less than $2$ GeV/$c$. The $\Delta\eta$-dependent part may be attributed to nonflow correlations, and is found to be $5\% \pm 2\% (sys.)$ relative to the flow of the measured second harmonic cumulant at $|\Delta\eta| > 0.7$.

Event-by-event fluctuations of the ratio of inclusive charged to photon multiplicities at forward rapidity in Au+Au collision at $\sqrt{s_{NN}}$=200 GeV have been studied. Dominant contribution to such fluctuations is expected to come from correlated production of charged and neutral pions. We search for evidences of dynamical fluctuations of different physical origins. Observables constructed out of moments of multiplicities are used as measures of fluctuations. Mixed events and model calculations are used as baselines. Results are compared to the dynamical net-charge fluctuations measured in the same acceptance. A non-zero statistically significant signal of dynamical fluctuations is observed in excess to the model prediction when charged particles and photons are measured in the same acceptance. We find that, unlike dynamical net-charge fluctuation, charge-neutral fluctuation is not dominated by correlation due to particle decay. Results are compared to the expectations based on the generic production mechanism of pions due to isospin symmetry, for which no significant (<1%) deviation is observed.

We report measurements of single- and double- spin asymmetries for $W^{\pm}$ and $Z/\gamma^*$ boson production in longitudinally polarized $p+p$ collisions at $\sqrt{s} = 510$ GeV by the STAR experiment at RHIC. The asymmetries for $W^{\pm}$ were measured as a function of the decay lepton pseudorapidity, which provides a theoretically clean probe of the proton's polarized quark distributions at the scale of the $W$ mass. The results are compared to theoretical predictions, constrained by recent polarized deep inelastic scattering measurements, and show a preference for a sizable, positive up antiquark polarization in the range $0.05<x<0.2$.