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.

Transverse momentum spectra and yields of hadrons are measured by the PHENIX collaboration in Au + Au collisions at sqrt(s_NN) = 130 GeV at the Relativistic Heavy Ion Collider (RHIC). The time-of-flight resolution allows identification of pions to transverse momenta of 2 GeV/c and protons and antiprotons to 4 GeV/c. The yield of pions rises approximately linearly with the number of nucleons participating in the collision, while the number of kaons, protons, and antiprotons increases more rapidly. The shape of the momentum distribution changes between peripheral and central collisions. Simultaneous analysis of all the p_T spectra indicates radial collective expansion, consistent with predictions of hydrodynamic models. Hydrodynamic analysis of the spectra shows that the expansion velocity increases with collision centrality and collision energy. This expansion boosts the particle momenta, causing the yield from soft processes to exceed that for hard to large transverse momentum, perhaps as large as 3 GeV/c.

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

Transverse momentum spectra of charged hadrons with p_T < 8 GeV/c and neutral pions with p_T < 10 GeV/c have been measured at mid-rapidity by the PHENIX experiment at RHIC in d+Au collisions at sqrt(s_NN) = 200 GeV. The measured yields are compared to those in p+p collisions at the same sqrt(s_NN) scaled up by the number of underlying nucleon-nucleon collisions in d+Au. The yield ratio does not show the suppression observed in central Au+Au collisions at RHIC. Instead, there is a small enhancement in the yield of high momentum particles.

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.

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.

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.

Longitudinal density correlations of produced matter in Au+Au collisions at sqrt(s_NN)=200 GeV have been measured from the inclusive charged particle distributions as a function of pseudorapidity window sizes. The extracted \alpha \xi parameter, related to the susceptibility of the density fluctuations in the long wavelength limit, exhibits a non-monotonic behavior as a function of the number of participant nucleons, N_part. A local maximum is seen at N_part ~ 90, with corresponding energy density based on the Bjorken picture of \epsilon_Bj \tau ~ 2.4 GeV/(fm^2 c) with a transverse area size of 60 fm^2. This behavior may suggest a critical phase boundary based on the Ginzburg-Landau framework.

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.

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.