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.

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.

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

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.

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.

Asymmetric nuclear collisions of $p$$+$Al, $p$$+$Au, $d$$+$Au, and $^{3}$He$+$Au at $\sqrt{s_{_{NN}}}=200$ GeV provide an excellent laboratory for understanding particle production, as well as exploring interactions among these particles after their initial creation in the collision. We present measurements of charged hadron production $dN_{\rm ch}/d\eta$ in all such collision systems over a broad pseudorapidity range and as a function of collision multiplicity. A simple wounded quark model is remarkably successful at describing the full data set. We also measure the elliptic flow $v_{2}$ over a similarly broad pseudorapidity range. These measurements provide key constraints on models of particle emission and their translation into flow.

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.

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.

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.

We report on event structure and double helicity asymmetry ($A_LL$) of jet production in longitudinally polarized p+p collisions at $\sqrt{s}$=200 GeV. Photons and charged particles were measured at midrapidity $|\eta| < 0.35$ with the requirement of a high-momentum ($>2$ GeV/$c$) photon in each event. Measured event structure is compared with {\sc pythia} and {\sc geant} simulations. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet $A_{LL}$, photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster $p_T$ sum ($p_T^{\rm reco}$). The effect of detector response and the underlying events on $p_T^{\rm reco}$ was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the NLO pQCD jet production cross section. For $4 < p_T^{\rm reco} < 12$ GeV/$c$ with an average beam polarization of $< P > = 49%$ we measured $A_{LL} = -0.0014 \pm 0.0037^{\rm stat}$ at the lowest $p_T^{\rm reco}$ bin (4-5 GeV/$c$) and $-0.0181 \pm 0.0282^{\rm stat}$ at the highest $p_T^{\rm reco}$ bin (10-12 GeV/$c$) with a beam polarization scale error of 9.4% and a $\pT$ scale error of 10%. Jets in the measured $p_T^{\rm reco}$ range arise primarily from hard-scattered gluons with momentum fraction $0.02 < x < 0.3$ according to {\sc pythia}. The measured $A_{LL}$ is compared with predictions that assume various $\Delta G(x)$ distributions based on the GRSV parameterization. The present result imposes the limit $-1.1 < \int_{0.02}^{0.3}dx \Delta G(x, \mu^2 = 1 {\rm GeV}^2) < 0.4$ at 95% confidence level or $\int_{0.02}^{0.3}dx \Delta G(x, \mu^2 = 1 {\rm GeV}^2) < 0.5$ at 99% confidence level.