Sensitivity study with a <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>D</mml:mi></mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>B</mml:mi></mml:math> mesons modular simulation code of heavy flavor <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub></mml:math> and azimuthal anisotropies based on beam energy, initial conditions, hadronization, and suppression mechanisms

Heavy flavor probes provide important information about the in-medium properties of the quark-gluon plasma produced in heavy-ion collisions. In this work, we investigate the effects of $2\text{D}+1$ event-by-event fluctuating hydrodynamic backgrounds on the nuclear suppression factor and momentum anisotropies of heavy flavor mesons and nonphotonic electrons. Using the state-of-the-art $D$ and $B$ mesons modular simulation code (called ``DAB-MOD''), we perform a systematic comparison of different transport equations in the same background, including a few energy-loss models---with and without energy-loss fluctuations---and a relativistic Langevin model with different drag parametrizations. We present the resulting $D$ and $B$ mesons ${R}_{AA}$, ${v}_{2}$, ${v}_{3}$, and ${v}_{4}$ as well as multiparticle cumulants, in $\text{AuAu}$ collisions at $\sqrt{{s}_{NN}}=200\phantom{\rule{0.16em}{0ex}}\mathrm{GeV}$ and $\text{PbPb}$ collisions at $\sqrt{{s}_{NN}}=2.76\phantom{\rule{0.16em}{0ex}}\mathrm{TeV}$ and $\sqrt{{s}_{NN}}=5.02\phantom{\rule{0.16em}{0ex}}\mathrm{TeV}$, and compare them to the available experimental data. The ${v}_{2}{4}/{v}_{2}{2}$ ratio, which is known to be a powerful probe of the initial conditions and flow fluctuations in the soft sector, is also studied in the context of heavy flavor. We also investigate the correlations between the transverse anisotropies of heavy mesons and all charged particles to better understand how heavy quarks couple to the hydrodynamically expanding quark-gluon plasma. We study the influence that different initial conditions and the implementation of heavy-light quark coalescence has on our results.