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.

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 experiment has studied nuclear effects in $p$$+$Al and $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV on charged hadron production at forward rapidity ($1.4<\eta<2.4$, $p$-going direction) and backward rapidity ($-2.2<\eta<-1.2$, $A$-going direction). Such effects are quantified by measuring nuclear modification factors as a function of transverse momentum and pseudorapidity in various collision multiplicity selections. In central $p$$+$Al and $p$$+$Au collisions, a suppression (enhancement) is observed at forward (backward) rapidity compared to the binary scaled yields in $p$+$p$ collisions. The magnitude of enhancement at backward rapidity is larger in $p$$+$Au collisions than in $p$$+$Al collisions, which have a smaller number of participating nucleons. However, the results at forward rapidity show a similar suppression within uncertainties. The results in the integrated centrality are compared with calculations using nuclear parton distribution functions, which show a reasonable agreement at the forward rapidity but fail to describe the backward rapidity enhancement.

A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV, and Cu+Cu collisions sqrt(s_NN) = 22.5, 62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for non-dynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. An analysis of Negative Binomial Distribution fits to the multiplicity distributions demonstrates that the heavy ion data exhibit weak clustering properties.

PHENIX has measured the centrality dependence of charged hadron p_T spectra from central Au+Au collisions at sqrt(s_NN)=130 GeV. The truncated mean p_T decreases with centrality for p_T > 2 GeV/c, indicating an apparent reduction of the contribution from hard scattering to high p_T hadron production. For central collisions the yield at high p_T is shown to be suppressed compared to binary nucleon-nucleon collision scaling of p+p data. This suppression is monotonically increasing with centrality, but most of the change occurs below 30% centrality, i.e. for collisions with less than about 140 participating nucleons. The observed p_T and centrality dependence is consistent with the particle production predicted by models including hard scattering and subsequent energy loss of the scattered partons in the dense matter created in the collisions.

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.

We report the measurement of direct photons at midrapidity in Au+Au collisions at sqrt{s_NN} = 200 GeV. The direct photon signal was extracted for the transverse-momentum range of 4 GeV/c < p_T < 22 GeV/c, using a statistical method to subtract decay photons from the inclusive-photon sample. The direct-photon nuclear-modification factor R_AA was calculated as a function of p_T for different Au+Au collision centralities using the measured p+p direct-photon spectrum and compared to theoretical predictions. R_AA was found to be consistent with unity for all centralities over the entire measured p_T range. Theoretical models that account for modifications of initial-direct-photon production due to modified-parton-distribution functions in Au and the different isospin composition of the nuclei, predict a modest change of R_AA from unity and are consistent with the data. Models with compensating effects of the quark-gluon plasma on high-energy photons, such as suppression of jet-fragmentation photons and induced-photon bremsstrahlung from partons traversing the medium, are also consistent with this measurement.

Transverse momentum spectra of electrons from Au+Au collisions at sqrt(s_NN) = 130 GeV have been measured by the PHENIX experiment at RHIC. The spectra show an excess above the background from photon conversions and light hadron decays. The electron signal is consistent with that expected from semi-leptonic decays of charm. The yield of the electron signal dN_e/dy for p_T > 0.8 GeV/c is 0.025 +/- 0.004 (stat.) +/- 0.010 (sys.) in central collisions, and the corresponding charm cross section is 380 +/- 60 (stat.) +/- 200 (sys.) micro barns per binary nucleon-nucleon collision.

We present results on the measurement of lambda and lambda^bar production in Au+Au collisions at sqrt(s_NN)=130 GeV with the PHENIX detector at RHIC. The transverse momentum spectra were measured for minimum bias and for the 5% most central events. The lambda^bar/lambda ratios are constant as a function of p_T and the number of participants. The measured net lambda density is significantly larger than predicted by models based on hadronic strings (e.g. HIJING) but in approximate agreement with models which include the gluon junction mechanism.