Measurements of the midrapidity transverse energy distribution, $d\Et/d\eta$, are presented for $p$$+$$p$, $d$$+$Au, and Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV and additionally for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ and 130 GeV. The $d\Et/d\eta$ distributions are first compared with the number of nucleon participants $N_{\rm part}$, number of binary collisions $N_{\rm coll}$, and number of constituent-quark participants $N_{qp}$ calculated from a Glauber model based on the nuclear geometry. For Au$+$Au, $\mean{d\Et/d\eta}/N_{\rm part}$ increases with $N_{\rm part}$, while $\mean{d\Et/d\eta}/N_{qp}$ is approximately constant for all three energies. This indicates that the two component ansatz, $dE_{T}/d\eta \propto (1-x) N_{\rm part}/2 + x N_{\rm coll}$, which has been used to represent $E_T$ distributions, is simply a proxy for $N_{qp}$, and that the $N_{\rm coll}$ term does not represent a hard-scattering component in $E_T$ distributions. The $dE_{T}/d\eta$ distributions of Au$+$Au and $d$$+$Au are then calculated from the measured $p$$+$$p$ $E_T$ distribution using two models that both reproduce the Au$+$Au data. However, while the number-of-constituent-quark-participant model agrees well with the $d$$+$Au data, the additive-quark model does not.

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.

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.

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.

The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) reports measurements of azimuthal dihadron correlations near midrapidity in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. These measurements complement recent analyses by experiments at the Large Hadron Collider (LHC) involving central $p$$+$Pb collisions at $\sqrt{s_{_{NN}}}$=5.02 TeV, which have indicated strong anisotropic long-range correlations in angular distributions of hadron pairs. The origin of these anisotropies is currently unknown. Various competing explanations include parton saturation and hydrodynamic flow. We observe qualitatively similar, but larger, anisotropies in $d$$+$Au collisions compared to those seen in $p$$+$Pb collisions at the LHC. The larger extracted $v_2$ values in $d$$+$Au collisions at RHIC are consistent with expectations from hydrodynamic calculations owing to the larger expected initial-state eccentricity compared with that from $p$$+$Pb collisions. When both are divided by an estimate of the initial-state eccentricity the scaled anisotropies follow a common trend with multiplicity that may extend to heavy ion data at RHIC and the LHC, where the anisotropies are widely thought to arise from hydrodynamic flow.

Muon production at forward rapidity (1.5 < |\eta| < 1.8) has been measured by the PHENIX experiment over the transverse momentum range 1 < p_T \le 3 GeV/c in sqrt(s) = 200 GeV p+p collisions at the Relativistic Heavy Ion Collider. After statistically subtracting contributions from light hadron decays an excess remains which is attributed to the semileptonic decays of hadrons carrying heavy flavor, i.e. charm quarks or, at high p_T, bottom quarks. The resulting muon spectrum from heavy flavor decays is compared to PYTHIA and a next-to-leading order perturbative QCD calculation. PYTHIA is used to determine the charm quark spectrum that would produce the observed muon excess. The corresponding differential cross section for charm quark production at forward rapidity is determined to be d\sigmac c^bar)/dy|_(y=1.6)=0.243 +/- 0.013 (stat.) +/- 0.105 (data syst.) ^(+0.049(-0.087) (PYTHIA syst.) mb.

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 experiement has measured the electron-positron pair mass spectrum from 0 to 8 GeV/c^2 in p+p collisions at sqrt(s)=200 GeV. The contributions from light meson decays to e^+e^- pairs have been determined based on measurements of hadron production cross sections by PHENIX. They account for nearly all e^+e^- pairs in the mass region below 1 GeV/c^2. The e^+e^- pair yield remaining after subtracting these contributions is dominated by semileptonic decays of charmed hadrons correlated through flavor conservation. Using the spectral shape predicted by PYTHIA, we estimate the charm production cross section to be 544 +/- 39(stat) +/- 142(syst) +/- 200(model) \mu b, which is consistent with QCD calculations and measurements of single leptons by PHENIX.

We present the midrapidity charged pion invariant cross sections and the ratio of $\pi^-$-to-$\pi^+$ production ($5<p_T<13$ GeV/$c$), together with the double-helicity asymmetries ($5<p_T<12$ GeV/$c$) in polarized $p$$+$$p$ collisions at $\sqrt{s} = 200$ GeV. The cross section measurements are consistent with perturbative calculations in quantum chromodynamics within large uncertainties in the calculation due to the choice of factorization, renormalization, and fragmentation scales. However, the theoretical calculation of the ratio of $\pi^-$-to-$\pi^+$ production when considering these scale uncertainties overestimates the measured value, suggesting further investigation of the uncertainties on the charge-separated pion fragmentation functions is needed. Due to cancellations of uncertainties in the charge ratio, direct inclusion of these ratio data in future parameterizations should improve constraints on the flavor dependence of quark fragmentation functions to pions. By measuring charge-separated pion asymmetries, one can gain sensitivity to the sign of $\Delta G$ through the opposite sign of the up and down quark helicity distributions in conjunction with preferential fragmentation of positive pions from up quarks and negative pions from down quarks. The double-helicity asymmetries presented are sensitive to the gluon helicity distribution over an $x$ range of $\sim$0.03--0.16.

We present measurements from the PHENIX experiment of large parity-violating single spin asymmetries of high transverse momentum electrons and positrons from $W^\pm/Z$ decays, produced in longitudinally polarized $p$$+$$p$ collisions at center of mass energies of $\sqrt{s}$=500 and 510~GeV. These asymmetries allow direct access to the anti-quark polarized parton distribution functions due to the parity-violating nature of the $W$-boson coupling to quarks and anti-quarks. The results presented are based on data collected in 2011, 2012, and 2013 with an integrated luminosity of 240 pb$^{-1}$, which exceeds previous PHENIX published results by a factor of more than 27. These high $Q^2$ data provide an important addition to our understanding of anti-quark parton helicity distribution functions.