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.

Two particle azimuthal correlation functions are presented for charged hadrons produced in Au + Au collisions at RHIC sqrt(s_NN) = 130 GeV. The measurements permit determination of elliptic flow without event-by-event estimation of the reaction plane. The extracted elliptic flow values v_2 show significant sensitivity to both the collision centrality and the transverse momenta of emitted hadrons, suggesting rapid thermalization and relatively strong velocity fields. When scaled by the eccentricity of the collision zone, epsilon, the scaled elliptic flow shows little or no dependence on centrality for charged hadrons with relatively low p_T. A breakdown of this epsilon scaling is observed for charged hadrons with p_T > 1.0 GeV/c for the most central collisions.

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.

Fast parton probes produced by hard scattering and embedded within collisions of large nuclei have shown that partons suffer large energy loss and that the produced medium may respond collectively to the lost energy. We present measurements of neutral pion trigger particles at transverse momenta p^t_T = 4-12 GeV/c and associated charged hadrons (p^a_T = 0.5-7 GeV/c) as a function of relative azimuthal angle Delta Phi at midrapidity in Au+Au and p+p collisions at sqrt(s_NN) = 200 GeV. These data lead to two major observations. First, the relative angular distribution of low momentum hadrons, whose shape modification has been interpreted as a medium response to parton energy loss, is found to be modified only for p^t_T &lt; 7 GeV/c. At higher p^t_T, the data are consistent with unmodified or very weakly modified shapes, even for the lowest measured p^a_T. This observation presents a quantitative challenge to medium response scenarios. Second, the associated yield of hadrons opposite to the trigger particle in Au+Au relative to that in p+p (I_AA) is found to be suppressed at large momentum (IAA ~ 0.35-0.5), but less than the single particle nuclear modification factor (R_AA ~0.2).

Emission source functions are extracted from correlation functions constructed from charged pions produced at mid-rapidity in Au+Au collisions at sqrt(s_NN)=200 GeV. The source parameters extracted from these functions at low k_T, give first indications of a long tail for the pion emission source. The source extension cannot be explained solely by simple kinematic considerations. The possible role of a halo of secondary pions from resonance emissions is explored.

The azimuthal distribution of identified pi^0 and inclusive photons has been measured in sqrt{s_{NN}} = 200 GeV Au+Au collisions with the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). The second harmonic parameter (v_2) was measured to describe the observed anisotropy of the azimuthal distribution. The measured inclusive photon v_2 is consistent with the value expected for the photons from hadron decay and is also consistent with the lack of direct photon signal over the measured p_T range 1-6 GeV/c. An attempt is made to extract v_2 of direct photons.

We present azimuthal angle correlations of intermediate transverse momentum (1-4 GeV/c) hadrons from {dijets} in Cu+Cu and Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV. The away-side dijet induced azimuthal correlation is broadened, non-Gaussian, and peaked away from \Delta\phi=\pi in central and semi-central collisions in all the systems. The broadening and peak location are found to depend upon the number of participants in the collision, but not on the collision energy or beam nuclei. These results are consistent with sound or shock wave models, but pose challenges to Cherenkov gluon radiation models.

We have studied the dependence of azimuthal anisotropy $v_2$ for inclusive and identified charged hadrons in Au$+$Au and Cu$+$Cu collisions on collision energy, species, and centrality. The values of $v_2$ as a function of transverse momentum $p_T$ and centrality in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV and 62.4 GeV are the same within uncertainties. However, in Cu$+$Cu collisions we observe a decrease in $v_2$ values as the collision energy is reduced from 200 to 62.4 GeV. The decrease is larger in the more peripheral collisions. By examining both Au$+$Au and Cu$+$Cu collisions we find that $v_2$ depends both on eccentricity and the number of participants, $N_{\rm part}$. We observe that $v_2$ divided by eccentricity ($\varepsilon$) monotonically increases with $N_{\rm part}$ and scales as ${N_{\rm part}^{1/3}}$. The Cu$+$Cu data at 62.4 GeV falls below the other scaled $v_{2}$ data. For identified hadrons, $v_2$ divided by the number of constituent quarks $n_q$ is independent of hadron species as a function of transverse kinetic energy $KE_T=m_T-m$ between $0.1<KE_T/n_q<1$ GeV. Combining all of the above scaling and normalizations, we observe a near-universal scaling, with the exception of the Cu$+$Cu data at 62.4 GeV, of $v_2/(n_q\cdot\varepsilon\cdot N^{1/3}_{\rm part})$ vs $KE_T/n_q$ for all measured particles.

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.