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.

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 PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured $\phi$ meson production and its nuclear modification in asymmetric Cu$+$Au heavy-ion collisions at $\sqrt{s_{NN}}=200$ GeV at both forward Cu-going direction ($1.2<y<2.2$) and backward Au-going direction ($-2.2<y<-1.2$), rapidities. The measurements are performed via the dimuon decay channel and reported as a function of the number of participating nucleons, rapidity, and transverse momentum. In the most central events, 0\%--20\% centrality, the $\phi$ meson yield integrated over $1<p_T<5$ GeV/$c$ prefers a smaller value, which means a larger nuclear modification, in the Cu-going direction compared to the Au-going direction. Additionally, the nuclear-modification factor in Cu$+$Au collisions averaged over all centrality is measured to be similar to the previous PHENIX result in $d$$+$Au collisions for these rapidities.

We report the first measurement of the full angular distribution for inclusive $J/\psi\rightarrow\mu^{+}\mu^{-}$ decays in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. The measurements are made for $J/\psi$ transverse momentum $2<p_{T}<10$ GeV/$c$ and rapidity $1.2<y<2.2$ in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient $\lambda_{\theta}$ is strongly negative at low $p_{T}$ and becomes close to zero at high $p_{T}$, while the azimuthal coefficient $\lambda_{\phi}$ is close to zero at low $p_{T}$, and becomes slightly negative at higher $p_{T}$. The frame-independent coefficient $\tilde{\lambda}$ is strongly negative at all $p_{T}$ in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

We have measured the azimuthal anisotropy of pi0's for 1 < pT < 18 GeV/c for Au+Au collisions at sqrt s_NN = 200 GeV. The observed anisotropy shows a gradual decrease in 3 < pT < 7 - 10 GeV/c, but remains positive beyond 10 GeV/c. The magnitude of this anisotropy is under-predicted, up to at least 10 GeV/c, by current perturbative QCD (pQCD) energy-loss model calculations. An estimate of the increase in anisotropy expected from initial-geometry modification due to gluon saturation effects and initial-geometry fluctuations is insufficient to account for this discrepancy. Calculations which implement a path length dependence steeper than what is implied by current pQCD energy-loss models, show reasonable agreement with the data.

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.

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.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured 2nd and 3rd order Fourier coefficients of the azimuthal distributions of direct photons emitted at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for various collision centralities. Combining two different analysis techniques, results were obtained in the transverse momentum range of $0.4<p_{T}<4.0$ GeV/$c$. At low $p_T$ the second-order coefficients, $v_2$, are similar to the ones observed in hadrons. Third order coefficients, $v_3$, are nonzero and almost independent of centrality. These new results on $v_2$ and $v_3$, combined with previously published results on yields, are compared to model calculations that provide yields and asymmetries in the same framework. Those models are challenged to explain simultaneously the observed large yield and large azimuthal anisotropies.

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 collaboration presents first measurements of low-momentum ($0.4<p_T<3$ GeV/$c$) direct-photon yields from Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=39 and 62.4 GeV. For both beam energies the direct-photon yields are substantially enhanced with respect to expectations from prompt processes, similar to the yields observed in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200. Analyzing the photon yield as a function of the experimental observable $dN_{\rm ch}/d\eta$ reveals that the low-momentum ($>$1\,GeV/$c$) direct-photon yield $dN_{\gamma}^{\rm dir}/d\eta$ is a smooth function of $dN_{\rm ch}/d\eta$ and can be well described as proportional to $(dN_{\rm ch}/d\eta)^\alpha$ with $\alpha{\sim}$1.25. This new scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and Large Hadron Collider, for centrality selected samples, as well as for different, $A$$+$$A$ collision systems. At a given beam energy the scaling also holds for high $p_T$ ($>5$\,GeV/$c$) but when results from different collision energies are compared, an additional $\sqrt{s_{_{NN}}}$-dependent multiplicative factor is needed to describe the integrated-direct-photon yield.

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.

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

