We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.

New measurements by the PHENIX experiment at RHIC for eta production at midrapidity as a function of transverse momentum (p_T) and collision centrality in sqrt(s_NN) = 200 GeV Au+Au and p+p collisions are presented. They indicate nuclear modification factors (R_AA) that are similar both in magnitude and trend to those found in earlier pi^0 measurements. Linear fits to R_AA in the 5--20 GeV/c p_T region show that the slope is consistent with zero within two standard deviations at all centralities although a slow rise cannot be excluded. Having different statistical and systematic uncertainties the pi^0 and eta measurements are complementary at high p_T/ thus, along with the extended p_T range of these data they can provide additional constraints for theoretical modeling and the extraction of transport properties.

We present inclusive charged hadron elliptic flow v_2 measured over the pseudorapidity range |\eta| < 0.35 in Au+Au collisions at sqrt(s_NN) = 200 GeV. Results for v_2 are presented over a broad range of transverse momentum (p_T = 0.2-8.0 GeV/c) and centrality (0-60%). In order to study non-flow effects that are not correlated with the reaction plane, as well as the fluctuations of v_2, we compare two different analysis methods: (1) event plane method from two independent sub-detectors at forward (|\eta| = 3.1-3.9) and beam (|\eta| > 6.5) pseudorapidities and (2) two-particle cumulant method extracted using correlations between particles detected at midrapidity. The two event-plane results are consistent within systematic uncertainties over the measured p_T and in centrality 0-40%. There is at most 20% difference of the v_2 between the two event plane methods in peripheral (40-60%) collisions. The comparisons between the two-particle cumulant results and the standard event plane measurements are discussed.

A measurement of direct photons in p+p collisions at sqrt(s)=200 GeV is presented. A photon excess above background from pi^0 --> gamma+gamma, eta --> gamma+gamma, and other decays is observed in the transverse momentum range 5.5 < p_T < 7 GeV/c. The result is compared to a next-to-leading-order perturbative QCD calculation. Within errors, good agreement is found between the QCD calculation and the measured result.

The properties of jets produced in p+p collisions at sqrt(s)=200 GeV are measured using the method of two particle correlations. The trigger particle is a leading particle from a large transverse momentum jet while the associated particle comes from either the same jet or the away-side jet. Analysis of the angular width of the near-side peak in the correlation function determines the jet fragmentation transverse momentum j_T . The extracted value, sqrt(<j_T^2>)= 585 +/- 6(stat) +/- 15(sys) MeV/c, is constant with respect to the trigger particle transverse momentum, and comparable to the previous lower sqrt(s) measurements. The width of the away-side peak is shown to be a convolution of j_T with the fragmentation variable, z, and the partonic transverse momentum, k_T . The <z> is determined through a combined analysis of the measured pi^0 inclusive and associated spectra using jet fragmentation functions measured in e^+e^-. collisions. The final extracted values of k_T are then determined to also be independent of the trigger particle transverse momentum, over the range measured, with value of sqrt(<k_T^2>) = 2.68 +/- 0.07(stat) +/- 0.34(sys) GeV/c.

The invariant differential cross section for inclusive electron production in $p + p$ collisions at $\sqrt{s} = 200$~GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range $0.4 \le p_T \le 5.0$~GeV/$c$ in the central rapidity region ($|\eta| \le 0.35$). The contribution to the inclusive electron spectrum from semileptonic decays of hadrons carrying heavy flavor, {\it i.e.} charm quarks or, at high $p_T$, bottom quarks, is determined via three independent methods. The resulting electron spectrum from heavy flavor decays is compared to recent leading and next-to-leading order perturbative QCD calculations. The total cross section of charm quark-antiquark pair production is determined to be $\sigma_{c\bar{c}} = 0.92 \pm 0.15 {\rm (stat.)} \pm 0.54 {\rm (sys.)}$~mb.

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.

