We present measurements of elliptic and triangular azimuthal anisotropy of charged particles detected at forward rapidity $1<|\eta|<3$ in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, as a function of centrality. The multiparticle cumulant technique is used to obtain the elliptic flow coefficients $v_2\{2\}$, $v_2\{4\}$, $v_2\{6\}$, and $v_2\{8\}$, and triangular flow coefficients $v_3\{2\}$ and $v_3\{4\}$. Using the small-variance limit, we estimate the mean and variance of the event-by-event $v_2$ distribution from $v_2\{2\}$ and $v_2\{4\}$. In a complementary analysis, we also use a folding procedure to study the distributions of $v_2$ and $v_3$ directly, extracting both the mean and variance. Implications for initial geometrical fluctuations and their translation into the final state momentum distributions are discussed.

Measurements of anisotropic flow Fourier coefficients ($v_n$) for inclusive charged particles and identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ produced at midrapidity in Cu+Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV are presented. The data were collected in 2012 by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). The particle azimuthal distributions with respect to different order symmetry planes $\Psi_n$, for $n$~=~1, 2, and 3 are studied as a function of transverse momentum $p_T$ over a broad range of collisions centralities. Mass ordering, as expected from hydrodynamic flow, is observed for all three harmonics. The charged-particle results are compared to hydrodynamical and transport model calculations. We also compare these Cu$+$Au results with those in Cu$+$Cu and Au$+$Au collisions at the same $\sqrt{s_{_{NN}}}$, and find that the $v_2$ and $v_3$, as a function of transverse momentum, follow a common scaling with $1/(\varepsilon_n N_{\rm part}^{1/3})$.

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 present the first measurements of long-range angular correlations and the transverse momentum dependence of elliptic flow $v_2$ in high-multiplicity $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. A comparison of these results with previous measurements in high-multiplicity $d$$+$Au and $^3{\rm He}$$+$Au collisions demonstrates a relation between $v_2$ and the initial collision eccentricity $\varepsilon_2$, suggesting that the observed momentum-space azimuthal anisotropies in these small systems have a collective origin and reflect the initial geometry. Good agreement is observed between the measured $v_2$ and hydrodynamic calculations for all systems, and an argument disfavoring theoretical explanations based on momentum-space domain correlations is presented. The set of measurements presented here allows us to leverage the distinct intrinsic geometry of each of these systems to distinguish between different theoretical descriptions of the long-range correlations observed in small collision systems.

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

PHENIX measurements are presented for the cross section and double-helicity asymmetry ($A_{LL}$) in inclusive $\pi^0$ production at midrapidity from $p$$+$$p$ collisions at $\sqrt{s}=510$~GeV from data taken in 2012 and 2013 at the Relativistic Heavy Ion Collider. The next-to-leading-order perturbative-quantum-chromodynamics theory calculation is in excellent agreement with the presented cross section results. The calculation utilized parton-to-pion fragmentation functions from the recent DSS14 global analysis, which prefer a smaller gluon-to-pion fragmentation function. The $\pi^{0}A_{LL}$ results follow an increasingly positive asymmetry trend with $p_T$ and $\sqrt{s}$ with respect to the predictions and are in excellent agreement with the latest global analysis results. This analysis incorporated earlier results on $\pi^0$ and jet $A_{LL}$, and suggested a positive contribution of gluon polarization to the spin of the proton $\Delta G$ for the gluon momentum fraction range $x>0.05$. The data presented here extend to a currently unexplored region, down to $x\sim0.01$, and thus provide additional constraints on the value of $\Delta G$. The results confirm the evidence for nonzero $\Delta G$ using a different production channel in a complementary kinematic region.

PHENIX reports differential cross sections of $\mu\mu$ pairs from semileptonic heavy-flavor decays and the Drell-Yan production mechanism measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV at forward and backward rapidity ($1.2<|\eta|<2.2$). The $\mu\mu$ pairs from $c\bar{c}$, $b\bar{b}$, and Drell-Yan are separated using a simultaneous fit to unlike- and like-sign muon pair spectra in mass and $p_T$. The azimuthal opening angle correlation between the muons from $c\bar{c}$ and $b\bar{b}$ decays and the pair-$p_T$ distributions are compared to distributions generated using {\sc pythia} and {\sc powheg} models, which both include next-to-leading order processes. The measured distributions for pairs from $c\bar{c}$ are consistent with {\sc pythia} calculations. The $c\bar{c}$ data presents narrower azimuthal correlations and softer $p_T$ distributions compared to distributions generated from {\sc powheg}. The $b\bar{b}$ data are well described by both models. The extrapolated total cross section for bottom production is $3.75{\pm}0.24({\rm stat}){\pm}^{0.35}_{0.50}({\rm syst}){\pm}0.45({\rm global})$[$\mu$b], which is consistent with previous measurements at the Relativistic Heavy Ion Collider in the same system at the same collision energy, and is approximately a factor of two higher than the central value calculated with theoretical models. The measured Drell-Yan cross section is in good agreement with next-to-leading-order quantum-chromodynamics calculations.

