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.

Two-particle angular correlations between trigger particles in the forward pseudorapidity range ($2.5 < |\eta| < 4.0$) and associated particles in the central range ($|\eta| < 1.0$) are measured with the ALICE detector in p-Pb collisions at a nucleon-nucleon centre-of-mass energy of 5.02 TeV. The trigger particles are reconstructed using the muon spectrometer, and the associated particles by the central barrel tracking detectors. In high-multiplicity events, the double-ridge structure, previously discovered in two-particle angular correlations at midrapidity, is found to persist to the pseudorapidity ranges studied in this Letter. The second-order Fourier coefficients for muons in high-multiplicity events are extracted after jet-like correlations from low-multiplicity events have been subtracted. The coefficients are found to have a similar transverse momentum ($p_{\rm T}$) dependence in p-going (p-Pb) and Pb-going (Pb-p) configurations, with the Pb-going coefficients larger by about $16\pm6$%, rather independent of $p_{\rm T}$ within the uncertainties of the measurement. The data are compared with calculations using the AMPT model, which predicts a different $p_{\rm T}$ and $\eta$ dependence than observed in the data. The results are sensitive to the parent particle $v_2$ and composition of reconstructed muon tracks, where the contribution from heavy flavour decays are expected to dominate at $p_{\rm T}>2$ GeV/$c$.

We report the first results of elliptic ($v_2$), triangular ($v_3$) and quadrangular flow ($v_4$) of charged particles in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE detector at the CERN Large Hadron Collider. The measurements are performed in the central pseudorapidity region $|\eta|<0.8$ and for the transverse momentum range $0.2<p_{\rm T}<5$ GeV/$c$. The anisotropic flow is measured using two-particle correlations with a pseudorapidity gap greater than one unit and with the multi-particle cumulant method. Compared to results from Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV, the anisotropic flow coefficients $v_{2}$, $v_{3}$ and $v_{4}$ are found to increase by ($3.0\pm0.6$)%, ($4.3\pm1.4$)% and ($10.2\pm3.8$)%, respectively, in the centrality range 0-50%. This increase can be attributed mostly to an increase of the average transverse momentum between the two energies. The measurements are found to be compatible with hydrodynamic model calculations. This comparison provides a unique opportunity to test the validity of the hydrodynamic picture and the power to further discriminate between various possibilities for the temperature dependence of shear viscosity to entropy density ratio of the produced matter in heavy-ion collisions at the highest energies.

The inclusive J/$\psi$ elliptic ($v_2$) and triangular ($v_3$) flow coefficients measured at forward rapidity (2.5 $<y<$ 4) and the $v_2$ measured at midrapidity ($|y|<$ 0.9) in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV using the ALICE detector at the LHC are reported. The entire Pb-Pb data sample collected during Run 2 is employed, amounting to an integrated luminosity of 750 $\mu$b$^{-1}$ at forward rapidity and 93 $\mu$b$^{-1}$ at midrapidity. The results are obtained using the scalar product method and are reported as a function of transverse momentum $p_{\rm T}$ and collision centrality. At midrapidity, the J/$\psi$ $v_2$ is in agreement with the forward rapidity measurement. The centrality averaged results indicate a positive J/$\psi$ $v_3$ with a significance of more than 5$\sigma$ at forward rapidity in the $p_{\rm T}$ range $2<p_{\rm T}<5$ GeV/$c$. The forward rapidity $v_2$, $v_3$, and $v_3$/$v_2$ results at low and intermediate $p_{\rm T}$ ($p_{\rm T} \lesssim 8$ GeV/$c$) exhibit a mass hierarchy when compared to pions and D mesons, while converging into a species-independent curve at higher $p_{\rm T}$. At low and intermediate $p_{\rm T}$, the results could be interpreted in terms of a later thermalization of charm quarks compared to light quarks, while at high $p_{\rm T}$, path-length dependent effects seem to dominate. The J/$\psi$ $v_2$ measurements are further compared to a microscopic transport model calculation. Using a simplified extension of the quark scaling approach involving both light and charm quark flow components, it is shown that the D-meson $v_{\rm n}$ measurements can be described based on those for charged pions and J/$\psi$ flow.

The azimuthal anisotropy coefficient $v_2$ of prompt D$^0$, D$^+$, D$^{*+}$ and D$_s^+$ mesons was measured in mid-central (30-50% centrality class) Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV, with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decays at mid-rapidity, $|y|<0.8$, in the transverse momentum interval $1<p_{\rm T}<24$ GeV/$c$. The measured D-meson $v_2$ has similar values as that of charged pions. The D$_s^+$ $v_2$, measured for the first time, is found to be compatible with that of non-strange D mesons. The measurements are compared with theoretical calculations of charm-quark transport in a hydrodynamically expanding medium and have the potential to constrain medium parameters.

