Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStreaming Instability with Multiple Dust Species: I. Favourable Conditions for the Linear Growth
80
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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 timescale 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 $[T_{s,min}, T_{s,max}]$ and the total solid-to-gas mass ratio $varepsilon$ 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 $varepsilon$-$T_{s,max}$ space, while this boundary is not appreciably sensitive to $q$ or $T_{s,min}$. With a wide range of dust sizes present in the disc (e.g. $T_{s,min}lesssim10^{-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. $T_{s,max}/T_{s,min}=5$), on the other hand, the growth rate in most of the $varepsilon$-$T_{s,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 $varepsilongtrsim 1$ or $T_{s,max}gtrsim 1$, which may represent the favourable conditions for planetesimal formation.
rate research
Read More
Planet formation via core accretion requires the production of km-sized planetesimals from cosmic dust. This process must overcome barriers to simple collisional growth, for which the Streaming Instability (SI) is often invoked. Dust evolution is still required to create particles large enough to undergo vigorous instability. The SI has been studied primarily with single size dust, and the role of the full evolved dust distribution is largely unexplored. We survey the Polydispserse Streaming Instability (PSI) with physical parameters corresponding to plausible conditions in protoplanetary discs. We consider a full range of particle stopping times, generalized dust size distributions, and the effect of turbulence. We find that, while the PSI grows in many cases more slowly with a interstellar power-law dust distribution than with a single size, reasonable collisional dust evolution, producing an enhancement of the largest dust sizes, produces instability behaviour similar to the monodisperse case. Considering turbulent diffusion the trend is similar. We conclude that if fast linear growth of PSI is required for planet formation, then dust evolution producing a distribution with peak stopping times on the order of 0.1 orbits and an enhancement of the largest dust significantly above the single power-law distribution produced by a fragmentation cascade is sufficient, along with local enhancement of the dust to gas volume mass density ratio to order unity.
Occurring in protoplanetary discs composed of dust and gas, streaming instabilities are a favoured mechanism to drive the formation of planetesimals. The Polydispserse Streaming Instability is a generalisation of the Streaming Instability to a continuum of dust sizes. This second paper in the series provides a more in-depth derivation of the governing equations and presents novel numerical methods for solving the associated linear stability problem. In addition to the direct discretisation of the eigenproblem at second order introduced in the previous paper, a new technique based on numerically reducing the system of integral equations to a complex polynomial combined with root finding is found to yield accurate results at much lower computational cost. A related method for counting roots of the dispersion relation inside a contour without locating those roots is also demonstrated. Applications of these methods show they can reproduce and exceed the accuracy of previous results in the literature, and new benchmark results are provided. Implementations of the methods described are made available in an accompanying Python package psitools.
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.
The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfven waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($lambda ll rm pc$), low cosmic ray number density ($n_{rm CR}/n_{rm ISM} sim 10^{-9}$), and weak CR anisotropy ($sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $delta f$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfven modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics.
The recently discovered resonant drag instability (RDI) of dust streaming in protoplanetary disc is considered as the mode coupling of subsonic gas-dust mixture perturbations. This mode coupling is coalescence of two modes with nearly equal phase velocities: inertial wave (IW) having positive energy and a streaming dust wave (SDW) having negative energy as measured in the frame of gas environment being at rest in vertical hydrostatic equilibrium. SDW is a trivial mode produced by the bulk streaming of dust, which transports perturbations of dust density. In this way, settling combined with radial drift of the dust makes possible coupling of SDW with IW and the onset of the instability. In accordance with the concept of the mode coupling, RDI growth rate is proportional to the square root of the coupling term of the dispersion equation, which itself is proportional to mass fraction of dust, $fll 1$. This clarifies why RDI growth rate $propto f^{1/2}$. When SDW has positive energy, its resonance with IW provides an avoided crossing instead of the mode coupling. In the high wavenumber limit RDI with unbounded growth rate $propto f^{1/3}$ is explained by the triple mode coupling, which is coupling of SDW with two IW. It coexists with a new quasi-resonant instability accompanied by bonding of two oppositely propagating low-frequency IW. The mode coupling does not exist for dust streaming only radially in a disc. In this case RDI is provided by the obscured mechanism associated with the inertia of solids.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