The invariant yields for $J/\psi$ production at forward rapidity $(1.2<|y|<2.2)$ in U$+$U collisions at $\sqrt{s_{_{NN}}}$=193 GeV have been measured as a function of collision centrality. The invariant yields and nuclear-modification factor $R_{AA}$ are presented and compared with those from Au$+$Au collisions in the same rapidity range. Additionally, the direct ratio of the invariant yields from U$+$U and Au$+$Au collisions within the same centrality class is presented, and used to investigate the role of $c\bar{c}$ coalescence. Two different parameterizations of the deformed Woods-Saxon distribution were used in Glauber calculations to determine the values of the number of nucleon-nucleon collisions in each centrality class, $N_{\rm coll}$, and these were found to give significantly different $N_{\rm coll}$ values. Results using $N_{\rm coll}$ values from both deformed Woods-Saxon distributions are presented. The measured ratios show that the $J/\psi$ suppression, relative to binary collision scaling, is similar in U$+$U and Au$+$Au for peripheral and midcentral collisions, but that $J/\psi$ show less suppression for the most central U$+$U collisions. The results are consistent with a picture in which, for central collisions, increase in the $J/\psi$ yield due to $c\bar{c}$ coalescence becomes more important than the decrease in yield due to increased energy density. For midcentral collisions, the conclusions about the balance between $c\bar{c}$ coalescence and suppression depend on which deformed Woods-Saxon distribution is used to determine $N_{\rm coll}$.

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.

We present a systematic study of charged pion and kaon interferometry in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. The kaon mean source radii are found to be larger than pion radii in the outward and longitudinal directions for the same transverse mass; this difference increases for more central collisions. The azimuthal-angle dependence of the radii was measured with respect to the second-order event plane and similar oscillations of the source radii were found for pions and kaons. Hydrodynamic models qualitatively describe the similar oscillations of the mean source radii for pions and kaons, but they do not fully describe the transverse-mass dependence of the oscillations.

We have measured direct photons for $p_T<5~$GeV/$c$ in minimum bias and 0\%--40\% most central events at midrapidity for Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The $e^{+}e^{-}$ contribution from quasi-real direct virtual photons has been determined as an excess over the known hadronic contributions in the $e^{+}e^{-}$ mass distribution. A clear enhancement of photons over the binary scaled $p$$+$$p$ fit is observed for $p_T<4$ GeV/$c$ in Cu$+$Cu data. The $p_T$ spectra are consistent with the Au$+$Au data covering a similar number of participants. The inverse slopes of the exponential fits to the excess after subtraction of the $p$$+$$p$ baseline are 285$\pm$53(stat)$\pm$57(syst)~MeV/$c$ and 333$\pm$72(stat)$\pm$45(syst)~MeV/$c$ for minimum bias and 0\%--40\% most central events, respectively. The rapidity density, $dN/dy$, of photons demonstrates the same power law as a function of $dN_{\rm ch}/d\eta$ observed in Au$+$Au at the same collision energy.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. Correlations of charged hadrons of $0.7<p_T<10$ GeV/$c$ with $\pi^0$ mesons of $4<p_T<15$ GeV/$c$ or isolated direct photons of $7<p_T<15$ GeV/$c$ are used to study nonperturbative effects generated by initial-state partonic transverse momentum and final-state transverse momentum from fragmentation. The nonperturbative behavior is characterized by measuring the out-of-plane transverse momentum component $p_{\rm out}$ perpendicular to the axis of the trigger particle, which is the high-$p_T$ direct photon or $\pi^0$. Nonperturbative evolution effects are extracted from Gaussian fits to the away-side inclusive-charged-hadron yields for different trigger-particle transverse momenta ($p_T^{\rm trig}$). The Gaussian widths and root mean square of $p_{\rm out}$ are reported as a function of the interaction hard scale $p_T^{\rm trig}$ to investigate possible transverse-momentum-dependent evolution differences between the $\pi^0$-h$^\pm$ and direct photon-h$^\pm$ correlations and factorization breaking effects. The widths are found to decrease with $p_T^{\rm trig}$, which indicates that the Collins-Soper-Sterman soft factor is not driving the evolution with the hard scale in nearly back-to-back dihadron and direct photon-hadron production in $p$$+$$p$ collisions. This behavior is in contrast to Drell-Yan and semi-inclusive deep-inelastic scattering measurements.

