No Arabic abstract
The streaming instability is often invoked as solution to the fragmentation and drift barriers in planetesimal formation, catalyzing the aggregation of dust on kyr timescales to grow km-sized cores. However there remains a lack of consensus on the physical mechanism(s) responsible for initiating it. One potential avenue is disc photoevaporation, wherein the preferential removal of relatively dust-free gas increases the disc metallicity. Late in the disc lifetime, photoevaporation dominates viscous accretion, creating a gradient in the depleted gas surface density near the location of the gap. This induces a local pressure maximum that collects drifting dust particles, which may then become susceptible to the streaming instability. Using a one-dimensional viscous evolution model of a disc subject to internal X-ray photoevaporation, we explore the efficacy of this process to build planetestimals. Over a range of parameters we find that the amount of dust mass converted into planetesimals is often < 1 M_Earth and at most a few M_Earth spread across tens of AU. We conclude that photoevaporation may at best be relevant for the formation of debris discs, rather than a common mechanism for the formation of planetary cores. Our results are in contrast to a recent, similar investigation that considered an FUV-driven photoevaporation model and reported the formation of tens of M_Earth at large (> 100 AU) disc radii. The discrepancies are primarily a consequence of the different photoevaporation profiles assumed. Until observations more tightly constrain photoevaporation models, the relevance of this process to the formation of planets remains uncertain.
Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_oplus$) planetesimal belts beyond 100 au and up to $sim 20 M_oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.
A critical step toward the emergence of planets in a protoplanetary disk consists in accretion of planetesimals, bodies 1-1000 km in size, from smaller disk constituents. This process is poorly understood partly because we lack good observational constraints on the complex physical processes that contribute to planetesimal formation. In the outer solar system, the best place to look for clues is the Kuiper belt, where icy planetesimals survived to this day. Here we report evidence that Kuiper belt planetesimals formed by the streaming instability, a process in which aerodynamically concentrated clumps of pebbles gravitationally collapse into 100-km-class bodies. Gravitational collapse was previously suggested to explain the ubiquity of equal-size binaries in the Kuiper belt. We analyze new hydrodynamical simulations of the streaming instability to determine the model expectations for the spatial orientation of binary orbits. The predicted broad inclination distribution with 80% of prograde binary orbits matches the observations of trans-Neptunian binaries. The formation models which imply predominantly retrograde binary orbits can be ruled out. Given its applicability over a broad range of protoplanetary disk conditions, it is expected that the streaming instability seeded planetesimal formation also elsewhere in the solar system, and beyond.
Comets and small planetesimals are believed to contain primordial building blocks in the form of millimeter to centimeter sized pebbles. One of the viable growing mechanisms to form these small bodies is through the streaming instability (SI) in which pebbles cluster and gravitationally collapse towards a planetesimal or comet in the presence of gas drag. However, most SI simulations are global and lack the resolution to follow the final collapse stage of a pebble cloud within its Hill radius. We aim to track the collapse of a gravitationally bound pebble cloud subject to mutual collisions and gas drag with the representative particle approach. We determine the radial pebble size distribution of the collapsed core and the impact of mutual pebble collisions on the pebble size distribution. We find that virial equilibrium is never reached during the cloud evolution and that, in general, pebbles with given Stokes number (St) collapse towards an optically thick core in a sequence from aerodynamically largest to aerodynamically smallest. We show that at the location for which the core becomes optically thick, the terminal velocity is well below the fragmentation threshold velocity. While collisional processing is negligible during cloud evolution, the collisions that do occur are sticking. These results support the observations that comets and small planetary bodies are composed of primordial pebbles in the milimeter to centimeter size range
The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.
Damping of the previously discovered resonant drag instability (RDI) of dust streaming in protoplanetary disc is studied using the local approach to dynamics of gas-dust perturbations in the limit of the small dust fraction. Turbulence in a disc is represented by the effective viscosity and diffusivity in equations of motion for gas and dust, respectively. In the standard case of the Schmidt number (ratio of the effective viscosity to diffusivity) Sc = 1, the reduced description of RDI in terms of the inertial wave (IW) and the streaming dust wave (SDW) falling in resonance with each other reveals that damping solution differs from the inviscid solution simply by adding the characteristic damping frequency to its growth rate. RDI is fully suppressed at the threshold viscosity, which is estimated analytically, first, for radial drift, next, for vertical settling of dust, and at last, in the case of settling combined with radial drift of the dust. In the last case, RDI survives up to the highest threshold viscosity, with a greater excess for smaller solids. Once Sc eq 1, a new instability specific for dissipative perturbations on the dust settling background emerges. This instability of the quasi-resonant nature is referred to as settling viscous instability (SVI). The mode akin to SDW (IW) becomes growing in a region of long waves provided that Sc > 1 (Sc < 1). SVI leads to an additional increase of the threshold viscosity.