The elliptic and triangular flow coefficients $v_2$ and $v_3$ of prompt D$^{0}$, D$^{+}$, and D$^{*+}$ mesons were measured at midrapidity ($|y|<0.8$) in Pb-Pb collisions at the centre-of-mass energy per nucleon pair of $\sqrt{s_{NN}} = 5.02$ TeV with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decays in the transverse momentum interval $1< p_{\rm T}<36$ GeV/$c$ in central (0-10%) and semi-central (30-50%) collisions. Compared to pions, protons, and J/$\psi$ mesons, the average D-meson $v_{n}$ harmonics are compatible within uncertainties with a mass hierarchy for $p_{\rm T} \lesssim 3$ GeV/$c$, and are similar to those of charged pions for higher $p_{\rm T}$. The coupling of the charm quark to the light quarks in the underlying medium is further investigated with the application of the event-shape engineering (ESE) technique to the D-meson $v_2$ and $p_{\rm T}$-differential yields. The D-meson $v_2$ is correlated with average bulk elliptic flow in both central and semi-central collisions. Within the current precision, the ratios of per-event D-meson yields in the ESE-selected and unbiased samples are found to be compatible with unity. All the measurements are found to be reasonably well described by theoretical calculations including the effects of charm-quark transport and the recombination of charm quarks with light quarks in a hydrodynamically expanding medium.

The elliptic flow ($v_2$) of $\rm D^0$ mesons from beauty-hadron decays (non-prompt $\rm D^0$) was measured in midcentral (30-50%) Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC. The $\rm D^0$ mesons were reconstructed at midrapidity ($|y|<0.8$) from their hadronic decay $\mathrm{D^0 \to K^-\pi^+}$, in the transverse momentum interval $2 < p_{\rm T} < 12$ GeV/$c$. The result indicates a positive $v_2$ for non-prompt $\rm D^0$ mesons with a significance of 2.7$\sigma$. The non-prompt $\rm D^0$-meson $v_2$ is lower than that of prompt non-strange D mesons with 3.2$\sigma$ significance in $2 < p_{\rm T} < 8$ GeV/$c$, and compatible with the $v_2$ of beauty-decay electrons. Theoretical calculations of beauty-quark transport in a hydrodynamically expanding medium describe the measurement within uncertainties.

The production yield of prompt D mesons and their elliptic flow coefficient $v_2$ were measured with the Event-Shape Engineering (ESE) technique applied to mid-central (10-30% and 30-50% centrality classes) Pb-Pb collisions at the centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} =5.02$ TeV, with the ALICE detector at the LHC. The ESE technique allows the classification of events, belonging to the same centrality, according to the azimuthal anisotropy of soft particle production in the collision. The reported measurements give the opportunity to investigate the dynamics of charm quarks in the Quark-Gluon Plasma and provide information on their participation in the collective expansion of the medium. D mesons were reconstructed via their hadronic decays at mid-rapidity, $|\eta|<0.8$, in the transverse momentum interval $1<p_{\rm T}<24$ GeV/$c$. The $v_2$ coefficient is found to be sensitive to the event-shape selection confirming a correlation between the D-meson azimuthal anisotropy and the collective expansion of the bulk matter, while the per-event D-meson yields do not show any significant modification within the current uncertainties.

The production yield and angular anisotropy of prompt ${D_s^+}$ mesons were measured as a function of transverse momentum ($p_{ T}$) in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{ NN}} = 5.02$ TeV collected with the ALICE detector at the LHC. ${D_s^+}$ mesons and their charge conjugates were reconstructed at midrapidity ($|y|<0.5$) from their hadronic decay channel ${D_s^+ \to \phi \pi^+}$, with ${\phi \to K^-K^+}$, in the $p_{ T}$ intervals $2<p_{ T}<50$ GeV/$c$ and $2<p_{ T}<36$ GeV/$c$ for the 0-10% and 30-50% centrality intervals. For $p_{ T}>10$ GeV/$c$, the measured ${D_s^+}$-meson nuclear modification factor $R_{ AA}$ is consistent with the one of non-strange D mesons within uncertainties, while at lower $p_{ T}$ a hint for a ${D_s^+}$-meson $R_{ AA}$ larger than that of non-strange D mesons is seen. The enhanced production of ${D_s^+}$ relative to non-strange D mesons is also studied by comparing the $p_{ T}$-dependent ${D_s^+/D^0}$ production yield ratios in Pb-Pb and in pp collisions. The ratio measured in Pb-Pb collisions is found to be on average higher than that in pp collisions in the interval $2<p_{ T} <8$ GeV/$c$ with a significance of 2.3$\sigma$ and 2.4$\sigma$ for the 0-10% and 30-50% centrality intervals. The azimuthal anisotropy coefficient $v_2$ of prompt ${D_s^+}$ mesons was measured in Pb-Pb collisions in the 30-50% centrality interval and is found to be compatible with that of non-strange D mesons. The main features of the measured $R_{ AA}$, ${D_s^+/D^0}$ ratio, and $v_2$ as a function of $p_{ T}$ are described by theoretical calculations of charm-quark transport in a hydrodynamically expanding quark-gluon plasma including hadronisation via charm-quark recombination with light quarks from the medium. The $p_{ T}$-integrated production yield of ${D_s^+}$ mesons is compatible with the prediction of the statistical hadronisation model.

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 $