The size of dust particles used in our model is determined by a two-populations algorithm derived from [1]}. This algorithm computes the representative size of particles through the estimate of the limiting Stokes number in various dynamical regimes. In our model, dust is initially present in the form of particles of sizes \(a_0 = 10^{-7}\)  m, and grow through mutual collisions. This growth is limited by the maximum sizes imposed by fragmentation or by the drift velocity of the grains (see [2]} for details). The dust surface density is the sum over all surface densities of available solids at given time and location, assuming a protosolar ice-to-rock ratio (\(\sim \) 1.58) [3]}, and a bulk density of 1000 kg m\(^{-3}\) .
