Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGlobal drag-induced instabilities in protoplanetary disks
93
0
0.0
(
0
)
Authors
Mir Abbas Jalali
Ask ChatGPT about the research
No Arabic abstract
We use the Fokker-Planck equation and model the dispersive dynamics of solid particles in annular protoplanetary disks whose gas component is more massive than the particle phase. We model particle--gas interactions as hard sphere collisions, determine the functional form of diffusion coefficients, and show the existence of two global unstable modes in the particle phase. These modes have spiral patterns with the azimuthal wavenumber $m=1$ and rotate slowly. We show that in ring-shaped disks, the phase space density of solid particles increases linearly in time towards an accumulation point near the location of pressure maximum, while instabilities grow exponentially. Therefore, planetesimals and planetary cores can be efficiently produced near the peaks of unstable density waves. In this mechanism, particles migrating towards the accumulation point will not participate in the formation of planets, and should eventually form a debris ring like the main asteroid belt or classical Kuiper belt objects. We present the implications of global instabilities to the formation of ice giants and terrestrial planets in the solar system.
rate research
Read More
The recently discovered resonant drag instability of dust settling 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: the first mode is inertial wave having positive energy, while the second mode is a settling 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 settling of dust, which transports perturbations of dust density. The phase velocity of SDW is equal to the bulk settling velocity times the cosine of the angle formed by the wave vector and the rotation axis. In this way, the bulk settling of dust makes possible the coupling of SDW with the inertial wave and the onset of the instability. In accordance with the concept of the mode coupling, the instability growth rate is proportional to the square root of the dispersion equation coupling term, which itself contains the small mass fraction of dust in gas-dust mixture, the squared radial wavenumber of the modes, and the squared bulk settling velocity. Thus, the higher is the bulk settling velocity, the heavier clumps of dust can be aggregated by the instability of the same rate.
Under the right conditions, the streaming instability between imperfectly coupled dust and gas is a powerful mechanism for planetesimal formation as it can concentrate dust grains to the point of gravitational collapse. In its simplest form, the streaming instability can be captured by analyzing the linear stability of unstratified disk models, which represent the midplane of protoplanetary disks. We extend such studies by carrying out vertically-global linear stability analyses of dust layers in protoplanetary disks. We find the dominant form of instability in stratified dust layers is one driven by the vertical gradient in the rotation velocity of the dust-gas mixture, but also requires partial dust-gas coupling. These vertically-shearing streaming instabilities grow on orbital timescales and occur on radial length scales $sim10^{-3}H_mathrm{g}$, where $H_mathrm{g}$ is the local pressure scale height. The classic streaming instability, associated with the relative radial drift between dust and gas, occur on radial length scales $sim10^{-2}H_mathrm{g}$, but have much smaller growth rates than vertically-shearing streaming instabilities. Including gas viscosity is strongly stabilizing and leads to vertically-elongated disturbances. We briefly discuss the potential effects of vertically-shearing streaming instabilities on planetesimal formation.
A global evolution picture of protoplanetary disks (PPDs) is key to understanding almost every aspect of planet formation, where standard alpha-disk models have been constantly employed for its simplicity. In the mean time, disk mass loss has been conventionally attributed to photoevaporation, which controls disk dispersal. However, a paradigm shift towards accretion driven by magnetized disk winds has been realized in the recent years, thanks to studies of non-ideal magneto-hydrodynamic effects in PPDs. I present a framework of global PPD evolution aiming to incorporate these advances, highlighting the role of wind-driven accretion and wind mass loss. Disk evolution is found to be largely dominated by wind-driven processes, and viscous spreading is suppressed. The timescale of disk evolution is controlled primarily by the amount of external magnetic flux threading the disks, and how rapidly the disk loses the flux. Rapid disk dispersal can be achieved if the disk is able to hold most of its magnetic flux during the evolution. In addition, because wind launching requires sufficient level of ionization at disk surface (mainly via external far-UV radiation), wind kinematics is also affected by far-UV penetration depth and disk geometry. For typical disk lifetime of a few Myrs, the disk loses approximately the same amount of mass through the wind as through accretion onto the protostar, and most of the wind mass loss proceeds from the outer disk via a slow wind. Fractional wind mass loss increases with increasing disk lifetime. Significant wind mass loss likely substantially enhances the dust to gas mass ratio, and promotes planet formation.
In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two components there. We address this problem by means of numerical 2.5D simulations in order to study the gas and dust interaction in fully global stratified disks. To this purpose, we employ the recently developed dust grain module in the PLUTO code. Our model focuses on a typical T Tauri disk model, simulating a short patch of the disk at 10 au which includes grains of constant Stokes number of $St=0.01$ and $St=0.1$, corresponding to grains with sizes of 0.9 cm and 0.9 mm, respectively, for the given disk model. By injecting a constant pebble flux at the outer domain, the system reaches a quasi steady state of turbulence and dust concentrations driven by the streaming instability. For our given setup and using resolutions up to 2500 cells per scale height we resolve the streaming instability, leading to local dust clumping and concentrations. Our results show dust density values of around 10-100 times the gas density with a steady state pebble flux between $3.5 times 10^{-4}$ and $2.5 times 10^{-3} M_{rm Earth}/mathit{year}$ for the models with $mathit{St}=0.01$ and $mathit{St}=0.1$. The grain size and pebble flux for model $mathit{St}=0.01$ compares well with dust evolution models of the first million years of disk evolution. For those grains the scatter opacity dominates the extinction coefficient at mm wavelengths. These types of global dust and gas simulations are a promising tool for studies of the gas and dust evolution at pressure bumps in protoplanetary disks.
In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $sim 1.0 {rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($sim 4.0$ mm) being a factor of $sim10$ larger than other vortex regions. Both large ($sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
