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.
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
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.
Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in high-multiplicity 0%-5% central $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.
Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in high-multiplicity 0%-5% central $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.
Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.
The experimental study of the collisions of heavy nuclei at relativistic energies has established the properties of the quark-gluon plasma (QGP), a state of hot, dense nuclear matter in which quarks and gluons are not bound into hadrons. In this state, matter behaves as a nearly inviscid fluid that efficiently translates initial spatial anisotropies into correlated momentum anisotropies among the produced particles, producing a common velocity field pattern known as collective flow. In recent years, comparable momentum anisotropies have been measured in small-system proton-proton ($p$$+$$p$) and proton-nucleus ($p$$+$$A$) collisions, despite expectations that the volume and lifetime of the medium produced would be too small to form a QGP. Here, we report on the observation of elliptic and triangular flow patterns of charged particles produced in proton-gold ($p$$+$Au), deuteron-gold ($d$$+$Au), and helium-gold ($^3$He$+$Au) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_{NN}}}$~=~200 GeV. The unique combination of three distinct initial geometries and two flow patterns provides unprecedented model discrimination. Hydrodynamical models, which include the formation of a short-lived QGP droplet, provide a simultaneous description of these measurements.
$v_2$for 0-5% central p+Au collisions
$v_2$for 0-5% central d+Au collisions
$v_2$for 0-5% central $^3$He+Au collisions
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.
Fourier fit coefficients for CNT-MPCS (Au-going) correlations, as a function of collision system and $\pi^0$ $p_T$: (a) the negative of the dipole coefficient, $-c_1$; (b) the quadrupole coefficient $c_2$; (c) the ratio ${-c_2}/{c_1}$.
Fourier fit coefficients for CNT-MPCS (Au-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au.
Fourier fit coefficients for CNT-MPCS (Au-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $p$+$p$.
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.
Values of $v_2$($p_T$) for kaons in central 0-5% $p$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.
Values of $v_2$($p_T$) for kaons in central 0-5% $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.
Transverse momentum dependence of $v_2$ for identified pions and protons within $|\eta|$ < 0.35 in 0-5% central $p$+Au 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).
The $\phi$-meson-production cross section d$\sigma_{\phi}$/dy in $p$ + $p$ collisions at $\sqrt{s}$ = 510 GeV integrated in the transverse-momentum range 2 < $p_T$ < 7 GeV/$c$.
The $\phi$-meson-production cross section d$\sigma_{\phi}$/dy in $p$ + $p$ collisions at $\sqrt{s}$ = 510 GeV integrated in the transverse-momentum range 2 < $p_T$ < 7 GeV/$c$.
The $\phi$-meson-differential-production cross section d${}^{2}$$\sigma_{\phi}/dp_T dy$ for 1.2 < |y| < 2.2 in $p$ + $p$ collisions at $\sqrt{s}$ = 510 GeV.
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.
Example fits of Bose-Einstein correlation functions of (a) $\pi^{-}\pi^{-}$ pair with $m_{T}$ between 0.331 and 0.349 GeV/$c^2$ and of (b) $\pi^{+}\pi^{+}$ pair with $m_T$ between 0.655 and 0.675 GeV/$c^2$, as a function $Q$ ≡ |$q_{LCMS}$|, defined in Eq. (26). Both fits show the measured correlation function and the complete fit function (described in VI A), while a Bose-Einstein fit function $C^{(0)}_{2} (Q)$ is also shown, with the Coulomb-corrected data, i.e. the raw data multiplied by $C^{(0)}_{2} (Q)/C_{2}(Q)$. In this analysis we measured 62 such correlation functions (for ++ and -- pairs, in 31 $m_T$ bins), and fitted all of them with the method described in VIA. The first visible point on both panels corresponds to $Q$ values below the accessible range (based on an evaluation of the two-track cuts), these were not taken into account in the fitting.
Example fits of Bose-Einstein correlation functions of (a) $\pi^{-}\pi^{-}$ pair with $m_{T}$ between 0.331 and 0.349 GeV/$c^2$ and of (b) $\pi^{+}\pi^{+}$ pair with $m_T$ between 0.655 and 0.675 GeV/$c^2$, as a function $Q$ ≡ |$q_{LCMS}$|, defined in Eq. (26). Both fits show the measured correlation function and the complete fit function (described in VI A), while a Bose-Einstein fit function $C^{(0)}_{2} (Q)$ is also shown, with the Coulomb-corrected data, i.e. the raw data multiplied by $C^{(0)}_{2} (Q)/C_{2}(Q)$. In this analysis we measured 62 such correlation functions (for ++ and -- pairs, in 31 $m_T$ bins), and fitted all of them with the method described in VIA. The first visible point on both panels corresponds to $Q$ values below the accessible range (based on an evaluation of the two-track cuts), these were not taken into account in the fitting.
Correlation strength parameter $\lambda$ versus average $m_T$ of the pair, for 0%-30% centrality collisions. Statistical and systematic uncertainties are shown as bars and boxes.
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.
Forward neutron $A_N$ in $p$+$A$ collisions for $A$ = 1 ($p$), 27 (Al), and 197 (Au), for ZDC inclusive, ZDC $\otimes$ BBC-tag, and ZDC $\otimes$ BBC-veto triggered samples. The 3% scale uncertainty is from the polarization normalization uncertaintiy.
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 ratio of the two harmonics ($\Sigma Q^{BBC-S}$)$_{p+p}$/($\Sigma Q^{BBC-S}$)$_{pAu}$.
$v_2$ of charged hadrons within $|\eta|$ < 0.35 in 0-5% $p$+Au central collisions, compared to hydrodynamic calculations using the SONIC model, matched to the same multiplicity as the data.