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.
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.
$A_N$ of negatively-charged muons from open heavy-flavor decays as a function of $p_T$ in the backward ($x_F$ < 0) and forward ($x_F$ > 0) regions.
$A_N$ of positively-charged muons from open heavy-flavor decays as a function of $p_T$ in the backward ($x_F$ < 0) and forward ($x_F$ > 0) regions.
$A_N$ of negatively-charged and positively-charged muons from open heavy-flavor decays as a function of $x_F$, where $x_F$ > 0 is along the direction of the polarized proton.
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.
Differential cross section for $p$+$p$ $\rightarrow$ $b\bar{b}$ at $\sqrt{s}$ = 200 GeV.
Fraction F$_{B \rightarrow J/\psi}$ of $B$-meson decays in the inclusive $J/\psi$ sample in $p$+$p$ and Cu+Au collisions at $sqrt{s_{NN}}$ = 200 GeV versus rapidity along with a theoretical estimation based on fixed-order plus next-to-leading logs (FONLL) for the $B \rightarrow J/\psi$ cross section and Color-Evaporation-Model (CEM) for the prompt $J/\psi$.
Fraction F$_{B \rightarrow J/\psi}$ of $B$-meson decays in the inclusive $J/\psi$ sample in $p$+$p$ and Cu+Au collisions at $sqrt{s_{NN}}$ = 200 GeV versus rapidity along with a theoretical estimation based on fixed-order plus next-to-leading logs (FONLL) for the $B \rightarrow J/\psi$ cross section and Color-Evaporation-Model (CEM) for the prompt $J/\psi$.
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.
Step by step extrapolation from the number of $e^+e^-$ pairs for $m_{e^+e^-}$ > 1.16 GeV/$c^2$ from $c\bar{c}$ in the PHENIX acceptance to the number of $c\bar{c}$ pairs in 4$\pi$ for PYTHIA, MC@NLO, and POWHEG. Numbers are in units of pairs per event using the $c\bar{c}$ cross sections determined in this paper.
Step by step extrapolation from the number of $e^+e^-$ pairs for $m_{e^+e^-}$ > 1.16 GeV/$c^2$ from $b\bar{b}$ in the PHENIX acceptance to the number of $b\bar{b}$ pairs in 4$\pi$ for PYTHIA, MC@NLO, and POWHEG. Numbers are in units of pairs per event using the $b\bar{b}$ cross sections determined in this paper.
Summary of $c\bar{c}$ and $b\bar{b}$ cross sections measured in $p$+$p$ collisions using three different generators, PYTHIA, MC@NLO, and POWHEG.
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.
Fraction of $B$-meson decays in $J/\psi$ samples obtained in $p$+$p$ collisions at $\sqrt{S}$ = 510 GeV.
Fraction of $B$-meson decays in $J/\psi$ samples obtained in $p$+$p$ collisions at $\sqrt{S}$ = 510 GeV.
The $b\bar{b}$ cross section per unit rapidity ($d\sigma/dy(pp→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.
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.
Polar angular decay coefficient $\lambda_{\theta}$ as a function of transverse momentum for four reference frames and three $p_T$ bins. The numbers in the CS frame for the $p_T$ = 2-3 GeV/$c$ bin are 90% confidence level upper limits.
"Mixed" angular decay coefficient $\lambda_{\theta \phi}$ as a function of transverse momentum for four reference frames and three $p_T$ bins.
Azimuthal angular decay coefficient $\lambda_{\phi}$ as a function of transverse momentum for four reference frames and three $p_T$ bins.
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.
Summary of the measured ratios of $\Psi$(2S)/$\Psi$(1S) mesons.
Summary of the measured ratios of $\Psi$(2S)/$\Psi$(1S) mesons.
Summary of the measured ratios of $\Psi$(2S)/$\Psi$(1S) mesons.