Summary
By utilizing a soft-hard factorized model, which combines a thermal perturbative description of soft scatterings and a perturbative QCD-based calculation for hard collisions, we study the energy and temperature dependence of the heavy quark diffusion coefficients in Langevin dynamics. The adjustable parameters are fixed from the comprehensive model-data comparison. We find that a small value of the spatial diffusion coefficient at transition temperature is preferred by data 2 pi TDs(Tc)similar or equal to 6. With the parameter-optimized model, we are able to describe simultaneously the prompt D0RAA and v2 data at pT <= 8 GeV in Pb-Pb collisions at <mml:msqrt>sNN</mml:msqrt>=2.76 and <mml:msqrt>sNN</mml:msqrt>=5.02 TeV. We also make predictions for non-prompt D0 meson for future experimental tests down to the low momentum region.