Transverse momentum spectra of $\pi^{\pm}$, $p$ and $\bar{p}$ up to 12 GeV/c at mid-rapidity in centrality selected Au+Au collisions at $\sqrt{s_{_{NN}}} = 200$ GeV are presented. In central Au+Au collisions, both $\pi^{\pm}$ and $p(\bar{p})$ show significant suppression with respect to binary scaling at $p_T > $ 4 GeV/c. Protons and anti-protons are less suppressed than $\pi^{\pm}$, in the range 1.5 $< p_{T} <$6 GeV/c. The $\pi^-/\pi^+$ and $\bar{p}/p$ ratios show at most a weak $p_T$ dependence and no significant centrality dependence. The $p/\pi$ ratios in central Au+Au collisions approach the values in p+p and d+Au collisions at $p_T >$ 5 GeV/c. The results at high $p_T$ indicate that the partonic sources of $\pi^{\pm}$, $p$ and $\bar{p}$ have similar energy loss when traversing the nuclear medium.

We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.

New measurements by the PHENIX experiment at RHIC for eta production at midrapidity as a function of transverse momentum (p_T) and collision centrality in sqrt(s_NN) = 200 GeV Au+Au and p+p collisions are presented. They indicate nuclear modification factors (R_AA) that are similar both in magnitude and trend to those found in earlier pi^0 measurements. Linear fits to R_AA in the 5--20 GeV/c p_T region show that the slope is consistent with zero within two standard deviations at all centralities although a slow rise cannot be excluded. Having different statistical and systematic uncertainties the pi^0 and eta measurements are complementary at high p_T/ thus, along with the extended p_T range of these data they can provide additional constraints for theoretical modeling and the extraction of transport properties.

Differential measurements of the elliptic (v_2) and hexadecapole (v_4) Fourier flow coefficients are reported for charged hadrons as a function of transverse momentum (p_T) and collision centrality or the number of participant nucleons (N_part) for Au+Au collisions at sqrt(s_NN)=200 GeV. The v_{2,4} measurements at pseudorapidity |\eta|<=0.35 obtained with four separate reaction plane detectors positioned in the range 1.0<|\eta|<3.9 show good agreement, indicating the absence of significant \eta-dependent nonflow perturbations. Sizable values for v_4(p_T) are observed with a ratio v_4(p_T,N_part)/v_2^2(p_T,N_part)~0.8 for 50<N_part<200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N_part>200 this ratio increases up to 1.7 in the most central collisions.

We present measurements of the J/psi invariant yields in sqrt(s_NN)=39 and 62.4 GeV Au+Au collisions at forward rapidity (1.2<|y|<2.2). Invariant yields are presented as a function of both collision centrality and transverse momentum. Nuclear modifications are obtained for central relative to peripheral Au+Au collisions (R_CP) and for various centrality selections in Au+Au relative to scaled p+p cross sections obtained from other measurements (R_AA). The observed suppression patterns at 39 and 62.4 GeV are quite similar to those previously measured at 200 GeV. This similar suppression presents a challenge to theoretical models that contain various competing mechanisms with different energy dependencies, some of which cause suppression and others enhancement.

We report the measurement of the transverse momentum dependence of inclusive J/psi polarization in p+p collisions at sqrt(s)=200 GeV performed by the PHENIX Experiment at RHIC. The polarization is studied in the helicity, Gottfried-Jackson, and Collins-Soper frames for p_T < 5 GeV/c and |y| < 0.35. The J/psi polarization in the helicity and Gottfried-Jackson frames is consistent with zero for all transverse momenta, with a slight (1.8 sigma) trend towards longitudinal polarization for transverse momenta above 2 GeV/c. No conclusion is allowed due to the limited acceptance in the Collins-Soper frame and the uncertainties of the current data. The results are compared to observations for other collision systems and center of mass energies and to different quarkonia production models.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured 2nd and 3rd order Fourier coefficients of the azimuthal distributions of direct photons emitted at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for various collision centralities. Combining two different analysis techniques, results were obtained in the transverse momentum range of $0.4<p_{T}<4.0$ GeV/$c$. At low $p_T$ the second-order coefficients, $v_2$, are similar to the ones observed in hadrons. Third order coefficients, $v_3$, are nonzero and almost independent of centrality. These new results on $v_2$ and $v_3$, combined with previously published results on yields, are compared to model calculations that provide yields and asymmetries in the same framework. Those models are challenged to explain simultaneously the observed large yield and large azimuthal anisotropies.

It has been postulated that partonic orbital angular momentum can lead to a significant double-helicity dependence in the net transverse momentum of Drell-Yan dileptons produced in longitudinally polarized p+p collisions. Analogous effects are also expected for dijet production. If confirmed by experiment, this hypothesis, which is based on semi-classical arguments, could lead to a new approach for studying the contributions of orbital angular momentum to the proton spin. We report the first measurement of the double-helicity dependence of the dijet transverse momentum in longitudinally polarized p+p collisions at sqrt(s) = 200 GeV from data taken by the PHENIX experiment in 2005 and 2006. The analysis deduces the transverse momentum of the dijet from the widths of the near- and far-side peaks in the azimuthal correlation of the dihadrons. When averaged over the transverse momentum of the triggered particle, the difference of the root-mean-square of the dijet transverse momentum between like- and unlike-helicity collisions is found to be -37 +/- 88(stat) +/- 14(syst) MeV/c.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K^0_S , \omega, \eta prime, and \phi mesons in p + p collisions at = 200 GeV. Measurements \omega and \phi production in different decay channels give consistent results. New results for the \phi are in agreement with previously published data and extend the measured pT coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-pT and characterizing the low-pT regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

Neutral pion transverse momentum (pT) spectra at mid-rapidity (|y| < 0.35) were measured in Cu+Cu collisions at \sqrt s_NN = 22.4, 62.4, and 200 GeV. Relative to pi -zero yields in p+p collisions scaled by the number of inelastic nucleon-nucleon collisions (Ncoll) at the respective energies, the pi-zero yields for pT \ge 2 GeV/c in central Cu+Cu collisions at 62.4 and 200 GeV are suppressed, whereas an enhancement is observed at 22.4 GeV. A comparison with a jet quenching model suggests that final state parton energy loss dominates in central Cu+Cu collisions at 62.4 GeV and 200 GeV, while the enhancement at 22.4 GeV is consistent with nuclear modifications in the initial state alone.