Dust Diffusion in Protoplanetary Disks by Magnetorotational Turbulence

We measure the turbulent diffusion coefficient of dust grains embedded in magnetorotational turbulence in a protoplanetary disc directly from numerical simulations and compare it to the turbulent viscosity of the flow. The simulations are done in a local coordinate frame comoving with the gas in Keplerian rotation. Periodic boundary conditions are used in all directions, and vertical gravity is not applied to the gas. Using a two-fluid approach, small dust grains of various sizes (with friction times up to ₀ ₅=0. 02) are allowed to move under the influence of friction with the turbulent gas. We measure the turbulent diffusion coefficient of the dust grains by applying an external sinusoidal force field acting in the vertical direction on the dust component only. This concentrates the dust around the mid-plane of the disc, and an equilibrium distribution of the dust density is achieved when the vertical settling is counteracted by the turbulent diffusion away from the mid-plane. Comparing with analytical expressions for the equilibrium concentration we deduce the vertical turbulent diffusion coefficient. The vertical diffusion coefficient is found to be lower than the turbulent viscosity and to have an associated vertical diffusion Prandtl number of about 1. 5. A similar radial force field also allows us to measure the radial turbulent diffusion coefficient. We find a radial diffusion Prandtl number of about 0. 85 and also find that the radial turbulent diffusion coefficient is around 70% higher than the vertical. We also find evidence for trapping of dust grains of intermediate friction time in turbulent eddies.