Dihadron correlations at high transverse momentum in d+Au collisions at sqrt(s_NN) = 200 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). From these correlations we extract several structural characteristics of jets; the root-mean-squared (RMS) transverse momentum of fragmenting hadrons with respect to the jet sqrt(<j_T^2>), the mean sine-squared angle between the scattered partons <sin^2(phi_jj)>, and the number of particles produced within the dijet that are associated with a high-p_T particle (dN/dx_E distributions). We observe that the fragmentation characteristics of jets in d+Au collisions are very similar to those in p+p collisions and that there is also little dependence on the centrality of the d+Au collision. This is consistent with the nuclear medium having little influence on the fragmentation process. Furthermore, there is no statistically significant increase in the value of <sin^2(phi_jj)> from p+p to d+Au collisions. This constrains the amount of multiple scattering that partons undergo in the cold nuclear medium before and after a hard-collision.

Production of $\pi^0$ and $\eta$ mesons has been measured at midrapidity in Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. Measurements were performed in $\pi^0(\eta)\rightarrow\gamma\gamma$ decay channel in the 1(2)-20 GeV/$c$ transverse momentum range. A strong suppression is observed for $\pi^0$ and $\eta$ meson production at high transverse momentum in central Cu$+$Au collisions relative to the $p$$+$$p$ results scaled by the number of nucleon-nucleon collisions. In central collisions the suppression is similar to Au$+$Au with comparable nuclear overlap. The $\eta/\pi^0$ ratio measured as a function of transverse momentum is consistent with $m_T$-scaling parameterization down to $p_T=$2 GeV/$c$, its asymptotic value is constant and consistent with Au$+$Au and $p$$+$$p$ and does not show any significant dependence on collision centrality. Similar results were obtained in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions as well as in $e^+e^-$ collisions in a range of collision energies $\sqrt{s_{_{NN}}}=$3--1800 GeV. This suggests that the quark-gluon-plasma medium produced in Cu$+$Cu collisions either does not affect the jet fragmentation into light mesons or it affects the $\pi^0$ and $\eta$ the same way.

We present measurements of the J/psi invariant yields in sqrt(s_NN)=39 and 62.4 GeV Au+Au collisions at forward rapidity (1.2<|y|<2.2). Invariant yields are presented as a function of both collision centrality and transverse momentum. Nuclear modifications are obtained for central relative to peripheral Au+Au collisions (R_CP) and for various centrality selections in Au+Au relative to scaled p+p cross sections obtained from other measurements (R_AA). The observed suppression patterns at 39 and 62.4 GeV are quite similar to those previously measured at 200 GeV. This similar suppression presents a challenge to theoretical models that contain various competing mechanisms with different energy dependencies, some of which cause suppression and others enhancement.

The PHENIX experiment at RHIC has measured charged hadron yields at mid-rapidity over a wide range of transverse momentum (0.5 < p_T < 10 GeV/c) in Au+Au collisions at sqrt(s_NN)=200 GeV. The data are compared to pi^zero measurements from the same experiment. For both charged hadrons and neutral pions, the yields per nucleon-nucleon collision are significantly suppressed in central compared to peripheral and nucleon-nucleon collisions. The suppression sets in gradually and increases with increasing centrality of the collisions. Above 4-5 GeV/c in p_T, a constant and almost identical suppression of charged hadrons and pi^zeroes is observed. The p_T spectra are compared to published spectra from Au+Au at sqrt(s_NN)=130 in terms of x_t scaling. Central and peripheral pi^zero as well as peripheral charged spectra exhibit the same x_t scaling as observed in p+p data.

We report the measurement of the transverse momentum dependence of inclusive J/psi polarization in p+p collisions at sqrt(s)=200 GeV performed by the PHENIX Experiment at RHIC. The polarization is studied in the helicity, Gottfried-Jackson, and Collins-Soper frames for p_T < 5 GeV/c and |y| < 0.35. The J/psi polarization in the helicity and Gottfried-Jackson frames is consistent with zero for all transverse momenta, with a slight (1.8 sigma) trend towards longitudinal polarization for transverse momenta above 2 GeV/c. No conclusion is allowed due to the limited acceptance in the Collins-Soper frame and the uncertainties of the current data. The results are compared to observations for other collision systems and center of mass energies and to different quarkonia production models.

Two-pion correlations in sqrt(s_NN)=130 GeV Au+Au collisions at RHIC have been measured over a broad range of pair transverse momentum k_T by the PHENIX experiment at RHIC. The k_T dependent transverse radii are similar to results from heavy ion collisions at sqrt(s_NN) = 4.1, 4.9, and 17.3 GeV, whereas the longitudinal radius increases monotonically with beam energy. The ratio of the outwards to sidewards transverse radii (R_out/R_side) is consistent with unity and independent of k_T.

