Heavy meson flow harmonics in event-by-event viscous relativistic hydrodynamics

Event-by-event viscous hydrodynamics is combined with heavy quark energy loss models to compute heavy flavor flow cumulants $v_2\{2\}$, $v_3\{2\}$, and $v_2\{4\}$ as well as the nuclear modification factors of $D^0$ and $B^0$ mesons in PbPb collisions at 2.76 TeV. Our results indicate that bottom quarks can flow as much as charm quarks in the $p_T$ range 8--30 GeV.