Transverse-momentum and event-shape dependence of D-meson flow harmonics in Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

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 found to follow a mass ordering for $p_{\rm T}<3$ GeV/$c$, and to be 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.