Neutral pion transverse momentum (pT) spectra at mid-rapidity (|y| < 0.35) were measured in Cu+Cu collisions at \sqrt s_NN = 22.4, 62.4, and 200 GeV. Relative to pi -zero yields in p+p collisions scaled by the number of inelastic nucleon-nucleon collisions (Ncoll) at the respective energies, the pi-zero yields for pT \ge 2 GeV/c in central Cu+Cu collisions at 62.4 and 200 GeV are suppressed, whereas an enhancement is observed at 22.4 GeV. A comparison with a jet quenching model suggests that final state parton energy loss dominates in central Cu+Cu collisions at 62.4 GeV and 200 GeV, while the enhancement at 22.4 GeV is consistent with nuclear modifications in the initial state alone.

Detailed differential measurements of the elliptic flow for particles produced in Au+Au and Cu+Cu collisions at sqrt(s_NN) = 200 GeV are presented. Predictions from perfect fluid hydrodynamics for the scaling of the elliptic flow coefficient v_2 with eccentricity, system size and transverse energy are tested and validated. For transverse kinetic energies KE_T ~ m_T-m up to ~1 GeV, scaling compatible with the hydrodynamic expansion of a thermalized fluid is observed for all produced particles. For large values of KE_T, the mesons and baryons scale separately. A universal scaling for the flow of both mesons and baryons is observed for the full transverse kinetic energy range of the data when quark number scaling is employed. In both cases the scaling is more pronounced in terms of KE_T rather than transverse momentum.

We report on event structure and double helicity asymmetry ($A_LL$) of jet production in longitudinally polarized p+p collisions at $\sqrt{s}$=200 GeV. Photons and charged particles were measured at midrapidity $|\eta| < 0.35$ with the requirement of a high-momentum ($>2$ GeV/$c$) photon in each event. Measured event structure is compared with {\sc pythia} and {\sc geant} simulations. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet $A_{LL}$, photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster $p_T$ sum ($p_T^{\rm reco}$). The effect of detector response and the underlying events on $p_T^{\rm reco}$ was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the NLO pQCD jet production cross section. For $4 < p_T^{\rm reco} < 12$ GeV/$c$ with an average beam polarization of $< P > = 49%$ we measured $A_{LL} = -0.0014 \pm 0.0037^{\rm stat}$ at the lowest $p_T^{\rm reco}$ bin (4-5 GeV/$c$) and $-0.0181 \pm 0.0282^{\rm stat}$ at the highest $p_T^{\rm reco}$ bin (10-12 GeV/$c$) with a beam polarization scale error of 9.4% and a $\pT$ scale error of 10%. Jets in the measured $p_T^{\rm reco}$ range arise primarily from hard-scattered gluons with momentum fraction $0.02 < x < 0.3$ according to {\sc pythia}. The measured $A_{LL}$ is compared with predictions that assume various $\Delta G(x)$ distributions based on the GRSV parameterization. The present result imposes the limit $-1.1 < \int_{0.02}^{0.3}dx \Delta G(x, \mu^2 = 1 {\rm GeV}^2) < 0.4$ at 95% confidence level or $\int_{0.02}^{0.3}dx \Delta G(x, \mu^2 = 1 {\rm GeV}^2) < 0.5$ at 99% confidence level.

Back-to-back hadron pair yields in d+Au and p+p collisions at sqrt(s_NN)=200 GeV were measured with the PHENIX detector at the Relativistic Heavy Ion Collider. Rapidity separated hadron pairs were detected with the trigger hadron at pseudorapidity |eta|<0.35 and the associated hadron at forward rapidity (deuteron direction, 3.0<eta<3.8). Pairs were also detected with both hadrons measured at forward rapidity; in this case the yield of back-to-back hadron pairs in d+Au collisions with small impact parameters is observed to be suppressed by a factor of 10 relative to p+p collisions. The kinematics of these pairs is expected to probe partons in the Au nucleus with low fraction x of the nucleon momenta, where the gluon densities rise sharply. The observed suppression as a function of nuclear thickness, p_T, and eta points to cold nuclear matter effects arising at high parton densities.

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 present the first measurement of elliptic ($v_2$) and triangular ($v_3$) flow in high-multiplicity $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Two-particle correlations, where the particles have a large separation in pseudorapidity, are compared in $^{3}$He$+$Au and in $p$$+$$p$ collisions and indicate that collective effects dominate the second and third Fourier components for the correlations observed in the $^{3}$He$+$Au system. The collective behavior is quantified in terms of elliptic $v_2$ and triangular $v_3$ anisotropy coefficients measured with respect to their corresponding event planes. The $v_2$ values are comparable to those previously measured in $d$$+$Au collisions at the same nucleon-nucleon center-of-mass energy. Comparison with various theoretical predictions are made, including to models where the hot spots created by the impact of the three $^{3}$He nucleons on the Au nucleus expand hydrodynamically to generate the triangular flow. The agreement of these models with data may indicate the formation of low-viscosity quark-gluon plasma even in these small collision systems.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured $\phi$ meson production and its nuclear modification in asymmetric Cu$+$Au heavy-ion collisions at $\sqrt{s_{NN}}=200$ GeV at both forward Cu-going direction ($1.2<y<2.2$) and backward Au-going direction ($-2.2<y<-1.2$), rapidities. The measurements are performed via the dimuon decay channel and reported as a function of the number of participating nucleons, rapidity, and transverse momentum. In the most central events, 0\%--20\% centrality, the $\phi$ meson yield integrated over $1<p_T<5$ GeV/$c$ prefers a smaller value, which means a larger nuclear modification, in the Cu-going direction compared to the Au-going direction. Additionally, the nuclear-modification factor in Cu$+$Au collisions averaged over all centrality is measured to be similar to the previous PHENIX result in $d$$+$Au collisions for these rapidities.

