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.

Transverse momentum distributions and yields for $\pi^{\pm}$, $K^{\pm}$, $p$ and $\bar{p}$ in $p+p$ collisions at $\sqrt{s}$=200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter $T_{\rm inv}$, mean transverse momentum $<p_T>$ and yield per unit rapidity $dN/dy$ at each energy, and compare them to other measurements at different $\sqrt{s}$ in $p+p$ and $p+\bar{p}$ collisions. We also present the scaling properties such as $m_T$ scaling, $x_T$ scaling on the $p_T$ spectra between different energies. To discuss the mechanism of the particle production in $p+p$ collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.

The differential cross section for the process $\gamma \gamma \to \eta \pi^0$ has been measured in the kinematic range $0.84 \GeV < W < 4.0 \GeV$, $|\cos \theta^*|<0.8$, where $W$ and $\theta^*$ are the energy and $\pi^0$ (or $\eta$) scattering angle, respectively, in the $\gamma\gamma$ center-of-mass system. The results are based on a 223 fb$^{-1}$ data sample collected with the Belle detector at the KEKB $e^+ e^-$ collider. Clear peaks due to the $a_0(980)$ and $a_2(1320)$ are visible. The differential cross sections are fitted in the energy region $0.9 \GeV < W < 1.46 \GeV$ to obtain the parameters of the $a_0(980)$. Its mass, width and $\Gamma_{\gamma \gamma} \B (\eta \pi^0)$ are measured to be $982.3 ^{+0.6}_{-0.7} ^{+3.1}_{-4.7} \MeV/c^2$, $75.6 \pm 1.6 ^{+17.4}_{-10.0} \MeV$ and $128 ^{+3}_{-2} ^{+502}_{-43} \eV$, respectively. The energy and angular dependences above 3.1 GeV are compared with those measured in the $\pi^0 \pi^0$ channel. The integrated cross section over $|\cos \theta^*|<0.8$ has a $W^{-n}$ dependence with $n = 10.5 \pm 1.2 \pm 0.5$, which is slightly larger than that for $\pi^0 \pi^0$. The differential cross sections show a $\sin^{-4} \theta^*$ dependence similar to $\gamma \gamma \to \pi^0 \pi^0$. The measured cross section ratio, $\sigma(\eta \pi^0)/\sigma(\pi^0 \pi^0) = 0.48 \pm 0.05 \pm 0.04$, is consistent with a QCD-based prediction.

We report a high-statistics measurement of differential cross sections for the process gamma gamma -> pi^0 pi^0 in the kinematic range 0.6 GeV <= W <= 4.0 GeV and |cos theta*| <= 0.8, where W and theta* are the energy and pion scattering angle, respectively, in the gamma gamma center-of-mass system. Differential cross sections are fitted to obtain information on S, D_0, D_2, G_0 and G_2 waves. The G waves are important above W ~= 1.6 GeV. For W <= 1.6 GeV the D_2 wave is dominated by the f_2(1270) resonance while the S wave requires at least one additional resonance besides the f_0(980), which may be the f_0(1370) or f_0(1500). The differential cross sections are fitted with a simple parameterization to determine the parameters (the mass, total width and Gamma_{gamma gamma}B(f_0 -> pi^0 pi^0)) of this scalar meson as well as the f_0(980). The helicity 0 fraction of the f_2(1270) meson, taking into account interference for the first time, is also obtained.

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

Using a secondary pion beam from the Argonne Zero Gradient Synchrotron we have studied the process π−p→φn in the region of the cross-section enhancement near kinematic threshold. For incident momenta between 1.6 and 2 GeV/c, we have determined production and decay angular distributions and extrapolated total cross sections from a sample of about 160 φ's above background. The production and decay distributions are consistent with isotropy over this entire incident-momentum range. The extrapolated total cross section varies between 19 and 25 μb.