We present measurements of azimuthal correlations of charged hadron pairs in $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions after subtracting an underlying event using a model that includes higher-order azimuthal anisotropy $v_2$, $v_3$, and $v_4$. After subtraction, the away-side ($\Delta\phi\sim\pi)$ of the highest transverse-momentum trigger ($p_T>4$ GeV/$c$) correlations is suppressed compared to that of correlations measured in $p$$+$$p$ collisions. At the lowest associated particle $p_T$, the away-side shape and yield are modified. These observations are consistent with the scenario of radiative-jet energy loss. For the lowest-$p_T$ trigger correlations, an away-side yield exists and we explore the dependence of the shape of the away-side within the context of an underlying-event model. Correlations are also studied differentially versus event-plane angle $\Psi_n$. The angular correlations show an asymmetry when selecting the sign of the trigger-particle azimuthal angle with respect to the $\Psi_2$ event plane. This asymmetry and the measured suppression of the pair yield out of plane is consistent with a path-length-dependent energy loss. No $\Psi_3$ dependence can be resolved within experimental uncertainties.

Differential elliptic flow (v_2) for phi mesons and (anti)deuterons (d^bar)d is measured for Au+Au collisions at sqrt(s_NN) = 200 GeV. The v_2 for phi mesons follows the trend of lighter pi^+/- and K^+/- mesons, suggesting that ordinary hadrons interacting with standard hadronic cross sections are not the primary driver for elliptic flow development. The v_2 values for (d^bar)d suggest that elliptic flow is additive for composite particles. This further validation of the universal scaling of v_2 per constituent quark for baryons and mesons suggests that partonic collectivity dominates the transverse expansion dynamics.

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.

The production of low mass e+e- pairs for m_{e+e-} < 300 MeV/c^2 and 1 < p_T <5 GeV/c is measured in p+p and Au+Au collisions at sqrt(s_NN)=200 GeV. Enhanced yield above hadronic sources is observed. Treating the excess as internal conversions, the invariant yield of direct photons is deduced. In central Au+Au collisions, the excess of direct photon yield over p+p is exponential in transverse momentum, with inverse slope T = 221 +/- 19 (stat) +/- 19 (syst) MeV. Hydrodynamical models with initial temperatures ranging from 300--600 MeV at times of ~ 0.6 - 0.15 fm/c after the collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to quark gluon plasma at ~ 170 MeV.

We report on charged hadron production in deuteron-gold reactions at sqrt(s_NN) = 200 GeV. Our measurements in the deuteron-direction cover 1.4 < eta < 2.2, referred to as forward rapidity, and in the gold-direction -2.0 < eta < -1.4, referred to as backward rapidity, and a transverse momentum range p_T = 0.5-4.0 GeV/c. We compare the relative yields for different deuteron-gold collision centrality classes. We observe a suppression relative to binary collision scaling at forward rapidity, sensitive to low momentum fraction (x) partons in the gold nucleus, and an enhancement at backward rapidity, sensitive to high momentum fraction partons in the gold nucleus.

Inclusive transverse momentum spectra of eta mesons have been measured within p_T = 2-10 GeV/c at mid-rapidity by the PHENIX experiment in Au+Au collisions at sqrt(s_NN) = 200 GeV. In central Au+Au the eta yields are significantly suppressed compared to peripheral Au+Au, d+Au and p+p yields scaled by the corresponding number of nucleon-nucleon collisions. The magnitude, centrality and p_T dependence of the suppression is common, within errors, for eta and pi^0. The ratio of eta to pi^0 spectra at high p_T amounts to 0.40 < R_eta/pi^0 < 0.48 for the three systems in agreement with the world average measured in hadronic and nuclear reactions and, at large scaled momentum, in e^+e^- collisions.

New measurements are presented for charged hadron azimuthal correlations at mid-rapidity in Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV. They are compared to earlier measurements obtained at sqrt(s_NN) = 130 GeV and in Pb+Pb collisions at sqrt(s_NN) = 17.2 GeV. Sizeable anisotropies are observed with centrality and transverse momentum (p_T) dependence characteristic of elliptic flow (v_2). For a broad range of centralities, the observed magnitudes and trends of the differential anisotropy, v_2(p_T), change very little over the collision energy range sqrt(s_NN) = 62-200 GeV, indicating saturation of the excitation function for v_2 at these energies. Such a saturation may be indicative of the dominance of a very soft equation of state for sqrt(s_NN) = 62-200 GeV.

PHENIX has measured the centrality dependence of mid-rapidity pion, kaon and proton transverse momentum distributions in d+Au and p+p collisions at sqrt(s_NN) = 200 GeV. The p+p data provide a reference for nuclear effects in d+Au and previously measured Au+Au collisions. Hadron production is enhanced in d+Au, relative to independent nucleon-nucleon scattering, as was observed in lower energy collisions. The nuclear modification factor for (anti) protons is larger than that for pions. The difference increases with centrality, but is not sufficient to account for the abundance of baryon production observed in central Au+Au collisions at RHIC. The centrality dependence in d+Au shows that the nuclear modification factor increases gradually with the number of collisions suffered by each participant nucleon. We also present comparisons with lower energy data as well as with parton recombination and other theoretical models of nuclear effects on particle production.

Bose-Einstein correlations of identically charged pion pairs were measured by the PHENIX experiment at mid-rapidity in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The Bertsch-Pratt radius parameters were determined as a function of the transverse momentum of the pair and as a function of the centrality of the collision. Using the \it{full} Coulomb correction, the ratio $R_{\rm out}/R_{\rm side}$ is smaller than unity for $<k_{\rm T}>$ from 0.25 to 1.2 GeV/c and for all measured centralities. However, using recently developed partial Coulomb correction methods, we find that $R_{\rm out}/R_{\rm side}$ is 0.8-1.1 for the measured $<k_{\rm T}>$ range, and approximately constant at unity with the number of participants.

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.