Yields for J/psi production in Cu+Cu collisions at sqrt (s_NN)= 200 GeV have been measured by the PHENIX experiment over the rapidity range |y| < 2.2 at transverse momenta from 0 to beyond 5 GeV/c. The invariant yield is obtained as a function of rapidity, transverse momentum and collision centrality, and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data provide greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, providing a key constraint that is needed for disentangling cold and hot nuclear matter effects.

J/psi production has been measured in proton-proton collisions at sqrt(s)= 200 GeV over a wide rapidity and transverse momentum range by the PHENIX experiment at RHIC. Distributions of the rapidity and transverse momentum, along with measurements of the mean transverse momentum and total production cross section are presented and compared to available theoretical calculations. The total J/psi cross section is 3.99 +/- 0.61(stat) +/- 0.58(sys) +/- 0.40(abs) micro barns. The mean transverse momentum is 1.80 +/- 0.23(stat) +/- 0.16(sys) GeV/c.

The invariant differential cross section for inclusive neutral pion production in p+p collisions at sqrt(s_NN) = 200 GeV has been measured at mid-rapidity |eta| < 0.35 over the range 1 < p_T <~ 14 GeV/c by the PHENIX experiment at RHIC. Predictions of next-to-leading order perturbative QCD calculations are consistent with these measurements. The precision of our result is sufficient to differentiate between prevailing gluon-to-pion fragmentation functions.

Transverse momentum distributions and yields for $\pi^{\pm}$, $K^{\pm}$, $p$ and $\bar{p}$ in $p+p$ collisions at $\sqrt{s}$=200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter $T_{\rm inv}$, mean transverse momentum $<p_T>$ and yield per unit rapidity $dN/dy$ at each energy, and compare them to other measurements at different $\sqrt{s}$ in $p+p$ and $p+\bar{p}$ collisions. We also present the scaling properties such as $m_T$ scaling, $x_T$ scaling on the $p_T$ spectra between different energies. To discuss the mechanism of the particle production in $p+p$ collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.

Flow coefficients v_n for n = 2, 3, 4, characterizing the anisotropic collective flow in Au+Au collisions at sqrt(s_NN) = 200 GeV, are measured relative to event planes Ψ_n determined at large rapidity. We report v_n as a function of transverse momentum and collision centrality, and study the correlations among the event planes of different order n. The v_n are well described by hydrodynamic models which employ a Glauber Monte Carlo initial state geometry with fluctuations, providing additional constraining power on the interplay between initial conditions and the effects of viscosity as the system evolves. This new constraint improves precision of the extracted viscosity to entropy density ratio eta/s.

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.

The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) reports measurements of azimuthal dihadron correlations near midrapidity in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. These measurements complement recent analyses by experiments at the Large Hadron Collider (LHC) involving central $p$$+$Pb collisions at $\sqrt{s_{_{NN}}}$=5.02 TeV, which have indicated strong anisotropic long-range correlations in angular distributions of hadron pairs. The origin of these anisotropies is currently unknown. Various competing explanations include parton saturation and hydrodynamic flow. We observe qualitatively similar, but larger, anisotropies in $d$$+$Au collisions compared to those seen in $p$$+$Pb collisions at the LHC. The larger extracted $v_2$ values in $d$$+$Au collisions at RHIC are consistent with expectations from hydrodynamic calculations owing to the larger expected initial-state eccentricity compared with that from $p$$+$Pb collisions. When both are divided by an estimate of the initial-state eccentricity the scaled anisotropies follow a common trend with multiplicity that may extend to heavy ion data at RHIC and the LHC, where the anisotropies are widely thought to arise from hydrodynamic flow.

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.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured electrons from heavy flavor (charm and bottom) decays for 0.3 < p_T < 9 GeV/c at midrapidity (|y| < 0.35) in Au+Au collisions at sqrt(s_NN) = 200 GeV. The nuclear modification factor R_AA relative to p+p collisions shows a strong suppression in central Au+Au collisions, indicating substantial energy loss of heavy quarks in the medium produced at RHIC. A large azimuthal anisotropy, v_2, with respect to the reaction plane is observed for 0.5 < p_T < 5 GeV/c indicating non-zero heavy flavor elliptic flow. Both R_AA and v_2 show a p_T dependence different from those of neutral pions. A comparison to transport models which simultaneously describe R_AA(p_T) and v_2(p_T) suggests that the viscosity to entropy density ratio is close to the conjectured quantum lower bound, i.e., near a perfect fluid.

