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.

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.

The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < p_T < 5.0 GeV/c) of electrons as a function of centrality in Au+Au collisions at sqrt(s_NN)=200 GeV. Contributions from photon conversions and from light hadron decays, mainly Dalitz decays of pi^0 and eta mesons, were removed. The resulting non-photonic electron spectra are primarily due to the semi-leptonic decays of hadrons carrying heavy quarks. Nuclear modification factors were determined by comparison to non-photonic electrons in p+p collisions. A significant suppression of electrons at high p_T is observed in central Au+Au collisions, indicating substantial energy loss of heavy quarks.

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.

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 momentum distribution of electrons from decays of heavy flavor (charm and beauty) for midrapidity |y| < 0.35 in p+p collisions at sqrt(s) = 200 GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over the transverse momentum range 0.3 < p_T < 9 GeV/c. Two independent methods have been used to determine the heavy flavor yields, and the results are in good agreement with each other. A fixed-order-plus-next-to-leading-log pQCD calculation agrees with the data within the theoretical and experimental uncertainties, with the data/theory ratio of 1.72 +/- 0.02^stat +/- 0.19^sys for 0.3 < p_T < 9 GeV/c. The total charm production cross section at this energy has also been deduced to be sigma_(c c^bar) = 567 +/- 57^stat +/- 224^sys micro barns.

We present the first measurement of photoproduction of J/psi and of two-photon production of high-mass e+e- pairs in electromagnetic (or ultra-peripheral) nucleus-nucleus interactions, using Au+Au data at sqrt(s_NN) = 200 GeV. The events are tagged with forward neutrons emitted following Coulomb excitation of one or both Au^{star} nuclei. The event sample consists of 28 events with m_{e+e-} > 2 GeV/c^2 with zero like-sign background. The measured cross sections at midrapidity of d\sigma / dy (J/psi + Xn, y=0) = 76 +/- 33 (stat) +/- 11 (syst) micro b and d^2\sigma/dm dy (e^+e^- + Xn, y=0) = 86 +/- 23 (stat) +/- 16 (syst) micro b/(GeV/c^2) for m_{e+e-} \in [2.0,2.8] GeV/c^2 are consistent with various theoretical predictions.

First results on charm quarkonia production in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) are presented. The yield of J/Psi's measured in the PHENIX experiment via electron-positron decay pairs at mid-rapidity for Au-Au reactions at sqrt(s_NN) = 200 GeV are analyzed as a function of collision centrality. For this analysis we have studied 49.3 million minimum bias Au-Au reactions. We present the J/Psi invariant yield dN/dy for peripheral and mid-central reactions. For the most central collisions where we observe no signal above background, we quote 90% confidence level upper limits. We compare these results with our J/Psi measurement from proton-proton reactions at the same energy. We find that our measurements are not consistent with models that predict strong enhancement relative to binary collision scaling.

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 PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured J/psi production for rapidities 2.2 < y < 2.2 in Au+Au collisions at sqrt(s_NN) = 200 GeV. The J/psi invariant yield and nuclear modification factor R_AA as a function of centrality, transverse momentum and rapidity are reported. A suppression of J/psi relative to binary collision scaling of proton-proton reaction yields is observed. Models which describe the lower energy J/Psi data at the Super Proton Synchrotron (SPS) invoking only J/psi destruction based on the local medium density would predict a significantly larger suppression at RHIC and more suppression at mid rapidity than at forward rapidity. Both trends are contradicted by our data.