Heavy quark diffusion coefficient in heavy-ion collisions via kinetic theory

We compute the heavy quark momentum diffusion coefficient κ using QCD kinetic theory for a system going through bottom-up isotropization in the initial stages of a heavy ion collision. We find that the values of κ are within 30% from a thermal system at the same energy density. When matching for other quantities we observe considerably larger deviations. We also observe that the diffusion coefficient in the transverse direction is larger at high occupation numbers, whereas for an underoccupied system the longitudinal diffusion coefficient dominates. The behavior of the diffusion coefficient can be understood on a qualitative level based on the Debye mass mD and the effective temperature of soft modes T∗. Our findings for the kinetic evolution of κ in different directions can be used in phenomenological descriptions of heavy quark diffusion and quarkonium dynamics to include the impact of preequilibrium stages.