Charmonium is a valuable probe in heavy-ion collisions to study the properties of the quark gluon plasma, and is also an interesting probe in small collision systems to study cold nuclear matter effects, which are also present in large collision systems. With the recent observations of collective behavior of produced particles in small system collisions, measurements of the modification of charmonium in small systems have become increasingly relevant. We present the results of J/ψ measurements at forward and backward rapidity in various small collision systems, p+p, p+Al, p+Au and 3He+Au, at √sNN =200 GeV. The results are presented in the form of the observable RAB, the nuclear modification factor, a measure of the ratio of the J/ψ invariant yield compared to the scaled yield in p+p collisions. We examine the rapidity, transverse momentum, and collision centrality dependence of nuclear effects on J/ψ production with different projectile sizes p and 3He, and different target sizes Al and Au. The modification is found to be strongly dependent on the target size, but to be very similar for p+Au and 3He+Au. However, for 0%–20% central collisions at backward rapidity, the modification for 3He+Au is found to be smaller than that for p+Au, with a mean fit to the ratio of 0.89±0.03(stat)±0.08(syst), possibly indicating final state effects due to the larger projectile size.

The PHENIX experiment has studied nuclear effects in $p$$+$Al and $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV on charged hadron production at forward rapidity ($1.4<\eta<2.4$, $p$-going direction) and backward rapidity ($-2.2<\eta<-1.2$, $A$-going direction). Such effects are quantified by measuring nuclear modification factors as a function of transverse momentum and pseudorapidity in various collision multiplicity selections. In central $p$$+$Al and $p$$+$Au collisions, a suppression (enhancement) is observed at forward (backward) rapidity compared to the binary scaled yields in $p$+$p$ collisions. The magnitude of enhancement at backward rapidity is larger in $p$$+$Au collisions than in $p$$+$Al collisions, which have a smaller number of participating nucleons. However, the results at forward rapidity show a similar suppression within uncertainties. The results in the integrated centrality are compared with calculations using nuclear parton distribution functions, which show a reasonable agreement at the forward rapidity but fail to describe the backward rapidity enhancement.

A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV, and Cu+Cu collisions sqrt(s_NN) = 22.5, 62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for non-dynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. An analysis of Negative Binomial Distribution fits to the multiplicity distributions demonstrates that the heavy ion data exhibit weak clustering properties.

We report the measurement of direct photons at midrapidity in Au+Au collisions at sqrt{s_NN} = 200 GeV. The direct photon signal was extracted for the transverse-momentum range of 4 GeV/c < p_T < 22 GeV/c, using a statistical method to subtract decay photons from the inclusive-photon sample. The direct-photon nuclear-modification factor R_AA was calculated as a function of p_T for different Au+Au collision centralities using the measured p+p direct-photon spectrum and compared to theoretical predictions. R_AA was found to be consistent with unity for all centralities over the entire measured p_T range. Theoretical models that account for modifications of initial-direct-photon production due to modified-parton-distribution functions in Au and the different isospin composition of the nuclei, predict a modest change of R_AA from unity and are consistent with the data. Models with compensating effects of the quark-gluon plasma on high-energy photons, such as suppression of jet-fragmentation photons and induced-photon bremsstrahlung from partons traversing the medium, are also consistent with this measurement.

During 2015 the Relativistic Heavy Ion Collider (RHIC) provided collisions of transversely polarized protons with Au and Al nuclei for the first time, enabling the exploration of transverse-single-spin asymmetries with heavy nuclei. Large single-spin asymmetries in very forward neutron production have been previously observed in transversely polarized $p$$+$$p$ collisions at RHIC, and the existing theoretical framework that was successful in describing the single-spin asymmetry in $p$$+$$p$ collisions predicts only a moderate atomic-mass-number ($A$) dependence. In contrast, the asymmetries observed at RHIC in $p$$+$$A$ collisions showed a surprisingly strong $A$ dependence in inclusive forward neutron production. The observed asymmetry in $p$$+$Al collisions is much smaller, while the asymmetry in $p$$+$Au collisions is a factor of three larger in absolute value and of opposite sign. The interplay of different neutron production mechanisms is discussed as a possible explanation of the observed $A$ dependence.

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.

Measurements of double-helicity asymmetries for inclusive hadron production in polarized p+p collisions are sensitive to helicity--dependent parton distribution functions, in particular to the gluon helicity distribution, Delta(g). This study focuses on the extraction of the double-helicity asymmetry in eta production: polarized p+p --> eta + X, the eta cross section, and the eta/pi^0 cross section ratio. The cross section and ratio measurements provide essential input for the extraction of fragmentation functions that are needed to access the helicity-dependent parton distribution functions.

Asymmetric nuclear collisions of $p$$+$Al, $p$$+$Au, $d$$+$Au, and $^{3}$He$+$Au at $\sqrt{s_{_{NN}}}=200$ GeV provide an excellent laboratory for understanding particle production, as well as exploring interactions among these particles after their initial creation in the collision. We present measurements of charged hadron production $dN_{\rm ch}/d\eta$ in all such collision systems over a broad pseudorapidity range and as a function of collision multiplicity. A simple wounded quark model is remarkably successful at describing the full data set. We also measure the elliptic flow $v_{2}$ over a similarly broad pseudorapidity range. These measurements provide key constraints on models of particle emission and their translation into flow.

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.