We present results for the charged-particle multiplicity distribution at mid-rapidity in Au - Au collisions at sqrt(s_NN)=130 GeV measured with the PHENIX detector at RHIC. For the 5% most central collisions we find $dN_{ch}/d\eta_{|\eta=0} = 622 \pm 1 (stat) \pm 41 (syst)$. The results, analyzed as a function of centrality, show a steady rise of the particle density per participating nucleon with centrality.

The PHENIX experiement has measured the electron-positron pair mass spectrum from 0 to 8 GeV/c^2 in p+p collisions at sqrt(s)=200 GeV. The contributions from light meson decays to e^+e^- pairs have been determined based on measurements of hadron production cross sections by PHENIX. They account for nearly all e^+e^- pairs in the mass region below 1 GeV/c^2. The e^+e^- pair yield remaining after subtracting these contributions is dominated by semileptonic decays of charmed hadrons correlated through flavor conservation. Using the spectral shape predicted by PYTHIA, we estimate the charm production cross section to be 544 +/- 39(stat) +/- 142(syst) +/- 200(model) \mu b, which is consistent with QCD calculations and measurements of single leptons by PHENIX.

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.

J/Psi production in p+p collisions at sqrt(s) = 200 GeV has been Measured in the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over a rapidity range of -2.2 < y < 2.2 and a transverse momentum range of 0 < pT < 9 GeV/c. The statistics available allow a detailed measurement of both the pT and rapidity distributions and are sufficient to constrain production models. The total cross section times branching ratio determined for J/Psi production is B_{ll} sigma_pp^J/psi = 178 +/- 3(stat) +/- 53(syst) +/- 18(norm) nb.

First results on charm quarkonia production in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) are presented. The yield of J/Psi's measured in the PHENIX experiment via electron-positron decay pairs at mid-rapidity for Au-Au reactions at sqrt(s_NN) = 200 GeV are analyzed as a function of collision centrality. For this analysis we have studied 49.3 million minimum bias Au-Au reactions. We present the J/Psi invariant yield dN/dy for peripheral and mid-central reactions. For the most central collisions where we observe no signal above background, we quote 90% confidence level upper limits. We compare these results with our J/Psi measurement from proton-proton reactions at the same energy. We find that our measurements are not consistent with models that predict strong enhancement relative to binary collision scaling.

The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < p_T < 4.0 GeV/c) of single electrons as a function of centrality in Au+Au collisions at sqrt(s_NN) = 200 GeV. Contributions to the raw spectra from photon conversions and Dalitz decays of light neutral mesons are measured by introducing a thin (1.7% X_0) converter into the PHENIX acceptance and are statistically removed. The subtracted ``non-photonic'' electron spectra are primarily due to the semi-leptonic decays of hadrons containing heavy quarks (charm and bottom). For all centralities, charm production is found to scale with the nuclear overlap function, T_AA. For minimum-bias collisions the charm cross section per binary collision is N_cc^bar/T_AA = 622 +/- 57 (stat.) +/- 160 (sys.) microbarns.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured J/psi production for rapidities 2.2 < y < 2.2 in Au+Au collisions at sqrt(s_NN) = 200 GeV. The J/psi invariant yield and nuclear modification factor R_AA as a function of centrality, transverse momentum and rapidity are reported. A suppression of J/psi relative to binary collision scaling of proton-proton reaction yields is observed. Models which describe the lower energy J/Psi data at the Super Proton Synchrotron (SPS) invoking only J/psi destruction based on the local medium density would predict a significantly larger suppression at RHIC and more suppression at mid rapidity than at forward rapidity. Both trends are contradicted by our data.

J/psi production in d+Au and p+p collisions at sqrt(s_NN) = 200 GeV has been measured by the PHENIX experiment at rapidities -2.2 < y < +2.4. The cross sections and nuclear dependence of J/\psi production versus rapidity, transverse momentum, and centrality are obtained and compared to lower energy p+A results and to theoretical models. The observed nuclear dependence in d+Au collisions is found to be modest, suggesting that the absorption in the final state is weak and the shadowing of the gluon distributions is small and consistent with Dokshitzer-Gribov-Lipatov-Altarelli-Parisi-based parameterizations that fit deep-inelastic scattering and Drell-Yan data at lower energies.

