Measurements of anisotropic flow coefficients with two- and multi-particle cumulants for inclusive charged particles in Pb-Pb collisions at $\sqrt{{\textit s}_\text{NN}} = 5.02$ and 2.76 TeV are reported in the pseudorapidity range $|\eta| < 0.8$ and transverse momentum $0.2 < p_\text{T} < 50$ GeV/$c$. The full data sample collected by the ALICE detector in 2015 (2010), corresponding to an integrated luminosity of 12.7 (2.0) $\mu$b$^{-1}$ in the centrality range 0-80%, is analysed. Flow coefficients up to the sixth flow harmonic ($v_6$) are reported and a detailed comparison among results at the two energies is carried out. The $p_\text{T}$ dependence of anisotropic flow coefficients and its evolution with respect to centrality and harmonic number $n$ are investigated. An approximate power-law scaling of the form $v_n(p_\text{T}) \sim p_\text{T}^{n/3}$ is observed for all flow harmonics at low $p_\text{T}$ ($0.2 < p_\text{T} < 3$ GeV/$c$). At the same time, the ratios $v_n/v_m^{n/m}$ are observed to be essentially independent of $p_\text{T}$ for most centralities up to about $p_\text{T} = 10$ GeV/$c$. Analysing the differences among higher-order cumulants of elliptic flow ($v_2$), which have different sensitivities to flow fluctuations, a measurement of the standardised skewness of the event-by-event $v_2$ distribution $P(v_2)$ is reported and constraints on its higher moments are provided. The Elliptic Power distribution is used to parametrise $P(v_2)$, extracting its parameters from fits to cumulants. The measurements are compared to different model predictions in order to discriminate among initial-state models and to constrain the temperature dependence of the shear viscosity to entropy-density ratio.

Charged-particle spectra at midrapidity are measured in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV and presented in centrality classes ranging from most central (0-5%) to most peripheral (95-100%) collisions. Possible medium effects are quantified using the nuclear modification factor ($R_{\rm AA}$) by comparing the measured spectra with those from proton-proton collisions, scaled by the number of independent nucleon-nucleon collisions obtained from a Glauber model. At large transverse momenta ($8<p_{\rm T}<20$ GeV/$c$), the average $R_{\rm AA}$ is found to increase from about $0.15$ in 0-5% central to a maximum value of about $0.8$ in 75-85% peripheral collisions, beyond which it falls off strongly to below $0.2$ for the most peripheral collisions. Furthermore, $R_{\rm AA}$ initially exhibits a positive slope as a function of $p_{\rm T}$ in the $8$-$20$ GeV/$c$ interval, while for collisions beyond the 80% class the slope is negative. To reduce uncertainties related to event selection and normalization, we also provide the ratio of $R_{\rm AA}$ in adjacent centrality intervals. Our results in peripheral collisions are consistent with a PYTHIA-based model without nuclear modification, demonstrating that biases caused by the event selection and collision geometry can lead to the apparent suppression in peripheral collisions. This explains the unintuitive observation that $R_{\rm AA}$ is below unity in peripheral Pb-Pb, but equal to unity in minimum-bias p-Pb collisions despite similar charged-particle multiplicities.

The azimuthal variation of jet yields in heavy-ion collisions provides information about the path-length dependence of the energy loss experienced by partons passing through the hot, dense nuclear matter known as the quark-gluon plasma. This paper presents the azimuthal anisotropy coefficients $v_2$, $v_3$, and $v_4$ measured for jets in Pb+Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV using the ATLAS detector at the LHC. The measurement uses data collected in 2015 and 2018, corresponding to an integrated luminosity of 2.2 nb$^{-1}$. The $v_n$ values are measured as a function of the transverse momentum of the jets between 71 GeV and 398 GeV and the event centrality. A nonzero value of $v_2$ is observed in all but the most central collisions. The value of $v_2$ is largest for jets with lower transverse momentum, with values up to 0.05 in mid-central collisions. A smaller, nonzero value of $v_3$ of approximately 0.01 is measured with no significant dependence on jet $p_T$ or centrality, suggesting that fluctuations in the initial state play a small but distinct role in jet energy loss. No significant deviation of $v_4$ from zero is observed in the measured kinematic region.

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.

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.