The pseudorapidity density of charged particles with minimum transverse momentum ($p_{\rm T}$) thresholds of 0.15, 0.5, 1, and 2 GeV$/c$ is measured in pp collisions at the centre of mass energies of $\sqrt{s} =$ 5.02 and 13 TeV with the ALICE detector. The study is carried out for inelastic collisions with at least one primary charged particle having a pseudorapidity ($\eta$) within $\pm0.8$ and $p_{\rm T}$ larger than the corresponding threshold. In addition, measurements without $p_{\rm T}$-thresholds are performed for inelastic and non-single-diffractive events as well as for inelastic events with at least one charged particle having $|\eta|<1$ in pp collisions at $\sqrt{s} =$ 5.02 TeV for the first time at the LHC. These measurements are compared to the PYTHIA 6, PYTHIA 8, and EPOS-LHC models. In general, the models describe the $\eta$ dependence of particle production well. However, discrepancies are observed for the highest transverse momentum threshold ($p_{\rm T}>2 {\rm\ GeV}/c$), highlighting the importance of such measurements for tuning event generators. The new measurements agree within uncertainties with results from the ATLAS and CMS experiments obtained at $\sqrt{s} = 13$ TeV.

The first measurements of the production of muons and electrons from heavy-flavour hadron decays in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV, using the ALICE detector at the LHC, are reported. The measurement of the nuclear modification factor $R_{\rm AA}$ is performed as a function of transverse momentum $p_{\rm T}$ in several centrality classes at forward rapidity ($2.5 < y <4$) and midrapidity ($\vert y \vert < 0.8$) for muons and electrons from heavy-flavour hadron decays, respectively. A suppression by a factor up to about 2.5 compared to the binary-scaled pp reference is observed in central collisions at both central and forward rapidities. The $R_{\rm AA}$ of muons from heavy-flavour hadron decays is compared to previous measurements in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. When the nuclear modification factors are compared in the centrality classes 0-10% for Xe-Xe collisions and 10-20% for Pb-Pb collisions, which have similar charged-particle multiplicity density, a similar suppression, with $R_{\rm AA} \sim 0.4$ in the $p_{\rm T}$ interval $4 < p_{\rm T} < 8$ GeV/$c$, is observed. The comparison of the measured $R_{\rm AA}$ values in the two collision systems brings new insights on the properties of the quark-gluon plasma by investigating the system-size and path-length dependence of medium-induced parton energy loss. The results of muons and electrons from heavy-flavour hadron decays provide new constraints to model calculations.

The first measurements of elliptic flow of $\pi^\pm$, ${\rm K}^\pm$, p+$\overline{\rm p}$, ${\rm K_{S}^0}$, $\Lambda$+$\overline{\Lambda}$, $\phi$, $\Xi^-$+$\Xi^+$, and $\Omega^-$+$\Omega^+$ using multiparticle cumulants in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented. Results obtained with two- ($v_2\{2\}$) and four-particle cumulants ($v_2\{4\}$) are shown as a function of transverse momentum, $p_{\rm T}$, for various collision centrality intervals. Combining the data for both $v_2\{2\}$ and $v_2\{4\}$ also allows us to report the first measurements of the mean elliptic flow, elliptic flow fluctuations, and relative elliptic flow fluctuations for various hadron species. These observables probe the event-by-event eccentricity fluctuations in the initial state and the contributions from the dynamic evolution of the expanding quark-gluon plasma. The characteristic features observed in previous $p_{\rm T}$-differential anisotropic flow measurements for identified hadrons with two-particle correlations, namely the mass ordering at low $p_{\rm T}$ and the approximate scaling with the number of constituent quarks at intermediate $p_{\rm T}$, are similarly present in the four-particle correlations and the combinations of $v_2\{2\}$ and $v_2\{4\}$. In addition, a particle species dependence of flow fluctuations is observed that could indicate a significant contribution from final state hadronic interactions. The comparison between experimental measurements and CoLBT model calculations, which combine the various physics processes of hydrodynamics, quark coalescence, and jet fragmentation, illustrates their importance over a wide $p_{\rm T}$ range.