Measurements of the midrapidity transverse energy distribution, $d\Et/d\eta$, are presented for $p$$+$$p$, $d$$+$Au, and Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV and additionally for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ and 130 GeV. The $d\Et/d\eta$ distributions are first compared with the number of nucleon participants $N_{\rm part}$, number of binary collisions $N_{\rm coll}$, and number of constituent-quark participants $N_{qp}$ calculated from a Glauber model based on the nuclear geometry. For Au$+$Au, $\mean{d\Et/d\eta}/N_{\rm part}$ increases with $N_{\rm part}$, while $\mean{d\Et/d\eta}/N_{qp}$ is approximately constant for all three energies. This indicates that the two component ansatz, $dE_{T}/d\eta \propto (1-x) N_{\rm part}/2 + x N_{\rm coll}$, which has been used to represent $E_T$ distributions, is simply a proxy for $N_{qp}$, and that the $N_{\rm coll}$ term does not represent a hard-scattering component in $E_T$ distributions. The $dE_{T}/d\eta$ distributions of Au$+$Au and $d$$+$Au are then calculated from the measured $p$$+$$p$ $E_T$ distribution using two models that both reproduce the Au$+$Au data. However, while the number-of-constituent-quark-participant model agrees well with the $d$$+$Au data, the additive-quark model does not.

Direct photons have been measured in sqrt(s_NN)=200 GeV d+Au collisions at midrapidity. A wide p_T range is covered by measurements of nearly-real virtual photons (1<p_T<6 GeV/c) and real photons (5<p_T<16 GeV/c). The invariant yield of the direct photons in d+Au collisions over the scaled p+p cross section is consistent with unity. Theoretical calculations assuming standard cold nuclear matter effects describe the data well for the entire p_T range. This indicates that the large enhancement of direct photons observed in Au+Au collisions for 1.0<p_T<2.5 GeV/c is due to a source other than the initial-state nuclear effects.

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.

We present transverse momentum (p_T) spectra of charged hadrons measured in deuteron-gold and nucleon-gold collisions at \sqrts = 200 GeV for four centrality classes. Nucleon-gold collisions were selected by tagging events in which a spectator nucleon was observed in one of two forward rapidity detectors. The spectra and yields were investigated as a function of the number of binary nucleon-nucleon collisions, \nu, suffered by deuteron nucleons. A comparison of charged particle yields to those in p+p collisions show that the yield per nucleon-nucleon collision saturates with \nu for high momentum particles. We also present the charged hadron to neutral pion ratios as a function of p_T.

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.

Inclusive transverse momentum spectra of eta mesons in the range p_T~2-12 GeV/c have been measured at mid-rapidity (|\eta| < 0,35) by the PHENIX experiment at RHIC in p+p, d+Au and Au+Au collisions at sqrt(s_NN) = 200 GeV. The eta mesons are reconstructed through their eta--> \gamma\gamma channel for the three colliding systems as well as through the eta-->pi^0 pi+ pi- decay mode in p+p and d+Au collisions. The nuclear modification factor in d+Au collisions, R_dAu(p_T~1.0-1.1, suggests at most only modest p_T broadening (Cronin enhancement). In central Au+Au reactions, the eta yields are significantly suppressed, with R_AuAu(pT)~0.2. The ratio of eta to pi^0 yields is approximately constant as a function of p_T for the three colliding systems in agreement with the high-p_T world average of R_eta/pi^0 \approx 0.5 in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions for a wide range of center-of-mass energies [sqrt(s_NN)~3-1800 GeV] as well as, for high scaled momentum x_p, in e+e- annihilations at sqrt(s)=91.2 GeV. These results are consistent with a scenario where high-p_T eta production in nuclear collisions at RHIC is largely unaffected by initial-state effects, but where light-quark mesons (pi^0:eta) are equally suppressed due to final-state interactions of the parent partons in the dense medium produced in Au+Au reactions.

Measurements of neutral pion production at midrapidity in sqrt(s_NN) = 200 GeV Au+Au collisions as a function of transverse momentum, p_T, collision centrality, and angle with respect to reaction plane are presented. The data represent the final pi^0 results from the PHENIX experiment for the first RHIC Au+Au run at design center-of-mass-energy. They include additional data obtained using the PHENIX Level-2 trigger with more than a factor of three increase in statistics over previously published results for p_T > 6 GeV/c. We evaluate the suppression in the yield of high-p_T pi^0's relative to point-like scaling expectations using the nuclear modification factor R_AA. We present the p_T dependence of R_AA for nine bins in collision centrality. We separately integrate R_AA over larger p_T bins to show more precisely the centrality dependence of the high-p_T suppression. We then evaluate the dependence of the high-p_T suppression on the emission angle \Delta\phi of the pions with respect to event reaction plane for 7 bins in collision centrality. We show that the yields of high-p_T pi^0's vary strongly with \Delta\phi, consistent with prior measurements. We show that this variation persists in the most peripheral bin accessible in this analysis. For the peripheral bins we observe no suppression for neutral pions produced aligned with the reaction plane while the yield of pi^0's produced perpendicular to the reaction plane is suppressed by more than a factor of 2. We analyze the combined centrality and \Delta\phi dependence of the pi^0 suppression in different p_T bins using different possible descriptions of parton energy loss dependence on jet path-length averages to determine whether a single geometric picture can explain the observed suppression pattern.

The PHENIX experiment at RHIC has measured the invariant cross section for omega-meson production at mid-rapidity in the transverse momentum range 2.5 < p_T < 9.25 GeV/c in p+p and d+Au collisions at sqrt(s_NN) = 200 GeV. Measurements in two decay channels (omega --> pi^0 pi^+ pi^- and omega --> pi^0 gamma) yield consistent results, and the reconstructed omega mass agrees with the accepted value within the p_T range of the measurements. The omega/pi^0 ratio is found to be 0.85 +/- 0.05(stat) +/- 0.09(sys) and 0.94 +/- 0.08(stat) +/- 0.12(sys) in p+p and d+Au collisions respectively, independent of p_T . The nuclear modification factor R_dA is 1.03 +/- 0.12(stat) +/- 0.21(sys) and 0.83 +/- 0.21(stat) +/- 0.17(sys) in minimum bias and central (0-20%) d+Au collisions, respectively.

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.

The dependence of transverse momentum spectra of neutral pions and eta mesons with p_T <16 GeV/c and p_T < 12 GeV/c, respectively, on the centrality of the collision has 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 by the number of underlying nucleon-nucleon collisions in d+Au. At all centralities the yield ratios show no suppression, in contrast to the strong suppression seen for central Au+Au collisions at RHIC. Only a weak p_T and centrality dependence can be observed.