The PHENIX experiment has measured the spin alignment for inclusive $J/\psi\rightarrow e^{+}e^{-}$ decays in $p$+$p$ collisions at $\sqrt{s}=510$ GeV at midrapidity. The angular distributions have been measured in three different polarization frames, and the three decay angular coefficients have been extracted in a full two-dimensional analysis. Previously, PHENIX saw large longitudinal net polarization at forward rapidity at the same collision energy. This analysis at midrapidity, complementary to the previous PHENIX results, sees no sizable polarization in the measured transverse momentum range of $0.0<p_T<10.0$ GeV/$c$. The results are consistent with a previous one-dimensional analysis at midrapidity at $\sqrt{s}=200$ GeV. The transverse-momentum-dependent cross section for midrapidity $J/\psi$ production has additionally been measured, and after comparison to world data we find a simple logarithmic dependence of the cross section on $\sqrt{s}$.

We report the first measurement of the fraction of $J/\psi$ mesons coming from $B$-meson decay ($F_{B{\rightarrow}J/\psi}$) in $p$+$p$ collisions at $\sqrt{s}=$ 510 GeV. The measurement is performed using the forward silicon vertex detector and central vertex detector at PHENIX, which provide precise tracking and distance-of-closest-approach determinations, enabling the statistical separation of $J/\psi$ due to $B$-meson decays from prompt $J/\psi$. The measured value of $F_{B{\rightarrow}J/\psi}$ is 8.1\%$\pm$2.3\% (stat)$\pm$1.9\% (syst) for $J/\psi$ with transverse momenta $0<p_T<5$ GeV/$c$ and rapidity $1.2<|y|<2.2$. The measured fraction $F_{B{\rightarrow}J/\psi}$ at PHENIX is compared to values measured by other experiments at higher center of mass energies and to fixed-order-next-to-leading-logarithm and color-evaporation-model predictions. The $b\bar{b}$ cross section per unit rapidity ($d\sigma/dy(pp{\rightarrow}b\bar{b})$) extracted from the obtained $F_{B{\rightarrow}J/\psi}$ and the PHENIX inclusive $J/\psi$ cross section measured at 200 GeV scaled with color-evaporation-model calculations, at the mean $B$ hadron rapidity $y={\pm}1.7$ in 510 GeV $p$$+$$p$ collisions, is $3.63^{+1.92}_{-1.70}\mu$b, and it is consistent with the fixed-order-next-to-leading-logarithm calculations.

The PHENIX experiment at the Relativistic Heavy Ion Collider measured $\pi^0$ and $\eta$ mesons at midrapidity in U$+$U collisions at $\sqrt{s_{_{NN}}}=192$ GeV in a wide transverse momentum range. Measurements were performed in the $\pi^0(\eta)\rightarrow\gamma\gamma$ decay modes. A strong suppression of $\pi^0$ and $\eta$ meson production at high transverse momentum was observed in central U$+$U collisions relative to binary scaled $p$$+$$p$ results. Yields of $\pi^0$ and $\eta$ mesons measured in U$+$U collisions show similar suppression pattern to the ones measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for similar numbers of participant nucleons. The $\eta$/$\pi^0$ ratios do not show dependence on centrality or transverse momentum, and are consistent with previously measured values in hadron-hadron, hadron-nucleus, nucleus-nucleus, and $e^+e^-$ collisions.

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.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured 2nd and 3rd order Fourier coefficients of the azimuthal distributions of direct photons emitted at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for various collision centralities. Combining two different analysis techniques, results were obtained in the transverse momentum range of $0.4<p_{T}<4.0$ GeV/$c$. At low $p_T$ the second-order coefficients, $v_2$, are similar to the ones observed in hadrons. Third order coefficients, $v_3$, are nonzero and almost independent of centrality. These new results on $v_2$ and $v_3$, combined with previously published results on yields, are compared to model calculations that provide yields and asymmetries in the same framework. Those models are challenged to explain simultaneously the observed large yield and large azimuthal anisotropies.

