Stochastic Dynamics of Heavy Quarks in Strongly Coupled Plasma

We study the stochastic dynamics of heavy quarks propagating through the strongly coupled plasma of $\mathcal{N}=4$ supersymmetric Yang-Mills (SYM) theory at nonzero temperature in terms of the corresponding Kolmogorov equation, which correctly describes their kinetic equilibration and the non-Gaussian fluctuations in their momenta without having to restrict their velocity to the non-relativistic regime. Leveraging the heavy quark limit, we show that the evolution of the momentum space distribution function can be reformulated as a Hamilton-Jacobi problem, and therefore can be solved in terms of first-order ordinary differential equations. We solve these evolution equations in an infinite thermal plasma with a constant temperature for initial conditions specified by spherically symmetric heavy quark momentum distributions with a phenomenologically motivated shape that is steeply falling at large momentum. To highlight the distinctive features of the kinetic equilibration process, we compare their solutions with Fokker-Planck dynamics with the same drag coefficient and fluctuations chosen by hand to guarantee equilibration. We find qualitatively similar dynamics at small momentum, and very different dynamics at large momentum, where, much like in jet quenching phenomena, the steepness of the momentum distribution gives a larger relevance to unlikely events in which a heavy quark loses little momentum in a given time step. Such events are much less unlikely in the Kolmogorov evolution than in Fokker-Planck evolution with the same mean energy loss, meaning that equilibration at large momentum is significantly delayed. Our results provide a systematic description of heavy quarks propagating through strongly coupled plasma from the ultra-relativistic to the non-relativistic regime and point the way towards implementation in phenomenological studies.