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.

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.

Measurements of midrapidity charged particle multiplicity distributions, $dN_{\rm ch}/d\eta$, and midrapidity transverse-energy distributions, $dE_T/d\eta$, are presented for a variety of collision systems and energies. Included are distributions for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV, Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ and 62.4 GeV, Cu$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, U$+$U collisions at $\sqrt{s_{_{NN}}}=193$ GeV, $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Centrality-dependent distributions at midrapidity are presented in terms of the number of nucleon participants, $N_{\rm part}$, and the number of constituent quark participants, $N_{q{\rm p}}$. For all $A$$+$$A$ collisions down to $\sqrt{s_{_{NN}}}=7.7$ GeV, it is observed that the midrapidity data are better described by scaling with $N_{q{\rm p}}$ than scaling with $N_{\rm part}$. Also presented are estimates of the Bjorken energy density, $\varepsilon_{\rm BJ}$, and the ratio of $dE_T/d\eta$ to $dN_{\rm ch}/d\eta$, the latter of which is seen to be constant as a function of centrality for all systems.

Neutral pion transverse momentum (pT) spectra at mid-rapidity (|y| < 0.35) were measured in Cu+Cu collisions at \sqrt s_NN = 22.4, 62.4, and 200 GeV. Relative to pi -zero yields in p+p collisions scaled by the number of inelastic nucleon-nucleon collisions (Ncoll) at the respective energies, the pi-zero yields for pT \ge 2 GeV/c in central Cu+Cu collisions at 62.4 and 200 GeV are suppressed, whereas an enhancement is observed at 22.4 GeV. A comparison with a jet quenching model suggests that final state parton energy loss dominates in central Cu+Cu collisions at 62.4 GeV and 200 GeV, while the enhancement at 22.4 GeV is consistent with nuclear modifications in the initial state alone.

Measurements of transverse-single-spin asymmetries ($A_{N}$) in $p$$+$$p$ collisions at $\sqrt{s}=$62.4 and 200 GeV with the PHENIX detector at RHIC are presented. At midrapidity, $A_{N}$ is measured for neutral pion and eta mesons reconstructed from diphoton decay, and at forward rapidities, neutral pions are measured using both diphotons and electromagnetic clusters. The neutral-pion measurement of $A_{N}$ at midrapidity is consistent with zero with uncertainties a factor of 20 smaller than previous publications, which will lead to improved constraints on the gluon Sivers function. At higher rapidities, where the valence quark distributions are probed, the data exhibit sizable asymmetries. In comparison with previous measurements in this kinematic region, the new data extend the kinematic coverage in $\sqrt{s}$ and $p_T$, and it is found that the asymmetries depend only weakly on $\sqrt{s}$. The origin of the forward $A_{N}$ is presently not understood quantitatively. The extended reach to higher $p_T$ probes the transition between transverse momentum dependent effects at low $p_T$ and multi-parton dynamics at high $p_T$.

Measurements of the fractional momentum loss ($S_{\rm loss}\equiv{\delta}p_T/p_T$) of high-transverse-momentum-identified hadrons in heavy ion collisions are presented. Using $\pi^0$ in Au$+$Au and Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=62.4$ and 200 GeV measured by the PHENIX experiment at the Relativistic Heavy Ion Collider and and charged hadrons in Pb$+$Pb collisions measured by the ALICE experiment at the Large Hadron Collider, we studied the scaling properties of $S_{\rm loss}$ as a function of a number of variables: the number of participants, $N_{\rm part}$, the number of quark participants, $N_{\rm qp}$, the charged-particle density, $dN_{\rm ch}/d\eta$, and the Bjorken energy density times the equilibration time, $\varepsilon_{\rm Bj}\tau_{0}$. We find that the $p_T$ where $S_{\rm loss}$ has its maximum, varies both with centrality and collision energy. Above the maximum, $S_{\rm loss}$ tends to follow a power-law function with all four scaling variables. The data at $\sqrt{s_{_{NN}}}$=200 GeV and 2.76 TeV, for sufficiently high particle densities, have a common scaling of $S_{\rm loss}$ with $dN_{\rm ch}/d\eta$ and $\varepsilon_{\rm Bj}\tau_{0}$, lending insight on the physics of parton energy loss.

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.