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
PHENIX presents a simultaneous measurement of the production of direct $\gamma$ and $\pi^0$ in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV over a $p_T$ range of 7.5 to 18 GeV/$c$ for different event samples selected by event activity, i.e. charged-particle multiplicity detected at forward rapidity. Direct-photon yields are used to empirically estimate the contribution of hard-scattering processes in the different event samples. Using this estimate, the average nuclear-modification factor $R_{d\rm Au,EXP}^{\gamma^{\rm dir}}$ is $0.925{\pm}0.023({\rm stat}){\pm}0.15^{\rm (scale)}$, consistent with unity for minimum-bias (MB) $d$$+$Au events. For event classes with moderate event activity, $R_{d\rm Au,EXP}^{\gamma^{\rm dir}}$ is consistent with the MB value within 5% uncertainty. These results confirm that the previously observed enhancement of high-$p_T$$\pi^0$ production found in small-system collisions with low event activity is a result of a bias in interpreting event activity within the Glauber framework. In contrast, for the top 5% of events with the highest event activity, $R_{d\rm Au,EXP}^{\gamma^{\rm dir}}$ is suppressed by 20% relative to the MB value with a significance of $4.5\sigma$, which may be due to final-state effects.
The $p_{T}$ distribution at high $p_{T}$ of neutral pions (panel (a)) for different d+Au event activity classes compared to those from p+p collisions (Phys. Rev. C 105, 064902).
The $p_{T}$ distribution at high $p_{T}$ of direct photons (panel (b)) for different d+Au event activity classes compared to those from p+p collisions (Phys. Rev. D 86, 072008).
Panel (c) shows the ratio $\gamma^{dir}/\pi^{0}$. For better visibility the points are slightly shifted in $p_{T}$ .
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.
Data from Figure 2 (a) of the $\pi^{0}$ transverse single-spin asymmetry in $\sqrt{s_{NN}}=200$ GeV $p^{\uparrow}+$Au and $p^{\uparrow}+$Al collisions as a function of $p_{T}$.
Data from Figure 2 (b) of the $\eta$ transverse single-spin asymmetry in $\sqrt{s_{NN}}=200$ GeV $p^{\uparrow}+$Au and $p^{\uparrow}+$Al collisions as a function of $p_{T}$.
We present the first measurements of transverse momentum spectra of $\pi^{\pm}$, $K^{\pm}$, $p(\bar{p})$ at midrapidity ($|y| < 0.1$) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The centrality dependence of particle yields, average transverse momenta, particle ratios and kinetic freeze-out parameters are discussed. The results are compared with the published results from Au+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV in STAR. The results are also compared to those from A Multi Phase Transport (AMPT) model.
'Identified transverse momentum spectra of $\pi^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'
'Identified transverse momentum spectra of $K^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'
'Identified transverse momentum spectra of p at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'
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.
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.
Data from Figure 1 of open heavy flavor $e^{\pm}$ transverse single-spin asymmetries in transversely polarized p+p collisions as a function of $p_{T}$.
Notwithstanding decades of progress since Yukawa first developed a description of the force between nucleons in terms of meson exchange, a full understanding of the strong interaction remains a major challenge in modern science. One remaining difficulty arises from the non-perturbative nature of the strong force, which leads to the phenomenon of quark confinement at distances on the order of the size of the proton. Here we show that in relativistic heavy-ion collisions, where quarks and gluons are set free over an extended volume, two species of produced vector (spin-1) mesons, namely $\phi$ and $K^{*0}$, emerge with a surprising pattern of global spin alignment. In particular, the global spin alignment for $\phi$ is unexpectedly large, while that for $K^{*0}$ is consistent with zero. The observed spin-alignment pattern and magnitude for the $\phi$ cannot be explained by conventional mechanisms, while a model with a connection to strong force fields, i.e. an effective proxy description within the Standard Model and Quantum Chromodynamics, accommodates the current data. This connection, if fully established, will open a potential new avenue for studying the behaviour of strong force fields.
Global spin alignment of $\phi$ and $K^{*0}$ vector mesons in heavy-ion collisions. The measured matrix element $\rho_{00}$ as a function of beam energy for the $\phi$ and $K^{*0}$ vector mesons within the indicated windows of centrality, transverse momentum ($p_T$) and rapidity ($y$). The open symbols indicate ALICE results for Pb+Pb collisions at 2.76 TeV at $p_{T}$ values of 2.0 and 1.4 GeV/c for the $\phi$ and $K^{*0}$ mesons, respectively, corresponding to the $p_{T}$ bin nearest to the mean $p_{T}$ for the 1.0 – 5.0 GeV/$c$ range assumed for each meson in the present analysis. The red solid curve is a fit to data in the range of $\sqrt{s_{NN}} = 19.6$ to 200 GeV, based on a theoretical calculation with a $\phi$-meson field. Parameter sensitivity of $\rho_{00}$ to the $\phi$-meson field is shown in Ref.5. The red dashed line is an extension of the solid curve with the fitted parameter $G_s^{(y)}$. The black dashed line represents $\rho_{00}=1/3.$
Global spin alignment of $\phi$ and $K^{*0}$ vector mesons in heavy-ion collisions. The measured matrix element $\rho_{00}$ as a function of beam energy for the $\phi$ and $K^{*0}$ vector mesons within the indicated windows of centrality, transverse momentum ($p_T$) and rapidity ($y$). The open symbols indicate ALICE results for Pb+Pb collisions at 2.76 TeV at $p_{T}$ values of 2.0 and 1.4 GeV/c for the $\phi$ and $K^{*0}$ mesons, respectively, corresponding to the $p_{T}$ bin nearest to the mean $p_{T}$ for the 1.0 – 5.0 GeV/$c$ range assumed for each meson in the present analysis. The red solid curve is a fit to data in the range of $\sqrt{s_{NN}} = 19.6$ to 200 GeV, based on a theoretical calculation with a $\phi$-meson field. Parameter sensitivity of $\rho_{00}$ to the $\phi$-meson field is shown in Ref.5. The red dashed line is an extension of the solid curve with the fitted parameter $G_s^{(y)}$. The black dashed line represents $\rho_{00}=1/3.$
Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.
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<p_T<6.0$ GeV/$c$ and suggest that the spectrum has an increasing inverse slope from ${\approx}0.2$ to 0.4 GeV/$c$ with increasing $p_T$, which indicates a possible sensitivity of the measurement to photons from earlier stages of the evolution of the collision. In addition, like the direct-photon production, the $p_T$-integrated nonprompt direct-photon yields also follow a power-law scaling behavior as a function of collision-system size. The exponent, $\alpha$, for the nonprompt component is found to be consistent with 1.1 with no apparent $p_T$ dependence.
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<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.
$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.