The first measurement of energy produced transverse to the beam direction at RHIC is presented. The mid-rapidity transverse energy density per participating nucleon rises steadily with the number of participants, closely paralleling the rise in charged-particle density, such that E_T / N_ch remains relatively constant as a function of centrality. The energy density calculated via Bjorken's prescription for the 2% most central Au+Au collisions at sqrt(s_NN)=130 GeV is at least epsilon_Bj = 4.6 GeV/fm^3 which is a factor of 1.6 larger than found at sqrt(s_NN)=17.2 GeV (Pb+Pb at CERN).

The momentum distribution of electrons from decays of heavy flavor (charm and beauty) for midrapidity |y| < 0.35 in p+p collisions at sqrt(s) = 200 GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over the transverse momentum range 0.3 < p_T < 9 GeV/c. Two independent methods have been used to determine the heavy flavor yields, and the results are in good agreement with each other. A fixed-order-plus-next-to-leading-log pQCD calculation agrees with the data within the theoretical and experimental uncertainties, with the data/theory ratio of 1.72 +/- 0.02^stat +/- 0.19^sys for 0.3 < p_T < 9 GeV/c. The total charm production cross section at this energy has also been deduced to be sigma_(c c^bar) = 567 +/- 57^stat +/- 224^sys micro barns.

We present the first measurement of photoproduction of J/psi and of two-photon production of high-mass e+e- pairs in electromagnetic (or ultra-peripheral) nucleus-nucleus interactions, using Au+Au data at sqrt(s_NN) = 200 GeV. The events are tagged with forward neutrons emitted following Coulomb excitation of one or both Au^{star} nuclei. The event sample consists of 28 events with m_{e+e-} > 2 GeV/c^2 with zero like-sign background. The measured cross sections at midrapidity of d\sigma / dy (J/psi + Xn, y=0) = 76 +/- 33 (stat) +/- 11 (syst) micro b and d^2\sigma/dm dy (e^+e^- + Xn, y=0) = 86 +/- 23 (stat) +/- 16 (syst) micro b/(GeV/c^2) for m_{e+e-} \in [2.0,2.8] GeV/c^2 are consistent with various theoretical predictions.

Results are presented from data recorded in 2009 by the PHENIX experiment at the Relativistic Heavy Ion Collider for the double-longitudinal spin asymmetry, $A_{LL}$, for $\pi^0$ and $\eta$ production in $\sqrt{s} = 200$ GeV polarized $p$$+$$p$ collisions. Comparison of the $\pi^0$ results with different theory expectations based on fits of other published data showed a preference for small positive values of gluon polarization, $\Delta G$, in the proton in the probed Bjorken $x$ range. The effect of adding the new 2009 \pz data to a recent global analysis of polarized scattering data is also shown, resulting in a best fit value $\Delta G^{[0.05,0.2]}_{\mbox{DSSV}} = 0.06^{+0.11}_{-0.15}$ in the range $0.05<x<0.2$, with the uncertainty at $\Delta \chi^2 = 9$ when considering only statistical experimental uncertainties. Shifting the PHENIX data points by their systematic uncertainty leads to a variation of the best-fit value of $\Delta G^{[0.05,0.2]}_{\mbox{DSSV}}$ between $0.02$ and $0.12$, demonstrating the need for full treatment of the experimental systematic uncertainties in future global analyses.