Two particle azimuthal correlation functions are presented for charged hadrons produced in Au + Au collisions at RHIC sqrt(s_NN) = 130 GeV. The measurements permit determination of elliptic flow without event-by-event estimation of the reaction plane. The extracted elliptic flow values v_2 show significant sensitivity to both the collision centrality and the transverse momenta of emitted hadrons, suggesting rapid thermalization and relatively strong velocity fields. When scaled by the eccentricity of the collision zone, epsilon, the scaled elliptic flow shows little or no dependence on centrality for charged hadrons with relatively low p_T. A breakdown of this epsilon scaling is observed for charged hadrons with p_T > 1.0 GeV/c for the most central collisions.

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.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section, mean transverse momentum, mean transverse momentum squared of inclusive $J/\psi$ and cross-section ratio of $\psi(2S)$ to $J/\psi$ at forward rapidity in \pp collisions at \sqrts = 510 GeV via the dimuon decay channel. Comparison is made to inclusive $J/\psi$ cross sections measured at \sqrts = 200 GeV and 2.76--13 TeV. The result is also compared to leading-order nonrelativistic QCD calculations coupled to a color-glass-condensate description of the low-$x$ gluons in the proton at low transverse momentum ($p_T$) and to next-to-leading order nonrelativistic QCD calculations for the rest of the $p_T$ range. These calculations overestimate the data at low $p_T$. While consistent with the data within uncertainties above $\approx3$ GeV/$c$, the calculations are systematically below the data. The total cross section times the branching ratio is BR $d\sigma^{J/\psi}_{pp}/dy (1.2<|y|<2.2, 0<p_T<10~\mbox{GeV/$c$}) =$ 54.3 $\pm$ 0.5 (stat) $\pm$ 5.5 (syst) nb.

Transverse momentum spectra for charged hadrons and for neutral pions in the range 1 GeV/c $< p_T <$ 5 GeV/c have been measured by the PHENIX experiment at RHIC in Au+Au collisions at $\sqrt{s_{_{NN}}}=130$ GeV. At high $p_T$ the spectra from peripheral nuclear collisions are consistent with the naive expectation of scaling the spectra from p+p collisions by the average number of binary nucleon- nucleon collisions. The spectra from central collisions are significantly suppressed when compared to the binary- scaled p+p expectation, and also when compared to similarly binary-scaled peripheral collisions, indicating a novel nuclear effect in central nuclear collisions at RHIC energies.

A three-dimensional (3D) correlation function obtained from mid-rapidity, low pT pion pairs in central Au+Au collisions at sqrt(s_NN)=200 GeV is studied. The extracted model-independent source function indicates a long range tail in the directions of the pion pair transverse momentum (out) and the beam (long). Model comparisons to these distensions indicate a proper breakup time \tau_0 ~ 9 fm/c and a mean proper emission duration \Delta\tau ~ 2 fm/c, leading to sizable emission time differences (<|\Delta \tau_LCM |> ~ 12 fm/c), partly due to resonance decays. They also suggest an outside-in 'burning' of the emission source reminiscent of many hydrodynamical models.

It has been postulated that partonic orbital angular momentum can lead to a significant double-helicity dependence in the net transverse momentum of Drell-Yan dileptons produced in longitudinally polarized p+p collisions. Analogous effects are also expected for dijet production. If confirmed by experiment, this hypothesis, which is based on semi-classical arguments, could lead to a new approach for studying the contributions of orbital angular momentum to the proton spin. We report the first measurement of the double-helicity dependence of the dijet transverse momentum in longitudinally polarized p+p collisions at sqrt(s) = 200 GeV from data taken by the PHENIX experiment in 2005 and 2006. The analysis deduces the transverse momentum of the dijet from the widths of the near- and far-side peaks in the azimuthal correlation of the dihadrons. When averaged over the transverse momentum of the triggered particle, the difference of the root-mean-square of the dijet transverse momentum between like- and unlike-helicity collisions is found to be -37 +/- 88(stat) +/- 14(syst) MeV/c.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K^0_S , \omega, \eta prime, and \phi mesons in p + p collisions at = 200 GeV. Measurements \omega and \phi production in different decay channels give consistent results. New results for the \phi are in agreement with previously published data and extend the measured pT coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-pT and characterizing the low-pT regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

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.

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.