We report the first measurement of the full angular distribution for inclusive $J/\psi\rightarrow\mu^{+}\mu^{-}$ decays in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. The measurements are made for $J/\psi$ transverse momentum $2<p_{T}<10$ GeV/$c$ and rapidity $1.2<y<2.2$ in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient $\lambda_{\theta}$ is strongly negative at low $p_{T}$ and becomes close to zero at high $p_{T}$, while the azimuthal coefficient $\lambda_{\phi}$ is close to zero at low $p_{T}$, and becomes slightly negative at higher $p_{T}$. The frame-independent coefficient $\tilde{\lambda}$ is strongly negative at all $p_{T}$ in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The correlations are sensitive to nonperturbative initial-state and final-state transverse momentum $k_T$ and $j_T$ in the azimuthal nearly back-to-back region $\Delta\phi\sim\pi$. In this region, transverse-momentum-dependent evolution can be studied when several different hard scales are measured. To have sensitivity to small transverse momentum scales, nonperturbative momentum widths of $p_{\rm out}$, the out-of-plane transverse momentum component perpendicular to the trigger particle, are measured. These widths are used to investigate possible effects from transverse-momentum-dependent factorization breaking. When accounting for the longitudinal momentum fraction of the away-side hadron with respect to the near-side trigger particle, the widths are found to increase with the hard scale; this is qualitatively similar to the observed behavior in Drell-Yan and semi-inclusive deep-inelastic scattering interactions. The momentum widths are also studied as a function of center-of-mass energy by comparing to previous measurements at $\sqrt{s}=510$ GeV. The nonperturbative jet widths also appear to increase with $\sqrt{s}$ at a similar $x_T$, which is qualitatively consistent to similar measurements in Drell-Yan interactions. To quantify the magnitude of any transverse-momentum-dependent factorization breaking effects, calculations will need to be performed to compare to these measurements.

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.

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.

The cross section and transverse single-spin asymmetries of $\mu^{-}$ and $\mu^{+}$ from open heavy-flavor decays in polarized $p$+$p$ collisions at $\sqrt{s}=200$ GeV were measured by the PHENIX experiment during 2012 at the Relativistic Heavy Ion Collider. Because heavy-flavor production is dominated by gluon-gluon interactions at $\sqrt{s}=200$ GeV, these measurements offer a unique opportunity to obtain information on the trigluon correlation functions. The measurements are performed at forward and backward rapidity ($1.4<|y|<2.0$) over the transverse momentum range of $1.25<p_T<7$ GeV/$c$ for the cross section and $1.25<p_T<5$ GeV/$c$ for the asymmetry measurements. The obtained cross section is compared to a fixed-order-plus-next-to-leading-log perturbative-quantum-chromodynamics calculation. The asymmetry results are consistent with zero within uncertainties, and a model calculation based on twist-3 three-gluon correlations agrees with the data.

Measurements of the fractional momentum loss ($S_{\rm loss}\equiv{\delta}p_T/p_T$) of high-transverse-momentum-identified hadrons in heavy ion collisions are presented. Using $\pi^0$ in Au$+$Au and Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=62.4$ and 200 GeV measured by the PHENIX experiment at the Relativistic Heavy Ion Collider and and charged hadrons in Pb$+$Pb collisions measured by the ALICE experiment at the Large Hadron Collider, we studied the scaling properties of $S_{\rm loss}$ as a function of a number of variables: the number of participants, $N_{\rm part}$, the number of quark participants, $N_{\rm qp}$, the charged-particle density, $dN_{\rm ch}/d\eta$, and the Bjorken energy density times the equilibration time, $\varepsilon_{\rm Bj}\tau_{0}$. We find that the $p_T$ where $S_{\rm loss}$ has its maximum, varies both with centrality and collision energy. Above the maximum, $S_{\rm loss}$ tends to follow a power-law function with all four scaling variables. The data at $\sqrt{s_{_{NN}}}$=200 GeV and 2.76 TeV, for sufficiently high particle densities, have a common scaling of $S_{\rm loss}$ with $dN_{\rm ch}/d\eta$ and $\varepsilon_{\rm Bj}\tau_{0}$, lending insight on the physics of parton energy loss.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured $\phi$ meson production and its nuclear modification in asymmetric Cu$+$Au heavy-ion collisions at $\sqrt{s_{NN}}=200$ GeV at both forward Cu-going direction ($1.2<y<2.2$) and backward Au-going direction ($-2.2<y<-1.2$), rapidities. The measurements are performed via the dimuon decay channel and reported as a function of the number of participating nucleons, rapidity, and transverse momentum. In the most central events, 0\%--20\% centrality, the $\phi$ meson yield integrated over $1<p_T<5$ GeV/$c$ prefers a smaller value, which means a larger nuclear modification, in the Cu-going direction compared to the Au-going direction. Additionally, the nuclear-modification factor in Cu$+$Au collisions averaged over all centrality is measured to be similar to the previous PHENIX result in $d$$+$Au collisions for these rapidities.

