Event-by-event fluctuations of the average transverse momentum of produced particles near mid-rapidity have been measured by the PHENIX Collaboration in sqrt(s_NN)=200 GeV Au+Au and p+p collisions at the Relativistic Heavy Ion Collider. The fluctuations are observed to be in excess of the expectation for statistically independent particle emission for all centralities. The excess fluctuations exhibit a dependence on both the centrality of the collision and on the transverse momentum window over which the average is calculated. Both the centrality and p_T dependence can be well reproduced by a simulation of random particle production with the addition of contributions from hard scattering processes.

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.

The PHENIX experiment at RHIC has measured transverse energy and charged particle multiplicity at mid-rapidity in Au+Au collisions at sqrt(s_NN) = 19.6, 130 and 200 GeV as a function of centrality. The presented results are compared to measurements from other RHIC experiments, and experiments at lower energies. The sqrt(s_NN) dependence of dE_T/deta and dN_ch/deta per pair of participants is consistent with logarithmic scaling for the most central events. The centrality dependence of dE_T/deta and dN_ch/deta is similar at all measured incident energies. At RHIC energies the ratio of transverse energy per charged particle was found independent of centrality and growing slowly with sqrt(s_NN). A survey of comparisons between the data and available theoretical models is also presented.

We present inclusive charged hadron elliptic flow v_2 measured over the pseudorapidity range |\eta| < 0.35 in Au+Au collisions at sqrt(s_NN) = 200 GeV. Results for v_2 are presented over a broad range of transverse momentum (p_T = 0.2-8.0 GeV/c) and centrality (0-60%). In order to study non-flow effects that are not correlated with the reaction plane, as well as the fluctuations of v_2, we compare two different analysis methods: (1) event plane method from two independent sub-detectors at forward (|\eta| = 3.1-3.9) and beam (|\eta| > 6.5) pseudorapidities and (2) two-particle cumulant method extracted using correlations between particles detected at midrapidity. The two event-plane results are consistent within systematic uncertainties over the measured p_T and in centrality 0-40%. There is at most 20% difference of the v_2 between the two event plane methods in peripheral (40-60%) collisions. The comparisons between the two-particle cumulant results and the standard event plane measurements are discussed.

Bose-Einstein correlations of charged kaons are measured for Au+Au collisions at sqrt(s_NN) = 200 GeV and are compared to charged pion probes, which have a larger hadronic scattering cross section. Three dimensional Gaussian source radii are extracted, along with a one-dimensional kaon emission source function. The centrality dependences of the three Gaussian radii are well described by a single linear function if N_part^1/3 with zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r >~ 10 fm, although the bulk emission at lower radius is well-described by a Gaussian. The presence of a non-Gaussian tail in the kaon source reaffirms that the particle emission region in a heavy ion collision is extended, and that similar measurements with pions are not solely due to the decay of long-lived resonances.

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.

Correlations between p and pbar's at transverse momenta typical of enhanced baryon production in Au+Au collisions are reported. The PHENIX experiment measures same and opposite sign baryon pairs in Au+Au collisions at sqrt(s_NN) = 200 GeV. Correlated production of p and p^bar with the trigger particle from the range 2.5 < p_T < 4.0 GeV/c and the associated particle with 1.8 < p_T < 2.5 GeV/c is observed to be nearly independent of the centrality of the collisions. Same sign pairs show no correlation at any centrality. The conditional yield of mesons triggered by baryons (and anti-baryons) and mesons in the same pT range rises with increasing centrality, except for the most central collisions, where baryons show a significantly smaller number of associated mesons. These data are consistent with a picture in which hard scattered partons produce correlated p and p^bar in the p_T region of the baryon excess.

A three-dimensional (3D) correlation function obtained from mid-rapidity, low pT pion pairs in central Au+Au collisions at sqrt(s_NN)=200 GeV is studied. The extracted model-independent source function indicates a long range tail in the directions of the pion pair transverse momentum (out) and the beam (long). Model comparisons to these distensions indicate a proper breakup time \tau_0 ~ 9 fm/c and a mean proper emission duration \Delta\tau ~ 2 fm/c, leading to sizable emission time differences (<|\Delta \tau_LCM |> ~ 12 fm/c), partly due to resonance decays. They also suggest an outside-in 'burning' of the emission source reminiscent of many hydrodynamical models.

Detailed differential measurements of the elliptic flow for particles produced in Au+Au and Cu+Cu collisions at sqrt(s_NN) = 200 GeV are presented. Predictions from perfect fluid hydrodynamics for the scaling of the elliptic flow coefficient v_2 with eccentricity, system size and transverse energy are tested and validated. For transverse kinetic energies KE_T ~ m_T-m up to ~1 GeV, scaling compatible with the hydrodynamic expansion of a thermalized fluid is observed for all produced particles. For large values of KE_T, the mesons and baryons scale separately. A universal scaling for the flow of both mesons and baryons is observed for the full transverse kinetic energy range of the data when quark number scaling is employed. In both cases the scaling is more pronounced in terms of KE_T rather than transverse momentum.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has performed systematic measurements of phi meson production in the K+K- decay channel at midrapidity in p+p, d+Au, Cu+Cu and Au+Au collisions at sqrt(S_NN)=200 GeV. Results are presented on the phi invariant yield and the nuclear modification factor R_AA for Au+Au and Cu+Cu, and R_dA for d+Au collisions, studied as a function of transverse momentum (1&lt;p_T&lt;7 GeV/c) and centrality. In central and mid-central Au+Au collisions, the R_AA of phi exhibits a suppression relative to expectations from binary scaled p+p results. The amount of suppression is smaller than that of the neutral pion and the eta meson in the intermediate p_T range (2--5 GeV/c); whereas at higher p_T the phi, pi^0, and eta show similar suppression. The baryon (protons and anti-protons) excess observed in central Au+Au collisions at intermediate p_T is not observed for the phi meson despite the similar mass of the proton and the phi. This suggests that the excess is linked to the number of constituent quarks rather than the hadron mass. The difference gradually disappears with decreasing centrality and for peripheral collisions the R_AA values for both particles are consistent with binary scaling. Cu+Cu collisions show the same yield and suppression as Au+Au collisions for the same number of N_part. The R_dA of phi shows no evidence for cold nuclear effects within uncertainties.