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.

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.

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.

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.

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.

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

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.

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

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.