We report the first measurement of the full angular distribution for inclusive J/ψ→μ+μ- decays in p+p collisions at s=510 GeV. The measurements are made for J/ψ transverse momentum 2<pT<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 λθ is strongly negative at low pT and becomes close to zero at high pT, while the azimuthal coefficient λφ is close to zero at low pT, and becomes slightly negative at higher pT. The frame-independent coefficient λ is strongly negative at all pT in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

We report the first measurement of the fraction of J/ψ mesons coming from B-meson decay (FB→J/ψ) in p+p collisions at s=510 GeV. The measurement is performed using the forward silicon vertex detector and central vertex detector at PHENIX, which provide precise tracking and distance-of-closest-approach determinations, enabling the statistical separation of J/ψ due to B-meson decays from prompt J/ψ. The measured value of FB→J/ψ is 8.1%±2.3%(stat)±1.9%(syst) for J/ψ with transverse momenta 0<pT<5 GeV/c and rapidity 1.2<|y|<2.2. The measured fraction FB→J/ψ at PHENIX is compared to values measured by other experiments at higher center of mass energies and to fixed-order-next-to-leading-logarithm and color-evaporation-model predictions. The bb cross section per unit rapidity [dσ/dy(pp→bb)] extracted from the obtained FB→J/ψ and the PHENIX inclusive J/ψ cross section measured at 200 GeV scaled with color-evaporation-model calculations, at the mean B hadron rapidity y=±1.7 in 510 GeV p+p collisions, is 3.63-1.70+1.92 μb. It is consistent with the fixed-order-next-to-leading-logarithm calculations.

The PHENIX Collaboration has measured the ratio of the yields of ψ(2S) to ψ(1S) mesons produced in p+p, p+Al, p+Au, and He3+Au collisions at sNN=200 GeV over the forward and backward rapidity intervals 1.2<|y|<2.2. We find that the ratio in p+p collisions is consistent with measurements at other collision energies. In collisions with nuclei, we find that in the forward (p-going or He3-going) direction, the relative yield of ψ(2S) mesons to ψ(1S) mesons is consistent with the value measured in p+p collisions. However, in the backward (nucleus-going) direction, the ψ(2S) meson is preferentially suppressed by a factor of ∼2. This suppression is attributed in some models to the breakup of the weakly bound ψ(2S) meson through final-state interactions with comoving particles, which have a higher density in the nucleus-going direction. These breakup effects may compete with color screening in a deconfined quark-gluon plasma to produce sequential suppression of excited quarkonia states.

PHENIX measurements are presented for the cross section and double-helicity asymmetry (A(LL)) in inclusive pi(0) production at midrapidity from p + p collisions at root s = 510 GeV from data taken in 2012 and 2013 at the Relativistic Heavy Ion Collider. The next-to-leading-order perturbative-quantum-chromodynamics theory calculation is in excellent agreement with the presented cross section results. The calculation utilized parton-to-pion fragmentation functions from the recent DSS14 global analysis, which prefer a smaller gluon-to-pion fragmentation function. The pi(0)A(LL) results follow an increasingly positive asymmetry trend with p(T) and root s with respect to the predictions and are in excellent agreement with the latest global analysis results. This analysis incorporated earlier results on pi(0) and jet A(LL) and suggested a positive contribution of gluon polarization to the spin of the proton Delta G for the gluon momentum fraction range x > 0.05. The data presented here extend to a currently unexplored region, down to x similar to 0.01, and thus provide additional constraints on the value of Delta G.

We report on J/psi production from asymmetric Cu + Au heavy-ion collisions at root S-NN = 200 GeV at the Relativistic Heavy Ion Collider at both forward (Cu-going direction) and backward (Au-going direction) rapidities. The nuclear modification of J/psi yields in Cu + Au collisions in the Au-going direction is found to be comparable to that inAu + Au collisions when plotted as a function of the number of participating nucleons. In the Cu-going direction, J/psi production shows a stronger suppression. This difference is comparable in magnitude and has the same sign as the difference expected from shadowing effects due to stronger low-x gluon suppression in the larger Au nucleus.

Measurements of electrons from the decay of open-heavy-flavor mesons have shown that the yields are suppressed in Au+Au collisions compared to expectations from binary-scaled p+p collisions. These measurements indicate that charm and bottom quarks interact with the hot dense matter produced in heavy-ion collisions much more than expected. Here we extend these studies to two-particle correlations where one particle is an electron from the decay of a heavy-flavor meson and the other is a charged hadron from either the decay of the heavy meson or from jet fragmentation. These measurements provide more detailed information about the interactions between heavy quarks and the matter, such as whether the modification of the away-side-jet shape seen in hadron-hadron correlations is present when the trigger particle is from heavy-meson decay and whether the overall level of away-side-jet suppression is consistent. We statistically subtract correlations of electrons arising from background sources from the inclusive electron-hadron correlations and obtain two-particle azimuthal correlations at root s(NN) = 200 GeV between electrons from heavy-flavor decay with charged hadrons in p+p and also first results in Au+Au collisions. We find the away-side-jet shape and yield to be modified in Au+Au collisions compared to p+p collisions.

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 root 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 Rp(d)+A, for which there is a bias in the measured centrality-dependent yields owing to auto correlations between the process of interest and the backward-rapidity multiplicity. We determine the bias-correction factors within this framework. This method is further tested using the HIJING Monte Carlo generator. We find that for d + Au collisions at root s(NN) = 200 GeV, these bias corrections are small and vary by less than 5% (10%) up to p(T) = 10 (20) GeV/c. In contrast, for p + Pb collisions at v root s(NN) = 5.02 TeV we find that these bias factors are an order of magnitude larger and strongly pT dependent, likely attributable to the larger effect of multiparton interactions.

The PHENIX experiment has measured open heavy-flavor production via semileptonic decay over the transverse momentum range 1 < p(T) < 6 GeV/c at forward and backward rapidity (1.4 < vertical bar y vertical bar < 2.0) in d + Au and p + p collisions at root s(NN) = 200 GeV. In central d + Au collisions, relative to the yield in p + p collisions scaled by the number of binary nucleon-nucleon collisions, a suppression is observed at forward rapidity (in the d-going direction) and an enhancement at backward rapidity (in the Au-going direction). Predictions using nuclear-modified-parton-distribution functions, even with additional nuclear-p(T) broadening, cannot simultaneously reproduce the data at both rapidity ranges, which implies that these models are incomplete and suggests the possible importance of final-state interactions in the asymmetric d + Au collision system. These results can be used to probe cold-nuclear-matter effects, which may significantly affect heavy-quark production, in addition to helping constrain the magnitude of charmonia-breakup effects in nuclear matter.

Measurements of double-helicity asymmetries in 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 ((p) over right arrow + (p) over right arrow -> 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.

Transverse momentum distributions and yields for pi(+/-), K-+/-, p, and (p) over bar in p + p collisions at root s = 200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter T-inv, mean transverse momentum < p(T)>, and yield per unit rapidity dN/dy at each energy, and compare them to other measurements at different root s in p + p and p + (p) over bar collisions. We also present the scaling properties such as m(T) scaling and x(T) scaling on the p(T) spectra between different energies. To discuss the mechanism of the particle production in p + p collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.