Azimuthal angle (Delta phi) correlations are presented for a broad range of transverse momentum (0.4 < pT < 10 GeV/c) and centrality (0-92%) selections for charged hadrons from di-jets in Au+Au collisions at sqrt(s_NN) = 200 GeV. With increasing pT, the away-side Delta phi distribution evolves from a broad and relatively flat shape to a concave shape, then to a convex shape. Comparisons to p+p data suggest that the away-side distribution can be divided into a partially suppressed head region centered at Delta phi ~ \pi, and an enhanced shoulder region centered at Delta phi ~ \pi \pm 1:1. The pT spectrum for the associated hadrons in the head region softens toward central collisions. The spectral slope for the shoulder region is independent of centrality and trigger pT . The properties of the near-side distributions are also modified relative to those in p + p collisions, reflected by the broadening of the jet shape in Delta phi and Delta eta, and an enhancement of the per-trigger yield. However, these modifications seem to be limited to pT < 4 GeV/c, above which both the dihadron pair shape and per-trigger yield become similar to p + p collisions. These observations suggest that both the away- and near-side distributions contain a jet fragmentation component which dominates for pT \ge 5GeV and a medium-induced component which is important for pT \le 4 GeV/c. We also quantify the role of jets at intermediate and low pT through the yield of jet-induced pairs in comparison to binary scaled p + p pair yield. The yield of jet-induced pairs is suppressed at high pair proxy energy (sum of the pT magnitudes of the two hadrons) and is enhanced at low pair proxy energy. The former is consistent with jet quenching/ the latter is consistent with the enhancement of soft hadron pairs due to transport of lost energy to lower pT.

High-energy proton- and deuteron-nucleus collisions provide an excellent tool for studying a wide array of physics effects, including modifications of parton distribution functions in nuclei, gluon saturation, and color neutralization and hadronization in a nuclear environment, among others. All of these effects are expected to have a significant dependence on the size of the nuclear target and the impact parameter of the collision, also known as the collision centrality. In this article, we detail a method for determining centrality classes in p(d)+A collisions via cuts on the multiplicity at backward rapidity (i.e., the nucleus-going direction) and for determining systematic uncertainties in this procedure. For d+Au collisions at sqrt(s_NN) = 200 GeV we find that the connection to geometry is confirmed by measuring the fraction of events in which a neutron from the deuteron does not interact with the nucleus. As an application, we consider the nuclear modification factors R_{p(d)+A}, for which there is a potential bias in the measured centrality dependent yields due to auto-correlations between the process of interest and the backward rapidity multiplicity. We determine the bias correction factor within this framework. This method is further tested using the HIJING Monte Carlo generator. We find that for d+Au collisions at sqrt(s_NN)=200 GeV, these bias corrections are small and vary by less than 5% (10%) up to p_T = 10 (20) GeV. In contrast, for p+Pb collisions at sqrt(s_NN) = 5.02 TeV we find these bias factors are an order of magnitude larger and strongly p_T dependent, likely due to the larger effect of multi-parton interactions.

Yields for J/psi production in Cu+Cu collisions at sqrt (s_NN)= 200 GeV have been measured by the PHENIX experiment over the rapidity range |y| < 2.2 at transverse momenta from 0 to beyond 5 GeV/c. The invariant yield is obtained as a function of rapidity, transverse momentum and collision centrality, and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data provide greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, providing a key constraint that is needed for disentangling cold and hot nuclear matter effects.

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.

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.

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.

The second Fourier component v_2 of the azimuthal anisotropy with respect to the reaction plane was measured for direct photons at midrapidity and transverse momentum (p_T) of 1--13 GeV/c in Au+Au collisions at sqr(s_NN)=200 GeV. Previous measurements of this quantity for hadrons with p_T < 6 GeV/c indicate that the medium behaves like a nearly perfect fluid, while for p_T > 6 GeV/c a reduced anisotropy is interpreted in terms of a path-length dependence for parton energy loss. In this measurement with the PHENIX detector at the Relativistic Heavy Ion Collider we find that for p_T > 4 GeV/c the anisotropy for direct photons is consistent with zero, as expected if the dominant source of direct photons is initial hard scattering. However, in the p_T < 4 GeV/c region dominated by thermal photons, we find a substantial direct photon v_2 comparable to that of hadrons, whereas model calculations for thermal photons in this kinematic region significantly underpredict the observed v_2.