Muon production at forward rapidity (1.5 < |\eta| < 1.8) has been measured by the PHENIX experiment over the transverse momentum range 1 < p_T \le 3 GeV/c in sqrt(s) = 200 GeV p+p collisions at the Relativistic Heavy Ion Collider. After statistically subtracting contributions from light hadron decays an excess remains which is attributed to the semileptonic decays of hadrons carrying heavy flavor, i.e. charm quarks or, at high p_T, bottom quarks. The resulting muon spectrum from heavy flavor decays is compared to PYTHIA and a next-to-leading order perturbative QCD calculation. PYTHIA is used to determine the charm quark spectrum that would produce the observed muon excess. The corresponding differential cross section for charm quark production at forward rapidity is determined to be d\sigmac c^bar)/dy|_(y=1.6)=0.243 +/- 0.013 (stat.) +/- 0.105 (data syst.) ^(+0.049(-0.087) (PYTHIA syst.) mb.

Identified pi^[+/-] K^[+/-], p and p-bar transverse momentum spectra at mid-rapidity in sqrt(s_NN)=130 GeV Au-Au collisions were measured by the PHENIX experiment at RHIC as a function of collision centrality. Average transverse momenta increase with the number of participating nucleons in a similar way for all particle species. The multiplicity densities scale faster than the number of participating nucleons. Kaon and nucleon yields per participant increase faster than the pion yields. In central collisions at high transverse momenta (p_T greater than 2 GeV/c), anti-proton and proton yields are comparable to the pion yields.

We present a measurement of the cross section and transverse single-spin asymmetry ($A_N$) for $\eta$ mesons at large pseudorapidity from $\sqrt{s}=200$~GeV $p^{\uparrow}+p$ collisions. The measured cross section for $0.5<p_T<5.0$~GeV/$c$ and $3.0<|\eta|<3.8$ is well described by a next-to-leading-order perturbative-quantum-chromodynamics calculation. The asymmetries $A_N$ have been measured as a function of Feynman-$x$ ($x_F$) from $0.2<|x_{F}|<0.7$, as well as transverse momentum ($p_T$) from $1.0<p_T<4.5$~GeV/$c$. The asymmetry averaged over positive $x_F$ is $\langle{A_{N}}\rangle=0.061{\pm}0.014$. The results are consistent with prior transverse single-spin measurements of forward $\eta$ and $\pi^{0}$ mesons at various energies in overlapping $x_F$ ranges. Comparison of different particle species can help to determine the origin of the large observed asymmetries in $p^{\uparrow}+p$ collisions.

We present the midrapidity charged pion invariant cross sections and the ratio of $\pi^-$-to-$\pi^+$ production ($5<p_T<13$ GeV/$c$), together with the double-helicity asymmetries ($5<p_T<12$ GeV/$c$) in polarized $p$$+$$p$ collisions at $\sqrt{s} = 200$ GeV. The cross section measurements are consistent with perturbative calculations in quantum chromodynamics within large uncertainties in the calculation due to the choice of factorization, renormalization, and fragmentation scales. However, the theoretical calculation of the ratio of $\pi^-$-to-$\pi^+$ production when considering these scale uncertainties overestimates the measured value, suggesting further investigation of the uncertainties on the charge-separated pion fragmentation functions is needed. Due to cancellations of uncertainties in the charge ratio, direct inclusion of these ratio data in future parameterizations should improve constraints on the flavor dependence of quark fragmentation functions to pions. By measuring charge-separated pion asymmetries, one can gain sensitivity to the sign of $\Delta G$ through the opposite sign of the up and down quark helicity distributions in conjunction with preferential fragmentation of positive pions from up quarks and negative pions from down quarks. The double-helicity asymmetries presented are sensitive to the gluon helicity distribution over an $x$ range of $\sim$0.03--0.16.

The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < p_T < 5.0 GeV/c) of electrons as a function of centrality in Au+Au collisions at sqrt(s_NN)=200 GeV. Contributions from photon conversions and from light hadron decays, mainly Dalitz decays of pi^0 and eta mesons, were removed. The resulting non-photonic electron spectra are primarily due to the semi-leptonic decays of hadrons carrying heavy quarks. Nuclear modification factors were determined by comparison to non-photonic electrons in p+p collisions. A significant suppression of electrons at high p_T is observed in central Au+Au collisions, indicating substantial energy loss of heavy quarks.

The momentum distribution of electrons from semi-leptonic decays of charm and bottom for mid-rapidity |y|<0.35 in p+p collisions at sqrt(s)=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over the transverse momentum range 2 < p_T < 7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D^bar --> e^{+/-} K^{-/+} X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in p_T. A fixed-order-plus-next-to-leading-log (FONLL) perturbative quantum chromodynamics (pQCD) calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is \sigma_{b\b^bar}= 3.2 ^{+1.2}_{-1.1}(stat) ^{+1.4}_{-1.3}(syst) micro b.

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.