We report the transverse single-spin asymmetries of $J/\psi$ production at forward and backward rapidity, $1.2<|y|<2.2$, as a function of $J/\psi$ transverse momentum ($p_T$) and Feynman-$x$ ($x_F$). The data analyzed were recorded by the PHENIX experiment at the Relativistic Heavy Ion Collider in 2015 from $p$$+$$p$, $p$$+$Al, and $p$$+$Au collisions with transversely polarized proton beams at $\sqrt{s_{_{NN}}}=200$ GeV. At this collision energy, single-spin asymmetries for heavy-flavor particle production of $p$$+$$p$ collisions provide access to the spin-dependent gluon distribution and higher-twist correlation functions inside the nucleon, such as the gluon Qiu-Sterman and trigluon correlation functions. Proton+nucleus collisions offer an excellent opportunity to study nuclear effects on the correlation functions. The data indicate negative asymmetries at the two-standard-deviation level in the $p$$+$Au data for $p_T<2$ GeV/$c$ at both forward and backward rapidity, while in $p$$+$$p$ and $p$$+$Al collisions the asymmetries are consistent with zero within the range of experimental uncertainties.

We report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $p$+$p$ collisions at $\sqrt{s_{NN}}=200$~GeV. The $e^+e^-$ pair yield from $b\bar{b}$ and $c\bar{c}$ is separated by exploiting a double differential fit done simultaneously in dielectron invariant mass and $p_T$. We used three different event generators, {\sc pythia}, {\sc mc@nlo}, and {\sc powheg}, to simulate the $e^+e^-$ spectra from $c\bar{c}$ and $b\bar{b}$ production. The data can be well described by all three generators within the detector acceptance. However, when using the generators to extrapolate to $4\pi$, significant differences are observed for the total cross section. These difference are less pronounced for $b\bar{b}$ than for $c\bar{c}$. The same model dependence was observed in already published $d$+$A$ data. The $p$+$p$ data are also directly compared with $d$+$A$ data in mass and $p_T$, and within the statistical accuracy no nuclear modification is seen.

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.

Measurements of the differential production of electrons from open-heavy-flavor hadrons with charm- and bottom-quark content in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV are presented. The measurements proceed through displaced-vertex analyses of electron tracks from the semileptonic decay of charm and bottom hadrons using the PHENIX silicon-vertex detector. The relative contribution of electrons from bottom decays to inclusive heavy-flavor-electron production is found to be consistent with fixed-order-plus-next-to-leading-log perturbative-QCD calculations within experimental and theoretical uncertainties. These new measurements in $p$$+$$p$ collisions provide a precision baseline for comparable forthcoming measurements in A$+$A collisions.

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.

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.

We present measurements of the transverse-momentum dependence of elliptic flow $v_2$ for identified pions and (anti)protons at midrapidity ($|\eta|<0.35$), in 0%--5% central $p$$+$Au and $^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. When taken together with previously published measurements in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, the results cover a broad range of small-collision-system multiplicities and intrinsic initial geometries. We observe a clear mass-dependent splitting of $v_2(p_{T})$ in $d$$+$Au and $^3$He$+$Au collisions, just as in large nucleus-nucleus ($A$$+$$A$) collisions, and a smaller splitting in $p$$+$Au collisions. Both hydrodynamic and transport model calculations successfully describe the data at low $p_T$ ($< 1.5$ GeV/$c$), but fail to describe various features at higher $p_T$. In all systems, the $v_2$ values follow an approximate quark-number scaling as a function of the hadron transverse kinetic energy per constituent quark($KE_T/n_q$), which was also seen previously in $A$$+$$A$ collisions.

We report on $J/\psi$ production from asymmetric Cu+Au heavy-ion collisions at $\sqrt{s_{_{NN}}}$=200 GeV at the Relativistic Heavy Ion Collider at both forward (Cu-going direction) and backward (Au-going direction) rapidities. The nuclear modification of $J/\psi$ yields in Cu$+$Au collisions in the Au-going direction is found to be comparable to that in Au$+$Au collisions when plotted as a function of the number of participating nucleons. In the Cu-going direction, $J/\psi$ production shows a stronger suppression. This difference is comparable in magnitude and has the same sign as the difference expected from shadowing effects due to stronger low-$x$ gluon suppression in the larger Au nucleus. The relative suppression is opposite to that expected from hot nuclear matter dissociation, since a higher energy density is expected in the Au-going direction.