The invariant yields for $J/\psi$ production at forward rapidity $(1.2<|y|<2.2)$ in U$+$U collisions at $\sqrt{s_{_{NN}}}$=193 GeV have been measured as a function of collision centrality. The invariant yields and nuclear-modification factor $R_{AA}$ are presented and compared with those from Au$+$Au collisions in the same rapidity range. Additionally, the direct ratio of the invariant yields from U$+$U and Au$+$Au collisions within the same centrality class is presented, and used to investigate the role of $c\bar{c}$ coalescence. Two different parameterizations of the deformed Woods-Saxon distribution were used in Glauber calculations to determine the values of the number of nucleon-nucleon collisions in each centrality class, $N_{\rm coll}$, and these were found to give significantly different $N_{\rm coll}$ values. Results using $N_{\rm coll}$ values from both deformed Woods-Saxon distributions are presented. The measured ratios show that the $J/\psi$ suppression, relative to binary collision scaling, is similar in U$+$U and Au$+$Au for peripheral and midcentral collisions, but that $J/\psi$ show less suppression for the most central U$+$U collisions. The results are consistent with a picture in which, for central collisions, increase in the $J/\psi$ yield due to $c\bar{c}$ coalescence becomes more important than the decrease in yield due to increased energy density. For midcentral collisions, the conclusions about the balance between $c\bar{c}$ coalescence and suppression depend on which deformed Woods-Saxon distribution is used to determine $N_{\rm coll}$.

The measurement of the direct-photon spectrum from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV is presented by the PHENIX collaboration using the external-photon-conversion technique for 0%--93% central collisions in a transverse-momentum ($p_T$) range of 0.8--10 GeV/$c$. An excess of direct photons, above prompt-photon production from hard-scattering processes, is observed for $p_T<6$ GeV/$c$. Nonprompt direct photons are measured by subtracting the prompt component, which is estimated as $N_{\rm coll}$-scaled direct photons from $p$ $+$ $p$ collisions at 200 GeV, from the direct-photon spectrum. Results are obtained for $0.8<p_T<6.0$ GeV/$c$ and suggest that the spectrum has an increasing inverse slope from ${\approx}0.2$ to 0.4 GeV/$c$ with increasing $p_T$, which indicates a possible sensitivity of the measurement to photons from earlier stages of the evolution of the collision. In addition, like the direct-photon production, the $p_T$-integrated nonprompt direct-photon yields also follow a power-law scaling behavior as a function of collision-system size. The exponent, $\alpha$, for the nonprompt component is found to be consistent with 1.1 with no apparent $p_T$ dependence.

The differential cross section for the production of direct photons in p+p collisions at sqrt(s)=200 GeV at midrapidity was measured in the PHENIX detector at the Relativistic Heavy Ion Collider. Inclusive-direct photons were measured in the transverse-momentum range from 5.5--25 GeV/c, extending the range beyond previous measurements. Event structure was studied with an isolation criterion. Next-to-leading-order perturbative-quantum-chromodynamics calculations give a good description of the spectrum. When the cross section is expressed versus x_T, the PHENIX data are seen to be in agreement with measurements from other experiments at different center-of-mass energies.

Heavy quarkonia are observed to be suppressed in relativistic heavy ion collisions relative to their production in p+p collisions scaled by the number of binary collisions. In order to determine if this suppression is related to color screening of these states in the produced medium, one needs to account for other nuclear modifications including those in cold nuclear matter. In this paper, we present new measurements from the PHENIX 2007 data set of J/psi yields at forward rapidity (1.2<|y|<2.2) in Au+Au collisions at sqrt(s_NN)=200 GeV. The data confirm the earlier finding that the suppression of J/psi at forward rapidity is stronger than at midrapidity, while also extending the measurement to finer bins in collision centrality and higher transverse momentum (pT). We compare the experimental data to the most recent theoretical calculations that incorporate a variety of physics mechanisms including gluon saturation, gluon shadowing, initial-state parton energy loss, cold nuclear matter breakup, color screening, and charm recombination. We find J/psi suppression beyond cold-nuclear-matter effects. However, the current level of disagreement between models and d+Au data precludes using these models to quantify the hot-nuclear-matter suppression.

Flow coefficients v_n for n = 2, 3, 4, characterizing the anisotropic collective flow in Au+Au collisions at sqrt(s_NN) = 200 GeV, are measured relative to event planes Ψ_n determined at large rapidity. We report v_n as a function of transverse momentum and collision centrality, and study the correlations among the event planes of different order n. The v_n are well described by hydrodynamic models which employ a Glauber Monte Carlo initial state geometry with fluctuations, providing additional constraining power on the interplay between initial conditions and the effects of viscosity as the system evolves. This new constraint improves precision of the extracted viscosity to entropy density ratio eta/s.

The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) reports measurements of azimuthal dihadron correlations near midrapidity in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. These measurements complement recent analyses by experiments at the Large Hadron Collider (LHC) involving central $p$$+$Pb collisions at $\sqrt{s_{_{NN}}}$=5.02 TeV, which have indicated strong anisotropic long-range correlations in angular distributions of hadron pairs. The origin of these anisotropies is currently unknown. Various competing explanations include parton saturation and hydrodynamic flow. We observe qualitatively similar, but larger, anisotropies in $d$$+$Au collisions compared to those seen in $p$$+$Pb collisions at the LHC. The larger extracted $v_2$ values in $d$$+$Au collisions at RHIC are consistent with expectations from hydrodynamic calculations owing to the larger expected initial-state eccentricity compared with that from $p$$+$Pb collisions. When both are divided by an estimate of the initial-state eccentricity the scaled anisotropies follow a common trend with multiplicity that may extend to heavy ion data at RHIC and the LHC, where the anisotropies are widely thought to arise from hydrodynamic flow.

We have measured direct photons for $p_T<5~$GeV/$c$ in minimum bias and 0\%--40\% most central events at midrapidity for Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The $e^{+}e^{-}$ contribution from quasi-real direct virtual photons has been determined as an excess over the known hadronic contributions in the $e^{+}e^{-}$ mass distribution. A clear enhancement of photons over the binary scaled $p$$+$$p$ fit is observed for $p_T<4$ GeV/$c$ in Cu$+$Cu data. The $p_T$ spectra are consistent with the Au$+$Au data covering a similar number of participants. The inverse slopes of the exponential fits to the excess after subtraction of the $p$$+$$p$ baseline are 285$\pm$53(stat)$\pm$57(syst)~MeV/$c$ and 333$\pm$72(stat)$\pm$45(syst)~MeV/$c$ for minimum bias and 0\%--40\% most central events, respectively. The rapidity density, $dN/dy$, of photons demonstrates the same power law as a function of $dN_{\rm ch}/d\eta$ observed in Au$+$Au at the same collision energy.

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.

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

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.

