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 distributions of event-by-event harmonic flow coefficients v_n for n=2-4 are measured in sqrt(s_NN)=2.76 TeV Pb+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using charged particles with transverse momentum pT> 0.5 GeV and in the pseudorapidity range |eta|<2.5 in a dataset of approximately 7 ub^-1 recorded in 2010. The shapes of the v_n distributions are described by a two-dimensional Gaussian function for the underlying flow vector in central collisions for v_2 and over most of the measured centrality range for v_3 and v_4. Significant deviations from this function are observed for v_2 in mid-central and peripheral collisions, and a small deviation is observed for v_3 in mid-central collisions. It is shown that the commonly used multi-particle cumulants are insensitive to the deviations for v_2. The v_n distributions are also measured independently for charged particles with 0.5<pT<1 GeV and pT>1 GeV. When these distributions are rescaled to the same mean values, the adjusted shapes are found to be nearly the same for these two pT ranges. The v_n distributions are compared with the eccentricity distributions from two models for the initial collision geometry: a Glauber model and a model that includes corrections to the initial geometry due to gluon saturation effects. Both models fail to describe the experimental data consistently over most of the measured centrality range.

The integrated elliptic flow of charged particles produced in Pb+Pb collisions at sqrt(s_NN)=2.76 TeV has been measured with the ATLAS detector using data collected at the Large Hadron Collider. The anisotropy parameter, v_2, was measured in the pseudorapidity range |eta| <= 2.5 with the event-plane method. In order to include tracks with very low transverse momentum pT, thus reducing the uncertainty in v_2 integrated over pT, a 1 mu b-1 data sample without a magnetic field in the tracking detectors is used. The centrality dependence of the integrated v_2 is compared to other measurements obtained with higher pT thresholds. A weak pseudorapidity dependence of the integrated elliptic flow is observed for central collisions, and a small decrease when moving away from mid-rapidity is observed only in peripheral collisions. The integrated v2 transformed to the rest frame of one of the colliding nuclei is compared to the lower-energy RHIC data.

Correlations between the elliptic or triangular flow coefficients $v_m$ ($m$=2 or 3) and other flow harmonics $v_n$ ($n$=2 to 5) are measured using $\sqrt{s_{NN}}=2.76$ TeV Pb+Pb collision data collected in 2010 by the ATLAS experiment at the LHC, corresponding to an integrated lumonisity of 7 $\mu$b$^{-1}$. The $v_m$-$v_n$ correlations are measured in midrapidity as a function of centrality, and, for events within the same centrality interval, as a function of event ellipticity or triangularity defined in a forward rapidity region. For events within the same centrality interval, $v_3$ is found to be anticorrelated with $v_2$ and this anticorrelation is consistent with similar anticorrelations between the corresponding eccentricities $\epsilon_2$ and $\epsilon_3$. On the other hand, it is observed that $v_4$ increases strongly with $v_2$, and $v_5$ increases strongly with both $v_2$ and $v_3$. The trend and strength of the $v_m$-$v_n$ correlations for $n$=4 and 5 are found to disagree with $\epsilon_m$-$\epsilon_n$ correlations predicted by initial-geometry models. Instead, these correlations are found to be consistent with the combined effects of a linear contribution to $v_n$ and a nonlinear term that is a function of $v_2^2$ or of $v_2v_3$, as predicted by hydrodynamic models. A simple two-component fit is used to separate these two contributions. The extracted linear and nonlinear contributions to $v_4$ and $v_5$ are found to be consistent with previously measured event-plane correlations.

Charged-particle spectra obtained in 0.15 nb${}^{-1}$ of Pb+Pb interactions at $\sqrt{{s}_\mathsf{{NN}}}=2.76$TeV and 4.2 pb${}^{-1}$ of pp interactions at $\sqrt{s}=2.76$ TeV with the ATLAS detector at the LHC are presented in a wide transverse momentum ($0.5 < p_{\mathrm{T}} < 150$ GeV) and pseudorapidity ($|\eta|<2$) range. For Pb+Pb collisions, the spectra are presented as a function of collision centrality, which is determined by the response of the forward calorimeter located on both sides of the interaction point. The nuclear modification factors $R_{\mathrm{AA}}$ and $R_{\mathrm{CP}}$ are presented in detail as function of centrality, $p_{\mathrm{T}}$ and $\eta$. They show a distinct $p_{\mathrm{T}}$-dependence with a pronounced minimum at about 7 GeV. Above 60 GeV, $R_{\mathrm{AA}}$ is consistent with a plateau at a centrality-dependent value, within the uncertainties. The value is $0.55\pm0.01(stat.)\pm0.04(syst.)$ in the most central collisions. The $R_{\mathrm{AA}}$ distribution is consistent with flat $|\eta|$ dependence over the whole transverse momentum range in all centrality classes.

Measurements of the azimuthal anisotropy in lead-lead collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV are presented using a data sample corresponding to 0.49 $\mathrm{nb}^{-1}$ integrated luminosity collected by the ATLAS experiment at the LHC in 2015. The recorded minimum-bias sample is enhanced by triggers for "ultra-central" collisions, providing an opportunity to perform detailed study of flow harmonics in the regime where the initial state is dominated by fluctuations. The anisotropy of the charged-particle azimuthal angle distributions is characterized by the Fourier coefficients, $v_{2}-v_{7}$, which are measured using the two-particle correlation, scalar-product and event-plane methods. The goal of the paper is to provide measurements of the differential as well as integrated flow harmonics $v_{n}$ over wide ranges of the transverse momentum, 0.5 $

Measurements of two-particle correlation functions and the first five azimuthal harmonics, $v_1$ to $v_5$, are presented, using 28 $\mathrm{nb}^{-1}$ of $p$+Pb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV measured with the ATLAS detector at the LHC. Significant long-range "ridge-like" correlations are observed for pairs with small relative azimuthal angle ($|\Delta\phi|<\pi/3$) and back-to-back pairs ($|\Delta\phi| > 2\pi/3$) over the transverse momentum range $0.4 < p_{\rm T} < 12$ GeV and in different intervals of event activity. The event activity is defined by either the number of reconstructed tracks or the total transverse energy on the Pb-fragmentation side. The azimuthal structure of such long-range correlations is Fourier decomposed to obtain the harmonics $v_n$ as a function of $p_{\rm T}$ and event activity. The extracted $v_n$ values for $n=2$ to 5 decrease with $n$. The $v_2$ and $v_3$ values are found to be positive in the measured $p_{\rm T}$ range. The $v_1$ is also measured as a function of $p_{\rm T}$ and is observed to change sign around $p_{\rm T}\approx 1.5$-2.0 GeV and then increase to about 0.1 for $p_{\rm T}>4$ GeV. The $v_2(p_{\rm T})$, $v_3(p_{\rm T})$ and $v_4(p_{\rm T})$ are compared to the $v_n$ coefficients in Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}} =2.76$ TeV with similar event multiplicities. Reasonable agreement is observed after accounting for the difference in the average $p_{\rm T}$ of particles produced in the two collision systems.

ATLAS measurements of the azimuthal anisotropy in lead-lead collisions at $\sqrt{s_{NN}}=2.76$ TeV are shown using a dataset of approximately 7 $\mu$b$^{-1}$ collected at the LHC in 2010. The measurements are performed for charged particles with transverse momenta $0.5<p_T<20$ GeV and in the pseudorapidity range $|\eta|<2.5$. The anisotropy is characterized by the Fourier coefficients, $v_n$, of the charged-particle azimuthal angle distribution for n = 2-4. The Fourier coefficients are evaluated using multi-particle cumulants calculated with the generating function method. Results on the transverse momentum, pseudorapidity and centrality dependence of the $v_n$ coefficients are presented. The elliptic flow, $v_2$, is obtained from the two-, four-, six- and eight-particle cumulants while higher-order coefficients, $v_3$ and $v_4$, are determined with two- and four-particle cumulants. Flow harmonics $v_n$ measured with four-particle cumulants are significantly reduced compared to the measurement involving two-particle cumulants. A comparison to $v_n$ measurements obtained using different analysis methods and previously reported by the LHC experiments is also shown. Results of measurements of flow fluctuations evaluated with multi-particle cumulants are shown as a function of transverse momentum and the collision centrality. Models of the initial spatial geometry and its fluctuations fail to describe the flow fluctuations measurements.

The anisotropy of the azimuthal distributions of charged particles produced in PbPb collisions with a nucleon-nucleon center-of-mass energy of 2.76 TeV is studied with the CMS experiment at the LHC. The elliptic anisotropy parameter defined as the second coefficient in a Fourier expansion of the particle invariant yields, is extracted using the event-plane method, two- and four-particle cumulants, and Lee--Yang zeros. The anisotropy is presented as a function of transverse momentum (pt), pseudorapidity (eta) over a broad kinematic range: 0.3 < pt < 20 GeV, abs(eta) < 2.4, and in 12 classes of collision centrality from 0 to 80%. The results are compared to those obtained at lower center-of-mass energies, and various scaling behaviors are examined. When scaled by the geometric eccentricity of the collision zone, the elliptic anisotropy is found to obey a universal scaling with the transverse particle density for different collision systems and center-of-mass energies.

Differential measurements of charged particle azimuthal anisotropy are presented for lead-lead collisions at sqrt(s_NN) = 2.76 TeV with the ATLAS detector at the LHC, based on an integrated luminosity of approximately 8 mb^-1. This anisotropy is characterized via a Fourier expansion of the distribution of charged particles in azimuthal angle (phi), with the coefficients v_n denoting the magnitude of the anisotropy. Significant v_2-v_6 values are obtained as a function of transverse momentum (0.5<pT<20 GeV), pseudorapidity (|eta|<2.5) and centrality using an event plane method. The v_n values for n>=3 are found to vary weakly with both eta and centrality, and their pT dependencies are found to follow an approximate scaling relation, v_n^{1/n}(pT) \propto v_2^{1/2}(pT). A Fourier analysis of the charged particle pair distribution in relative azimuthal angle (Dphi=phi_a-phi_b) is performed to extract the coefficients v_{n,n}=<cos (n Dphi)>. For pairs of charged particles with a large pseudorapidity gap (|Deta=eta_a-eta_b|>2) and one particle with pT<3 GeV, the v_{2,2}-v_{6,6} values are found to factorize as v_{n,n}(pT^a,pT^b) ~ v_n(pT^a)v_n(pT^b) in central and mid-central events. Such factorization suggests that these values of v_{2,2}-v_{6,6} are primarily due to the response of the created matter to the fluctuations in the geometry of the initial state. A detailed study shows that the v_{1,1}(pT^a,pT^b) data are consistent with the combined contributions from a rapidity-even v_1 and global momentum conservation. A two-component fit is used to extract the v_1 contribution. The extracted v_1 is observed to cross zero at pT\sim1.0 GeV, reaches a maximum at 4-5 GeV with a value comparable to that for v_3, and decreases at higher pT.