The PHENIX collaboration has measured high-$p_T$ dihadron correlations in $p$$+$$p$, $p$$+$Al, and $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The correlations arise from inter- and intra-jet correlations and thus have sensitivity to nonperturbative effects in both the initial and final states. The distributions of $p_{\rm out}$, the transverse momentum component of the associated hadron perpendicular to the trigger hadron, are sensitive to initial and final state transverse momenta. These distributions are measured multi-differentially as a function of $x_E$, the longitudinal momentum fraction of the associated hadron with respect to the trigger hadron. The near-side $p_{\rm out}$ widths, sensitive to fragmentation transverse momentum, show no significant broadening between $p$$+$Au, $p$$+$Al, and $p$$+$$p$. The away-side nonperturbative $p_{\rm out}$ widths are found to be broadened in $p$$+$Au when compared to $p$$+$$p$; however, there is no significant broadening in $p$$+$Al compared to $p$$+$$p$ collisions. The data also suggest that the away-side $p_{\rm out}$ broadening is a function of $N_{\rm coll}$, the number of binary nucleon-nucleon collisions, in the interaction. The potential implications of these results with regard to initial and final state transverse momentum broadening and energy loss of partons in a nucleus, among other nuclear effects, are discussed.

We present the first measurement of elliptic ($v_2$) and triangular ($v_3$) flow in high-multiplicity $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Two-particle correlations, where the particles have a large separation in pseudorapidity, are compared in $^{3}$He$+$Au and in $p$$+$$p$ collisions and indicate that collective effects dominate the second and third Fourier components for the correlations observed in the $^{3}$He$+$Au system. The collective behavior is quantified in terms of elliptic $v_2$ and triangular $v_3$ anisotropy coefficients measured with respect to their corresponding event planes. The $v_2$ values are comparable to those previously measured in $d$$+$Au collisions at the same nucleon-nucleon center-of-mass energy. Comparison with various theoretical predictions are made, including to models where the hot spots created by the impact of the three $^{3}$He nucleons on the Au nucleus expand hydrodynamically to generate the triangular flow. The agreement of these models with data may indicate the formation of low-viscosity quark-gluon plasma even in these small collision systems.

The fraction of $J/\psi$ mesons which come from B-meson decay, $\textrm{F}_{B{\rightarrow}J/\psi}$, is measured for J/$\psi$ rapidity \mbox{$1.2<|y|<2.2$} and $p_T>0$ in $p$+$p$ and Cu+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV with the PHENIX detector. The extracted fraction is $\textrm{F}_{B{\rightarrow}J/\psi}$ = 0.025 $\pm$ 0.006(stat) $\pm$ 0.010(syst) for $p$+$p$ collisions. For Cu+Au collisions, $\textrm{F}_{B{\rightarrow}J/\psi}$ is 0.094 $\pm$ 0.028(stat) $\pm$ 0.037(syst) in the Au-going direction ($-2.2<y<-1.2$) and 0.089 $\pm$ 0.026(stat) $\pm$ 0.040(syst) in the Cu-going direction ($1.2<y<2.2$). The nuclear modification factor, $R_{\rm CuAu}$, of B mesons in Cu+Au collisions is consistent with binary scaling of measured yields in $p$+$p$ at both forward and backward rapidity.

We have measured the cross section and single spin asymmetries from forward $W^{\pm}\rightarrow\mu^{\pm}\nu$ production in longitudinally polarized $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV using the PHENIX detector at the Relativistic Heavy Ion Collider. The cross sections are consistent with previous measurements at this collision energy, while the most forward and backward longitudinal single spin asymmetries provide new insights into the sea quark helicities in the proton. The charge of the W bosons provides a natural flavor separation of the participating partons.

The PHENIX experiment has measured the spin alignment for inclusive $J/\psi\rightarrow e^{+}e^{-}$ decays in $p$+$p$ collisions at $\sqrt{s}=510$ GeV at midrapidity. The angular distributions have been measured in three different polarization frames, and the three decay angular coefficients have been extracted in a full two-dimensional analysis. Previously, PHENIX saw large longitudinal net polarization at forward rapidity at the same collision energy. This analysis at midrapidity, complementary to the previous PHENIX results, sees no sizable polarization in the measured transverse momentum range of $0.0<p_T<10.0$ GeV/$c$. The results are consistent with a previous one-dimensional analysis at midrapidity at $\sqrt{s}=200$ GeV. The transverse-momentum-dependent cross section for midrapidity $J/\psi$ production has additionally been measured, and after comparison to world data we find a simple logarithmic dependence of the cross section on $\sqrt{s}$.

We present measurements of two-particle angular correlations between high-transverse-momentum ($2<p_T<11$ GeV/$c$) $\pi^0$ observed at midrapidity ($|\eta|<0.35$) and particles produced either at forward ($3.1<\eta<3.9$) or backward ($-3.7<\eta<-3.1$) rapidity in $d$$+$Au and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The azimuthal angle correlations for particle pairs with this large rapidity gap in the Au-going direction exhibit a ridge-like structure that persists up to $p_T{\approx}6$ GeV/$c$ and which strongly depends on collision centrality, which is a similar characteristic to the hydrodynamical particle flow in A+A collisions. The ridge-like structure is absent in the $d$-going direction as well as in $p$$+$$p$ collisions, in the transverse-momentum range studied. The results indicate that the ridge-like structure is shifted in the Au-going direction toward more central collisions, similar to the charged-particle pseudorapidity distributions.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section, mean transverse momentum, mean transverse momentum squared of inclusive $J/\psi$ and cross-section ratio of $\psi(2S)$ to $J/\psi$ at forward rapidity in \pp collisions at \sqrts = 510 GeV via the dimuon decay channel. Comparison is made to inclusive $J/\psi$ cross sections measured at \sqrts = 200 GeV and 2.76--13 TeV. The result is also compared to leading-order nonrelativistic QCD calculations coupled to a color-glass-condensate description of the low-$x$ gluons in the proton at low transverse momentum ($p_T$) and to next-to-leading order nonrelativistic QCD calculations for the rest of the $p_T$ range. These calculations overestimate the data at low $p_T$. While consistent with the data within uncertainties above $\approx3$ GeV/$c$, the calculations are systematically below the data. The total cross section times the branching ratio is BR $d\sigma^{J/\psi}_{pp}/dy (1.2<|y|<2.2, 0<p_T<10~\mbox{GeV/$c$}) =$ 54.3 $\pm$ 0.5 (stat) $\pm$ 5.5 (syst) nb.

We report the first measurement of the fraction of $J/\psi$ mesons coming from $B$-meson decay ($F_{B{\rightarrow}J/\psi}$) in $p$+$p$ collisions at $\sqrt{s}=$ 510 GeV. The measurement is performed using the forward silicon vertex detector and central vertex detector at PHENIX, which provide precise tracking and distance-of-closest-approach determinations, enabling the statistical separation of $J/\psi$ due to $B$-meson decays from prompt $J/\psi$. The measured value of $F_{B{\rightarrow}J/\psi}$ is 8.1\%$\pm$2.3\% (stat)$\pm$1.9\% (syst) for $J/\psi$ with transverse momenta $0<p_T<5$ GeV/$c$ and rapidity $1.2<|y|<2.2$. The measured fraction $F_{B{\rightarrow}J/\psi}$ at PHENIX is compared to values measured by other experiments at higher center of mass energies and to fixed-order-next-to-leading-logarithm and color-evaporation-model predictions. The $b\bar{b}$ cross section per unit rapidity ($d\sigma/dy(pp{\rightarrow}b\bar{b})$) extracted from the obtained $F_{B{\rightarrow}J/\psi}$ and the PHENIX inclusive $J/\psi$ cross section measured at 200 GeV scaled with color-evaporation-model calculations, at the mean $B$ hadron rapidity $y={\pm}1.7$ in 510 GeV $p$$+$$p$ collisions, is $3.63^{+1.92}_{-1.70}\mu$b, and it is consistent with the fixed-order-next-to-leading-logarithm calculations.

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.

Presented are the first measurements of the transverse single-spin asymmetries ($A_N$) for neutral pions and eta mesons in $p$+Au and $p$+Al collisions at $\sqrt{s_{_{NN}}}=200$ GeV in the pseudorapidity range $|\eta|<$0.35 with the PHENIX detector at the Relativistic Heavy Ion Collider. The asymmetries are consistent with zero, similar to those for midrapidity neutral pions and eta mesons produced in $p$+$p$ collisions. These measurements show no evidence of additional effects that could potentially arise from the more complex partonic environment present in proton-nucleus collisions.

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.

Recently, the PHENIX Collaboration has published second- and third-harmonic Fourier coefficients $v_2$ and $v_3$ for midrapidity ($|\eta|<0.35$) charged hadrons in 0%--5% central $p$$+$Au, $d$ $+$Au, and $^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV utilizing three sets of two-particle correlations for two detector combinations with different pseudorapidity acceptance [Phys. Rev. C {\bf 105}, 024901 (2022)]. This paper extends these measurements of $v_2$ to all centralities in $p$ $+$Au, $d$ $+$Au, and $^3$He$+$Au collisions, as well as $p$$+$$p$ collisions, as a function of transverse momentum ($p_T$) and event multiplicity. The kinematic dependence of $v_2$ is quantified as the ratio $R$ of $v_2$ between the two detector combinations as a function of event multiplicity for $0.5$ $<$ $p_T$ $<$ $1$ and $2$ $<$ $p_T$ $<$ $2.5$ GeV/$c$. A multiphase-transport (AMPT) model can reproduce the observed $v_2$ in most-central to midcentral $d$$+$Au and $^3$He$+$Au collisions. However, the AMPT model systematically overestimates the measurements in $p$ $+$ $p$, $p$ $+$Au, and peripheral $d$$+$Au and $^3$He$+$Au collisions, indicating a higher nonflow contribution in AMPT than in the experimental data. The AMPT model fails to describe the observed $R$ for $0.5$ $<$ $p_T$$<$ $1$ GeV/$c$, but there is qualitative agreement with the measurements for $2$ $<$ $p_T$ $<$ $2.5$ GeV/$c$.

