Thousands of planetesimals: Simulating the streaming Instability in very large computational domains

TL;DR

This study uses large-scale simulations to analyze the streaming instability in protoplanetary disks, revealing filament structures limited to about one gas scale height and forming thousands of planetesimals, thereby advancing understanding of planetesimal formation.

Contribution

The paper presents the first large-scale simulations of the streaming instability in domains 32 times larger than previous studies, showing filament structures and forming thousands of planetesimals.

Findings

Filaments are limited to approximately one gas scale height.

The simulations form up to 4,000 planetesimals, enabling better statistical analysis.

The mass distribution follows a steep exponential taper at high masses.

Abstract

The streaming instability is a mechanism whereby pebble-sized particles in protoplanetary discs spontaneously come together in dense filaments, which collapse gravitationally to form planetesimals upon reaching the Roche density. The extent of the filaments along the orbital direction is nevertheless poorly characterised, due to a focus in the literature on small simulation domains where the behaviour of the streaming instability on large scales cannot be determined. We present here computer simulations of the streaming instability in boxes with side lengths up to 6.4 scale heights in the plane. This is 32 times larger than typically considered simulation domains and nearly a factor 1,000 times the volume. We show that the azimuthal extent of filaments in the non-linear state of the streaming instability is limited to approximately one gas scale height. The streaming instability will therefore not transform the pebble density field into axisymmetric rings; rather the non-linear state of the streaming instability appears as a complex structure of loosely connected filaments. Including the self-gravity of the pebbles, our simulations form up to 4,000 planetesimals. This allows us to probe the high-mass end of the initial mass function of planetesimals with much higher statistical confidence than previously. We find that this end is well-described by a steep exponential tapering. Since the resolution of our simulations is moderate -- a necessary trade-off given the large domains -- the mass distribution is incomplete at the low-mass end. When putting comparatively less weight on the numbers at low masses, at intermediate masses we nevertheless reproduce the power-law shape of the distribution established in previous studies.