The PHENIX experiment has measured the suppression of semi-inclusive single high transverse momentum pi^0's in Au+Au collisions at sqrt(s_NN) = 200 GeV. The present understanding of this suppression is in terms of energy-loss of the parent (fragmenting) parton in a dense color-charge medium. We have performed a quantitative comparison between various parton energy-loss models and our experimental data. The statistical point-to-point uncorrelated as well as correlated systematic uncertainties are taken into account in the comparison. We detail this methodology and the resulting constraint on the model parameters, such as the initial color-charge density dN^g/dy, the medium transport coefficient <q^hat>, or the initial energy-loss parameter epsilon_0. We find that high transverse momentum pi^0 suppression in Au+Au collisions has sufficient precision to constrain these model dependent parameters at the +/1 20%-25% (one standard deviation) level. These constraints include only the experimental uncertainties, and further studies are needed to compute the corresponding theoretical uncertainties.

A three-dimensional (3D) correlation function obtained from mid-rapidity, low pT pion pairs in central Au+Au collisions at sqrt(s_NN)=200 GeV is studied. The extracted model-independent source function indicates a long range tail in the directions of the pion pair transverse momentum (out) and the beam (long). Model comparisons to these distensions indicate a proper breakup time \tau_0 ~ 9 fm/c and a mean proper emission duration \Delta\tau ~ 2 fm/c, leading to sizable emission time differences (<|\Delta \tau_LCM |> ~ 12 fm/c), partly due to resonance decays. They also suggest an outside-in 'burning' of the emission source reminiscent of many hydrodynamical models.

Yields for J/psi production in Cu+Cu collisions at sqrt (s_NN)= 200 GeV have been measured by the PHENIX experiment over the rapidity range |y| < 2.2 at transverse momenta from 0 to beyond 5 GeV/c. The invariant yield is obtained as a function of rapidity, transverse momentum and collision centrality, and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data provide greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, providing a key constraint that is needed for disentangling cold and hot nuclear matter effects.

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 differential cross sections for the process $\gamma \gamma \to \pi^0 \pi^0$ have been measured in the kinematic range 0.6 GeV $< W < 4.1$ GeV, $|\cos \theta^*|<0.8$ in energy and pion 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. The differential cross sections are fitted in the energy region 1.7 GeV $< W <$ 2.5 GeV to confirm the two-photon production of two pions in the G wave. In the higher energy region, we observe production of the $\chi_{c0}$ charmonium state and obtain the product of its two-photon decay width and branching fraction to $\pi^0\pi^0$. We also compare the observed angular dependence and ratios of cross sections for neutral-pion and charged-pion pair production to QCD models. The energy and angular dependence above 3.1 GeV are compatible with those measured in the $\pi^+\pi^-$ channel, and in addition we find that the cross section ratio, $\sigma(\pi^0\pi^0)/\sigma(\pi^+\pi^-)$, is $0.32 \pm 0.03 \pm 0.05$ on average in the 3.1-4.1 GeV region.

The PHENIX Collaboration has measured the ratio of the yields of $\psi(2S)$ to $\psi(1S)$ mesons produced in $p$$+$$p$, $p$$+$Al, $p$$+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV over the forward and backward rapidity intervals $1.2<|y|<2.2$. We find that the ratio in $p$$+$$p$ collisions is consistent with measurements at other collision energies. In collisions with nuclei, we find that in the forward ($p$-going or $^{3}$He-going) direction, the relative yield of $\psi(2S)$ mesons to $\psi(1S)$ mesons is consistent with the value measured in \pp collisions. However, in the backward (nucleus-going) direction, the $\psi(2S)$ is preferentially suppressed by a factor of $\sim$2. This suppression is attributed in some models to breakup of the weakly-bound $\psi(2S)$ through final state interactions with comoving particles, which have a higher density in the nucleus-going direction. These breakup effects may compete with color screening in a deconfined quark-gluon plasma to produce sequential suppression of excited quarkonia states.

We report on the nuclear dependence of transverse single-spin asymmetries (TSSAs) in the production of positively-charged hadrons in polarized $p^{\uparrow}+p$, $p^{\uparrow}+$Al and $p^{\uparrow}+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The measurements have been performed at forward rapidity ($1.4<\eta<2.4$) over the range of $1.8<p_{T}<7.0$ GeV$/c$ and $0.1<x_{F}<0.2$. We observed a positive asymmetry $A_{N}$ for positively-charged hadrons in \polpp collisions, and a significantly reduced asymmetry in $p^{\uparrow}$+$A$ collisions. These results reveal a nuclear dependence of charged hadron $A_N$ in a regime where perturbative techniques are relevant. These results provide new opportunities to use \polpA collisions as a tool to investigate the rich phenomena behind TSSAs in hadronic collisions and to use TSSA as a new handle in studying small-system collisions.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section of $\phi$(1020) meson production at forward rapidity in $p$$+$$p$ collisions at $\sqrt{s}=$510 GeV via the dimuon decay channel. The integrated cross section in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $2<p_T<7$ GeV/$c$ is $\sigma_\phi=2.79 \pm 0.20\,{\rm (stat)} \pm 0.17\,{\rm (syst)} \pm 0.34\, {\rm (norm)} \times 10^{-2}$~mb. The energy dependence of $\sigma_\phi$ ($1.2<|y|<2.2$; $2<p_T<5$ GeV/$c$) is studied using the PHENIX measurements at $\sqrt{s}=$200 and 510 GeV and the Large-Hadron-Collider measurements at $\sqrt{s}=$2.76 and 7 TeV. The experimental results are compared to various event generator predictions (pythia6, pythia8, phojet, ampt, epos3, and epos-lhc).

Charged-pion-interferometry measurements were made with respect to the 2$^{\rm nd}$- and 3$^{\rm rd}$-order event plane for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. A strong azimuthal-angle dependence of the extracted Gaussian-source radii was observed with respect to both the 2$^{\rm nd}$- and 3$^{\rm rd}$-order event planes. The results for the 2$^{\rm nd}$-order dependence indicate that the initial eccentricity is reduced during the medium evolution, but not reversed in the final state, which is consistent with previous results. In contrast, the results for the 3$^{\rm rd}$-order dependence indicate that the initial triangular shape is significantly reduced and potentially reversed by the end of the medium evolution, and that the 3$^{\rm rd}$-order oscillations are largely dominated by the dynamical effects from triangular flow.