The $J/\psi$ and $\psi(2S)$ charmonium states, composed of $c\bar{c}$ quark pairs and known since the 1970s, are widely believed to serve as ideal probes to test quantum chromodynamics in high-energy hadronic interactions. However, there is not yet a complete understanding of the charmonium-production mechanism. Recent measurements of $J/\psi$ production as a function of event charged-particle multiplicity at the collision energies of both the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) show enhanced $J/\psi$ production yields with increasing multiplicity. One potential explanation for this type of dependence is multiparton interactions (MPI). We carry out the first measurements of self-normalized $J/\psi$ yields and the $\psi(2S)$ to $J/\psi$ ratio at both forward and backward rapidities as a function of self-normalized charged-particle multiplicity in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. In addition, detailed {\sc pythia} studies tuned to RHIC energies were performed to investigate the MPI impacts. We find that the PHENIX data at RHIC are consistent with recent LHC measurements and can only be described by {\sc pythia} calculations that include MPI effects. The forward and backward $\psi(2S)$ to $J/\psi$ ratio, which serves as a unique and powerful approach to study final-state effects on charmonium production, is found to be less dependent on the charged-particle multiplicity.
Self-normalized $J/\psi$ yields as a function of self-normalized $N_{ch}$ for the same arm before subtraction
Self-normalized $J/\psi$ yields as a function of self-normalized $N_{ch}$ for the same arm after subtraction
Self-normalized $J/\psi$ yields as a function of self-normalized $N_{ch}$ for opposite arms
The jet cross-section and jet-substructure observables in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV were measured by the PHENIX Collaboration at the Relativistic Heavy Ion Collider (RHIC). Jets are reconstructed from charged-particle tracks and electromagnetic-calorimeter clusters using the anti-$k_{t}$ algorithm with a jet radius $R=0.3$ for jets with transverse momentum within $8.0
The jet differential cross section as a function of jet $p_T$. Statistical uncertainties are typically smaller than the data points while systematic uncertainties are shown with boxes. An overall normalization systematic of 7% is not included in the point-by-point systematic uncertainties.
Distribution of the SoftDrop groomed momentum fraction $z_g$ for different jet $p_T$ bins. Standard SoftDrop parameters were used ($z_{cut}<0.1$ and $\beta=0$).
$\xi$ distributions for different jet $p_T$ bins.
The PHENIX experiment reports systematic measurements at the Relativistic Heavy Ion Collider of $\phi$-meson production in asymmetric Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV and in U$+$U collisions at $\sqrt{s_{_{NN}}}$=193 GeV. Measurements were performed via the $\phi\rightarrow K^{+}K^{-}$ decay channel at midrapidity $|\eta|<0.35$. Features of $\phi$-meson production measured in Cu$+$Cu, Cu$+$Au, Au$+$Au, and U$+$U collisions were found to not depend on the collision geometry, which was expected because the yields are averaged over the azimuthal angle and follow the expected scaling with nuclear-overlap size. The elliptic flow of the $\phi$ meson in Cu$+$Au, Au$+$Au, and U$+$U collisions scales with second-order-participant eccentricity and the length scale of the nuclear-overlap region (estimated with the number of participating nucleons). At moderate $p_T$, $\phi$-meson production measured in Cu$+$Au and U$+$U collisions is consistent with coalescence-model predictions, whereas at high $p_T$ the production is in agreement with expectations for in-medium energy loss of parent partons prior to their fragmentation. The elliptic flow for $\phi$ mesons measured in Cu$+$Au and U$+$U collisions is well described by a (2+1)D viscous-hydrodynamic model with specific-shear viscosity $\eta/s=1/4\pi$.
Invariant transverse momentum spectra measured for $\phi$ mesons in (a) Cu+Au and (b) U+U collisions at $\sqrt{s}$ = 200 GeV at midrapidity
Invariant transverse momentum spectra measured for $\phi$ mesons in (c) Cu+Au and (d) U+U collisions at $\sqrt{s}$ = 200 GeV at midrapidity. Data-to- Levy fit ratios.
The $\phi$-meson nuclear modification factors $R_{AB}$ measured as a function of $p_T$ in different centrality intervals of (a) to (d) Cu+Au collisions at $\sqrt{s}$ = 200 GeV and (e) to (h) U+U collisions at $\sqrt{s}$ = 193 GeV at midrapidity $|\eta|<0.35$. The normalization uncertainty from p+p of about $\sim9.7\%$ is not shown.
The measurement of the direct-photon spectrum from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV is presented by the PHENIX collaboration using the external-photon-conversion technique for 0%--93% central collisions in a transverse-momentum ($p_T$) range of 0.8--10 GeV/$c$. An excess of direct photons, above prompt-photon production from hard-scattering processes, is observed for $p_T<6$ GeV/$c$. Nonprompt direct photons are measured by subtracting the prompt component, which is estimated as $N_{\rm coll}$-scaled direct photons from $p$ $+$ $p$ collisions at 200 GeV, from the direct-photon spectrum. Results are obtained for $0.8
Direct photon $R_{\gamma}$, every 20% centrality
Direct photon $R_{\gamma}$, every 10% centrality
Invariant yield of direct photons, every 10% centrality
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
$R_{\gamma}$ for minimum bias (0-86%) Au+Au collision at $\sqrt{s_{NN}} = 62.4$ GeV (a) and $39$ GeV (b). For $62.4$ GeV also centrality bins of 0-20% (c) and 20-40% (d) are shown. Data points are shown with statistical (bar) and systematic uncertainties (box)
$R_{\gamma}$ for minimum bias (0-86%) Au+Au collision at $\sqrt{s_{NN}} = 62.4$ GeV (a) and $39$ GeV (b). For $62.4$ GeV also centrality bins of 0-20% (c) and 20-40% (d) are shown. Data points are shown with statistical (bar) and systematic uncertainties (box)
Direct photon spectra for minimum bias (0-86%) Au+Au collision at $\sqrt{s_{NN}} = 62.4$ GeV (a) and $39$ GeV (b). For $62.4$ GeV also centrality bins of 0-20% (c) and 20-40% (d) are shown. Data points are shown with statistical and systematic uncertainties, unless the central value is negative (arrows) or is consistent with zero within the statistical uncertainties (arrows with data point). In these cases upper limit with CL = 95$%$ are given.
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$.
Azimuthal anisotropy $v_2\{BB\}$ as a function of transverse momentum $p_T$ in $p$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.
Azimuthal anisotropy $v_2\{BF\}$ as a function of transverse momentum $p_T$ in $p$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.
Azimuthal anisotropy $v_2\{BB\}$ as a function of transverse momentum $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.
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.
Invariant transverse momentum spectra measured for $\phi$ mesons in (a) $p$+Al, (b) $p$+Au, and (c) $^{3}$He+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV at midrapidity.
Comparison of $\phi$-meson nuclear-modification factors in $p$+Al, $p$+Au, $d$+Au [2], and $^{3}$He+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV at midrapidity. The normalization uncertainty from $p$+$p$ of about $9.7 \%$ is not shown [28].
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
Cross sections for inclusive and isolated direct photons as a function of $p_T$. Not shown are 10% absolute luminosity uncertainties.
Double helicity asymmetry $A_{LL}$ $vs$ $p_{T}$ for isolated direct-photon production in polarized $\vec{p}+\vec{p}$ collisions at $\sqrt{s}=510$ GeV at midrapidity. Not shown are $3.9 \times 10^{-4}$ shift uncertainty from relative luminosity and 6.6% scale uncertainty from polarization.
The PHENIX collaboration presents a systematic study of $\pi^0$ production from $p$ $+$ $p$, $p$ $+$Al, $p$ $+$Au, $d$ $+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Measurements were performed with different centrality selections as well as the total inelastic, 0%--100%, selection for all collision systems. For 0%--100% collisions, the nuclear modification factors, $R_{xA}$, are consistent with unity for $p_T$ above 8 GeV/$c$, but exhibit an enhancement in peripheral collisions and a suppression in central collisions. The enhancement and suppression characteristics are similar for all systems for the same centrality class. It is shown that for high-$p_T$-$\pi^0$ production, the nucleons in the $d$ and $^3$He interact mostly independently with the Au nucleus and that the counter intuitive centrality dependence is likely due to a physical correlation between multiplicity and the presence of a hard scattering process. These observations disfavor models where parton energy loss has a significant contribution to nuclear modifications in small systems. Nuclear modifications at lower $p_T$ resemble the Cronin effect -- an increase followed by a peak in central or inelastic collisions and a plateau in peripheral collisions. The peak height has a characteristic ordering by system size as $p$ $+$Au $>$ $d$ $+$Au $>$ $^{3}$He$+$Au $>$ $p$ $+$Al. For collisions with Au ions, current calculations based on initial state cold nuclear matter effects result in the opposite order, suggesting the presence of other contributions to nuclear modifications, in particular at lower $p_T$.
Differential cross section of $\pi^0$ in p+p collisions at $\sqrt{s}$ = 200 GeV
Invariant yield of $\pi^0$ from (a) p+Al, (b) p+Au, (c) d+Au, and (d) $^{3}$HeAu in different centrality selections at $\sqrt{s}$ = 200 GeV
Nuclear modification factors from inelastic (a) p+Al, (b) p+Au, (c) d+Au, and (d) $^{3}$HeAu collisions at $\sqrt{s}$ = 200 GeV. The right boxes are the $N_{coll}$ uncertainties from the Glauber model, while the left box represents the overall normalization uncertainty from p+p collisions
There is strong evidence for the formation of small droplets of quark-gluon plasma in $p/d/^{3}$He+Au collisions at the Relativistic Heavy Ion Collider (RHIC) and in $p$+$p$/Pb collisions at the Large Hadron Collider. In particular, the analysis of data at RHIC for different geometries obtained by varying the projectile size and shape has proven insightful. In the present analysis, we find excellent agreement with the previously published PHENIX at RHIC results on elliptical and triangular flow with an independent analysis via the two-particle correlation method, which has quite different systematic uncertainties and an independent code base. In addition, the results are extended to other detector combinations with different kinematic (pseudorapidity) coverage. These results provide additional constraints on contributions from nonflow and longitudinal decorrelations.
$v_2$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination
$v_2$ vs $p_T$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination
$v_2$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination