We report on the first measurement of the triangular $v_3$, quadrangular $v_4$, and pentagonal $v_5$ charged particle flow in Pb-Pb collisions at 2.76 TeV measured with the ALICE detector at the CERN Large Hadron Collider. We show that the triangular flow can be described in terms of the initial spatial anisotropy and its fluctuations, which provides strong constraints on its origin. In the most central events, where the elliptic flow $v_2$ and $v_3$ have similar magnitude, a double peaked structure in the two-particle azimuthal correlations is observed, which is often interpreted as a Mach cone response to fast partons. We show that this structure can be naturally explained from the measured anisotropic flow Fourier coefficients.

A measurement of the transverse momentum spectra of jets in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV is reported. Jets are reconstructed from charged particles using the anti-$k_{\rm T}$ jet algorithm with jet resolution parameters $R$ of $0.2$ and $0.3$ in pseudo-rapidity $|\eta|<0.5$. The transverse momentum $p_{\rm T}$ of charged particles is measured down to $0.15$ GeV/$c$ which gives access to the low $p_{\rm T}$ fragments of the jet. Jets found in heavy-ion collisions are corrected event-by-event for average background density and on an inclusive basis (via unfolding) for residual background fluctuations and detector effects. A strong suppression of jet production in central events with respect to peripheral events is observed. The suppression is found to be similar to the suppression of charged hadrons, which suggests that substantial energy is radiated at angles larger than the jet resolution parameter $R=0.3$ considered in the analysis. The fragmentation bias introduced by selecting jets with a high $p_{\rm T}$ leading particle, which rejects jets with a soft fragmentation pattern, has a similar effect on the jet yield for central and peripheral events. The ratio of jet spectra with $R=0.2$ and $R=0.3$ is found to be similar in Pb-Pb and simulated PYTHIA pp events, indicating no strong broadening of the radial jet structure in the reconstructed jets with $R<0.3$.

Measurements of multi-particle azimuthal correlations (cumulants) for charged particles in p-Pb and Pb-Pb collisions are presented. They help address the question of whether there is evidence for global, flow-like, azimuthal correlations in the p-Pb system. Comparisons are made to measurements from the larger Pb-Pb system, where such evidence is established. In particular, the second harmonic two-particle cumulants are found to decrease with multiplicity, characteristic of a dominance of few-particle correlations in p-Pb collisions. However, when a $|\Delta \eta|$ gap is placed to suppress such correlations, the two-particle cumulants begin to rise at high-multiplicity, indicating the presence of global azimuthal correlations. The Pb-Pb values are higher than the p-Pb values at similar multiplicities. In both systems, the second harmonic four-particle cumulants exhibit a transition from positive to negative values when the multiplicity increases. The negative values allow for a measurement of $v_{2}\{4\}$ to be made, which is found to be higher in Pb-Pb collisions at similar multiplicities. The second harmonic six-particle cumulants are also found to be higher in Pb-Pb collisions. In Pb-Pb collisions, we generally find $v_{2}\{4\} \simeq v_{2}\{6\}\neq 0$ which is indicative of a Bessel-Gaussian function for the $v_{2}$ distribution. For very high-multiplicity Pb-Pb collisions, we observe that the four- and six-particle cumulants become consistent with 0. Finally, third harmonic two-particle cumulants in p-Pb and Pb-Pb are measured. These are found to be similar for overlapping multiplicities, when a $|\Delta\eta| > 1.4$ gap is placed.

We report on results obtained with the Event Shape Engineering technique applied to Pb-Pb collisions at $\sqrt{s_\rm{NN}}=2.76$ TeV. By selecting events in the same centrality interval, but with very different average flow, different initial state conditions can be studied. We find the effect of the event-shape selection on the elliptic flow coefficient $v_2$ to be almost independent of transverse momentum $p_\rm{T}$, as expected if this effect is due to fluctuations in the initial geometry of the system. Charged hadron, pion, kaon, and proton transverse momentum distributions are found to be harder in events with higher-than-average elliptic flow, indicating an interplay between radial and elliptic flow.

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.

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 $

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.

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.

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.

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.