The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section, mean transverse momentum, mean transverse momentum squared of inclusive $J/\psi$ and cross-section ratio of $\psi(2S)$ to $J/\psi$ at forward rapidity in \pp collisions at \sqrts = 510 GeV via the dimuon decay channel. Comparison is made to inclusive $J/\psi$ cross sections measured at \sqrts = 200 GeV and 2.76--13 TeV. The result is also compared to leading-order nonrelativistic QCD calculations coupled to a color-glass-condensate description of the low-$x$ gluons in the proton at low transverse momentum ($p_T$) and to next-to-leading order nonrelativistic QCD calculations for the rest of the $p_T$ range. These calculations overestimate the data at low $p_T$. While consistent with the data within uncertainties above $\approx3$ GeV/$c$, the calculations are systematically below the data. The total cross section times the branching ratio is BR $d\sigma^{J/\psi}_{pp}/dy (1.2<|y|<2.2, 0<p_T<10~\mbox{GeV/$c$}) =$ 54.3 $\pm$ 0.5 (stat) $\pm$ 5.5 (syst) nb.

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

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured $\phi$ meson production and its nuclear modification in asymmetric Cu$+$Au heavy-ion collisions at $\sqrt{s_{NN}}=200$ GeV at both forward Cu-going direction ($1.2<y<2.2$) and backward Au-going direction ($-2.2<y<-1.2$), rapidities. The measurements are performed via the dimuon decay channel and reported as a function of the number of participating nucleons, rapidity, and transverse momentum. In the most central events, 0\%--20\% centrality, the $\phi$ meson yield integrated over $1<p_T<5$ GeV/$c$ prefers a smaller value, which means a larger nuclear modification, in the Cu-going direction compared to the Au-going direction. Additionally, the nuclear-modification factor in Cu$+$Au collisions averaged over all centrality is measured to be similar to the previous PHENIX result in $d$$+$Au collisions for these rapidities.

Vector mesons may be photoproduced in relativistic heavy-ion collisions when a virtual photon emitted by one nucleus scatters from the other nucleus, emerging as a vector meson. The STAR Collaboration has previously presented measurements of coherent $\rho^0$ photoproduction at center of mass energies of 130 GeV and 200 GeV in AuAu collisions. Here, we present a measurement of the cross section at 62.4 GeV; we find that the cross section for coherent $\rho^0$ photoproduction with nuclear breakup is $10.5\pm1.5\pm 1.6$ mb at 62.4 GeV. The cross-section ratio between 200 GeV and 62.4 GeV is $2.8\pm0.6$, less than is predicted by most theoretical models. It is, however, proportionally much larger than the previously observed $15\pm 55$% increase between 130 GeV and 200 GeV.

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 report the first measurement of the full angular distribution for inclusive $J/\psi\rightarrow\mu^{+}\mu^{-}$ decays in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. The measurements are made for $J/\psi$ transverse momentum $2<p_{T}<10$ GeV/$c$ and rapidity $1.2<y<2.2$ in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient $\lambda_{\theta}$ is strongly negative at low $p_{T}$ and becomes close to zero at high $p_{T}$, while the azimuthal coefficient $\lambda_{\phi}$ is close to zero at low $p_{T}$, and becomes slightly negative at higher $p_{T}$. The frame-independent coefficient $\tilde{\lambda}$ is strongly negative at all $p_{T}$ in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

We present two-dimensional (2D) two-particle angular correlations on relative pseudorapidity $\eta$ and azimuth $\phi$ for charged particles from Au-Au collisions at $\sqrt{s_{\rm NN}} = 62$ and 200 GeV with transverse momentum $p_t \geq 0.15$ GeV/$c$, $|\eta| \leq 1$ and $2\pi$ azimuth. Observed correlations include a {same-side} (relative azimuth $< \pi/2$) 2D peak, a closely-related away-side azimuth dipole, and an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton-proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling in Au-Au collisions until mid-centrality where a transition to a qualitatively different centrality trend occurs within a small centrality interval. Above the transition point the number of same-side and away-side correlated pairs increases rapidly {relative to} binary-collision scaling, the $\eta$ width of the same-side 2D peak also increases rapidly ($\eta$ elongation) and the $\phi$ width actually decreases significantly. Those centrality trends are more remarkable when contrasted with expectations of jet quenching in a dense medium. Observed centrality trends are compared to {\sc hijing} predictions and to the expected trends for semihard parton scattering and fragmentation in a thermalized opaque medium. We are unable to reconcile a semihard parton scattering and fragmentation origin for the observed correlation structure and centrality trends with heavy ion collision scenarios which invoke rapid parton thermalization. On the other hand, if the collision system is effectively opaque to few-GeV partons the observations reported here would be inconsistent with a minijet picture.