Issue |
A&A
Volume 695, March 2025
|
|
---|---|---|
Article Number | A158 | |
Number of page(s) | 16 | |
Section | Planets, planetary systems, and small bodies | |
DOI | https://doi.org/10.1051/0004-6361/202453072 | |
Published online | 18 March 2025 |
Polydisperse formation of planetesimals
The dust size distribution in clumps
Planetary Exploration, Technical University Delft,
Kluyverweg 1 2629 HS Delft,
The Netherlands
★ Corresponding author; j.p.matthijsse@tudelft.nl
Received:
19
November
2024
Accepted:
31
January
2025
Aims. To form kilometre-sized planetesimals, the streaming instability is an efficient method for overcoming the barriers to planet formation in protoplanetary discs. The streaming instability has been extensively modelled by hydrodynamic simulations of gas and a single dust size. However, recent studies considering a more realistic case of a particle size distribution have shown that this will significantly decrease the growth rate of the instability. We follow up on these studies by evaluating the polydisperse streaming instability in the non-linear regime to see if clumping can occur in the same manner as the monodisperse streaming instability and determine the size distribution in the densest dust structures.
Methods. We employ 2D hydrodynamic simulations in an unstratified shearing box with multiple dust species representing an underlying continuous dust size spectrum using FARGO3D. We use the Gauss-Legendre quadrature in dust size space to calculate the drag force on the gas due to a continuous dust size distribution. These simulations are compared to previous analytical results of the polydisperse streaming instability in the linear phase. We then look at the saturated non-linear phase of the instability at the highest density regions and investigate the dust size distribution in the densest dust structures.
Results. When sampling the size distribution, the error in the growth rate converges significantly faster with the number of dust sizes using the Gauss-Legendre quadrature method than the usual uniform sampling method. In the non-linear regime, the maximum dust density reached in the polydisperse case is reduced compared to the monodisperse case. Larger dust particles are most abundant in the densest dust structure because they are less coupled to the gas and can therefore clump together more than the smaller dust grains. Contrary to expectations based solely on dust-gas coupling, our results reveal a distinct peak in the size distribution that arises from the size-dependent spatial segregation of the highest-density regions, where particles with the largest Stokes numbers are located just outside the densest areas of the combined dust species.
Conclusions. The 2D unstratified polydisperse streaming instability is less efficient than its monodisperse counterpart at generating dense clumps that may collapse into planetesimals, and in the densest regions, the distinct dust size distribution could be related to the size distribution that ends up in the planetesimal and can mimic the size distribution of dust growth.
Key words: hydrodynamics / instabilities / methods: numerical / planets and satellites: formation / protoplanetary disks
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.