The PHENIX experiment at the Relativistic Heavy Ion Collider has measured low mass vector meson, $\omega$, $\rho$, and $\phi$, production through the dimuon decay channel at forward rapidity ($1.2<|y|<2.2$) in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The differential cross sections for these mesons are measured as a function of both $p_T$ and rapidity. We also report the integrated differential cross sections over $1<p_T<7$ GeV/$c$ and $1.2<|y|<2.2$: $d\sigma/dy(\omega+\rho\rightarrow\mu\mu) = 80 \pm 6 \mbox{(stat)} \pm 12 \mbox{(syst)}$ nb and $d\sigma/dy(\phi\rightarrow\mu\mu) = 27 \pm 3 \mbox{(stat)} \pm 4 \mbox{(syst)}$ nb. These results are compared with midrapidity measurements and calculations.

The PHENIX experiment at RHIC has measured the centrality dependence of the direct photon yield from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV down to $p_T=0.4$ GeV/$c$. Photons are detected via photon conversions to $e^+e^-$ pairs and an improved technique is applied that minimizes the systematic uncertainties that usually limit direct photon measurements, in particular at low $p_T$. We find an excess of direct photons above the $N_{\rm coll}$-scaled yield measured in $p$$+$$p$ collisions. This excess yield is well described by an exponential distribution with an inverse slope of about 240 MeV/$c$ in the $p_T$ range from 0.6--2.0 GeV/$c$. While the shape of the $p_T$ distribution is independent of centrality within the experimental uncertainties, the yield increases rapidly with increasing centrality, scaling approximately with $N_{\rm part}^\alpha$, where $\alpha=1.48{\pm}0.08({\rm stat}){\pm}0.04({\rm syst})$.

We report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Exploring the mass and transverse-momentum dependence of the yield, the bottom decay contribution can be isolated from charm, and quantified by comparison to {\sc pythia} and {\sc mc@nlo} simulations. The resulting $b\bar{b}$-production cross section is $\sigma^{d{\rm Au}}_{b\bar{b}}=1.37{\pm}0.28({\rm stat}){\pm}0.46({\rm syst})$~mb, which is equivalent to a nucleon-nucleon cross section of $\sigma^{NN}_{bb}=3.4\pm0.8({\rm stat}){\pm}1.1({\rm syst})\ \mu$b.

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.

Measurements of bottomonium production in heavy ion and $p$$+$$p$ collisions at the Relativistic Heavy Ion Collider (RHIC) are presented. The inclusive yield of the three $\Upsilon$ states, $\Upsilon(1S+2S+3S)$, was measured in the PHENIX experiment via electron-positron decay pairs at midrapidity for Au$+$Au and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The $\Upsilon(1S+2S+3S)\rightarrow e^+e^-$ differential cross section at midrapidity was found to be $B_{\rm ee} d\sigma/dy =$ 108 $\pm$ 38 (stat) $\pm$ 15(syst) $\pm$ 11 (luminosity) pb in $p$$+$$p$ collisions. The nuclear modification factor in the 30\% most central Au$+$Au collisions indicates a suppression of the total $\Upsilon$ state yield relative to the extrapolation from $p$$+$$p$ collision data. The suppression is consistent with measurements made by STAR at RHIC and at higher energies by the CMS experiment at the Large Hadron Collider.

Results are presented from data recorded in 2009 by the PHENIX experiment at the Relativistic Heavy Ion Collider for the double-longitudinal spin asymmetry, $A_{LL}$, for $\pi^0$ and $\eta$ production in $\sqrt{s} = 200$ GeV polarized $p$$+$$p$ collisions. Comparison of the $\pi^0$ results with different theory expectations based on fits of other published data showed a preference for small positive values of gluon polarization, $\Delta G$, in the proton in the probed Bjorken $x$ range. The effect of adding the new 2009 \pz data to a recent global analysis of polarized scattering data is also shown, resulting in a best fit value $\Delta G^{[0.05,0.2]}_{\mbox{DSSV}} = 0.06^{+0.11}_{-0.15}$ in the range $0.05<x<0.2$, with the uncertainty at $\Delta \chi^2 = 9$ when considering only statistical experimental uncertainties. Shifting the PHENIX data points by their systematic uncertainty leads to a variation of the best-fit value of $\Delta G^{[0.05,0.2]}_{\mbox{DSSV}}$ between $0.02$ and $0.12$, demonstrating the need for full treatment of the experimental systematic uncertainties in future global analyses.

Charged-pion-interferometry measurements were made with respect to the 2$^{\rm nd}$- and 3$^{\rm rd}$-order event plane for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. A strong azimuthal-angle dependence of the extracted Gaussian-source radii was observed with respect to both the 2$^{\rm nd}$- and 3$^{\rm rd}$-order event planes. The results for the 2$^{\rm nd}$-order dependence indicate that the initial eccentricity is reduced during the medium evolution, but not reversed in the final state, which is consistent with previous results. In contrast, the results for the 3$^{\rm rd}$-order dependence indicate that the initial triangular shape is significantly reduced and potentially reversed by the end of the medium evolution, and that the 3$^{\rm rd}$-order oscillations are largely dominated by the dynamical effects from triangular flow.

Measurements of transverse-single-spin asymmetries ($A_{N}$) in $p$$+$$p$ collisions at $\sqrt{s}=$62.4 and 200 GeV with the PHENIX detector at RHIC are presented. At midrapidity, $A_{N}$ is measured for neutral pion and eta mesons reconstructed from diphoton decay, and at forward rapidities, neutral pions are measured using both diphotons and electromagnetic clusters. The neutral-pion measurement of $A_{N}$ at midrapidity is consistent with zero with uncertainties a factor of 20 smaller than previous publications, which will lead to improved constraints on the gluon Sivers function. At higher rapidities, where the valence quark distributions are probed, the data exhibit sizable asymmetries. In comparison with previous measurements in this kinematic region, the new data extend the kinematic coverage in $\sqrt{s}$ and $p_T$, and it is found that the asymmetries depend only weakly on $\sqrt{s}$. The origin of the forward $A_{N}$ is presently not understood quantitatively. The extended reach to higher $p_T$ probes the transition between transverse momentum dependent effects at low $p_T$ and multi-parton dynamics at high $p_T$.

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.

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.