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 measured the centrality dependence of two-pion Bose-Einstein correlation functions in $\sqrt{s_{_{NN}}}=200$~GeV Au$+$Au collisions at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. The data are well represented by Lévy-stable source distributions. The extracted source parameters are the correlation-strength parameter $\lambda$, the Lévy index of stability $\alpha$, and the Lévy-scale parameter $R$ as a function of transverse mass $m_T$ and centrality. The $\lambda(m_T)$ parameter is constant at larger values of $m_T$, but decreases as $m_T$ decreases. The Lévy scale parameter $R(m_T)$ decreases with $m_T$ and exhibits proportionality to the length scale of the nuclear overlap region. The Lévy exponent $\alpha(m_T)$ is independent of $m_T$ within uncertainties in each investigated centrality bin, but shows a clear centrality dependence. At all centralities, the Lévy exponent $\alpha$ is significantly different from that of Gaussian ($\alpha=2$) or Cauchy ($\alpha=1$) source distributions. Comparisons to the predictions of Monte-Carlo simulations of resonance-decay chains show that in all but the most peripheral centrality class (50%-60%), the obtained results are inconsistent with the measurements, unless a significant reduction of the in-medium mass of the $\eta'$ meson is included. In each centrality class, the best value of the in-medium $\eta'$ mass is compared to the mass of the $\eta$ meson, as well as to several theoretical predictions that consider restoration of $U_A(1)$ symmetry in hot hadronic matter.
The transverse-mass dependence of the correlation-strength parameter $\lambda$ in 0-10% centrality bin obtained from Lévy fits with Eq. (9).
The transverse-mass dependence of the correlation-strength parameter $\lambda$ in 10-20% centrality bin obtained from Lévy fits with Eq. (9).
The transverse-mass dependence of the correlation-strength parameter $\lambda$ in 20-30% centrality bin obtained from Lévy fits with Eq. (9).
Atomic nuclei are self-organized, many-body quantum systems bound by strong nuclear forces within femtometer-scale space. These complex systems manifest a variety of shapes, traditionally explored using non-invasive spectroscopic techniques at low energies. However, at these energies, their instantaneous shapes are obscured by long-timescale quantum fluctuations, making direct observation challenging. Here we introduce the ``collective flow assisted nuclear shape imaging'' method, which images the nuclear global shape by colliding them at ultrarelativistic speeds and analyzing the collective response of outgoing debris. This technique captures a collision-specific snapshot of the spatial matter distribution within the nuclei, which, through the hydrodynamic expansion, imprints patterns on the particle momentum distribution observed in detectors. We benchmark this method in collisions of ground state Uranium-238 nuclei, known for their elongated, axial-symmetric shape. Our findings show a large deformation with a slight deviation from axial symmetry in the nuclear ground state, aligning broadly with previous low-energy experiments. This approach offers a new method for imaging nuclear shapes, enhances our understanding of the initial conditions in high-energy collisions and addresses the important issue of nuclear structure evolution across energy scales.
Data from Figure 2, panel a, Au+Au, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c, systematics include non-flow difference in the two systems, but correlated non-flow systematics with the value of $\delta \left\langle v_{2}^{2}\right\rangle$ =+-3.2e-5 included
Data from Figure 2, panel a, U+U, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c, systematics include non-flow difference in the two systems, but correlated non-flow systematics with the value of $\delta \left\langle v_{2}^{2}\right\rangle$ =+-3.2e-5 included
Data from Figure 2, panel b, Au+Au, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c
For the search of the chiral magnetic effect (CME), STAR previously presented the results from isobar collisions (${^{96}_{44}\text{Ru}}+{^{96}_{44}\text{Ru}}$, ${^{96}_{40}\text{Zr}}+{^{96}_{40}\text{Zr}}$) obtained through a blind analysis. The ratio of results in Ru+Ru to Zr+Zr collisions for the CME-sensitive charge-dependent azimuthal correlator ($\Delta\gamma$), normalized by elliptic anisotropy ($v_{2}$), was observed to be close to but systematically larger than the inverse multiplicity ratio. The background baseline for the isobar ratio, $Y = \frac{(\Delta\gamma/v_{2})^{\text{Ru}}}{(\Delta\gamma/v_{2})^{\text{Zr}}}$, is naively expected to be $\frac{(1/N)^{\text{Ru}}}{(1/N)^{\text{Zr}}}$; however, genuine two- and three-particle correlations are expected to alter it. We estimate the contributions to $Y$ from those correlations, utilizing both the isobar data and HIJING simulations. After including those contributions, we arrive at a final background baseline for $Y$, which is consistent with the isobar data. We extract an upper limit for the CME fraction in the $\Delta\gamma$ measurement of approximately $10\%$ at a $95\%$ confidence level on in isobar collisions at $\sqrt{s_{\text{NN}}} = 200$ GeV, with an expected $15\%$ difference in their squared magnetic fields.
Figure 1a, upper panel, full-event
Figure 1a, lower panel, full-event
Figure 1b, upper panel, subevent
We report results on an elastic cross section measurement in proton-proton collisions at a center-of-mass energy $\sqrt{s}=510$ GeV, obtained with the Roman Pot setup of the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The elastic differential cross section is measured in the four-momentum transfer squared range $0.23 \leq -t \leq 0.67$ GeV$^2$. We find that a constant slope $B$ does not fit the data in the aforementioned $t$ range, and we obtain a much better fit using a second-order polynomial for $B(t)$. The $t$ dependence of $B$ is determined using six subintervals of $t$ in the STAR measured $t$ range, and is in good agreement with the phenomenological models. The measured elastic differential cross section $\mathrm{d}\sigma/\mathrm{dt}$ agrees well with the results obtained at $\sqrt{s} = 546$ GeV for proton--antiproton collisions by the UA4 experiment. We also determine that the integrated elastic cross section within the STAR $t$-range is $\sigma^\mathrm{fid}_\mathrm{el} = 462.1 \pm 0.9 (\mathrm{stat.}) \pm 1.1 (\mathrm {syst.}) \pm 11.6 (\mathrm {scale})$~$\mu\mathrm{b}$.
Top panel: The $pp$ elastic differential cross section $d\sigma/dt$ fitted with an exponential $A e^{-B(t)|t|}$. Bottom panel: Residuals (Data - Fit)/Error. Uncertainties on the data points are smaller than the symbol size. The vertical scale uncertainty of 2.5% is not included in in the full error.
Results of the exponential function $A e^{-B(t)|t|}$ fit to the elastic differential cross section data as well as the integrated fiducial cross section are listed. Also listed are the corresponding values of the statistical and systematic uncertainties. The scale (luminosity and trigger efficiency) uncertainty of 2.5% applicable to the fit parameter $A$ and fiducial cross section $\sigma^\mathrm{fid}_\mathrm{el}$ is not included in the full error.
The chiral magnetic effect (CME) is a phenomenon that arises from the QCD anomaly in the presence of an external magnetic field. The experimental search for its evidence has been one of the key goals of the physics program of the Relativistic Heavy-Ion Collider. The STAR collaboration has previously presented the results of a blind analysis of isobar collisions (${^{96}_{44}\text{Ru}}+{^{96}_{44}\text{Ru}}$, ${^{96}_{40}\text{Zr}}+{^{96}_{40}\text{Zr}}$) in the search for the CME. The isobar ratio ($Y$) of CME-sensitive observable, charge separation scaled by elliptic anisotropy, is close to but systematically larger than the inverse multiplicity ratio, the naive background baseline. This indicates the potential existence of a CME signal and the presence of remaining nonflow background due to two- and three-particle correlations, which are different between the isobars. In this post-blind analysis, we estimate the contributions from those nonflow correlations as a background baseline to $Y$, utilizing the isobar data as well as Heavy Ion Jet Interaction Generator simulations. This baseline is found consistent with the isobar ratio measurement, and an upper limit of 10% at 95% confidence level is extracted for the CME fraction in the charge separation measurement in isobar collisions at $\sqrt{s_{\rm NN}}=200$ GeV.
Figure 1a
Figure 1b
Figure 1c
The deconfined quark-gluon plasma (QGP) created in relativistic heavy-ion collisions enables the exploration of the fundamental properties of matter under extreme conditions. Non-central collisions can produce strong magnetic fields on the order of $10^{18}$ Gauss, which offers a probe into the electrical conductivity of the QGP. In particular, quarks and anti-quarks carry opposite charges and receive contrary electromagnetic forces that alter their momenta. This phenomenon can be manifested in the collective motion of final-state particles, specifically in the rapidity-odd directed flow, denoted as $v_1(\mathsf{y})$. Here we present the charge-dependent measurements of $dv_1/d\mathsf{y}$ near midrapidities for $\pi^{\pm}$, $K^{\pm}$, and $p(\bar{p})$ in Au+Au and isobar ($_{44}^{96}$Ru+$_{44}^{96}$Ru and $_{40}^{96}$Zr+$_{40}^{96}$Zr) collisions at $\sqrt{s_{\rm NN}}=$ 200 GeV, and in Au+Au collisions at 27 GeV, recorded by the STAR detector at the Relativistic Heavy Ion Collider. The combined dependence of the $v_1$ signal on collision system, particle species, and collision centrality can be qualitatively and semi-quantitatively understood as several effects on constituent quarks. While the results in central events can be explained by the $u$ and $d$ quarks transported from initial-state nuclei, those in peripheral events reveal the impacts of the electromagnetic field on the QGP. Our data put valuable constraints on the electrical conductivity of the QGP in theoretical calculations.
Directed flow of $p$ and $\bar{p}$ vs rapidity in Au+Au 200 GeV 50-80% centrality.
Directed flow of $p$ and $\bar{p}$ vs rapidity in Zr+Zr and Ru+Ru 200 GeV (combined) 50-80% centrality.
Directed flow of $p$ and $\bar{p}$ vs rapidity in Au+Au 27 GeV 50-80% centrality.
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}$.