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.

Measurements of neutral pion production at midrapidity in sqrt(s_NN) = 200 GeV Au+Au collisions as a function of transverse momentum, p_T, collision centrality, and angle with respect to reaction plane are presented. The data represent the final pi^0 results from the PHENIX experiment for the first RHIC Au+Au run at design center-of-mass-energy. They include additional data obtained using the PHENIX Level-2 trigger with more than a factor of three increase in statistics over previously published results for p_T > 6 GeV/c. We evaluate the suppression in the yield of high-p_T pi^0's relative to point-like scaling expectations using the nuclear modification factor R_AA. We present the p_T dependence of R_AA for nine bins in collision centrality. We separately integrate R_AA over larger p_T bins to show more precisely the centrality dependence of the high-p_T suppression. We then evaluate the dependence of the high-p_T suppression on the emission angle \Delta\phi of the pions with respect to event reaction plane for 7 bins in collision centrality. We show that the yields of high-p_T pi^0's vary strongly with \Delta\phi, consistent with prior measurements. We show that this variation persists in the most peripheral bin accessible in this analysis. For the peripheral bins we observe no suppression for neutral pions produced aligned with the reaction plane while the yield of pi^0's produced perpendicular to the reaction plane is suppressed by more than a factor of 2. We analyze the combined centrality and \Delta\phi dependence of the pi^0 suppression in different p_T bins using different possible descriptions of parton energy loss dependence on jet path-length averages to determine whether a single geometric picture can explain the observed suppression pattern.

Transverse momentum spectra of charged hadrons with p_T < 8 GeV/c and neutral pions with p_T < 10 GeV/c have been measured at mid-rapidity by the PHENIX experiment at RHIC in d+Au collisions at sqrt(s_NN) = 200 GeV. The measured yields are compared to those in p+p collisions at the same sqrt(s_NN) scaled up by the number of underlying nucleon-nucleon collisions in d+Au. The yield ratio does not show the suppression observed in central Au+Au collisions at RHIC. Instead, there is a small enhancement in the yield of high momentum particles.

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 report on two-particle azimuthal angle correlations between charged hadrons at forward/backward (deuteron/gold going direction) rapidity and charged hadrons at mid-rapidity in deuteron-gold (d+Au) and proton-proton (p+p) collisions at sqrt(s_NN) = 200 GeV. Jet structures are observed in the correlations which we quantify in terms of the conditional yield and angular width of away side partners. The kinematic region studied here samples partons in the gold nucleus carrying nucleon momentum fraction x~0.1 to x~0.01. Within this range, we find no x dependence of the jet structure in d+Au collisions.

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

Bose-Einstein correlations of identically charged pion pairs were measured by the PHENIX experiment at mid-rapidity in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The Bertsch-Pratt radius parameters were determined as a function of the transverse momentum of the pair and as a function of the centrality of the collision. Using the \it{full} Coulomb correction, the ratio $R_{\rm out}/R_{\rm side}$ is smaller than unity for $<k_{\rm T}>$ from 0.25 to 1.2 GeV/c and for all measured centralities. However, using recently developed partial Coulomb correction methods, we find that $R_{\rm out}/R_{\rm side}$ is 0.8-1.1 for the measured $<k_{\rm T}>$ range, and approximately constant at unity with the number of participants.

We present transverse momentum (p_T) spectra of charged hadrons measured in deuteron-gold and nucleon-gold collisions at \sqrts = 200 GeV for four centrality classes. Nucleon-gold collisions were selected by tagging events in which a spectator nucleon was observed in one of two forward rapidity detectors. The spectra and yields were investigated as a function of the number of binary nucleon-nucleon collisions, \nu, suffered by deuteron nucleons. A comparison of charged particle yields to those in p+p collisions show that the yield per nucleon-nucleon collision saturates with \nu for high momentum particles. We also present the charged hadron to neutral pion ratios as a function of p_T.

We present results for the charged-particle multiplicity distribution at mid-rapidity in Au - Au collisions at sqrt(s_NN)=130 GeV measured with the PHENIX detector at RHIC. For the 5% most central collisions we find $dN_{ch}/d\eta_{|\eta=0} = 622 \pm 1 (stat) \pm 41 (syst)$. The results, analyzed as a function of centrality, show a steady rise of the particle density per participating nucleon with centrality.

The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < p_T < 4.0 GeV/c) of single electrons as a function of centrality in Au+Au collisions at sqrt(s_NN) = 200 GeV. Contributions to the raw spectra from photon conversions and Dalitz decays of light neutral mesons are measured by introducing a thin (1.7% X_0) converter into the PHENIX acceptance and are statistically removed. The subtracted ``non-photonic'' electron spectra are primarily due to the semi-leptonic decays of hadrons containing heavy quarks (charm and bottom). For all centralities, charm production is found to scale with the nuclear overlap function, T_AA. For minimum-bias collisions the charm cross section per binary collision is N_cc^bar/T_AA = 622 +/- 57 (stat.) +/- 160 (sys.) microbarns.

