-
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
partial 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.28 \pm 0.09\,{\rm (stat)} \pm 0.14\,{\rm
(syst)} \pm 0.27\, {\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 ({\sc pythia6, pythia8, phojet,
ampt, epos3,} and {\sc epos-lhc}).
-
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.
-
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.
-
Recently, multiparticle-correlation measurements of relativistic
$p/d/^3$He$+$Au, $p$$+$Pb, and even $p$$+$$p$ collisions have shown surprising
collective signatures. Here we present beam-energy-scan measurements of 2-, 4-,
and 6-particle angular correlations in $d$$+$Au collisions at
$\sqrt{s_{_{NN}}}$=200, 62.4, 39, and 19.6 GeV. We also present measurements of
2- and 4-particle angular correlations in $p$$+$Au collisions at
$\sqrt{s_{_{NN}}}$=200 GeV. We find the 4-particle cumulant to be real-valued
for $d$$+$Au collisions at all four energies. We also find that the 4-particle
cumulant in $p$$+$Au has the opposite sign as that in $d$$+$Au, indicating that
the correlations are geometrical and therefore collective in origin. Further we
find that the 6-particle cumulant agrees with the 4-particle cumulant in
$d$$+$Au collisions at 200 GeV, which indicates that nonflow effects are
subdominant and provides strong evidence of collectivity.
-
We present measurements of the elliptic flow, $v_2$, as a function of
transverse momentum ($p_T$), pseudorapidity ($\eta$), and centrality in
$d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200, 62.4, 39, and 19.6 GeV. The
beam-energy scan of $d$$+$Au collisions provides a testing ground for the onset
of flow signatures in small collision systems. We measure a nonzero $v_2$
signal at all four collision energies, which, at midrapidity and low $p_T$, is
consistent with predictions from viscous hydrodynamic models. Comparisons with
calculations from parton transport models (based on the {\sc ampt} Monte Carlo
generator) show good agreement with the data at midrapidity to forward
($d$-going) rapidities and low $p_T$. At backward (Au-going) rapidities and
$p_T>1.5$ GeV/$c$, the data diverges from {\sc ampt} calculations of $v_2$
relative to the initial geometry, indicating the possible dominance of
nongeometry related correlations, referred to as nonflow. We also present
measurements of the charged-particle multiplicity ($dN_{\rm ch}/d\eta$) as a
function of $\eta$ in central $d$$+$Au collisions at the same energies. We find
that in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200 GeV the $v_2$ scales with
$dN_{\rm ch}/d\eta$ over all $\eta$ in the PHENIX acceptance. At
$\sqrt{s_{_{NN}}}=$62.4, and 39 GeV, $v_2$ scales with $dN_{\rm ch}/d\eta$ at
midrapidity and forward rapidity, but falls off at backward rapidity. This
departure from the $dN_{\rm ch}/d\eta$ scaling may be a further indication of
nonflow effects dominating at backward rapidity.
-
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 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 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.
-
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.
-
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.
-
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.
-
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.
-
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.
Recently, multiparticle-correlation measurements of relativistic $p/d/^3$He$+$Au, $p$$+$Pb, and even $p$$+$$p$ collisions have shown surprising collective signatures. Here we present beam-energy-scan measurements of 2-, 4-, and 6-particle angular correlations in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200, 62.4, 39, and 19.6 GeV. We also present measurements of 2- and 4-particle angular correlations in $p$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. We find the 4-particle cumulant to be real-valued for $d$$+$Au collisions at all four energies. We also find that the 4-particle cumulant in $p$$+$Au has the opposite sign as that in $d$$+$Au, indicating that the correlations are geometrical and therefore collective in origin. Further we find that the 6-particle cumulant agrees with the 4-particle cumulant in $d$$+$Au collisions at 200 GeV, which indicates that nonflow effects are subdominant and provides strong evidence of collectivity.
We present measurements of the elliptic flow, $v_2$, as a function of transverse momentum ($p_T$), pseudorapidity ($\eta$), and centrality in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200, 62.4, 39, and 19.6 GeV. The beam-energy scan of $d$$+$Au collisions provides a testing ground for the onset of flow signatures in small collision systems. We measure a nonzero $v_2$ signal at all four collision energies, which, at midrapidity and low $p_T$, is consistent with predictions from viscous hydrodynamic models. Comparisons with calculations from parton transport models (based on the {\sc ampt} Monte Carlo generator) show good agreement with the data at midrapidity to forward ($d$-going) rapidities and low $p_T$. At backward (Au-going) rapidities and $p_T>1.5$ GeV/$c$, the data diverges from {\sc ampt} calculations of $v_2$ relative to the initial geometry, indicating the possible dominance of nongeometry related correlations, referred to as nonflow. We also present measurements of the charged-particle multiplicity ($dN_{\rm ch}/d\eta$) as a function of $\eta$ in central $d$$+$Au collisions at the same energies. We find that in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200 GeV the $v_2$ scales with $dN_{\rm ch}/d\eta$ over all $\eta$ in the PHENIX acceptance. At $\sqrt{s_{_{NN}}}=$62.4, and 39 GeV, $v_2$ scales with $dN_{\rm ch}/d\eta$ at midrapidity and forward rapidity, but falls off at backward rapidity. This departure from the $dN_{\rm ch}/d\eta$ scaling may be a further indication of nonflow effects dominating at backward rapidity.
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 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 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.
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.
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.
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.
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.
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.