We present measurements of azimuthal correlations of charged hadron pairs in $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions after subtracting an underlying event using a model that includes higher-order azimuthal anisotropy $v_2$, $v_3$, and $v_4$. After subtraction, the away-side ($\Delta\phi\sim\pi)$ of the highest transverse-momentum trigger ($p_T>4$ GeV/$c$) correlations is suppressed compared to that of correlations measured in $p$$+$$p$ collisions. At the lowest associated particle $p_T$, the away-side shape and yield are modified. These observations are consistent with the scenario of radiative-jet energy loss. For the lowest-$p_T$ trigger correlations, an away-side yield exists and we explore the dependence of the shape of the away-side within the context of an underlying-event model. Correlations are also studied differentially versus event-plane angle $\Psi_n$. The angular correlations show an asymmetry when selecting the sign of the trigger-particle azimuthal angle with respect to the $\Psi_2$ event plane. This asymmetry and the measured suppression of the pair yield out of plane is consistent with a path-length-dependent energy loss. No $\Psi_3$ dependence can be resolved within experimental uncertainties.

The measurement of direct photons from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=39$ and 62.4 GeV in the transverse-momentum range $0.4<p_T<3$ Gev/$c$ is presented by the PHENIX collaboration at the Relativistic Heavy Ion Collider. A significant direct-photon yield is observed in both collision systems. A universal scaling is observed when the direct-photon $p_T$ spectra for different center-of-mass energies and for different centrality selections at $\sqrt{s_{_{NN}}}=62.4$ GeV is scaled with $(dN_{\rm ch}/d\eta)^{\alpha}$ for $\alpha=1.21{\pm}0.04$. This scaling also holds true for direct-photon spectra from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV measured earlier by PHENIX, as well as the spectra from Pb$+$Pb at $\sqrt{s_{_{NN}}}=2760$ GeV published by ALICE. The scaling power $\alpha$ seems to be independent of $p_T$, center of mass energy, and collision centrality. The spectra from different collision energies have a similar shape up to $p_T$ of 2 GeV/$c$. The spectra have a local inverse slope $T_{\rm eff}$ increasing with $p_T$ of $0.174\pm0.018$ GeV/$c$ in the range $0.4<p_T<1.3$ GeV/$c$ and increasing to $0.289\pm0.024$ GeV/$c$ for $0.9<p_T<2.1$ GeV/$c$. The observed similarity of low-$p_T$ direct-photon production from $\sqrt{s_{_{NN}}}= 39$ to 2760 GeV suggests a common source of direct photons for the different collision energies and event centrality selections, and suggests a comparable space-time evolution of direct-photon emission.

In 2015 the PHENIX collaboration at the Relativistic Heavy Ion Collider recorded $p+p$, $p+$Al, and $p+$Au collision data at center of mass energies of $\sqrt{s_{_{NN}}}=200$ GeV with the proton beam(s) transversely polarized. At very forward rapidities $\eta>6.8$ relative to the polarized proton beam, neutrons were detected either inclusively or in (anti)correlation with detector activity related to hard collisions. The resulting single spin asymmetries, that were previously reported, have now been extracted as a function of the transverse momentum of the neutron as well as its longitudinal momentum fraction $x_F$. The explicit kinematic dependence, combined with the correlation information allows for a closer look at the interplay of different mechanisms suggested to describe these asymmetries, such as hadronic interactions or electromagnetic interactions in ultra-peripheral collisions, UPC. Events that are correlated with a hard collision indeed display a mostly negative asymmetry that increases in magnitude as a function of transverse momentum with only little dependence on $x_F$. In contrast, events that are not likely to have emerged from a hard collision display positive asymmetries for the nuclear collisions with a kinematic dependence that resembles that of a UPC based model. Because the UPC interaction depends strongly on the charge of the nucleus, those effects are very small for $p+p$ collisions, moderate for $p+$Al collisions, and large for $p+$Au collisions.

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.

Studying spin-momentum correlations in hadronic collisions offers a glimpse into a three-dimensional picture of proton structure. The transverse single-spin asymmetry for midrapidity isolated direct photons in $p^\uparrow+p$ collisions at $\sqrt{s}=200$ GeV is measured with the PHENIX detector at the Relativistic Heavy Ion Collider (RHIC). Because direct photons in particular are produced from the hard scattering and do not interact via the strong force, this measurement is a clean probe of initial-state spin-momentum correlations inside the proton and is in particular sensitive to gluon interference effects within the proton. This is the first time direct photons have been used as a probe of spin-momentum correlations at RHIC. The uncertainties on the results are a fifty-fold improvement with respect to those of the one prior measurement for the same observable, from the Fermilab E704 experiment. These results constrain gluon spin-momentum correlations in transversely polarized protons.

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.

