The PHENIX experiment at the Relativistic Heavy Ion Collider has measured omega meson production via leptonic and hadronic decay channels in p + p, d + Au, Cu+ Cu, and Au + Au collisions at root s(NN) = 200 GeV. The invariant transverse momentum spectra measured in different decay modes give consistent results. Measurements in the hadronic decay channel in Cu Cu and Au + Au collisions show that. production has a suppression pattern at high transverse momentum, similar to that of pi(0) and eta in central collisions, but no suppression is observed in peripheral collisions. The nuclear modification factors, R-AA, are consistent in Cu + Cu and Au + Au collisions at similar numbers of participant nucleons.

The PHENIX experiment at the Relativistic Heavy Ion Collider 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 root 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 midcentral 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 pi(0) and the. in the intermediate p(T) range (2-5 GeV/c), whereas, at higher p(T), the phi, pi(0), and. show similar suppression. The baryon (proton and antiproton) excess observed in central Au + Au collisions at intermediate p(T) is not observed for the phi meson despite the similar masses of the proton and the phi. This suggests that the excess is linked to the number of valence quarks in the hadron rather than its mass. The difference gradually disappears with decreasing centrality, and, for peripheral collisions, the R-AA values for both particle species 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.

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 &lt;radical&gt;&lt;radicand&gt;s(NN)&lt;/radicand&gt;&lt;/radical&gt;=62.4 and 200 GeV, and Cu+Cu collisions at &lt;radical&gt;&lt;radicand&gt;s(NN)&lt;/radicand&gt;&lt;/radical&gt;=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 nondynamical 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. A comparison of the data with a model where hadrons are independently emitted from a number of hadron clusters suggests that the mean number of hadrons per cluster is small in heavy ion collisions.

We present azimuthal angle correlations of intermediate transverse momentum (1-4 GeV/c) hadrons from dijets in Cu+Cu and Au+Au collisions at root 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 semicentral 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 PHENIX experiment at the Relativistic Heavy Ion Collider has performed a systematic study of K-S(0) and K*(0) meson production at midrapidity in p + p, d + Au, and Cu + Cu collisions at root s(NN) = 200 GeV. The K-S(0) and K*(0) mesons are reconstructed via their K-S(0) -&gt; pi(0)(-&gt; gamma gamma) pi(0)(-&gt; gamma gamma) and K*(0) -&gt; K-+/-pi(-/+) decay modes, respectively. The measured transverse-momentum spectra are used to determine the nuclear modification factor of K-S(0) and K*(0) mesons in d + Au and Cu + Cu collisions at different centralities. In the d + Au collisions, the nuclear modification factor of K-S(0) and K*(0) mesons is almost constant as a function of transverse momentum and is consistent with unity, showing that cold-nuclear-matter effects do not play a significant role in the measured kinematic range. In Cu + Cu collisions, within the uncertainties no nuclear modification is registered in peripheral collisions. In central collisions, both mesons show suppression relative to the expectations from the p + p yield scaled by the number of binary nucleon-nucleon collisions in the Cu + Cu system. In the p(T) range 2-5 GeV/c, the strange mesons (K-S(0), K*(0)) similarly to the phi meson with hidden strangeness, showan intermediate suppression between the more suppressed light quark mesons (pi(0)) and the nonsuppressed baryons (p, (p) over bar). At higher transverse momentum, p(T) &gt; 5 GeV/c, production of all particles is similarly suppressed by a factor of approximate to 2.

Measurements of the midrapidity transverse-energy distribution, dE(T)/d eta, are presented for p + p, d+Au, and Au+Au collisions atv root s(NN) = 200 GeV and additionally for Au+Au collisions atv root s(NN) = 62.4 and 130 GeV. The dE(T)/d eta distributions are first compared with the number of nucleon participants N-part, number of binary collisions N-coll, and number of constituent-quark participants N-qp calculated from a Glauber model based on the nuclear geometry. For Au+Au, &lt; dE(T)/d eta &gt;/N-part increases with N-part, while &lt; dE(T)/d eta &gt;/N-qp is approximately constant for all three energies. This indicates that the two-component ansatz, dE(T)/d eta alpha (1 - x)N-part/2 + xN(coll), which was used to represent E-T distributions, is simply a proxy for N-qp, and that the N-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.

We present a new analysis of J/psi production yields in deuteron-gold collisions at root s(NN) =200 GeV using data taken from the PHENIX experiment in 2003 and previously published in S. S. Adler Phys. Rev. Lett 96, 012304 (2006). The high statistics proton-proton J/psi data taken in 2005 are used to improve the baseline measurement and thus construct updated cold nuclear matter modification factors (R-dAu). A suppression of J/psi in cold nuclear matter is observed as one goes forward in rapidity (in the deuteron-going direction), corresponding to a region more sensitive to initial-state low-x gluons in the gold nucleus. The measured nuclear modification factors are compared to theoretical calculations of nuclear shadowing to which a J/psi (or precursor) breakup cross section is added. Breakup cross sections of sigma(breakup)=2.8(-1.4)(+1.7) (2.2(-1.5)(+1.6)) mb are obtained by fitting these calculations to the data using two different models of nuclear shadowing. These breakup cross-section values are consistent within large uncertainties with the 4.2 +/- 0.5 mb determined at lower collision energies. Projecting this range of cold nuclear matter effects to copper-copper and gold-gold collisions reveals that the current constraints are not sufficient to firmly quantify the additional hot nuclear matter effect.

We present a measurement of the cross section and transverse single-spin asymmetry (AN) for. mesons at large pseudorapidity from root s = 200 GeV p up arrow + p collisions. The measured cross section for 0.5 &lt; p(T) &lt; 5.0 GeV/c and 3.0 &lt; vertical bar eta vertical bar &lt; 3.8 is well described by a next-to-leading-order perturbative-quantum-chromodynamics calculation. The asymmetries A(N) have been measured as a function of Feynman-x (x(F)) from 0.2 &lt; vertical bar x(F)vertical bar &lt; 0.7, as well as transverse momentum (p(T)) from 1.0 &lt; p(T) &lt; 4.5 GeV/c. The asymmetry averaged over positive x(F) is &lt; A(N)&gt; = 0.061 +/- 0.014. The results are consistent with prior transverse single-spin measurements of forward eta and pi(0) mesons at various energies in overlapping x(F) ranges. Comparison of different particle species can help to determine the origin of the large observed asymmetries in p up arrow + p collisions.

We present measurements of J/psi yields in d + Au collisions at root S-NN = 200 GeV recorded by the PHENIX experiment and compare them with yields in p + p collisions at the same energy per nucleon-nucleon collision. The measurements cover a large kinematic range in J/psi rapidity (-2.2 &lt; y &lt; 2.4) with high statistical precision and are compared with two theoretical models: one with nuclear shadowing combined with final state breakup and one with coherent gluon saturation effects. In order to remove model dependent systematic uncertainties we also compare the data to a simple geometric model. The forward rapidity data are inconsistent with nuclear modifications that are linear or exponential in the density weighted longitudinal thickness, such as those from the final state breakup of the bound state.