Dust dynamics in RAMSES -- II. Equilibrium drift velocity distributions of charged dust grains

Eric R. Moseley,Romain Teyssier

2024-05-30

Abstract:We investigate the gas-grain relative drift velocity distributions of charged astrophysical dust grains in MHD turbulence. We do this using a range of MHD-PIC simulations spanning different plasma-$\beta$, sonic/Alfvén Mach number, and with grains of varying size and charge-to-mass ratio. We find that the root-mean-square drift velocity is a strong function of the grain size, following a power law with a 1/2 slope. The r.m.s. value has only a very weak dependence on the charge-to-mass ratio. On the other hand, the shape of the distribution is a strong function of the grain charge-to-mass ratio, and in compressible turbulence, also the grain size. We then compare these results to simple analytic models based upon time-domain quasi-linear theory and solutions to the Fokker-Planck equation. These models explain qualitatively the r.m.s. drift velocity's lack of charge-to-mass ratio dependence, as well as why the shape of the distribution changes as the charge-to-mass ratio increases. Finally we scale our results to astrophysical conditions. As an example, at a length scale of one parsec in the cold neutral medium, 0.1 $\mu$m grains should be drifting at roughly 40% of the turbulent velocity dispersion. These findings may serve as a basis for a model for grain velocities in the context of grain-grain collisions, non-thermal sputtering, and accretion of metals. These findings also have implications for the transport of grains through the galaxy, suggesting that grains may have non-negligible random motions at length-scales that many modern galaxy simulations approach.

Astrophysics of Galaxies,Cosmology and Nongalactic Astrophysics,Plasma Physics

What problem does this paper attempt to address?

This paper discusses the distribution of relative drift velocities between charged dust particles and gas in magnetohydrodynamic (MHD) turbulence. The study employs a series of MHD particle-in-cell (PIC) simulations to analyze the effects of different parameters such as plasma beta, sound speed/Alfvén Mach number, and dust particle size and charge-to-mass ratio on the dust drift velocity. The results show that the root mean square drift velocity of dust is proportional to its size, following a power-law with an exponent of 1/2, and it has a weak dependence on the charge-to-mass ratio. The shape of the drift velocity distribution strongly depends on the charge-to-mass ratio and is also influenced by the particle size in compressible turbulence.

The paper also compares simple analytical models based on time-domain quasi-linear theory and Fokker-Planck equation, and finds that these models qualitatively explain the insensitivity of the root mean square drift velocity to the charge-to-mass ratio and the changes in the distribution shape as the charge-to-mass ratio increases. Applying these results to astrophysical conditions, for example, in a parsec-scale cold neutral medium, dust particles of 0.1 micrometer should disperse at approximately 40% of the turbulent velocity. This may affect processes such as dust particle collisions, non-thermal sputtering, and metal absorption, and has significant implications for the transport of dust in galaxies, indicating the presence of significant random motion of dust particles at multiple length scales.

The study highlights that previous work has mainly focused on neutral dust species and neglected the behavior of charged dust in MHD turbulence. Furthermore, by systematically studying the dust velocity distribution using a three-dimensional MHD-PIC simulation, this is the first attempt of its kind that covers scenarios with supersonic and super-Alfvén Mach numbers. The paper concludes by discussing how the results can be applied to real interstellar medium conditions and proposing future research directions.