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
Streaming instability with multiple dust species – I. Favourable conditions for the linear growth
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
- Award ID(s):
- 1753168
- PAR ID:
- 10273107
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 501
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 467 to 482
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
ABSTRACT We address the formation of giant clumps in violently unstable gas-rich disc galaxies at cosmic noon. While these are commonly thought to originate from gravitational Toomre instability, some cosmological simulations have indicated that clumps can form in Lagrangian proto-clump regions where the Toomre Q parameter is well above unity, which are linearly stable. Examining one of these cosmological simulations, we find that it exhibits an excess in compressive modes of turbulence with converging motions. The energy in converging motions within proto-clumps is $${\sim} 70~{{\ \rm per\ cent}}$$ of the total turbulent energy, compared to $${\sim} 17~{{\ \rm per\ cent}}$$ expected in equipartition. When averaged over the whole disc, $${\sim} 40~{{\ \rm per\ cent}}$$ of the turbulent energy is in compressive modes, mostly in converging motions, with the rest in solenoidal modes, compared to the $(1/3)-(2/3)$ division expected in equipartition. By contrast, we find that in an isolated-disc simulation with similar properties, resembling high-z star-forming galaxies, the different turbulence modes are in equipartition, both in proto-clumps and over the whole disc. We conclude that the origin of excessive converging motions in proto-clumps is external to the disc, and propose several mechanisms that can induce them. This is an additional mechanism for clump formation, complementary to and possibly preceding gravitational instability.more » « less
-
Abstract Millimeter emitting dust grains have sizes that make them susceptible to drift in protoplanetary disks due to the difference between their orbital speed and that of the gas. The characteristic drift timescale depends on the surface density of the gas. By comparing disk radius measurements from Atacama Large Millimeter/submillimeter Array CO and continuum observations at millimeter wavelengths, the gas surface density profile and dust drift time can be self-consistently determined. We find that profiles which match the measured dust mass have very short drift timescales, an order of magnitude or more shorter than the stellar age, whereas profiles for disks that are on the cusp of gravitational instability, defined via the minimum value of the Toomre parameter, , have drift timescales comparable to the stellar lifetime. This holds for disks with masses of dust ≳5M⊕across a range of absolute ages from less than 1 Myr to over 10 Myr. The inferred disk masses scale with stellar mass as . This interpretation of the gas and dust disk sizes simultaneously solves two long standing issues regarding the dust lifetime and exoplanet mass budget, and suggests that we consider millimeter wavelength observations as a window into an underlying population of particles with a wide size distribution in secular evolution with a massive planetesimal disk.more » « less
An official website of the United States government