A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV, and Cu+Cu collisions sqrt(s_NN) = 22.5, 62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for non-dynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. An analysis of Negative Binomial Distribution fits to the multiplicity distributions demonstrates that the heavy ion data exhibit weak clustering properties.

We present measurements of J/psi yields in d+Au collisions at sqrt(s_NN) = 200 GeV recorded by the PHENIX experiment and compare with yields in p+p collisions at the same energy per nucleon-nucleon collision. The measurements cover a large kinematic range in J/psi rapidity (-2.2 < y < 2.4) with high statistical precision and are compared with two theoretical models: one with nuclear shadowing combined with final state breakup and one with coherent gluon saturation effects. To remove model dependent systematic uncertainties we also compare the data to a simple geometric model. We find that calculations where the nuclear modification is linear or exponential in the density weighted longitudinal thickness are difficult to reconcile with the forward rapidity data.

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 PHENIX experiment has measured open heavy-flavor production via semileptonic decay muons over the transverse momentum range 1 < pT < 6 GeV/c at forward and backward rapidity (1.4 < |y| < 2.0) in d+Au and p+p collisions at ?sNN = 200 GeV. In central d+Au collisions an enhancement (suppression) of heavy-flavor muon production is observed at backward (forward) rapidity relative to the yield in p+p collisions scaled by the number of binary collisions. Modification of the gluon density distribution in the Au nucleus contributes in terms of anti-shadowing enhancement and shadowing suppression; however, the enhancement seen at backward rapidity exceeds expectations from this effect alone. These results, implying an important role for additional cold nuclear matter effects, serves as a key baseline for heavy-quark measurements in A+A collisions and in constraining the magnitude of charmonia breakup effects at the Relativistic Heavy Ion Collider and the Large Hadron Collider.

Correlations between p and pbar's at transverse momenta typical of enhanced baryon production in Au+Au collisions are reported. The PHENIX experiment measures same and opposite sign baryon pairs in Au+Au collisions at sqrt(s_NN) = 200 GeV. Correlated production of p and p^bar with the trigger particle from the range 2.5 < p_T < 4.0 GeV/c and the associated particle with 1.8 < p_T < 2.5 GeV/c is observed to be nearly independent of the centrality of the collisions. Same sign pairs show no correlation at any centrality. The conditional yield of mesons triggered by baryons (and anti-baryons) and mesons in the same pT range rises with increasing centrality, except for the most central collisions, where baryons show a significantly smaller number of associated mesons. These data are consistent with a picture in which hard scattered partons produce correlated p and p^bar in the p_T region of the baryon excess.

Large parity violating longitudinal single spin asymmetries A^{e^-}_L= -0.86^{+0.14}_{-0.30} and A^{e^+}_L= 0.88^{+0.12}_{-0.71} are observed for inclusive high transverse momentum electrons and positrons in polarized pp collisions at a center of mass energy of \sqrt{s}=500\ GeV with the PHENIX detector at RHIC. These e^{+/-} come mainly from the decay of W^{+/-} and Z^0 bosons, and the asymmetries directly demonstrate parity violation in the couplings of the W^{\pm} to the light quarks. The observed electron and positron yields were used to estimate W^\pm boson production cross sections equal to \sigma(pp \to W^+ X) \times BR(W^ \to \nu_e)= 144.1+/-21.2(stat)^{+3.4}_{-10.3}(syst) +/- 15%(norm) pb, and \sigma(pp \to W^{-}X) \times BR(W^\to e^-\bar{\nu_e}) = 31.7+/-12.1(stat)^{+10.1}_{-8.2}(syst)+/-15%(norm) pb.

