Heavy quarkonia are observed to be suppressed in relativistic heavy ion collisions relative to their production in p+p collisions scaled by the number of binary collisions. In order to determine if this suppression is related to color screening of these states in the produced medium, one needs to account for other nuclear modifications including those in cold nuclear matter. In this paper, we present new measurements from the PHENIX 2007 data set of J/psi yields at forward rapidity (1.2<|y|<2.2) in Au+Au collisions at sqrt(s_NN)=200 GeV. The data confirm the earlier finding that the suppression of J/psi at forward rapidity is stronger than at midrapidity, while also extending the measurement to finer bins in collision centrality and higher transverse momentum (pT). We compare the experimental data to the most recent theoretical calculations that incorporate a variety of physics mechanisms including gluon saturation, gluon shadowing, initial-state parton energy loss, cold nuclear matter breakup, color screening, and charm recombination. We find J/psi suppression beyond cold-nuclear-matter effects. However, the current level of disagreement between models and d+Au data precludes using these models to quantify the hot-nuclear-matter suppression.

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.

The PHENIX Collaboration at the Relativistic Heavy Ion Collider has measured open heavy-flavor production in minimum bias Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV via the yields of electrons from semileptonic decays of charm and bottom hadrons. Previous heavy-flavor electron measurements indicated substantial modification in the momentum distribution of the parent heavy quarks due to the quark-gluon plasma created in these collisions. For the first time, using the PHENIX silicon vertex detector to measure precision displaced tracking, the relative contributions from charm and bottom hadrons to these electrons as a function of transverse momentum are measured in Au$+$Au collisions. We compare the fraction of electrons from bottom hadrons to previously published results extracted from electron-hadron correlations in $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV and find the fractions to be similar within the large uncertainties on both measurements for $p_T>4$ GeV/$c$. We use the bottom electron fractions in Au$+$Au and $p$$+$$p$ along with the previously measured heavy flavor electron $R_{AA}$ to calculate the $R_{AA}$ for electrons from charm and bottom hadron decays separately. We find that electrons from bottom hadron decays are less suppressed than those from charm for the region $3<p_T<4$ GeV/$c$.

We report the first measurement of transverse single-spin asymmetries in $J/\psi$ production from transversely polarized $p+p$ collisions at $\sqrt{s} = 200$ GeV with data taken by the PHENIX experiment in 2006 and 2008. The measurement was performed over the rapidity ranges $1.2 &lt; |y| &lt; 2.2$ and $ |y| &lt; 0.35$ for transverse momenta up to 6 GeV/$c$. $J/\psi$ production at RHIC is dominated by processes involving initial-state gluons, and transverse single-spin asymmetries of the $J/\psi$ can provide access to gluon dynamics within the nucleon. Such asymmetries may also shed light on the long-standing question in QCD of the $J/\psi$ production mechanism. Asymmetries were obtained as a function of $J/\psi$ transverse momentum and Feynman-$x$, with a value of $-0.086 \pm 0.026^{\rm stat} \pm 0.003^{\rm syst}$ in the forward region. This result suggests possible nonzero trigluon correlation functions in transversely polarized protons and, if well defined in this reaction, a nonzero gluon Sivers distribution function.

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 has measured open heavy-flavor production via semileptonic decay muons over the transverse momentum range 1 < pT < 6 GeV/c at forward and backward rapidity (1.4 < |y| < 2.0) in d+Au and p+p collisions at ?sNN = 200 GeV. In central d+Au collisions an enhancement (suppression) of heavy-flavor muon production is observed at backward (forward) rapidity relative to the yield in p+p collisions scaled by the number of binary collisions. Modification of the gluon density distribution in the Au nucleus contributes in terms of anti-shadowing enhancement and shadowing suppression; however, the enhancement seen at backward rapidity exceeds expectations from this effect alone. These results, implying an important role for additional cold nuclear matter effects, serves as a key baseline for heavy-quark measurements in A+A collisions and in constraining the magnitude of charmonia breakup effects at the Relativistic Heavy Ion Collider and the Large Hadron Collider.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured low mass vector meson, $\omega$, $\rho$, and $\phi$, production through the dimuon decay channel at forward rapidity ($1.2<|y|<2.2$) in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The differential cross sections for these mesons are measured as a function of both $p_T$ and rapidity. We also report the integrated differential cross sections over $1<p_T<7$ GeV/$c$ and $1.2<|y|<2.2$: $d\sigma/dy(\omega+\rho\rightarrow\mu\mu) = 80 \pm 6 \mbox{(stat)} \pm 12 \mbox{(syst)}$ nb and $d\sigma/dy(\phi\rightarrow\mu\mu) = 27 \pm 3 \mbox{(stat)} \pm 4 \mbox{(syst)}$ nb. These results are compared with midrapidity measurements and calculations.

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 report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Exploring the mass and transverse-momentum dependence of the yield, the bottom decay contribution can be isolated from charm, and quantified by comparison to {\sc pythia} and {\sc mc@nlo} simulations. The resulting $b\bar{b}$-production cross section is $\sigma^{d{\rm Au}}_{b\bar{b}}=1.37{\pm}0.28({\rm stat}){\pm}0.46({\rm syst})$~mb, which is equivalent to a nucleon-nucleon cross section of $\sigma^{NN}_{bb}=3.4\pm0.8({\rm stat}){\pm}1.1({\rm syst})\ \mu$b.

Measurements of midrapidity charged particle multiplicity distributions, $dN_{\rm ch}/d\eta$, and midrapidity transverse-energy distributions, $dE_T/d\eta$, are presented for a variety of collision systems and energies. Included are distributions for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV, Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ and 62.4 GeV, Cu$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, U$+$U collisions at $\sqrt{s_{_{NN}}}=193$ GeV, $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Centrality-dependent distributions at midrapidity are presented in terms of the number of nucleon participants, $N_{\rm part}$, and the number of constituent quark participants, $N_{q{\rm p}}$. For all $A$$+$$A$ collisions down to $\sqrt{s_{_{NN}}}=7.7$ GeV, it is observed that the midrapidity data are better described by scaling with $N_{q{\rm p}}$ than scaling with $N_{\rm part}$. Also presented are estimates of the Bjorken energy density, $\varepsilon_{\rm BJ}$, and the ratio of $dE_T/d\eta$ to $dN_{\rm ch}/d\eta$, the latter of which is seen to be constant as a function of centrality for all systems.