Quasi-Particle Model (QPM) has provided a good description of the main features of lattice QCD (lQCD) and it play a key for studying heavy quark (HQ) dynamics in ultra-relativistic heavy-ion collisions. It can provide a phenomenological estimation of the HQ diffusion coefficient, $D_s(T)$. In this paper, we investigate the bottom quark dynamics within an event-by-event full Boltzmann transport approach coupled with QPM and followed by a hybrid hadronization via coalescence plus fragmentation. Within this approach we provide predictions for $R_{AA}$ and $v_{2}$ of B mesons.
Finally, we introduce an extension of QPM, where quark and gluon masses depend on their momentum inspired by Dyson-Schwinger studies and converging toward the current quark mass at high momenta. This extended model denoted as QPMp, successfully provides a simultaneous quantitative description of both the Equation of State (EoS) and quark susceptibilities.

