Initially Maxwellian distribution of test particles with thermal velocity v _th = c _s is simulated for t = 20 L/ c _s (many eddy turnover times) in our fiducial simulation (Table 1), with two different values of the imposed scattering rate ν. The final distribution functions are slightly non-Maxwellian but do not develop a significantly non-thermal tail. f( v) shown here is the one-dimensional distribution function, normalized such that the density is n ₀ = ∫ dvf( v).