J/psi production has been measured in proton-proton collisions at sqrt(s)= 200 GeV over a wide rapidity and transverse momentum range by the PHENIX experiment at RHIC. Distributions of the rapidity and transverse momentum, along with measurements of the mean transverse momentum and total production cross section are presented and compared to available theoretical calculations. The total J/psi cross section is 3.99 +/- 0.61(stat) +/- 0.58(sys) +/- 0.40(abs) micro barns. The mean transverse momentum is 1.80 +/- 0.23(stat) +/- 0.16(sys) GeV/c.

The PHENIX experiment at the Relativistic Heavy Ion Collider measured $\pi^0$ and $\eta$ mesons at midrapidity in U$+$U collisions at $\sqrt{s_{_{NN}}}=192$ GeV in a wide transverse momentum range. Measurements were performed in the $\pi^0(\eta)\rightarrow\gamma\gamma$ decay modes. A strong suppression of $\pi^0$ and $\eta$ meson production at high transverse momentum was observed in central U$+$U collisions relative to binary scaled $p$$+$$p$ results. Yields of $\pi^0$ and $\eta$ mesons measured in U$+$U collisions show similar suppression pattern to the ones measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for similar numbers of participant nucleons. The $\eta$/$\pi^0$ ratios do not show dependence on centrality or transverse momentum, and are consistent with previously measured values in hadron-hadron, hadron-nucleus, nucleus-nucleus, and $e^+e^-$ collisions.

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.

A measurement of direct photons in p+p collisions at sqrt(s)=200 GeV is presented. A photon excess above background from pi^0 --> gamma+gamma, eta --> gamma+gamma, and other decays is observed in the transverse momentum range 5.5 < p_T < 7 GeV/c. The result is compared to a next-to-leading-order perturbative QCD calculation. Within errors, good agreement is found between the QCD calculation and the measured result.

Production of $\pi^0$ and $\eta$ mesons has been measured at midrapidity in Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. Measurements were performed in $\pi^0(\eta)\rightarrow\gamma\gamma$ decay channel in the 1(2)-20 GeV/$c$ transverse momentum range. A strong suppression is observed for $\pi^0$ and $\eta$ meson production at high transverse momentum in central Cu$+$Au collisions relative to the $p$$+$$p$ results scaled by the number of nucleon-nucleon collisions. In central collisions the suppression is similar to Au$+$Au with comparable nuclear overlap. The $\eta/\pi^0$ ratio measured as a function of transverse momentum is consistent with $m_T$-scaling parameterization down to $p_T=$2 GeV/$c$, its asymptotic value is constant and consistent with Au$+$Au and $p$$+$$p$ and does not show any significant dependence on collision centrality. Similar results were obtained in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions as well as in $e^+e^-$ collisions in a range of collision energies $\sqrt{s_{_{NN}}}=$3--1800 GeV. This suggests that the quark-gluon-plasma medium produced in Cu$+$Cu collisions either does not affect the jet fragmentation into light mesons or it affects the $\pi^0$ and $\eta$ the same way.

The differential cross section for the production of direct photons in p+p collisions at sqrt(s)=200 GeV at midrapidity was measured in the PHENIX detector at the Relativistic Heavy Ion Collider. Inclusive-direct photons were measured in the transverse-momentum range from 5.5--25 GeV/c, extending the range beyond previous measurements. Event structure was studied with an isolation criterion. Next-to-leading-order perturbative-quantum-chromodynamics calculations give a good description of the spectrum. When the cross section is expressed versus x_T, the PHENIX data are seen to be in agreement with measurements from other experiments at different center-of-mass energies.

Cross sections for mid-rapidity production of direct photons in p+p collisions at the Relativistic Heavy Ion Collider (RHIC) are reported for 3 < p_T < 16 GeV/c. Next-to-leading order (NLO) perturbative QCD (pQCD) describes the data well for p_T > 5 GeV/c, where the uncertainties of the measurement and theory are comparable. We also report on the effect of requiring the photons to be isolated from parton jet energy. The observed fraction of isolated photons is well described by pQCD for p_T > 7 GeV/c.

The properties of jets produced in p+p collisions at sqrt(s)=200 GeV are measured using the method of two particle correlations. The trigger particle is a leading particle from a large transverse momentum jet while the associated particle comes from either the same jet or the away-side jet. Analysis of the angular width of the near-side peak in the correlation function determines the jet fragmentation transverse momentum j_T . The extracted value, sqrt(<j_T^2>)= 585 +/- 6(stat) +/- 15(sys) MeV/c, is constant with respect to the trigger particle transverse momentum, and comparable to the previous lower sqrt(s) measurements. The width of the away-side peak is shown to be a convolution of j_T with the fragmentation variable, z, and the partonic transverse momentum, k_T . The <z> is determined through a combined analysis of the measured pi^0 inclusive and associated spectra using jet fragmentation functions measured in e^+e^-. collisions. The final extracted values of k_T are then determined to also be independent of the trigger particle transverse momentum, over the range measured, with value of sqrt(<k_T^2>) = 2.68 +/- 0.07(stat) +/- 0.34(sys) GeV/c.

The cross section of bottom quark-antiquark ($b\bar{b}$) production in $p$+$p$ collisions at $\sqrt{s}=510$ GeV is measured with the PHENIX detector at the Relativistic Heavy Ion Collider. The results are based on the yield of high mass, like-sign muon pairs measured within the PHENIX muon arm acceptance ($1.2<|y|<2.2$). The $b\bar{b}$ signal is extracted from like-sign dimuons by utilizing the unique properties of neutral $B$ meson oscillation. We report a differential cross section of $d\sigma_{b\bar{b}\rightarrow \mu^\pm\mu^\pm}/dy = 0.16 \pm 0.01~(\mbox{stat}) \pm 0.02~(\mbox{syst}) \pm 0.02~(\mbox{global})$ nb for like-sign muons in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $p_T>1$ GeV/$c$, and dimuon mass of 5--10 GeV/$c^2$. The extrapolated total cross section at this energy for $b\bar{b}$ production is $13.1 \pm 0.6~(\mbox{stat}) \pm 1.5~(\mbox{syst}) \pm 2.7~(\mbox{global})~\mu$b. The total cross section is compared to a perturbative quantum chromodynamics calculation and is consistent within uncertainties. The azimuthal opening angle between muon pairs from $b\bar{b}$ decays and their $p_T$ distributions are compared to distributions generated using {\sc ps pythia 6}, which includes next-to-leading order processes. The azimuthal correlations and pair $p_T$ distribution are not very well described by {\sc pythia} calculations, but are still consistent within uncertainties. Flavor creation and flavor excitation subprocesses are favored over gluon splitting.

We report on charmonium measurements [J/psi(1S), psi'(2S), and chi_c(1P)] in p+p collisions at sqrt(s)=200 GeV. We find that the fraction of J/psi coming from the feed-down decay of psi' and chi_c in the midrapidity region ($|\eta|<0.35$) is 9.6+/-2.4% and 32+/-9%, respectively. We also report new, higher statistics p_T and rapidity dependencies of the J/psi yield via dielectron decay in the same midrapidity range and at forward rapidity (1.2<|eta|<2.4) via dimuon decay. These results are compared with measurements from other experiments and discussed in the context of current charmonium production models.