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.

Longitudinal density correlations of produced matter in Au+Au collisions at sqrt(s_NN)=200 GeV have been measured from the inclusive charged particle distributions as a function of pseudorapidity window sizes. The extracted \alpha \xi parameter, related to the susceptibility of the density fluctuations in the long wavelength limit, exhibits a non-monotonic behavior as a function of the number of participant nucleons, N_part. A local maximum is seen at N_part ~ 90, with corresponding energy density based on the Bjorken picture of \epsilon_Bj \tau ~ 2.4 GeV/(fm^2 c) with a transverse area size of 60 fm^2. This behavior may suggest a critical phase boundary based on the Ginzburg-Landau framework.

Transverse momentum spectra of electrons from Au+Au collisions at sqrt(s_NN) = 130 GeV have been measured by the PHENIX experiment at RHIC. The spectra show an excess above the background from photon conversions and light hadron decays. The electron signal is consistent with that expected from semi-leptonic decays of charm. The yield of the electron signal dN_e/dy for p_T > 0.8 GeV/c is 0.025 +/- 0.004 (stat.) +/- 0.010 (sys.) in central collisions, and the corresponding charm cross section is 380 +/- 60 (stat.) +/- 200 (sys.) micro barns per binary nucleon-nucleon collision.

Transverse momentum spectra for charged hadrons and for neutral pions in the range 1 GeV/c $< p_T <$ 5 GeV/c have been measured by the PHENIX experiment at RHIC in Au+Au collisions at $\sqrt{s_{_{NN}}}=130$ GeV. At high $p_T$ the spectra from peripheral nuclear collisions are consistent with the naive expectation of scaling the spectra from p+p collisions by the average number of binary nucleon- nucleon collisions. The spectra from central collisions are significantly suppressed when compared to the binary- scaled p+p expectation, and also when compared to similarly binary-scaled peripheral collisions, indicating a novel nuclear effect in central nuclear collisions at RHIC energies.

Transverse momentum spectra of neutral pions in the range 1 < p_T < 10 GeV/c have been measured at mid-rapidity by the PHENIX experiment at RHIC in Au+Au collisions at sqrt(s_NN) = 200 GeV. The pi^0 multiplicity in central reactions is significantly below the yields measured at the same sqrt(s_NN) in peripheral Au+Au and p+p reactions scaled by the number of nucleon-nucleon collisions. For the most central bin, the suppression factor is ~2.5 at p_T = 2 GeV/c and increases to ~4-5 at p_T ~= 4 GeV/c. At larger p_T, the suppression remains constant within errors. The deficit is already apparent in semi-peripheral reactions and increases smoothly with centrality.

We present the first results of meson production in the K^+K^- decay channel from Au+Au collisions at sqrt(s_NN) = 200 GeV as measured at mid-rapidity by the PHENIX detector at RHIC. Precision resonance centroid and width values are extracted as a function of collision centrality. No significant variation from the PDG accepted values is observed. The transverse mass spectra are fitted with a linear exponential function for which the derived inverse slope parameter is seen to be constant as a function of centrality. These data are also fitted by a hydrodynamic model with the result that the freeze-out temperature and the expansion velocity values are consistent with the values previously derived from fitting single hadron inclusive data. As a function of transverse momentum the collisions scaled peripheral.to.central yield ratio RCP for the is comparable to that of pions rather than that of protons. This result lends support to theoretical models which distinguish between baryons and mesons instead of particle mass for explaining the anomalous proton yield.

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.

Two-pion correlations in sqrt(s_NN)=130 GeV Au+Au collisions at RHIC have been measured over a broad range of pair transverse momentum k_T by the PHENIX experiment at RHIC. The k_T dependent transverse radii are similar to results from heavy ion collisions at sqrt(s_NN) = 4.1, 4.9, and 17.3 GeV, whereas the longitudinal radius increases monotonically with beam energy. The ratio of the outwards to sidewards transverse radii (R_out/R_side) is consistent with unity and independent of k_T.

We present results on the measurement of lambda and lambda^bar production in Au+Au collisions at sqrt(s_NN)=130 GeV with the PHENIX detector at RHIC. The transverse momentum spectra were measured for minimum bias and for the 5% most central events. The lambda^bar/lambda ratios are constant as a function of p_T and the number of participants. The measured net lambda density is significantly larger than predicted by models based on hadronic strings (e.g. HIJING) but in approximate agreement with models which include the gluon junction mechanism.

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.