Polarized proton-proton collisions provide leading-order access to gluons, presenting an opportunity to constrain gluon spin-momentum correlations within transversely polarized protons and enhance our understanding of the three-dimensional structure of the proton. Midrapidity open-heavy-flavor production at $\sqrt{s}=200$ GeV is dominated by gluon-gluon fusion, providing heightened sensitivity to gluon dynamics relative to other production channels. Transverse single-spin asymmetries of positrons and electrons from heavy-flavor hadron decays are measured at midrapidity using the PHENIX detector at the Relativistic Heavy Ion Collider. These charge-separated measurements are sensitive to gluon correlators that can in principle be related to gluon orbital angular momentum via model calculations. Explicit constraints on gluon correlators are extracted for two separate models, one of which had not been constrained previously.

Small nuclear collisions are mainly sensitive to cold-nuclear-matter effects; however, the collective behavior observed in these collisions shows a hint of hot-nuclear-matter effects. The identified-particle spectra, especially the $\phi$ mesons which contain strange and antistrange quarks and have a relatively small hadronic-interaction cross section, are a good tool to study these effects. The PHENIX experiment has measured $\phi$ mesons in a specific set of small collision systems $p$$+$Al, $p$$+$Au, and $^3$He$+$Au, as well as $d$$+$Au [Phys. Rev. C {\bf 83}, 024909 (2011)], at $\sqrt{s_{_{NN}}}=200$ GeV. The transverse-momentum spectra and nuclear-modification factors are presented and compared to theoretical-model predictions. The comparisons with different calculations suggest that quark-gluon plasma may be formed in these small collision systems at $\sqrt{s_{_{NN}}}=200$ GeV. However, the volume and the lifetime of the produced medium may be insufficient for observing strangeness-enhancement and jet-quenching effects. Comparison with calculations suggests that the main production mechanisms of $\phi$ mesons at midrapidity may be different in $p$$+$Al versus $p/d/$$^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. While thermal quark recombination seems to dominate in $p/d/$$^3$He$+$Au collisions, fragmentation seems to be the main production mechanism in $p$$+$Al collisions.

Suppression of the $J/\psi$ nuclear-modification factor has been seen as a trademark signature of final-state effects in large collision systems for decades. In small systems, the nuclear modification was attributed to cold-nuclear-matter effects until the observation of strong differential suppression of the $\psi(2S)$ state in $p/d$ $+$ $A$ collisions suggested the presence of final-state effects. Results of $J/\psi$ and $\psi(2S)$ measurements in the dimuon decay channel are presented here for $p$ $+$ $p$, $p$ $+$Al, and $p$ $+$Au collision systems at $\sqrt{s_{_{NN}}}=200$ GeV. The results are predominantly shown in the form of the nuclear-modification factor, $R_{pA}$, the ratio of the $\psi(2S)$ invariant yield per nucleon-nucleon collision in collisions of proton on target nucleus to that in $p$ $+$ $p$ collisions. Measurements of the $J/\psi$ and $\psi(2S)$ nuclear-modification factor are compared with shadowing and transport-model predictions, as well as to complementary measurements at Large-Hadron-Collider energies.

In 2015, the PHENIX collaboration has measured single-spin asymmetries for charged pions in transversely polarized proton-proton collisions at the center of mass energy of $\sqrt{s}=200$ GeV. The pions were detected at central rapidities of $|\eta|<0.35$. The single-spin asymmetries are consistent with zero for each charge individually, as well as consistent with the previously published neutral-pion asymmetries in the same rapidity range. However, they show a slight indication of charge-dependent differences which may suggest a flavor dependence in the underlying mechanisms that create these asymmetries.

We present measurements of the cross section and double-helicity asymmetry $A_{LL}$ of direct-photon production in $\vec{p}+\vec{p}$ collisions at $\sqrt{s}=510$ GeV. The measurements have been performed at midrapidity ($|\eta|<0.25$) with the PHENIX detector at the Relativistic Heavy Ion Collider. At relativistic energies, direct photons are dominantly produced from the initial quark-gluon hard scattering and do not interact via the strong force at leading order. Therefore, at $\sqrt{s}=510$ GeV, where leading-order-effects dominate, these measurements provide clean and direct access to the gluon helicity in the polarized proton in the gluon-momentum-fraction range $0.02<x<0.08$, with direct sensitivity to the sign of the gluon contribution.