The first measurement of direct photons in Au+Au collisions at sqrt(s_NN) = 200 GeV is presented. The direct photon signal is extracted as a function of the Au+Au collision centrality and compared to NLO pQCD calculations. The direct photon yield is shown to scale with the number of nucleon-nucleon collisions for all centralities.

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

Identified pi^[+/-] K^[+/-], p and p-bar transverse momentum spectra at mid-rapidity in sqrt(s_NN)=130 GeV Au-Au collisions were measured by the PHENIX experiment at RHIC as a function of collision centrality. Average transverse momenta increase with the number of participating nucleons in a similar way for all particle species. The multiplicity densities scale faster than the number of participating nucleons. Kaon and nucleon yields per participant increase faster than the pion yields. In central collisions at high transverse momenta (p_T greater than 2 GeV/c), anti-proton and proton yields are comparable to the pion yields.

The dependence of transverse momentum spectra of neutral pions and eta mesons with p_T <16 GeV/c and p_T < 12 GeV/c, respectively, on the centrality of the collision has been measured at mid-rapidity by the PHENIX experiment at RHIC in d+Au collisions at sqrt(s_(NN)) = 200 GeV. The measured yields are compared to those in p + p collisions at the same sqrt(s_(NN)) scaled by the number of underlying nucleon-nucleon collisions in d+Au. At all centralities the yield ratios show no suppression, in contrast to the strong suppression seen for central Au+Au collisions at RHIC. Only a weak p_T and centrality dependence can be observed.

PHENIX has measured the centrality dependence of charged hadron p_T spectra from central Au+Au collisions at sqrt(s_NN)=130 GeV. The truncated mean p_T decreases with centrality for p_T > 2 GeV/c, indicating an apparent reduction of the contribution from hard scattering to high p_T hadron production. For central collisions the yield at high p_T is shown to be suppressed compared to binary nucleon-nucleon collision scaling of p+p data. This suppression is monotonically increasing with centrality, but most of the change occurs below 30% centrality, i.e. for collisions with less than about 140 participating nucleons. The observed p_T and centrality dependence is consistent with the particle production predicted by models including hard scattering and subsequent energy loss of the scattered partons in the dense matter created in the collisions.

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.

The standard model (SM) of particle physics is spectacularly successful, yet the measured value of the muon anomalous magnetic moment $(g-2)_\mu$ deviates from SM calculations by 3.6$\sigma$. Several theoretical models attribute this to the existence of a "dark photon," an additional U(1) gauge boson, which is weakly coupled to ordinary photons. The PHENIX experiment at the Relativistic Heavy Ion Collider has searched for a dark photon, $U$, in $\pi^0,\eta \rightarrow \gamma e^+e^-$ decays and obtained upper limits of $\mathcal{O}(2\times10^{-6})$ on $U$-$\gamma$ mixing at 90% CL for the mass range $30<m_U<90$ MeV/$c^2$. Combined with other experimental limits, the remaining region in the $U$-$\gamma$ mixing parameter space that can explain the $(g-2)_\mu$ deviation from its SM value is nearly completely excluded at the 90% confidence level, with only a small region of $29<m_U<32$ MeV/$c^2$ remaining.

All of the experimental data points presented in the original paper are correct and unchanged (including statistical and systematic uncertainties). However, herein we correct a comparison between the experimental data and a theoretical picture, because we discovered a mistake in the code used. All of the most probable sigma_breakup values differ by less than 0.4 mb from those originally presented. However, the one standard deviation uncertainties (that include contributions from both the statistical and systematic uncertainties on the experimental data points) are approximately 30-60% larger than originally reported. We give a table of the new comparison results and corrected versions of Figs. 8-11 of the original paper and we note that no correction is needed for results from the data-driven method in Fig. 13.

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.

Inclusive transverse momentum spectra of eta mesons have been measured within p_T = 2-10 GeV/c at mid-rapidity by the PHENIX experiment in Au+Au collisions at sqrt(s_NN) = 200 GeV. In central Au+Au the eta yields are significantly suppressed compared to peripheral Au+Au, d+Au and p+p yields scaled by the corresponding number of nucleon-nucleon collisions. The magnitude, centrality and p_T dependence of the suppression is common, within errors, for eta and pi^0. The ratio of eta to pi^0 spectra at high p_T amounts to 0.40 < R_eta/pi^0 < 0.48 for the three systems in agreement with the world average measured in hadronic and nuclear reactions and, at large scaled momentum, in e^+e^- collisions.