  The dynamics of heavy quarks within the hot QCD medium have been revisited,
considering the influence of anomalous diffusion. This study has been conducted
using the framework of the fractional Langevin equation involving the Caputo
fractional derivative. We introduce a numerical scheme for the fractional
Langevin equation and demonstrate that the mean square displacement of the
particle exhibits anomalous diffusion, deviating from a linear relationship
with time. Our analysis calculates various entities, such as mean squared
momentum, momentum spread, and the nuclear suppression factor, $R_{AA}$.
Notably, our findings indicate that superdiffusion strongly suppresses the
$R_{AA}$ compared to normal diffusion in the hot QCD medium. The possible
impacts on other parameters are also discussed.
