The pseudorapidity dependence of elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of charged particles measured in Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}=5.02$ TeV and in Xe-Xe collisions at $\sqrt{s_{\rm NN}}=5.44$ TeV with ALICE at the LHC are presented. The measurements are performed in the pseudorapidity range $-3.5 < \eta < 5$ for various centrality intervals using two- and multi-particle cumulants with the subevent method. The flow probability density function (p.d.f.) is studied with the ratio of flow coefficient $v_2$ calculated with four- and two-particle cumulant, and suggests that the variance of flow p.d.f. is independent of pseudorapidity. The decorrelation of the flow vector in the longitudinal direction is probed using two-particle correlations. The results measured with respect to different reference regions in pseudorapidity exhibit differences, argued to be a result of saturating decorrelation effect above a certain pseudorapidity separation, in contrast to previous publications which assign this observation to non-flow effects. The results are compared to $3+1$ dimensional hydrodynamic and the AMPT transport model calculations. Neither of the models is able to simultaneously describe the pseudorapidity dependence of measurements of anisotropic flow and its fluctuations. The results presented in this work highlight shortcomings in our current understanding of initial conditions and subsequent system expansion in the longitudinal direction. Therefore, they provide input for its improvement.

Measurements of the suppression and correlations of dijets is performed using 3 $\mu$b$^{-1}$ of Xe+Xe data at $\sqrt{s_{\mathrm{NN}}} = 5.44$ TeV collected with the ATLAS detector at the LHC. Dijets with jets reconstructed using the $R=0.4$ anti-$k_t$ algorithm are measured differentially in jet $p_{\text{T}}$ over the range of 32 GeV to 398 GeV and the centrality of the collisions. Significant dijet momentum imbalance is found in the most central Xe+Xe collisions, which decreases in more peripheral collisions. Results from the measurement of per-pair normalized and absolutely normalized dijet $p_{\text{T}}$ balance are compared with previous Pb+Pb measurements at $\sqrt{s_{\mathrm{NN}}} =5.02$ TeV. The differences between the dijet suppression in Xe+Xe and Pb+Pb are further quantified by the ratio of pair nuclear-modification factors. The results are found to be consistent with those measured in Pb+Pb data when compared in classes of the same event activity and when taking into account the difference between the center-of-mass energies of the initial parton scattering process in Xe+Xe and Pb+Pb collisions. These results should provide input for a better understanding of the role of energy density, system size, path length, and fluctuations in the parton energy loss.

This paper presents measurements of charged-hadron spectra obtained in $pp$, $p$+Pb, and Pb+Pb collisions at $\sqrt{s}$ or $\sqrt{s_{_\text{NN}}}=5.02$ TeV, and in Xe+Xe collisions at $\sqrt{s_{_\text{NN}}}=5.44$ TeV. The data recorded by the ATLAS detector at the LHC have total integrated luminosities of 25 pb${}^{-1}$, 28 nb${}^{-1}$, 0.50 nb${}^{-1}$, and 3 $\mu$b${}^{-1}$, respectively. The nuclear modification factors $R_{p\text{Pb}}$ and $R_\text{AA}$ are obtained by comparing the spectra in heavy-ion and $pp$ collisions in a wide range of charged-particle transverse momenta and pseudorapidity. The nuclear modification factor $R_{p\text{Pb}}$ shows a moderate enhancement above unity with a maximum at $p_{\mathrm{T}} \approx 3$ GeV; the enhancement is stronger in the Pb-going direction. The nuclear modification factors in both Pb+Pb and Xe+Xe collisions feature a significant, centrality-dependent suppression. They show a similar distinct $p_{\mathrm{T}}$-dependence with a local maximum at $p_{\mathrm{T}} \approx 2$ GeV and a local minimum at $p_{\mathrm{T}} \approx 7$ GeV. This dependence is more distinguishable in more central collisions. No significant $|\eta|$-dependence is found. A comprehensive comparison with several theoretical predictions is also provided. They typically describe $R_\text{AA}$ better in central collisions and in the $p_{\mathrm{T}}$ range from about 10 to 100 GeV.

