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-squared 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 RAA. Notably, our findings indicate that superdiffusion strongly suppresses the RAA compared to normal diffusion in the hot QCD medium. The possible impacts on other parameters are also discussed. Published by the American Physical Society 2024