Polydisperse Formation of Planetesimals: The dust size distribution in clumps

Abstract

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, more recent studies considering a more realistic case of a particle size distribution show that this will significantly decrease the growth rate of the instability. We follow up on these studies by evaluating the polydisperse streaming instability, looking at the non-linear phase of the instability at the highest density regions, and investigating the dust size distribution in the densest dust structures. We employ 2D hydrodynamic simulations in an unstratified shearing box with multiple dust species representing an underlying continuous dust size spectrum using FARGO3D. To calculate the drag force on the gas due to a continuous dust size distribution, we apply the Gauss-Legendre quadrature method in dust size space. This method converges faster with the number of dust species than the usual uniform sampling method. The polydisperse streaming instability is less efficient than its monodisperse counterpart in generating dense clumps that could collapse into planetesimals. In the densest dust structure, the larger dust sizes are more abundant because they are less coupled to the gas and, therefore, can clump together more than the smaller dust grains. This trend is broken at the largest dust size due to 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. This is observed as a peak in the size distribution at the densest regions, which could relate to the size distribution that ends up in the planetesimal after collapse and can mimic the size distribution of dust growth.

16 pages, 19 figures, 2 tables. Accepted in A&A

References in corpus (16)

• Array Programming with NumPy
• Dust size distributions in coagulation/fragmentation equilibrium: Numerical solutions and analytical fits
• Evidence for the formation of comet 67P/Churyumov-Gerasimenko through gravitational collapse of a bound clump of pebbles
• Thresholds for Particle Clumping by the Streaming Instability
• Streaming Instability for Particle-Size Distributions
• Predicting the observational signature of migrating Neptune-sized planets in low-viscosity disks
• Streaming Instability with Multiple Dust Species: II. Turbulence and Dust-Gas Dynamics at Nonlinear Saturation
• Polydisperse Streaming Instability III. Dust evolution encourages fast instability
• Stratified and vertically-shearing streaming instabilities in protoplanetary disks
• Streaming instability of multiple particle species II -- Numerical convergence with increasing particle number
• Polydisperse Streaming Instability II. Methods for solving the linear stability problem
• Radially resolved simulations of collapsing pebble clouds in protoplanetary discs
• Constraints on planetesimal accretion inferred from particle-size distribution in CO chondrites
• The radial structure of planetary bodies formed by the streaming instability
• Instabilities in dusty non-isothermal proto-planetary discs
• Planetesimal formation via the streaming instability with multiple grain sizes

Cited by in corpus (2)

• Nonlinear Evolution of the Unstratified Polydisperse Dust Settling Instability
• A multi-fluid approach for polydisperse pebble accretion: From particles to fluids, establishing the multifluid framework