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.

We report a precise measurement of the J/$\psi$ elliptic flow in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector at the LHC. The J/$\psi$ mesons are reconstructed at mid-rapidity ($|y| < 0.9$) in the dielectron decay channel and at forward rapidity ($2.5<y<4.0$) in the dimuon channel, both down to zero transverse momentum. At forward rapidity, the elliptic flow $v_2$ of the J/$\psi$ is studied as a function of transverse momentum and centrality. A positive $v_2$ is observed in the transverse momentum range $2 < p_{\rm T} < 8$ GeV/$c$ in the three centrality classes studied and confirms with higher statistics our earlier results at $\sqrt{s_{\rm NN}} = 2.76$ TeV in semi-central collisions. At mid-rapidity, the J/$\psi$ $v_2$ is investigated as a function of transverse momentum in semi-central collisions and found to be in agreement with the measurements at forward rapidity. These results are compared to transport model calculations. The comparison supports the idea that at low $p_{\rm T}$ the elliptic flow of the J/$\psi$ originates from the thermalization of charm quarks in the deconfined medium, but suggests that additional mechanisms might be missing in the models.

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

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 present the first measurement of the two-particle transverse momentum differential correlation function, $P_2\equiv\langle \Delta p_{\rm T} \Delta p_{\rm T} \rangle /\langle p_{\rm T} \rangle^2$, in Pb-Pb collisions at $\sqrt{s_{_{\rm NN}}} =$ 2.76 TeV. Results for $P_2$ are reported as a function of relative pseudorapidity ($\Delta \eta$) and azimuthal angle ($\Delta \varphi$) between two particles for different collision centralities. The $\Delta \phi$ dependence is found to be largely independent of $\Delta \eta$ for $|\Delta \eta| \geq$ 0.9. In 5% most central Pb-Pb collisions, the two-particle transverse momentum correlation function exhibits a clear double-hump structure around $\Delta \varphi = \pi$ (i.e., on the away side), which is not observed in number correlations in the same centrality range, and thus provides an indication of the dominance of triangular flow in this collision centrality. Fourier decompositions of $P_2$, studied as a function of collision centrality, show that correlations at $|\Delta \eta| \geq$ 0.9 can be well reproduced by a flow ansatz based on the notion that measured momentum correlations are strictly determined by the collective motion of the system.

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