Measurements of transverse-single-spin asymmetries ($A_{N}$) in $p$$+$$p$ collisions at $\sqrt{s}=$62.4 and 200 GeV with the PHENIX detector at RHIC are presented. At midrapidity, $A_{N}$ is measured for neutral pion and eta mesons reconstructed from diphoton decay, and at forward rapidities, neutral pions are measured using both diphotons and electromagnetic clusters. The neutral-pion measurement of $A_{N}$ at midrapidity is consistent with zero with uncertainties a factor of 20 smaller than previous publications, which will lead to improved constraints on the gluon Sivers function. At higher rapidities, where the valence quark distributions are probed, the data exhibit sizable asymmetries. In comparison with previous measurements in this kinematic region, the new data extend the kinematic coverage in $\sqrt{s}$ and $p_T$, and it is found that the asymmetries depend only weakly on $\sqrt{s}$. The origin of the forward $A_{N}$ is presently not understood quantitatively. The extended reach to higher $p_T$ probes the transition between transverse momentum dependent effects at low $p_T$ and multi-parton dynamics at high $p_T$.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured omega meson production via leptonic and hadronic decay channels in p+p, d+Au, Cu+Cu, and Au+Au collisions at sqrt(s_NN) = 200 GeV. The invariant transverse momentum spectra measured in different decay modes give consistent results. Measurements in the hadronic decay channel in Cu+Cu and Au+Au collisions show that omega production has a suppression pattern at high transverse momentum, similar to that of pi^0 and eta in central collisions, but no suppression is observed in peripheral collisions. The nuclear modification factors, R_AA, are consistent in Cu+Cu and Au+Au collisions at similar numbers of participant nucleons.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K^0_S , \omega, \eta prime, and \phi mesons in p + p collisions at = 200 GeV. Measurements \omega and \phi production in different decay channels give consistent results. New results for the \phi are in agreement with previously published data and extend the measured pT coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-pT and characterizing the low-pT regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

Back-to-back hadron pair yields in d+Au and p+p collisions at sqrt(s_NN)=200 GeV were measured with the PHENIX detector at the Relativistic Heavy Ion Collider. Rapidity separated hadron pairs were detected with the trigger hadron at pseudorapidity |eta|<0.35 and the associated hadron at forward rapidity (deuteron direction, 3.0<eta<3.8). Pairs were also detected with both hadrons measured at forward rapidity; in this case the yield of back-to-back hadron pairs in d+Au collisions with small impact parameters is observed to be suppressed by a factor of 10 relative to p+p collisions. The kinematics of these pairs is expected to probe partons in the Au nucleus with low fraction x of the nucleon momenta, where the gluon densities rise sharply. The observed suppression as a function of nuclear thickness, p_T, and eta points to cold nuclear matter effects arising at high parton densities.

The energy dependence of the single-transverse-spin asymmetry, A_N, and the cross section for neutron production at very forward angles were measured in the PHENIX experiment at RHIC for polarized p+p collisions at sqrt(s)=200 GeV. The neutrons were observed in forward detectors covering an angular range of up to 2.2 mrad. We report results for neutrons with momentum fraction of x_F=0.45 to 1.0. The energy dependence of the measured cross sections were consistent with x_F scaling, compared to measurements by an ISR experiment which measured neutron production in unpolarized p+p collisions at sqrt(s)=30.6--62.7 GeV. The cross sections for large x_F neutron production for p+p collisions, as well as those in e+p collisions measured at HERA, are described by a pion exchange mechanism. The observed forward neutron asymmetries were large, reaching A_N=-0.08+/-0.02 for x_F=0.8; the measured backward asymmetries, for negative x_F, were consistent with zero. The observed asymmetry for forward neutron production is discussed within the pion exchange framework, with interference between the spin-flip amplitude due to the pion exchange and nonflip amplitudes from all Reggeon exchanges. Within the pion exchange description, the measured neutron asymmetry is sensitive to the contribution of other Reggeon exchanges even for small amplitudes.

Direct photons have been measured in sqrt(s_NN)=200 GeV d+Au collisions at midrapidity. A wide p_T range is covered by measurements of nearly-real virtual photons (1<p_T<6 GeV/c) and real photons (5<p_T<16 GeV/c). The invariant yield of the direct photons in d+Au collisions over the scaled p+p cross section is consistent with unity. Theoretical calculations assuming standard cold nuclear matter effects describe the data well for the entire p_T range. This indicates that the large enhancement of direct photons observed in Au+Au collisions for 1.0<p_T<2.5 GeV/c is due to a source other than the initial-state nuclear effects.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured 2nd and 3rd order Fourier coefficients of the azimuthal distributions of direct photons emitted at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for various collision centralities. Combining two different analysis techniques, results were obtained in the transverse momentum range of $0.4<p_{T}<4.0$ GeV/$c$. At low $p_T$ the second-order coefficients, $v_2$, are similar to the ones observed in hadrons. Third order coefficients, $v_3$, are nonzero and almost independent of centrality. These new results on $v_2$ and $v_3$, combined with previously published results on yields, are compared to model calculations that provide yields and asymmetries in the same framework. Those models are challenged to explain simultaneously the observed large yield and large azimuthal anisotropies.

Differential measurements of the elliptic (v_2) and hexadecapole (v_4) Fourier flow coefficients are reported for charged hadrons as a function of transverse momentum (p_T) and collision centrality or the number of participant nucleons (N_part) for Au+Au collisions at sqrt(s_NN)=200 GeV. The v_{2,4} measurements at pseudorapidity |\eta|<=0.35 obtained with four separate reaction plane detectors positioned in the range 1.0<|\eta|<3.9 show good agreement, indicating the absence of significant \eta-dependent nonflow perturbations. Sizable values for v_4(p_T) are observed with a ratio v_4(p_T,N_part)/v_2^2(p_T,N_part)~0.8 for 50<N_part<200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N_part>200 this ratio increases up to 1.7 in the most central collisions.

We present measurements of electrons and positrons from the semileptonic decays of heavy-flavor hadrons at midrapidity ($|y|<$ 0.35) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. The data were collected in 2010 by the PHENIX experiment that included the new hadron-blind detector. The invariant yield of electrons from heavy-flavor decays is measured as a function of transverse momentum in the range $1<p_T^e<5$ GeV/$c$. The invariant yield per binary collision is slightly enhanced above the $p$$+$$p$ reference in Au$+$Au 0%--20%, 20%--40% and 40%--60% centralities at a comparable level. This may be a result of the interplay between initial-state Cronin effects, final-state flow, and energy loss for heavy-quark production at this low beam energy. The $v_2$ of electrons from heavy-flavor decays is nonzero when averaged between $1.3<p_T^e<2.5$ GeV/$c$ from $0<{\rm centrality}<40$% collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. For 20%--40% centrality collisions, the $v_2$ at $\sqrt{s_{_{NN}}}=62.4$ GeV is smaller than that for heavy flavor decays at $\sqrt{s_{_{NN}}}=200$ GeV. The $v_2$ of the electrons from heavy-flavor decay at the lower beam energy is also smaller than $v_2$ for pions. Both results indicate that the heavy-quarks interact with the medium formed in these collisions, but they may not be at the same level of thermalization with the medium as observed at $\sqrt{s_{_{NN}}}=200$ GeV.

We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.