Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userParticle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks
84
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low mass planet ($8 M_{oplus}$), two narrow gaps start to open in the gas on each side of the planet where the density waves shock. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet ($>0.1 M_{J}$) carves out a deeper gap with sharp edges, which are unstable to the vortex formation. Particles with a wide range of sizes ($0.02<Omega t_{s}<20$) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with $Omega t_{s}sim 1$. The dust concentration is highly elongated in the $phi$ direction, and can be as wide as 4 disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10$^2$ in a very non-axisymmetric fashion. For very big particles ($Omega t_{s}gg 1$) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap opening there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g. the dead zone) dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with $M_{p}<M_{J}$, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and EVLA.
rate research
Read More
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.
Recent sub-millimeter observations show non-axisymmetric brightness distributions with a horseshoe-like morphology for more than a dozen transition disks. The most accepted explanation for the observed asymmetries is the accumulation of dust in large-scale vortices. Protoplanetary disks vortices can form by the excitation of Rossby-wave instability in the vicinity of a steep pressure gradient, which can develop at the edges of a giant planet carved gap or at the edges of an accretionally inactive zone. We studied the formation and evolution of vortices formed in these two distinct scenarios by means of two-dimensional locally isothermal hydrodynamic simulations. We found that the vortex formed at the edge of a planetary gap is short-lived, unless the disk is nearly inviscid. In contrast, the vortex formed at the outer edge of a dead zone is long-lived. The vortex morphology can be significantly different in the two scenarios: the vortex radial and azimuthal extensions are ~1.5 and ~3.5 times larger for the dead zone edge compared to gap models. In some particular cases, the vortex aspect ratios can be similar in the two scenarios, however, the vortex azimuthal extensions can be used to distinguish the vortex formation mechanisms. We calculate predictions for vortex observability in the sub-millimeter continuum with ALMA. We found that the azimuthal and radial extent of brightness asymmetry correlates with vortex formation process, within the limitations of alpha-viscosity prescription.
Dust gaps and rings appear ubiquitous in bright protoplanetary disks. Disk-planet interaction with dust-trapping at the edges of planet-induced gaps is one plausible explanation. However, the sharpness of some observed dust rings indicate that sub-mm-sized dust grains have settled to a thin layer in some systems. We test whether or not such dust around gas gaps opened by planets can remain settled by performing three-dimensional, dust-plus-gas simulations of protoplanetary disks with an embedded planet. We find planets massive enough to open gas gaps stir small, sub-mm-sized dust grains to high disk elevations at the gap edges, where the dust scale-height can reach ~70% of the gas scale-height. We attribute this dust puff-up to the planet-induced meridional gas flows previously identified by Fung & Chiang and others. We thus emphasize the importance of explicit 3D simulations to obtain the vertical distribution of sub-mm-sized grains around gas gaps opened by massive planets. We caution that the gas-gap-opening planet interpretation of well-defined dust rings is only self-consistent with large grains exceeding mm in size.
The presence of rings and gaps in protoplanetary discs are often ascribed to planet-disc interactions, where dust and pebbles are trapped at the edges of planetary induced gas gaps. Recent work has shown that these are likely sites for planetesimal formation via the streaming instability. Given the large amount of planetesimals that potentially form at gap edges, we address the question of their fate and ability to radially transport solids in protoplanetary discs. We perform a series of N-body simulations of planetesimal orbits, taking into account the effect of gas drag and mass loss via ablation. We consider two planetary systems: one akin to the young Solar System, and another one inspired by HL Tau. In both systems, the close proximity to the gap-opening planets results in large orbital excitations, causing the planetesimals to leave their birth locations and spread out across the disc soon after formation. Planetesimals that end up on eccentric orbits interior of 10au experience efficient ablation, and lose all mass before they reach the innermost disc region. In our nominal Solar System simulation with $dot{M}_0=10^{-7}, M_{odot}, textrm{yr}^{-1}$ and $alpha=10^{-2}$, we find that 70% of the initial planetesimal mass has been ablated after 500kyr. The ablation rate in HL Tau is lower, and only 11% of the initial planetesimal mass has been ablated after 1Myr. The ablated material consist of a mixture of solid grains and vaporized ices, where a large fraction of the vaporized ices re-condense to form solid ice. Assuming that the solids grow to pebbles in the disc midplane, this results in a pebble flux of $sim 10-100,M_{oplus}textrm{Myr}^{-1}$ through the inner disc. Our results demonstrate that scattered planetesimals can carry a significant flux of solids past planetary-induced gaps in young and massive protoplanetary discs.
The Streaming Instability (SI) is a mechanism to concentrate solids in protoplanetary disks. Nonlinear particle clumping from the SI can trigger gravitational collapse into planetesimals. To better understand the numerical robustness of the SI, we perform a suite of vertically-stratified 3D simulations with fixed physical parameters known to produce strong clumping. We vary the numerical implementation, namely the computational domain size and the vertical boundary conditions (vBCs), comparing newly-implemented outflow vBCs to the previously-used periodic and reflecting vBCs. We find strong particle clumping by the SI is mostly independent of the vBCs. However, peak particle densities are higher in larger simulation domains due to a larger particle mass reservoir. We report SI-triggered zonal flows, i.e., azimuthally-banded radial variations of gas pressure. These structures have low amplitudes, insufficient to halt particle radial drift, confirming that particle trapping in gas pressure maxima is not the mechanism of the SI. We find that outflow vBCs produce artificially large gas outflow rates at vertical boundaries. However, the outflow vBCs reduce artificial reflections at vertical boundaries, allowing more particle sedimentation, and showing less temporal variation and better convergence with box size. The radial spacing of dense particle filaments is $sim0.15$ gas scale heights ($H$) for all vBCs, which sets the feeding zone for planetesimal growth in self-gravitating simulations. Our results validate the use of the outflow vBCs in SI simulations, even with vertical boundaries close ($leq 0.4H$) to the disk midplane. Overall, our study demonstrates the numerical robustness of nonlinear particle clumping by the SI.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