Multiplicity ($N_{\rm ch}$) distributions and transverse momentum ($p_{\rm T}$) spectra of inclusive primary charged particles in the kinematic range of $|\eta| < 0.8$ and 0.15 GeV/$c$$< p_{T} <$ 10 GeV/$c$ are reported for pp, p-Pb, Xe-Xe and Pb-Pb collisions at centre-of-mass energies per nucleon pair ranging from $\sqrt{s_{\rm NN}} = 2.76$ TeV up to $13$ TeV. A sequential two-dimensional unfolding procedure is used to extract the correlation between the transverse momentum of primary charged particles and the charged-particle multiplicity of the corresponding collision. This correlation sharply characterises important features of the final state of a collision and, therefore, can be used as a stringent test of theoretical models. The multiplicity distributions as well as the mean and standard deviation derived from the $p_{\rm T}$ spectra are compared to state-of-the-art model predictions. Providing these fundamental observables of bulk particle production consistently across a wide range of collision energies and system sizes can serve as an important input for tuning Monte Carlo event generators.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons $v_2${$\Psi_{\rm SP}$} in collisions of Pb ions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$=2.76 TeV and Xe ions at $\sqrt{s_{\rm NN}}$=5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum. The ratio between $v_2${$\Psi_{\rm SP}$} and the elliptic flow coefficient relative to the participant plane $v_2$4, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the TRENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models of initial state fluctuations. A significant transverse momentum dependence of the ratio $v_2${$\Psi_{\rm SP}$}/$v_2${4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of $v_2${$\Psi_{\rm SP}$} and $v_2${4} to the corresponding initial state eccentricities for Xe-Xe and Pb-Pb collisions at similar initial entropy density show a difference of $(7.0 \pm 0.9)$% with an additional variation of +1.8% when including RHIC data in the TRENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark-gluon plasma.

The production of J/$\psi$ is measured as a function of charged-particle multiplicity at forward rapidity in proton$-$proton (pp) collisions at center-of-mass energies $\sqrt{s} =$ 5.02 and 13 TeV. The J/$\psi$ mesons are reconstructed via their decay into dimuons in the rapidity interval (2.5 $< y <$ 4.0), whereas the charged-particle multiplicity density (${\rm d}N_{\rm{ch}}/{\rm d}\eta$) is measured at midrapidity $(|\eta| < 1)$. The production rate as a function of multiplicity is reported as the ratio of the yield in a given multiplicity interval to the multiplicity-integrated one. This observable shows a linear increase with charged-particle multiplicity normalized to the corresponding average value for inelastic events (${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$), at both the colliding energies. Measurements are compared with available ALICE results at midrapidity and theoretical model calculations. First measurement of the mean transverse momentum ($\langle p_{\mathrm{T}}\rangle$) of J/$\psi$ in pp collisions exhibits an increasing trend as a function of ${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$ showing a saturation towards high charged-particle multiplicities.

Correlations between mean transverse momentum $[p_{\rm T}]$ and anisotropic flow coefficients $v_{\rm 2}$ or $v_{\rm 3}$ are measured as a function of centrality in Pb$-$Pb and Xe$-$Xe collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and 5.44 TeV, respectively, with ALICE. In addition, the recently proposed higher-order correlation between $[p_{\rm T}]$, $v_{\rm 2}$, and $v_{\rm 3}$ is measured for the first time, which shows an anticorrelation for the presented centrality ranges. These measurements are compared with hydrodynamic calculations using IP-Glasma and $\rm T_{R}ENTo$ initial-state shapes, the former based on the Color Glass Condensate effective theory with gluon saturation, and the latter a parameterized model with nucleons as the relevant degrees of freedom. The data are better described by the IP-Glasma rather than the $\rm T_{R}ENTo$ based calculations. In particular, Trajectum and JETSCAPE predictions, both based on the $\rm T_{R}ENTo$ initial state model but with different parameter settings, fail to describe the measurements. As the correlations between $[p_{\rm T}]$ and $v_{\rm n}$ are mainly driven by the correlations of the size and the shape of the system in the initial state, these new studies pave a novel way to characterize the initial state and help pin down the uncertainty of the extracted properties of the quark$-$gluon plasma recreated in relativistic heavy-ion collisions.

Measurements of elliptic ($v_2$) and triangular ($v_3$) flow coefficients of $\pi^{\pm}$, K$^{\pm}$, p+$\rm \overline{p}$, K$^0_{\rm S}$, and $\Lambda + \overline{\Lambda}$ obtained with the scalar product method in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV are presented. The results are obtained in the rapidity range $\left | y \right |<0.5$ and reported as a function of transverse momentum, $p_{\rm T}$, for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for $p_{\rm T}<3$ GeV/$c$, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3< $p_{\rm T}$ <8 GeV/$c$). The magnitude of the baryon $v_{2}$ is larger than that of mesons up to $p_{\rm T}$ = 6 GeV/$c$. The centrality dependence of the shape evolution of the $p_{\rm T}$-differential $v_2$ is studied for the various hadron species. The $v_2$ coefficients of $\pi^{\pm}$, K$^{\pm}$, and p+$\rm \overline{p}$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for $p_{\rm T} <1$ GeV/$c$. A comparison with $v_{\rm n}$ measurements in the corresponding centrality intervals in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV yields an enhanced $v_2$ in central collisions and diminished value in semicentral collisions.

The first measurement of the production of pions, kaons, (anti-)protons and $\phi$ mesons at midrapidity in Xe-Xe collisions at $\sqrt{s_{\rm NN}} = 5.44$ TeV is presented. Transverse momentum ($p_{\rm T}$) spectra and $p_{\rm T}$-integrated yields are extracted in several centrality intervals bridging from p-Pb to mid-central Pb-Pb collisions in terms of final-state multiplicity. The study of Xe-Xe and Pb-Pb collisions allows systems at similar charged-particle multiplicities but with different initial geometrical eccentricities to be investigated. A detailed comparison of the spectral shapes in the two systems reveals an opposite behaviour for radial and elliptic flow. In particular, this study shows that the radial flow does not depend on the colliding system when compared at similar charged-particle multiplicity. In terms of hadron chemistry, the previously observed smooth evolution of particle ratios with multiplicity from small to large collision systems is also found to hold in Xe-Xe. In addition, our results confirm that two remarkable features of particle production at LHC energies are also valid in the collision of medium-sized nuclei: the lower proton-to-pion ratio with respect to the thermal model expectations and the increase of the $\phi$-to-pion ratio with increasing final-state multiplicity.