Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform simulations of small particles with significantly higher resolutions and longer simulation times than in previous investigations. We find that particles of dimensionless stopping time τ s = 10 −2 and 10 −3 -representing mm-and cm-sized particles interior of the water ice line -concentrate themselves via the streaming instability at a solid abundance of a few percent. We thus revise a previously published critical solid abundance curve for the regime of τ s 1. The solid density in the concentrated regions reaches values higher than the Roche density, indicating that direct collapse of particles down to mm sizes into planetesimals is possible. Our results hence bridge the gap in particle size between direct dust growth limited by bouncing and the streaming instability.

The streaming instability is a mechanism to concentrate solid particles into overdense filaments that undergo gravitational collapse and form planetesimals. However, it remains unclear how the initial mass function of these planetesimals depends on the box dimensions of numerical simulations. To resolve this, we perform simulations of planetesimal formation with the largest box dimensions to date, allowing planetesimals to form simultaneously in multiple filaments that can only emerge within such large simulation boxes. In our simulations, planetesimals with sizes between 80 km and several hundred kilometers form. We find that a power law with a rather shallow exponential cutoff at the high-mass end represents the cumulative birth mass function better than an integrated power law. The steepness of the exponential cutoff is largely independent of box dimensions and resolution, while the exponent of the power law is not constrained at the resolutions we employ. Moreover, we find that the characteristic mass scale of the exponential cutoff correlates with the mass budget in each filament. Together with previous studies of high-resolution simulations with small box domains, our results therefore imply that the cumulative birth mass function of planetesimals is consistent with an exponentially tapered power law with a power-law exponent of approximately −1.6 and a steepness of the exponential cutoff in the range of 0.3-0.4.

The streaming instability is a promising mechanism to overcome the barriers in direct dust growth and lead to the formation of planetesimals. Most previous studies of the streaming instability, however, were focused on a local region of a protoplanetary disk with a limited simulation domain such that only one filamentary concentration of solids has been observed. The characteristic separation between filaments is therefore not known. To address this, we conduct the largest-scale simulations of the streaming instability to date, with computational domains up to 1.6 gas scale heights both horizontally and vertically. The large dynamical range allows the effect of vertical gas stratification to become prominent. We observe more frequent merging and splitting of filaments in simulation boxes of high vertical extent. We find multiple filamentary concentrations of solids with an average separation of about 0.2 local gas scale heights, much higher than the most unstable wavelength from linear stability analysis. This measures the characteristic separation of planetesimal forming events driven by the streaming instability and thus the initial feeding zone of planetesimals.

The streaming instability is a promising mechanism to drive the formation of planetesimals in protoplanetary disks. To trigger this process, it has been argued that sedimentation of solids onto the mid-plane needs to be efficient and therefore that a quiescent gaseous environment is required. It is often suggested that dead-zone or disk-wind structure created by non-ideal magnetohydrodynamical (MHD) effects meets this requirement. However, simulations have shown that the midplane of a dead zone is not completely quiescent. In order to examine the concentration of solids in such an environment, we use the local-shearing-box approximation to simulate a particle-gas system with an Ohmic dead zone including mutual drag force between the gas and the solids. We systematically compare the evolution of the system with ideal or non-ideal MHD, with or without back-reaction drag force from particles on gas, and with varying solid abundances. Similar to previous investigations of dead zone dynamics, we find that particles of dimensionless stopping time τ s = 0.1 do not sediment appreciably more than those in ideal magneto-rotational turbulence, resulting in a vertical scale height an order of magnitude larger than in a laminar disk. Contrary to the expectation that this should curb the formation of planetesimals, we nevertheless find that strong clumping of solids still occurs in the dead zone when solid abundances are similar to the critical value for a laminar environment. This can be explained by the weak radial diffusion of particles near the mid-plane. The results imply that the sedimentation of particles to the mid-plane is not a necessary criterion for the formation of planetesimals by the streaming instability.

Observations show that radial metallicity gradients in disk galaxies are relatively shallow, if not flat, especially at large galactocentric distances and for galaxies in the high-redshift universe. Given that star formation and metal production are centrally concentrated, this requires a mechanism to redistribute metals. However, the nature of this mechanism is poorly understood, let alone quantified. To address this problem, we conduct magnetohydrodynamical simulations of a local shearing sheet of a thin, thermally unstable, gaseous disk driven by a background stellar spiral potential, including metals modeled as passive scalar fields. Contrary to what a simple α prescription for the gas disk would suggest, we find that turbulence driven by thermal instability is very efficient at mixing metals, regardless of the presence or absence of stellar spiral potentials or magnetic fields. The timescale for homogenizing randomly distributed metals is comparable to or less than the local orbital time in the disk. This implies that turbulent mixing of metals is a significant process in the history of chemical evolution of disk galaxies.