The second- and third-order azimuthal anisotropy Fourier harmonics of charged particles produced in pPb collisions, at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV, are studied over a wide range of event multiplicities. Multiparticle correlations are used to isolate global properties stemming from the collision overlap geometry. The second-order "elliptic" harmonic moment is obtained with high precision through four-, six-, and eight-particle correlations and, for the first time, the third-order "triangular" harmonic moment is studied using four-particle correlations. A sample of peripheral PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV that covers a similar range of event multiplicities as the pPb results is also analyzed. Model calculations of initial-state fluctuations in pPb and PbPb collisions can be directly compared to the high precision experimental results. This work provides new insight into the fluctuation-driven origin of the $v_3$ coefficients in pPb and PbPb collisions, and into the dominating overall collision geometry in PbPb collisions at the earliest stages of heavy ion interactions.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_\mathrm{NN}} =$ 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

The elliptic flow coefficient ($v_{2}$) of identified particles in Pb-Pb collisions at $\sqrt{s_\mathrm{{NN}}} = 2.76$ TeV was measured with the ALICE detector at the LHC. The results were obtained with the Scalar Product method, a two-particle correlation technique, using a pseudo-rapidity gap of $|\Delta\eta| > 0.9$ between the identified hadron under study and the reference particles. The $v_2$ is reported for $\pi^{\pm}$, $\mathrm{K}^{\pm}$, $\mathrm{K}^0_\mathrm{S}$, p+$\overline{\mathrm{p}}$, $\mathrm{\phi}$, $\Lambda$+$\overline{\mathrm{\Lambda}}$, $\Xi^-$+$\overline{\Xi}^+$ and $\Omega^-$+$\overline{\Omega}^+$ in several collision centralities. In the low transverse momentum ($p_{\mathrm{T}}$) region, $p_{\mathrm{T}} < 2 $GeV/$c$, $v_2(p_\mathrm{T})$ exhibits a particle mass dependence consistent with elliptic flow accompanied by the transverse radial expansion of the system with a common velocity field. The experimental data for $\pi^{\pm}$ and $\mathrm{K}$ are described fairly well by hydrodynamical calculations coupled to a hadronic cascade model (VISHNU) for central collisions. However, the same calculations fail to reproduce the $v_2(p_\mathrm{T})$ for p+$\overline{\mathrm{p}}$, $\mathrm{\phi}$, $\Lambda$+$\overline{\mathrm{\Lambda}}$ and $\Xi^-$+$\overline{\Xi}^+$. For transverse momentum values larger than about 3 GeV/$c$, particles tend to group according to their type, i.e. mesons and baryons. However, the experimental data at the LHC exhibit deviations from the number of constituent quark (NCQ) scaling at the level of $\pm$20$\%$ for $p_{\mathrm{T}} > 3 $GeV/$c$.

Event-by-event fluctuations of the mean transverse momentum of charged particles produced in pp collisions at $\sqrt{s}$ = 0.9, 2.76 and 7 TeV, and Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV are studied as a function of the charged-particle multiplicity using the ALICE detector at the LHC. Dynamical fluctuations indicative of correlated particle emission are observed in all systems. The results in pp collisions show little dependence on collision energy. The Monte Carlo event generators PYTHIA and PHOJET are in qualitative agreement with the data. Peripheral Pb-Pb data exhibit a similar multiplicity dependence as that observed in pp. In central Pb-Pb, the results deviate from this trend, featuring a significant reduction of the fluctuation strength. The results in Pb--Pb are in qualitative agreement with previous measurements in Au-Au at lower collision energies and with expectations from models that incorporate collective phenomena.

Multi-particle cumulants and corresponding Fourier harmonics are measured for azimuthal angle distributions of charged particles in $pp$ collisions at $\sqrt{s}$ = 5.02 and 13 TeV and in $p$+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, and compared to the results obtained for low-multiplicity Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. These measurements aim to assess the collective nature of particle production. The measurements of multi-particle cumulants confirm the evidence for collective phenomena in $p$+Pb and low-multiplicity Pb+Pb collisions. On the other hand, the $pp$ results for four-particle cumulants do not demonstrate collective behaviour, indicating that they may be biased by contributions from non-flow correlations. A comparison of multi-particle cumulants and derived Fourier harmonics across different collision systems is presented as a function of the charged-particle multiplicity. For a given multiplicity, the measured Fourier harmonics are largest in Pb+Pb, smaller in $p$+Pb and smallest in $pp$ collisions. The $pp$ results show no dependence on the collision energy, nor on the multiplicity.

The second and the third order anisotropic flow, $V_{2}$ and $V_3$, are mostly determined by the corresponding initial spatial anisotropy coefficients, $\varepsilon_{2}$ and $\varepsilon_{3}$, in the initial density distribution. In addition to their dependence on the same order initial anisotropy coefficient, higher order anisotropic flow, $V_n$ ($n > 3$), can also have a significant contribution from lower order initial anisotropy coefficients, which leads to mode-coupling effects. In this Letter we investigate the linear and non-linear modes in higher order anisotropic flow $V_n$ for $n=4$, $5$, $6$ with the ALICE detector at the Large Hadron Collider. The measurements are done for particles in the pseudorapidity range $|\eta| < 0.8$ and the transverse momentum range $0.2 < p_{\rm T} < 5.0$ GeV/$c$ as a function of collision centrality. The results are compared with theoretical calculations and provide important constraints on the initial conditions, including initial spatial geometry and its fluctuations, as well as the ratio of the shear viscosity to entropy density of the produced system.