We present a measurement of the cross section and transverse single-spin asymmetry ($A_N$) for $\eta$ mesons at large pseudorapidity from $\sqrt{s}=200$~GeV $p^{\uparrow}+p$ collisions. The measured cross section for $0.5<p_T<5.0$~GeV/$c$ and $3.0<|\eta|<3.8$ is well described by a next-to-leading-order perturbative-quantum-chromodynamics calculation. The asymmetries $A_N$ have been measured as a function of Feynman-$x$ ($x_F$) from $0.2<|x_{F}|<0.7$, as well as transverse momentum ($p_T$) from $1.0<p_T<4.5$~GeV/$c$. The asymmetry averaged over positive $x_F$ is $\langle{A_{N}}\rangle=0.061{\pm}0.014$. The results are consistent with prior transverse single-spin measurements of forward $\eta$ and $\pi^{0}$ mesons at various energies in overlapping $x_F$ ranges. Comparison of different particle species can help to determine the origin of the large observed asymmetries in $p^{\uparrow}+p$ collisions.

Measurements of double-helicity asymmetries for inclusive hadron production in polarized p+p collisions are sensitive to helicity--dependent parton distribution functions, in particular to the gluon helicity distribution, Delta(g). This study focuses on the extraction of the double-helicity asymmetry in eta production: polarized p+p --> eta + X, the eta cross section, and the eta/pi^0 cross section ratio. The cross section and ratio measurements provide essential input for the extraction of fragmentation functions that are needed to access the helicity-dependent parton distribution functions.

We have measured the cross section and single spin asymmetries from forward $W^{\pm}\rightarrow\mu^{\pm}\nu$ production in longitudinally polarized $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV using the PHENIX detector at the Relativistic Heavy Ion Collider. The cross sections are consistent with previous measurements at this collision energy, while the most forward and backward longitudinal single spin asymmetries provide new insights into the sea quark helicities in the proton. The charge of the W bosons provides a natural flavor separation of the participating partons.

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.

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.

Unpolarized cross sections and double-helicity asymmetries of single-inclusive positive and negative charged hadrons at midrapidity from p+p collisions at sqrt(s)=62.4 GeV are presented. The PHENIX measurements for 1.0 < p_T < 4.5 GeV/c are consistent with perturbative QCD calculations at next-to-leading order in the strong coupling constant, alpha_s. Resummed pQCD calculations including terms with next-to-leading-log accuracy, yielding reduced theoretical uncertainties, also agree with the data. The double-helicity asymmetry, sensitive at leading order to the gluon polarization in a momentum-fraction range of 0.05 ~< x_gluon ~< 0.2, is consistent with recent global parameterizations disfavoring large gluon polarization.

PHENIX has measured the e^+e^- pair continuum in sqrt(s_NN)=200 GeV Au+Au and p+p collisions over a wide range of mass and transverse momenta. The e^+e^- yield is compared to the expectations from hadronic sources, based on PHENIX measurements. In the intermediate mass region, between the masses of the phi and the J/psi meson, the yield is consistent with expectations from correlated c^bar-c production, though other mechanisms are not ruled out. In the low mass region (below the phi) the p+p inclusive mass spectrum is well described by known contributions from light meson decays. In contrast, the Au+Au minimum bias inclusive mass spectrum in this region shows an enhancement by a factor of 4.7+/-0.4(stat)+/-1.5(syst)+/-0.9(model) At low mass (m_ee<0.3 GeV/c^2) and high p_T (1<p_T<5 GeV/c) an enhanced e^+e^- pair yield is observed that is consistent with production of virtual direct photons. This excess is used to infer the yield of real direct photons. In central Au+Au collisions, the excess of the direct photon yield over the p+p is exponential in p_T, with inverse slope T=221+/-19(stat)+/-19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T_init ~=300--600 MeV at times of 0.6--0.15 fm/c after the collision are in qualitative agreement with the direct photon data in Au+Au. For low p_T<1 GeV/c the low mass region shows a further significant enhancement that increases with centrality and has an inverse slope of T ~=100 MeV. Theoretical models under predict the low mass, low p_T enhancement.