The first measurement of energy produced transverse to the beam direction at RHIC is presented. The mid-rapidity transverse energy density per participating nucleon rises steadily with the number of participants, closely paralleling the rise in charged-particle density, such that E_T / N_ch remains relatively constant as a function of centrality. The energy density calculated via Bjorken's prescription for the 2% most central Au+Au collisions at sqrt(s_NN)=130 GeV is at least epsilon_Bj = 4.6 GeV/fm^3 which is a factor of 1.6 larger than found at sqrt(s_NN)=17.2 GeV (Pb+Pb at CERN).

The invariant differential cross section for inclusive neutral pion production in p+p collisions at sqrt(s_NN) = 200 GeV has been measured at mid-rapidity |eta| < 0.35 over the range 1 < p_T <~ 14 GeV/c by the PHENIX experiment at RHIC. Predictions of next-to-leading order perturbative QCD calculations are consistent with these measurements. The precision of our result is sufficient to differentiate between prevailing gluon-to-pion fragmentation functions.

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.

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

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.

J/Psi production in p+p collisions at sqrt(s) = 200 GeV has been Measured in the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over a rapidity range of -2.2 < y < 2.2 and a transverse momentum range of 0 < pT < 9 GeV/c. The statistics available allow a detailed measurement of both the pT and rapidity distributions and are sufficient to constrain production models. The total cross section times branching ratio determined for J/Psi production is B_{ll} sigma_pp^J/psi = 178 +/- 3(stat) +/- 53(syst) +/- 18(norm) nb.

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.

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.