Suppression of the $J/\psi$ nuclear-modification factor has been seen as a trademark signature of final-state effects in large collision systems for decades. In small systems, the nuclear modification was attributed to cold-nuclear-matter effects until the observation of strong differential suppression of the $\psi(2S)$ state in $p/d$ $+$ $A$ collisions suggested the presence of final-state effects. Results of $J/\psi$ and $\psi(2S)$ measurements in the dimuon decay channel are presented here for $p$ $+$ $p$, $p$ $+$Al, and $p$ $+$Au collision systems at $\sqrt{s_{_{NN}}}=200$ GeV. The results are predominantly shown in the form of the nuclear-modification factor, $R_{pA}$, the ratio of the $\psi(2S)$ invariant yield per nucleon-nucleon collision in collisions of proton on target nucleus to that in $p$ $+$ $p$ collisions. Measurements of the $J/\psi$ and $\psi(2S)$ nuclear-modification factor are compared with shadowing and transport-model predictions, as well as to complementary measurements at Large-Hadron-Collider energies.

Charmonium is a valuable probe in heavy-ion collisions to study the properties of the quark gluon plasma, and is also an interesting probe in small collision systems to study cold nuclear matter effects, which are also present in large collision systems. With the recent observations of collective behavior of produced particles in small system collisions, measurements of the modification of charmonium in small systems have become increasingly relevant. We present the results of J/ψ measurements at forward and backward rapidity in various small collision systems, p+p, p+Al, p+Au and 3He+Au, at √sNN =200 GeV. The results are presented in the form of the observable RAB, the nuclear modification factor, a measure of the ratio of the J/ψ invariant yield compared to the scaled yield in p+p collisions. We examine the rapidity, transverse momentum, and collision centrality dependence of nuclear effects on J/ψ production with different projectile sizes p and 3He, and different target sizes Al and Au. The modification is found to be strongly dependent on the target size, but to be very similar for p+Au and 3He+Au. However, for 0%–20% central collisions at backward rapidity, the modification for 3He+Au is found to be smaller than that for p+Au, with a mean fit to the ratio of 0.89±0.03(stat)±0.08(syst), possibly indicating final state effects due to the larger projectile size.

We present an analysis of di-hadron correlations using recently developed methods for background subtraction which allow for higher precision measurements with fewer assumptions about the background. These studies indicate that low momentum jets interacting with the medium do not equilibrate with the medium, but rather that interactions with the medium lead to more subtle increases in their widths and fragmentation functions, consistent with observations from studies of higher momentum fully reconstructed jets. The away-side shape is not consistent with a Mach cone.

We present measurements of the near-side of triggered di-hadron correlations using neutral strange baryons ($\Lambda$, $\bar{\Lambda}$) and mesons ($K^0_S$) at intermediate transverse momentum (3 $<$ $p_T$ $<$ 6 GeV/$c$) to look for possible flavor and baryon/meson dependence. This study is performed in $d$+Au, Cu+Cu and Au+Au collisions at $\sqrt{s_{{NN}}}$ = 200 GeV measured by the STAR experiment at RHIC. The near-side di-hadron correlation contains two structures, a peak which is narrow in azimuth and pseudorapidity consistent with correlations due to jet fragmentation, and a correlation in azimuth which is broad in pseudorapidity. The particle composition of the jet-like correlation is determined using identified associated particles. The dependence of the conditional yield of the jet-like correlation on the trigger particle momentum, associated particle momentum, and centrality for correlations with unidentified trigger particles are presented. The neutral strange particle composition in jet-like correlations with unidentified charged particle triggers is not well described by PYTHIA. However, the yield of unidentified particles in jet-like correlations with neutral strange particle triggers is described reasonably well by the same model.

We report on the measurement of $\rm{J}/\psi$ production in the dielectron channel at mid-rapidity (|y|<1) in p+p and d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the STAR experiment at the Relativistic Heavy Ion Collider. The transverse momentum $p_{T}$ spectra in p+p for $p_{T}$ < 4 GeV/c and d+Au collisions for $p_{T}$ < 3 GeV/c are presented. These measurements extend the STAR coverage for $\rm{J}/\psi$ production in p+p collisions to low $p_{T}$. The $<p_{T}^{2}>$ from the measured $\rm{J}/\psi$ invariant cross section in p+p and d+Au collisions are evaluated and compared to similar measurements at other collision energies. The nuclear modification factor for $\rm{J}/\psi$ is extracted as a function of $p_{T}$ and collision centrality in d+Au and compared to model calculations using the modified nuclear Parton Distribution Function and a final-state $\rm{J}/\psi$ nuclear absorption cross section.

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 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}$.

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.

We present measurements from the PHENIX experiment of large parity-violating single spin asymmetries of high transverse momentum electrons and positrons from $W^\pm/Z$ decays, produced in longitudinally polarized $p$$+$$p$ collisions at center of mass energies of $\sqrt{s}$=500 and 510~GeV. These asymmetries allow direct access to the anti-quark polarized parton distribution functions due to the parity-violating nature of the $W$-boson coupling to quarks and anti-quarks. The results presented are based on data collected in 2011, 2012, and 2013 with an integrated luminosity of 240 pb$^{-1}$, which exceeds previous PHENIX published results by a factor of more than 27. These high $Q^2$ data provide an important addition to our understanding of anti-quark parton helicity distribution functions.

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.