skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Dust–gas dynamics driven by the streaming instability with various pressure gradients
ABSTRACT The streaming instability, a promising mechanism to drive planetesimal formation in dusty protoplanetary discs, relies on aerodynamic drag naturally induced by the background radial pressure gradient. This gradient should vary in discs, but its effect on the streaming instability has not been sufficiently explored. For this purpose, we use numerical simulations of an unstratified disc to study the non-linear saturation of the streaming instability with mono-disperse dust particles and survey a wide range of gradients for two distinct combinations of the particle stopping time and the dust-to-gas mass ratio. As the gradient increases, we find most kinematic and morphological properties increase but not always in linear proportion. The density distributions of tightly coupled particles are insensitive to the gradient whereas marginally coupled particles tend to concentrate by more than an order of magnitude as the gradient decreases. Moreover, dust–gas vortices for tightly coupled particles shrink as the gradient decreases, and we note higher resolutions are required to trigger the instability in this case. In addition, we find various properties at saturation that depend on the gradient may be observable and may help reconstruct models of observed discs dominated by streaming turbulence. In general, increased dust diffusion from stronger gradients can lower the concentration of dust filaments and can explain the higher solid abundances needed to trigger strong particle clumping and the reduced planetesimal formation efficiency previously found in vertically stratified simulations.  more » « less
Award ID(s):
2007422
PAR ID:
10621847
Author(s) / Creator(s):
; ;
Publisher / Repository:
Royal Astronomical Society
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
529
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
275 to 295
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The streaming instability is a fundamental process that can drive dust–gas dynamics and ultimately planetesimal formation in protoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classified into two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio ϵ. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into non-linear saturation. We find that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximum dimensionless stopping time τs,max = 0.1 and ϵ = 2 drives turbulent vertical dust–gas vortices, while the other with τs,max = 2 and ϵ = 0.2 leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former is flat in low densities, while the one for the latter has a low-end cut-off. By contrast, the one slow-growth case results in a virtually quiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particles moving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust species is appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived from multiwavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latter along with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs. 
    more » « less
  2. null (Ed.)
    ABSTRACT We develop simple, physically motivated models for drag-induced dust–gas streaming instabilities, which are thought to be crucial for clumping grains to form planetesimals in protoplanetary discs. The models explain, based on the physics of gaseous epicyclic motion and dust–gas drag forces, the most important features of the streaming instability and its simple generalization, the disc settling instability. Some of the key properties explained by our models include the sudden change in the growth rate of the streaming instability when the dust-to-gas mass ratio surpasses one, the slow growth rate of the streaming instability compared to the settling instability for smaller grains, and the main physical processes underlying the growth of the most unstable modes in different regimes. As well as providing helpful simplified pictures for understanding the operation of an interesting and fundamental astrophysical fluid instability, our models may prove useful for analysing simulations and developing non-linear theories of planetesimal growth in discs. 
    more » « less
  3. ABSTRACT We investigate the dynamics of dust concentration in actively accreting, substructured, non-ideal magnetohydrodynamic wind-launching discs using two-dimensional and three-dimensional (3D) simulations incorporating pressureless dust fluids of various grain sizes and their aerodynamic feedback on gas dynamics. Our results reveal that mm/cm-sized grains are preferentially concentrated within the inner 5–10 au of the disc, where the dust-to-gas surface density ratio (local metallicity Z) significantly exceeds the canonical 0.01, reaching values up to 0.25. This enhancement arises from the interplay of dust settling and complex gas flows in the meridional plane, including mid-plane accretion streams at early times, mid-plane expansion driven by magnetically braked surface accretion at later times, and vigorous meridional circulation in spontaneously formed gas rings. The resulting size-dependent dust distribution has a strong spatial variation, with large grains preferentially accumulating in dense rings, particularly in the inner disc, while being depleted in low-density gas gaps. In 3D, these rings and gaps are unstable to Rossby wave instability, generating arc-shaped vortices that stand out more prominently than their gas counterparts in the inner disc because of preferential dust concentration at small radii. The substantial local enhancement of the dust relative to the gas could promote planetesimal formation via streaming instability, potentially aided by the ‘azimuthal drift’ streaming instability that operates efficiently in accreting discs and a lower Toomre Q expected in younger discs. Our findings suggest that actively accreting young discs may provide favourable conditions for early planetesimal formation, which warrants further investigation. 
    more » « less
  4. null (Ed.)
    ABSTRACT A recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth time-scale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles to mutually validate and study the unstable modes of the instability in more detail. We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution [Ts, min, Ts, max] and the total solid-to-gas mass ratio ε and introduce a third dimension with the slope q of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the ε–Ts, max space, while this boundary is not appreciably sensitive to q or Ts, min. With a wide range of dust sizes present in the disc (e.g. Ts, min ≲ 10−3), the growth rate in the slow-growth regime decreases as more dust species are considered. With a narrow range of dust sizes (e.g. Ts, max/Ts, min = 5), on the other hand, the growth rate in most of the ε–Ts, max space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by ε ≳ 1 or Ts, max ≳ 1, which may represent the favourable conditions for planetesimal formation. 
    more » « less
  5. Abstract The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface density ratio (Z), above which SI-induced concentration triggers planetesimal formation. The thresholdZdepends on the dimensionless stopping time (τs, a proxy for dust size). However, these studies neglected both particle self-gravity and external turbulence. Here, we perform 3D stratified shearing box simulations with both particle self-gravity and turbulent forcing, which we characterize via a turbulent diffusion parameter,αD. We find that forced turbulence, at amplitudes plausibly present in some protoplanetary disks, can increase the thresholdZby up to an order of magnitude. For example, forτs= 0.01, planetesimal formation occurs whenZ≳ 0.06, ≳0.1, and ≳0.2 atαD= 10−4, 10−3.5, and 10−3, respectively. We provide a single fit to the criticalZrequired for the SI to work as a function ofαDandτs(although limited to the rangeτs= 0.01–0.1). Our simulations also show that planetesimal formation requires a mid-plane particle-to-gas density ratio that exceeds unity, with the critical value being largely insensitive toαD. Finally, we provide an estimation of particle scale height that accounts for both particle feedback and external turbulence. 
    more » « less