The elliptic, triangular, quadrangular and pentagonal anisotropic flow coefficients for $\pi^{\pm}$, $\mathrm{K}^{\pm}$ and p+$\overline{\mathrm{p}}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{{NN}}} = 2.76$ TeV were measured with the ALICE detector at the Large Hadron Collider. The results were obtained with the Scalar Product method, correlating the identified hadrons with reference particles from a different pseudorapidity region. Effects not related to the common event symmetry planes (non-flow) were estimated using correlations in pp collisions and were subtracted from the measurement. The obtained flow coefficients exhibit a clear mass ordering for transverse momentum ($p_{\mathrm{T}}$) values below $\approx$ 3 GeV/$c$. In the intermediate $p_{\mathrm{T}}$ region ($3 < p_{\mathrm{T}} < 6$ GeV/$c$), particles group at an approximate level according to the number of constituent quarks, suggesting that coalescence might be the relevant particle production mechanism in this region. The results for $p_{\mathrm{T}} < 3$ GeV/$c$ are described fairly well by a hydrodynamical model (iEBE-VISHNU) that uses initial conditions generated by A Multi-Phase Transport model (AMPT) and describes the expansion of the fireball using a value of 0.08 for the ratio of shear viscosity to entropy density ($\eta/s$), coupled to a hadronic cascade model (UrQMD). Finally, expectations from AMPT alone fail to quantitatively describe the measurements for all harmonics throughout the measured transverse momentum region. However, the comparison to the AMPT model highlights the importance of the late hadronic rescattering stage to the development of the observed mass ordering at low values of $p_{\mathrm{T}}$ and of coalescence as a particle production mechanism for the particle type grouping at intermediate values of $p_{\mathrm{T}}$ for all harmonics.

The measurement of azimuthal correlations of charged particles is presented for Pb-Pb collisions at $\sqrt{s_{\rm NN}}=$ 2.76 TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}=$ 5.02 TeV with the ALICE detector at the CERN Large Hadron Collider. These correlations are measured for the second, third and fourth order flow vector in the pseudorapidity region $|\eta|<0.8$ as a function of centrality and transverse momentum $p_{\rm T}$ using two observables, to search for evidence of $p_{\rm T}$-dependent flow vector fluctuations. For Pb-Pb collisions at 2.76 TeV, the measurements indicate that $p_{\rm T}$-dependent fluctuations are only present for the second order flow vector. Similar results have been found for p-Pb collisions at 5.02 TeV. These measurements are compared to hydrodynamic model calculations with event-by-event geometry fluctuations in the initial state to constrain the initial conditions and transport properties of the matter created in Pb-Pb and p-Pb collisions.

The elliptic flow of inclusive and direct photons was measured at mid-rapidity in two centrality classes 0-20% and 20-40% in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV by ALICE. Photons were detected with the highly segmented electromagnetic calorimeter PHOS and via conversions in the detector material with the $e^{+}e^{-}$ pairs reconstructed in the central tracking system. The results of the two methods were combined and the direct photon elliptic flow was extracted in the transverse momentum range $0.9 < p_{\rm T} < 6.2$ GeV/$c$. A comparison to RHIC data shows a similar magnitude of the measured direct-photon elliptic flow. Hydrodynamic and transport model calculations are systematically lower than the data, but are found to be compatible.

Measurements of anisotropic flow coefficients ($v_n$) and their cross-correlations using two- and multi-particle cumulant methods are reported in collisions of pp at $\sqrt{s} = 13$ TeV, p-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV, Xe-Xe at $\sqrt{s_{_{\rm NN}}} = 5.44$ TeV, and Pb-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV recorded with the ALICE detector. The multiplicity dependence of $v_n$ is studied in a very wide range from 20 to 3000 particles produced in the mid-rapidity region $|\eta|<0.8$ for the transverse momentum range $0.2 < p_{\rm T} < 3.0$ GeV/$c$. An ordering of the coefficients $v_2 > v_3 > v_4$ is found in pp and p-Pb collisions, similar to that seen in large collision systems, while a weak $v_2$ multiplicity dependence is observed relative to nucleus-nucleus collisions in the same multiplicity range. Using a novel subevent method, $v_{2}$ measured with four-particle cumulants is found to be compatible with that from six-particle cumulants in pp and p-Pb collisions. The magnitude of the correlation between $v_n^2$ and $v_m^2$, evaluated with the symmetric cumulants SC$(m,n)$ is observed to be positive at all multiplicities for $v_2$ and $v_4$, while for $v_2$ and $v_3$ it is negative and changes sign for multiplicities below 100, which may indicate a different $v_n$ fluctuation pattern in this multiplicity range. The observed long-range multi-particle azimuthal correlations in high multiplicity pp and p-Pb collisions can neither be described by PYTHIA 8 nor by IP-Glasma+MUSIC+UrQMD model calculations, and hence provide new insights into the understanding of collective effects in small collision systems.