Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDust settling and rings in the outer regions of protoplanetary discs subject to ambipolar diffusion
184
0
0.0
(
0
)
Added by
Antoine Riols-Fonclare
Publication date
2018
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
MHD turbulence plays a crucial role in the dust dynamics of protoplanetary discs. It affects planet formation, vertical settling and is one possible origin of the large scale axisymmetric structures, such as rings, recently imaged by ALMA and SPHERE. Among the variety of MHD processes, the magnetorotational instability (MRI) has raised particular interest since it provides a source of turbulence and potentially organizes the flow into radial structures. However, the weak ionization of discs prevents the MRI from being excited beyond 1 AU. The strong sedimentation of millimetre dust measured in T-Tauri discs is also in contradiction with predictions based on ideal MRI turbulence. In this paper, we study the effects of non-ideal MHD and winds on the dynamics and sedimentation of dust grains. We consider a weakly ionized plasma subject to ambipolar diffusion characterizing the disc outer regions (>1 AU). For that, we perform numerical MHD simulations in the stratified shearing box, using the PLUTO code. Our simulations show that the mm-cm dust is contained vertically in a very thin layer, with typical heightscale ~0.4 AU at 30 AU, compatible with recent ALMA observations. Horizontally, the grains are trapped within pressure maxima induced by ambipolar diffusion, leading to the formation of dust rings. For micrometer grains, dust and gas scaleheights are similar. In this regime, the settling cannot be explained by a simple 1D diffusion theory but results from a large scale 2D circulation induced by both MHD winds and zonal flows. Overall, our results show that non-ideal MHD effects and their related winds play a major role in shaping the radial and vertical distribution of dust in protoplanetary discs. Leading to substantial accretion efficiency, non-ideal effects also a promising avenue to reconcile the low turbulent activity measured in discs with their relatively high accretion rates.
rate research
Read More
Gravitational instability (GI) controls the dynamics of young massive protoplanetary discs. Apart from facilitating gas accretion on to the central protostar, it must also impact on the process of planet formation: directly through fragmentation, and indirectly through the turbulent concentration of small solids. To understand the latter process, it is essential to determine the dust dynamics in such a turbulent flow. For that purpose, we conduct a series of 3D shearing box simulations of coupled gas and dust, including the gass self-gravity and scanning a range of Stokes numbers, from 0.001 to ~0.2. First, we show that the vertical settling of dust in the midplane is significantly impeded by gravitoturbulence, with the dust scale-height roughly 0.6 times the gas scale height for centimetre grains. This is a result of the strong vertical diffusion issuing from (a) small-scale inertial-wave turbulence feeding off the GI spiral waves and (b) the larger-scale vertical circulations that naturally accompany the spirals. Second, we show that at R=50 AU concentration events involving sub-metre particles and yielding order 1 dust to gas ratios are rare and last for less than an orbit. Moreover, dust concentration is less efficient in 3D than in 2D simulations. We conclude that GI is not especially prone to the turbulent accumulation of dust grains. Finally, the large dust scale-height measured in simulations could be, in the future, compared with that of edge-on discs seen by ALMA, thus aiding detection and characterisation of GI in real systems.
The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could trigger particle clumping away from the midplane. This work presents a range of linear calculations and non-linear simulations, performed with FARGO3D, to assess conditions for DSI growth. We expand on previous linear analyses by including particle size distributions and performing a detailed study of the amount of background turbulence needed to stabilize the DSI. When including binned size distributions, the DSI often produces converged growth rates with fewer bins than the standard SI. With background turbulence, we find that the most favorable conditions for DSI growth are weak turbulence, characterized by $alpha lesssim 10^{-6}$ with intermediate-sized grains that settle from one gas scale-height. These conditions could arise during a sudden decrease in disc turbulence following an accretion outburst. Ignoring background turbulence, we performed a parameter survey of local 2D DSI simulations. Particle clumping was either weak or occurred slower than particles settle. Clumping was reduced by a factor of two in a comparison 3D simulation. Overall, our results strongly disfavor the hypothesis that the DSI significantly promotes planetesimal formation. Non-linear simulations of the DSI with different numerical methods could support or challenge these findings.
Enhancing the local dust-to-gas ratio in protoplanetary discs is a necessary first step to planetesimal formation. In laminar discs, dust settling is an efficient mechanism to raise the dust-to-gas ratio at the disc midplane. However, turbulence, if present, can stir and lift dust particles, which ultimately hinders planetesimal formation. In this work, we study dust settling in protoplanetary discs with hydrodynamic turbulence sustained by the vertical shear instability. We perform axisymmetric numerical simulations to investigate the effect of turbulence, particle size, and solid abundance or metallicity on dust settling. We highlight the positive role of drag forces exerted onto the gas by the dust for settling to overcome the vertical shear instability. In typical disc models we find particles with a Stokes number $sim 10^{-3}$ can sediment to $lesssim 10%$ of the gas scale-height, provided that $Sigma_mathrm{d}/Sigma_mathrm{g}gtrsim 0.02$-$0.05$, where $Sigma_mathrm{d,g}$ are the surface densities in dust and gas, respectively. This coincides with the metallicity condition for small particles to undergo clumping via the streaming instability. Super-solar metallicities, at least locally, are thus required for a self-consistent picture of planetesimal formation. Our results also imply that dust rings observed in protoplanetary discs should have smaller scale-heights than dust gaps, provided that the metallicity contrast between rings and gaps exceed the corresponding contrast in gas density.
Planetesimal formation is a crucial yet poorly understood process in planet formation. It is widely believed that planetesimal formation is the outcome of dust clumping by the streaming instability (SI). However, recent analytical and numerical studies have shown that the SI can be damped or suppressed by external turbulence, and at least the outer regions of protoplanetary disks are likely weakly turbulent due to magneto-rotational instability (MRI). We conduct high-resolution local shearing-box simulations of hybrid particle-gas magnetohydrodynamics (MHD), incorporating ambipolar diffusion as the dominant non-ideal MHD effect, applicable to outer disk regions. We first show that dust backreaction enhances dust settling towards the midplane by reducing turbulence correlation time. Under modest level of MRI turbulence, we find that dust clumping is in fact easier than the conventional SI case, in the sense that the threshold of solid abundance for clumping is lower. The key to dust clumping includes dust backreaction and the presence of local pressure maxima, which in our work is formed by the MRI zonal flows overcoming background pressure gradient. Overall, our results support planetesimal formation in the MRI-turbulent outer protoplanetary disks, especially in ring-like substructures.
The structure and evolution of protoplanetary disks (PPDs) are largely governed by disk angular momentum transport, mediated by magnetic fields. In the most observable outer disk, PPD gas dynamics is primarily controlled by ambipolar diffusion as the dominant non-ideal magnetohydrodynamic (MHD) effect. In this work, we study the gas dynamics in outer PPDs by conducting a set of global 3D non-ideal MHD simulations with ambipolar diffusion and net poloidal magnetic flux, using the Athena++ MHD code, with resolution comparable to local simulations. Our simulations demonstrate the co-existence of magnetized disk wind and turbulence driven by the magneto-rotational instability (MRI). While MHD wind dominates disk angular momentum transport, the MRI turbulence also contributes significantly. We observe that magnetic flux spontaneously concentrate into axisymmetric flux sheets, leading to radial variations in turbulence levels, stresses, and accretion rates. Annular substructures arise as a natural consequence of magnetic flux concentration. The flux concentration phenomena show diverse properties with different levels of disk magnetization and ambipolar diffusion. The disk generally loses magnetic flux over time, though flux sheets could prevent the leak of magnetic flux in some cases. Our results demonstrate the ubiquity of disk annular substructures in weakly MRI turbulent outer PPDs, and imply a stochastic nature of disk evolution.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