We report the first measurement of the full angular distribution for inclusive $J/\psi\rightarrow\mu^{+}\mu^{-}$ decays in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. The measurements are made for $J/\psi$ transverse momentum $2<p_{T}<10$ GeV/$c$ and rapidity $1.2<y<2.2$ in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient $\lambda_{\theta}$ is strongly negative at low $p_{T}$ and becomes close to zero at high $p_{T}$, while the azimuthal coefficient $\lambda_{\phi}$ is close to zero at low $p_{T}$, and becomes slightly negative at higher $p_{T}$. The frame-independent coefficient $\tilde{\lambda}$ is strongly negative at all $p_{T}$ in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

PHENIX has measured the centrality dependence of mid-rapidity pion, kaon and proton transverse momentum distributions in d+Au and p+p collisions at sqrt(s_NN) = 200 GeV. The p+p data provide a reference for nuclear effects in d+Au and previously measured Au+Au collisions. Hadron production is enhanced in d+Au, relative to independent nucleon-nucleon scattering, as was observed in lower energy collisions. The nuclear modification factor for (anti) protons is larger than that for pions. The difference increases with centrality, but is not sufficient to account for the abundance of baryon production observed in central Au+Au collisions at RHIC. The centrality dependence in d+Au shows that the nuclear modification factor increases gradually with the number of collisions suffered by each participant nucleon. We also present comparisons with lower energy data as well as with parton recombination and other theoretical models of nuclear effects on particle production.

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 the first measurement of transverse single-spin asymmetries in $J/\psi$ production from transversely polarized $p+p$ collisions at $\sqrt{s} = 200$ GeV with data taken by the PHENIX experiment in 2006 and 2008. The measurement was performed over the rapidity ranges $1.2 &lt; |y| &lt; 2.2$ and $ |y| &lt; 0.35$ for transverse momenta up to 6 GeV/$c$. $J/\psi$ production at RHIC is dominated by processes involving initial-state gluons, and transverse single-spin asymmetries of the $J/\psi$ can provide access to gluon dynamics within the nucleon. Such asymmetries may also shed light on the long-standing question in QCD of the $J/\psi$ production mechanism. Asymmetries were obtained as a function of $J/\psi$ transverse momentum and Feynman-$x$, with a value of $-0.086 \pm 0.026^{\rm stat} \pm 0.003^{\rm syst}$ in the forward region. This result suggests possible nonzero trigluon correlation functions in transversely polarized protons and, if well defined in this reaction, a nonzero gluon Sivers distribution function.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) 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 sqrt(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 mid-central 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 neutral pion and the eta meson in the intermediate p_T range (2--5 GeV/c); whereas at higher p_T the phi, pi^0, and eta show similar suppression. The baryon (protons and anti-protons) excess observed in central Au+Au collisions at intermediate p_T is not observed for the phi meson despite the similar mass of the proton and the phi. This suggests that the excess is linked to the number of constituent quarks rather than the hadron mass. The difference gradually disappears with decreasing centrality and for peripheral collisions the R_AA values for both particles 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.

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.

PHENIX has measured the centrality dependence of charged hadron p_T spectra from central Au+Au collisions at sqrt(s_NN)=130 GeV. The truncated mean p_T decreases with centrality for p_T > 2 GeV/c, indicating an apparent reduction of the contribution from hard scattering to high p_T hadron production. For central collisions the yield at high p_T is shown to be suppressed compared to binary nucleon-nucleon collision scaling of p+p data. This suppression is monotonically increasing with centrality, but most of the change occurs below 30% centrality, i.e. for collisions with less than about 140 participating nucleons. The observed p_T and centrality dependence is consistent with the particle production predicted by models including hard scattering and subsequent energy loss of the scattered partons in the dense matter created in the collisions.

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.

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.

We present the first measurements of long-range angular correlations and the transverse momentum dependence of elliptic flow $v_2$ in high-multiplicity $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. A comparison of these results with previous measurements in high-multiplicity $d$$+$Au and $^3{\rm He}$$+$Au collisions demonstrates a relation between $v_2$ and the initial collision eccentricity $\varepsilon_2$, suggesting that the observed momentum-space azimuthal anisotropies in these small systems have a collective origin and reflect the initial geometry. Good agreement is observed between the measured $v_2$ and hydrodynamic calculations for all systems, and an argument disfavoring theoretical explanations based on momentum-space domain correlations is presented. The set of measurements presented here allows us to leverage the distinct intrinsic geometry of each of these systems to distinguish between different theoretical descriptions of the long-range correlations observed in small collision systems.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The correlations are sensitive to nonperturbative initial-state and final-state transverse momentum $k_T$ and $j_T$ in the azimuthal nearly back-to-back region $\Delta\phi\sim\pi$. In this region, transverse-momentum-dependent evolution can be studied when several different hard scales are measured. To have sensitivity to small transverse momentum scales, nonperturbative momentum widths of $p_{\rm out}$, the out-of-plane transverse momentum component perpendicular to the trigger particle, are measured. These widths are used to investigate possible effects from transverse-momentum-dependent factorization breaking. When accounting for the longitudinal momentum fraction of the away-side hadron with respect to the near-side trigger particle, the widths are found to increase with the hard scale; this is qualitatively similar to the observed behavior in Drell-Yan and semi-inclusive deep-inelastic scattering interactions. The momentum widths are also studied as a function of center-of-mass energy by comparing to previous measurements at $\sqrt{s}=510$ GeV. The nonperturbative jet widths also appear to increase with $\sqrt{s}$ at a similar $x_T$, which is qualitatively consistent to similar measurements in Drell-Yan interactions. To quantify the magnitude of any transverse-momentum-dependent factorization breaking effects, calculations will need to be performed to compare to these measurements.

We present direct photon-hadron correlations in 200 GeV/A Au+Au, d+Au, and p+p collisions, for direct photon pT from 5–12 GeV/c, collected by the PHENIX Collaboration in the years from 2006 to 2011. We observe no significant modification of jet fragmentation in d+Au collisions, indicating that cold nuclear matter effects are small or absent. Hadrons carrying a large fraction of the quark's momentum are suppressed in Au+Au compared to p+p and d+Au. As the momentum fraction decreases, the yield of hadrons in Au+Au increases to an excess over the yield in p+p collisions. The excess is at large angles and at low hadron pT and is most pronounced for hadrons associated with lower momentum direct photons. Comparison to theoretical calculations suggests that the hadron excess arises from medium response to energy deposited by jets.

We present measurements of the transverse-momentum dependence of elliptic flow $v_2$ for identified pions and (anti)protons at midrapidity ($|\eta|<0.35$), in 0%--5% central $p$$+$Au and $^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. When taken together with previously published measurements in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, the results cover a broad range of small-collision-system multiplicities and intrinsic initial geometries. We observe a clear mass-dependent splitting of $v_2(p_{T})$ in $d$$+$Au and $^3$He$+$Au collisions, just as in large nucleus-nucleus ($A$$+$$A$) collisions, and a smaller splitting in $p$$+$Au collisions. Both hydrodynamic and transport model calculations successfully describe the data at low $p_T$ ($< 1.5$ GeV/$c$), but fail to describe various features at higher $p_T$. In all systems, the $v_2$ values follow an approximate quark-number scaling as a function of the hadron transverse kinetic energy per constituent quark($KE_T/n_q$), which was also seen previously in $A$$+$$A$ collisions.

We present measurements of azimuthal correlations of charged hadron pairs in $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions after subtracting an underlying event using a model that includes higher-order azimuthal anisotropy $v_2$, $v_3$, and $v_4$. After subtraction, the away-side ($\Delta\phi\sim\pi)$ of the highest transverse-momentum trigger ($p_T>4$ GeV/$c$) correlations is suppressed compared to that of correlations measured in $p$$+$$p$ collisions. At the lowest associated particle $p_T$, the away-side shape and yield are modified. These observations are consistent with the scenario of radiative-jet energy loss. For the lowest-$p_T$ trigger correlations, an away-side yield exists and we explore the dependence of the shape of the away-side within the context of an underlying-event model. Correlations are also studied differentially versus event-plane angle $\Psi_n$. The angular correlations show an asymmetry when selecting the sign of the trigger-particle azimuthal angle with respect to the $\Psi_2$ event plane. This asymmetry and the measured suppression of the pair yield out of plane is consistent with a path-length-dependent energy loss. No $\Psi_3$ dependence can be resolved within experimental uncertainties.

For Au + Au collisions at 200 GeV we measure neutral pion production with good statistics for transverse momentum, p_T, up to 20 GeV/c. A fivefold suppression is found, which is essentially constant for 5 < p_T < 20 GeV/c. Experimental uncertainties are small enough to constrain any model-dependent parameterization for the transport coefficient of the medium, e.g. \mean(q^hat) in the parton quenching model. The spectral shape is similar for all collision classes, and the suppression does not saturate in Au+Au collisions/ instead, it increases proportional to the number of participating nucleons, as N_part^2/3.

Measurements of the anisotropy parameter v_2 of identified hadrons (pions, kaons, and protons) as a function of centrality, transverse momentum p_T, and transverse kinetic energy KE_T at midrapidity (|\eta|<0.35) in Au+Au collisions at sqrt(s_NN) = 200 GeV are presented. Pions and protons are identified up to p_T = 6 GeV/c, and kaons up to p_T = 4 GeV/c, by combining information from time-of-flight and aerogel Cherenkov detectors in the PHENIX Experiment. The scaling of v_2 with the number of valence quarks (n_q) has been studied in different centrality bins as a function of transverse momentum and transverse kinetic energy. A deviation from previously observed quark-number scaling is observed at large values of KE_T/n_q in noncentral Au+Au collisions (20--60%), but this scaling remains valid in central collisions (0--10%).

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.

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 transverse single-spin asymmetries of neutral pions and non-identified charged hadrons have been measured at mid-rapidity in polarized proton-proton collisions at sqrt(s) = 200 GeV. The data cover a transverse momentum (p_T) range 0.5-5.0 GeV/c for charged hadrons and 1.0-5.0 GeV/c for neutral pions, at a Feynman-x (x_F) value of approximately zero. The asymmetries seen in this previously unexplored kinematic region are consistent with zero within statistical errors of a few percent. In addition, the inclusive charged hadron cross section at mid-rapidity from 0.5 < p_T < 7.0 GeV/c is presented and compared to NLO pQCD calculations. Successful description of the unpolarized cross section above ~2 GeV/c using NLO pQCD suggests that pQCD is applicable in the interpretation of the asymmetry results in the relevant kinematic range.

We present measurements of two-particle angular correlations between high-transverse-momentum ($2<p_T<11$ GeV/$c$) $\pi^0$ observed at midrapidity ($|\eta|<0.35$) and particles produced either at forward ($3.1<\eta<3.9$) or backward ($-3.7<\eta<-3.1$) rapidity in $d$$+$Au and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The azimuthal angle correlations for particle pairs with this large rapidity gap in the Au-going direction exhibit a ridge-like structure that persists up to $p_T{\approx}6$ GeV/$c$ and which strongly depends on collision centrality, which is a similar characteristic to the hydrodynamical particle flow in A+A collisions. The ridge-like structure is absent in the $d$-going direction as well as in $p$$+$$p$ collisions, in the transverse-momentum range studied. The results indicate that the ridge-like structure is shifted in the Au-going direction toward more central collisions, similar to the charged-particle pseudorapidity distributions.

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 report on the nuclear dependence of transverse single-spin asymmetries (TSSAs) in the production of positively-charged hadrons in polarized $p^{\uparrow}+p$, $p^{\uparrow}+$Al and $p^{\uparrow}+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The measurements have been performed at forward rapidity ($1.4<\eta<2.4$) over the range of $1.8<p_{T}<7.0$ GeV$/c$ and $0.1<x_{F}<0.2$. We observed a positive asymmetry $A_{N}$ for positively-charged hadrons in \polpp collisions, and a significantly reduced asymmetry in $p^{\uparrow}$+$A$ collisions. These results reveal a nuclear dependence of charged hadron $A_N$ in a regime where perturbative techniques are relevant. These results provide new opportunities to use \polpA collisions as a tool to investigate the rich phenomena behind TSSAs in hadronic collisions and to use TSSA as a new handle in studying small-system collisions.