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 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 angular correlations between charged hadrons and energy deposited in calorimeter towers in central $d$$+$Au and minimum bias $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The charged hadron is measured at midrapidity $|\eta|<0.35$, and the energy is measured at large rapidity ($-3.7<\eta<-3.1$, Au-going direction). An enhanced near-side angular correlation across $|\Delta\eta| >$ 2.75 is observed in $d$$+$Au collisions. Using the event plane method applied to the Au-going energy distribution, we extract the anisotropy strength $v_2$ for inclusive charged hadrons at midrapidity up to $p_T=4.5$ GeV/$c$. We also present the measurement of $v_2$ for identified $\pi^{\pm}$ and (anti)protons in central $d$$+$Au collisions, and observe a mass-ordering pattern similar to that seen in heavy ion collisions. These results are compared with viscous hydrodynamic calculations and measurements from $p$$+$Pb at $\sqrt{s_{_{NN}}}=5.02$ TeV. The magnitude of the mass-ordering in $d$$+$Au is found to be smaller than that in $p$$+$Pb collisions, which may indicate smaller radial flow in lower energy $d$$+$Au collisions.

The cross section and transverse single-spin asymmetries of $\mu^{-}$ and $\mu^{+}$ from open heavy-flavor decays in polarized $p$+$p$ collisions at $\sqrt{s}=200$ GeV were measured by the PHENIX experiment during 2012 at the Relativistic Heavy Ion Collider. Because heavy-flavor production is dominated by gluon-gluon interactions at $\sqrt{s}=200$ GeV, these measurements offer a unique opportunity to obtain information on the trigluon correlation functions. The measurements are performed at forward and backward rapidity ($1.4<|y|<2.0$) over the transverse momentum range of $1.25<p_T<7$ GeV/$c$ for the cross section and $1.25<p_T<5$ GeV/$c$ for the asymmetry measurements. The obtained cross section is compared to a fixed-order-plus-next-to-leading-log perturbative-quantum-chromodynamics calculation. The asymmetry results are consistent with zero within uncertainties, and a model calculation based on twist-3 three-gluon correlations agrees with the data.

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

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.

The fraction of $J/\psi$ mesons which come from B-meson decay, $\textrm{F}_{B{\rightarrow}J/\psi}$, is measured for J/$\psi$ rapidity \mbox{$1.2<|y|<2.2$} and $p_T>0$ in $p$+$p$ and Cu+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV with the PHENIX detector. The extracted fraction is $\textrm{F}_{B{\rightarrow}J/\psi}$ = 0.025 $\pm$ 0.006(stat) $\pm$ 0.010(syst) for $p$+$p$ collisions. For Cu+Au collisions, $\textrm{F}_{B{\rightarrow}J/\psi}$ is 0.094 $\pm$ 0.028(stat) $\pm$ 0.037(syst) in the Au-going direction ($-2.2<y<-1.2$) and 0.089 $\pm$ 0.026(stat) $\pm$ 0.040(syst) in the Cu-going direction ($1.2<y<2.2$). The nuclear modification factor, $R_{\rm CuAu}$, of B mesons in Cu+Au collisions is consistent with binary scaling of measured yields in $p$+$p$ at both forward and backward rapidity.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured omega meson production via leptonic and hadronic decay channels in p+p, d+Au, Cu+Cu, and Au+Au collisions at sqrt(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 omega 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 (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.