We report on charmonium measurements [J/psi(1S), psi'(2S), and chi_c(1P)] in p+p collisions at sqrt(s)=200 GeV. We find that the fraction of J/psi coming from the feed-down decay of psi' and chi_c in the midrapidity region ($|\eta|<0.35$) is 9.6+/-2.4% and 32+/-9%, respectively. We also report new, higher statistics p_T and rapidity dependencies of the J/psi yield via dielectron decay in the same midrapidity range and at forward rapidity (1.2<|eta|<2.4) via dimuon decay. These results are compared with measurements from other experiments and discussed in the context of current charmonium production models.

Transverse momentum distributions and yields for $\pi^{\pm}$, $K^{\pm}$, $p$ and $\bar{p}$ in $p+p$ collisions at $\sqrt{s}$=200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter $T_{\rm inv}$, mean transverse momentum $<p_T>$ and yield per unit rapidity $dN/dy$ at each energy, and compare them to other measurements at different $\sqrt{s}$ in $p+p$ and $p+\bar{p}$ collisions. We also present the scaling properties such as $m_T$ scaling, $x_T$ scaling on the $p_T$ spectra between different energies. To discuss the mechanism of the particle production in $p+p$ collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.

We present measurements of electrons and positrons from the semileptonic decays of heavy-flavor hadrons at midrapidity ($|y|<$ 0.35) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. The data were collected in 2010 by the PHENIX experiment that included the new hadron-blind detector. The invariant yield of electrons from heavy-flavor decays is measured as a function of transverse momentum in the range $1<p_T^e<5$ GeV/$c$. The invariant yield per binary collision is slightly enhanced above the $p$$+$$p$ reference in Au$+$Au 0%--20%, 20%--40% and 40%--60% centralities at a comparable level. This may be a result of the interplay between initial-state Cronin effects, final-state flow, and energy loss for heavy-quark production at this low beam energy. The $v_2$ of electrons from heavy-flavor decays is nonzero when averaged between $1.3<p_T^e<2.5$ GeV/$c$ from $0<{\rm centrality}<40$% collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. For 20%--40% centrality collisions, the $v_2$ at $\sqrt{s_{_{NN}}}=62.4$ GeV is smaller than that for heavy flavor decays at $\sqrt{s_{_{NN}}}=200$ GeV. The $v_2$ of the electrons from heavy-flavor decay at the lower beam energy is also smaller than $v_2$ for pions. Both results indicate that the heavy-quarks interact with the medium formed in these collisions, but they may not be at the same level of thermalization with the medium as observed at $\sqrt{s_{_{NN}}}=200$ GeV.

The PHENIX experiment has measured electrons and positrons at midrapidity from the decays of hadrons containing charm and bottom quarks produced in d+Au and p+p collisions at sqrt(s_NN)=200 GeV at the Relativistic Heavy Ion Collider, in the transverse-momentum range 0.85 < pT < 8.5 GeV/c. In central d+Au collisions, the nuclear modification factor R_dA at 1.5 < pT < 5 GeV/c displays evidence of enhancement of these electrons, relative to those produced in p+p collisions, and shows that the mass-dependent Cronin enhancement observed at RHIC extends to the heavy-D-meson family. A comparison with the neutral-pion data suggests that the difference in cold-nuclear-matter effects on light- and heavy-flavor mesons could contribute to the observed differences between the pi0 and heavy-flavor-electron nuclear modification factor R_AA.

The azimuthal anisotropy coefficients v_2 and v_4 of pi^0 and eta mesons are measured in Au+Au collisions at sqrt(s_NN)=200 GeV, as a function of transverse momentum p_T (1-14 GeV/c) and centrality. The extracted v_2 coefficients are found to be consistent between the two meson species over the measured p_T range. The ratio of v_4/v_2^2 for pi^0 mesons is found to be independent of p_T for 1-9 GeV/c, implying a lack of sensitivity of the ratio to the change of underlying physics with p_T. Furthermore, the ratio of v_4/v_2^2 is systematically larger in central collisions, which may reflect the combined effects of fluctuations in the initial collision geometry and finite viscosity in the evolving medium.

Neutral-pion, pi^0, spectra were measured at midrapidity (|y|<0.35) in Au+Au collisions at sqrt(s_NN) = 39 and 62.4 GeV and compared to earlier measurements at 200 GeV in the 1<p_T<10 GeV/c transverse-momentum (p_T) range. The high-p_T tail is well described by a power law in all cases and the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very similar to that observed in the corresponding p+p-collision spectra. The nuclear-modification factors (R_AA) show significant suppression and a distinct energy dependence at moderate p_T in central collisions. At high p_T, R_AA is similar for 62.4 and 200 GeV at all centralities. Perturbative-quantum-chromodynamics calculations that describe R_AA well at 200 GeV, fail to describe the 39 GeV data, raising the possibility that the relative importance of initial-state effects and soft processes increases at lower energies. A conclusion that the region where hard processes are dominant is reached only at higher p_T, is also supported by the x_T dependence of the x_T-scaling power-law exponent.

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.

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.