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.

Azimuthally anisotropic distributions of D$^0$, D$^+$ and D$^{*+}$ mesons were studied in the central rapidity region ($|y|<0.8$) in Pb-Pb collisions at a centre-of-mass energy $\sqrt{s_{\rm NN}} = 2.76$ TeV per nucleon-nucleon collision, with the ALICE detector at the LHC. The second Fourier coefficient $v_2$ (commonly denoted elliptic flow) was measured in the centrality class 30-50% as a function of the D meson transverse momentum $p_{\rm T}$, in the range 2-16 GeV/$c$. The measured $v_2$ of D mesons is comparable in magnitude to that of light-flavour hadrons. It is positive in the range $2 < p_{\rm T} < 6$ GeV/$c$ with $5.7\sigma$ significance, based on the combination of statistical and systematic uncertainties.

We report the first measurement of charged particle elliptic flow in Pb-Pb collisions at 2.76 TeV with the ALICE detector at the CERN Large Hadron Collider. The measurement is performed in the central pseudorapidity region (|$\eta$|<0.8) and transverse momentum range 0.2< $p_{\rm T}$< 5.0 GeV/$c$. The elliptic flow signal v$_2$, measured using the 4-particle correlation method, averaged over transverse momentum and pseudorapidity is 0.087 $\pm$ 0.002 (stat) $\pm$ 0.004 (syst) in the 40-50% centrality class. The differential elliptic flow v$_2(p_{\rm T})$ reaches a maximum of 0.2 near $p_{\rm T}$ = 3 GeV/$c$. Compared to RHIC Au-Au collisions at 200 GeV, the elliptic flow increases by about 30%. Some hydrodynamic model predictions which include viscous corrections are in agreement with the observed increase.

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons $v_2${$\Psi_{\rm SP}$} in collisions of Pb ions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$=2.76 TeV and Xe ions at $\sqrt{s_{\rm NN}}$=5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum. The ratio between $v_2${$\Psi_{\rm SP}$} and the elliptic flow coefficient relative to the participant plane $v_2$4, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the TRENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models of initial state fluctuations. A significant transverse momentum dependence of the ratio $v_2${$\Psi_{\rm SP}$}/$v_2${4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of $v_2${$\Psi_{\rm SP}$} and $v_2${4} to the corresponding initial state eccentricities for Xe-Xe and Pb-Pb collisions at similar initial entropy density show a difference of $(7.0 \pm 0.9)$% with an additional variation of +1.8% when including RHIC data in the TRENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark-gluon plasma.

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 the first measurement of inclusive J/$\psi$ elliptic flow, $v_2$, in heavy-ion collisions at the LHC. The measurement is performed with the ALICE detector in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV in the rapidity range $2.5 < y < 4.0$. The dependence of the J/$\psi$ $v_2$ on the collision centrality and on the J/$\psi$ transverse momentum is studied in the range $0 < p_{\rm T} < 10$ GeV/$c$. For semi-central Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV, an indication of non-zero $v_2$ is observed with a maximum value of $v_2 = 0.116 \pm 0.046 {\rm (stat.)} \pm 0.029 {\rm (syst.)}$ for J/$\psi$ in the transverse momentum range $2 < p_{\rm T} < 4$ GeV/$c$. The elliptic flow measurement complements the previously reported ALICE results on the inclusive J/$\psi$ nuclear modification factor and favors the scenario of a significant fraction of J/$\psi$ production from charm quarks in a deconfined partonic phase.

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 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.

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$.

The study of the azimuthal anisotropy of inclusive muons produced in p-Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV, using the ALICE detector at the LHC is reported. The measurement of the second-order Fourier coefficient of the particle azimuthal distribution, $v_2$, is performed as a function of transverse momentum $p_{\rm T}$ in the 0-20% high-multiplicity interval at both forward ($2.03 < y_{\rm CMS} < 3.53$) and backward ($-4.46 < y_{\rm CMS} < -2.96$) rapidities over a wide $p_{\rm T}$ range, $0.5 < p_{\rm T} < 10$ GeV/$c$, in which a dominant contribution of muons from heavy-flavour hadron decays is expected at $p_{\rm T} > 2$ GeV/$c$. The $v_2$ coefficient of inclusive muons is extracted using two different techniques, namely two-particle cumulants, used for the first time for heavy-flavour measurements, and forward-central two-particle correlations. Both techniques give compatible results. A positive $v_2$ is measured at both forward and backward rapidities with a significance larger than $4.7\sigma$ and $7.6\sigma$, respectively, in the interval $2 < p_{\rm T} < 6$ GeV/$c$. Comparisons with previous measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, and with AMPT and CGC-based theoretical calculations are discussed. The findings impose new constraints on the theoretical interpretations of the origin of the collective behaviour in small collision systems.