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.

Two particle correlations between identified meson and baryon trigger particles with 2.5 < p_T < 4.0 GeV/c and lower p_T charged hadrons have been measured at midrapidity by the PHENIX experiment at RHIC in p+p, d+Au and Au+Au collisions at sqrt(s_NN) = 200 GeV. The probability of finding a hadron near in azimuthal angle to the trigger particle is almost identical for leading mesons and baryons for non-central Au+Au. The yield for both trigger baryons and mesons is significantly higher in Au+Au than in p+p and d+Au, except for trigger baryons in central collisions. The baryon excess is likely to arise predominantly from hard scattering processes.

We report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $p$+$p$ collisions at $\sqrt{s_{NN}}=200$~GeV. The $e^+e^-$ pair yield from $b\bar{b}$ and $c\bar{c}$ is separated by exploiting a double differential fit done simultaneously in dielectron invariant mass and $p_T$. We used three different event generators, {\sc pythia}, {\sc mc@nlo}, and {\sc powheg}, to simulate the $e^+e^-$ spectra from $c\bar{c}$ and $b\bar{b}$ production. The data can be well described by all three generators within the detector acceptance. However, when using the generators to extrapolate to $4\pi$, significant differences are observed for the total cross section. These difference are less pronounced for $b\bar{b}$ than for $c\bar{c}$. The same model dependence was observed in already published $d$+$A$ data. The $p$+$p$ data are also directly compared with $d$+$A$ data in mass and $p_T$, and within the statistical accuracy no nuclear modification is seen.

We report $e^\pm-\mu^\mp$ pair yield from charm decay measured between midrapidity electrons ($|\eta|<0.35$ and $p_T>0.5$ GeV/$c$) and forward rapidity muons ($1.4<\eta<2.1$ and $p_T>1.0$ GeV/$c$) as a function of $\Delta\phi$ in both $p$$+$$p$ and in $d$+Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Comparing the $p$$+$$p$ results with several different models, we find the results are consistent with a total charm cross section $\sigma_{c\bar{c}} =$ 538 $\pm$ 46 (stat) $\pm$ 197 (data syst) $\pm$ 174 (model syst) $\mu$b. These generators also indicate that the back-to-back peak at $\Delta\phi = \pi$ is dominantly from the leading order contributions (gluon fusion), while higher order processes (flavor excitation and gluon splitting) contribute to the yield at all $\Delta\phi$. We observe a suppression in the pair yield per collision in $d$+Au. We find the pair yield suppression factor for $2.7<\Delta\phi<3.2$ rad is $J_{dA}$ = 0.433 $\pm$ 0.087 (stat) $\pm$ 0.135 (syst), indicating cold nuclear matter modification of $c\bar{c}$ pairs.

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.

We present results for three charmonia states (psi^prime, chi_c and J/psi) in d+Au collisions at |y|<0.35 and sqrt(s_NN)=200 GeV. We find that the modification of the psi^prime yield relative to that of the J/psi scales approximately with charged-particle multiplicity at midrapidity across p+A, d+Au, and A+A results from the Super Proton Synchrotron and the Relativistic Heavy Ion Collider. In large impact-parameter collisions we observe a similar suppression for the psi^prime and J/psi, while in small impact-parameter collisions the more weakly bound psi^prime is more strongly suppressed. Owing to the short time spent traversing the Au nucleus, the larger psi^prime suppression in central events is not explained by an increase of the nuclear absorption due to meson formation time effects.

Back-to-back hadron pair yields in d+Au and p+p collisions at sqrt(s_NN)=200 GeV were measured with the PHENIX detector at the Relativistic Heavy Ion Collider. Rapidity separated hadron pairs were detected with the trigger hadron at pseudorapidity |eta|<0.35 and the associated hadron at forward rapidity (deuteron direction, 3.0<eta<3.8). Pairs were also detected with both hadrons measured at forward rapidity; in this case the yield of back-to-back hadron pairs in d+Au collisions with small impact parameters is observed to be suppressed by a factor of 10 relative to p+p collisions. The kinematics of these pairs is expected to probe partons in the Au nucleus with low fraction x of the nucleon momenta, where the gluon densities rise sharply. The observed suppression as a function of nuclear thickness, p_T, and eta points to cold nuclear matter effects arising at high parton densities.

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.

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.

The three Upsilon states, Upsilon(1S+2S+3S), are measured in d+Au and p+p collisions at sqrt(s_NN)=200 GeV and rapidities 1.2<|y|<2.2 by the PHENIX experiment at the Relativistic Heavy-Ion Collider. Cross sections for the inclusive Upsilon(1S+2S+3S) production are obtained. The inclusive yields per binary collision for d+Au collisions relative to those in p+p collisions (R_dAu) are found to be 0.62 +/- 0.26 (stat) +/- 0.13 (syst) in the gold-going direction and 0.91 +/- 0.33 (stat) +/- 0.16 (syst) in the deuteron-going direction. The measured results are compared to a nuclear-shadowing model, EPS09 [JHEP 04, 065 (2009)], combined with a final-state breakup cross section, sigma_br, and compared to lower energy p+A results. We also compare the results to the PHENIX J/psi results [Phys. Rev. Lett. 107, 142301 (2011)]. The rapidity dependence of the observed Upsilon suppression is consistent with lower energy p+A measurements.