Time evolution of the mSI (run: mSI₁₀₂₄) and PSI (run: PSI10₁₀₂₄). The dot-dashed line shows the mSI run, the solid lines show the ten individual dust species of the PSI run, and the dotted line is the sum over all the dust species. Plot A shows the amplitude of the largest mode in the shearing box (A sinusoid with wavenumber K = (30, 0, 30)^T), plot B shows the maximum density in the shearing box, plot C shows the mean density at every snapshot's 99^th percentile and plot D shows the average probability density function of the normalised density distribution in the non-linear regime between (140 ≤ Ωt ≤ 160).