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.

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.

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.

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.

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.

Two particle correlations between identified meson and baryon trigger particles with 2.5 < p_T < 4.0 GeV/c and lower p_T charged hadrons have been measured at midrapidity by the PHENIX experiment at RHIC in p+p, d+Au and Au+Au collisions at sqrt(s_NN) = 200 GeV. The probability of finding a hadron near in azimuthal angle to the trigger particle is almost identical for leading mesons and baryons for non-central Au+Au. The yield for both trigger baryons and mesons is significantly higher in Au+Au than in p+p and d+Au, except for trigger baryons in central collisions. The baryon excess is likely to arise predominantly from hard scattering processes.

Azimuthal angle \Delta\phi correlations are presented for charged hadrons from dijets for 0.4 < p_T < 10 GeV/c in Au+Au collisions at sqrt(s_NN) = 200 GeV. With increasing p_T, the away-side distribution evolves from a broad to a concave shape, then to a convex shape. Comparisons to p+p data suggest that the away-side can be divided into a partially suppressed 'head' region centered at Delta\phi ~ \pi, and an enhanced 'shoulder' region centered at Delta\phi ~ \pi +/- 1.1. The p_T spectrum for the 'head' region softens toward central collisions, consistent with the onset of jet quenching. The spectral slope for the 'shoulder' region is independent of centrality and trigger p_T, which offers constraints on energy transport mechanisms and suggests that the 'shoulder' region contains the medium response to energetic jets.

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).

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.