skip to main content

Title: Gap opening in protoplanetary discs: gas dynamics from global axisymmetric non-ideal MHD simulations with consistent thermochemistry

Recent high angular resolution ALMA observations have revealed numerous gaps in protoplanetary discs. A popular interpretation has been that planets open them. Most previous investigations of planet gap-opening have concentrated on viscous discs. Here, we carry out 2D (axisymmetric) global simulations of gap opening by a planet in a wind-launching non-ideal MHD disc with consistent thermochemistry. We find a strong concentration of poloidal magnetic flux in the planet-opened gap, where the gas dynamics are magnetically dominated. The magnetic field also drives a fast (nearly sonic) meridional gas circulation in the denser disc regions near the inner and outer edges of the gap, which may be observable through high-resolution molecular line observations. The gap is more ionized than its denser surrounding regions, with a better magnetic field–matter coupling. In particular, it has a much higher abundance of molecular ion HCO+, consistent with ALMA observations of the well-studied AS 209 protoplanetary disc that has prominent gaps and fast meridional motions reaching the local sound speed. Finally, we provide fitting formulae for the ambipolar and Ohmic diffusivities as a function of the disc local density, which can be used for future 3D simulations of planet gap-opening in non-ideal MHD discs where thermochemistry is too computationally expensive to evolve self-consistently with the magneto-hydrodynamics.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 4883-4894
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    Rings and gaps are commonly observed in the dust continuum emission of young stellar discs. Previous studies have shown that substructures naturally develop in the weakly ionized gas of magnetized, non-ideal MHD discs. The gas rings are expected to trap large mm/cm-sized grains through pressure gradient-induced radial dust–gas drift. Using 2D (axisymmetric) MHD simulations that include ambipolar diffusion and dust grains of three representative sizes (1 mm, 3.3 mm, and 1 cm), we show that the grains indeed tend to drift radially relative to the gas towards the centres of the gas rings, at speeds much higher than in a smooth disc because of steeper pressure gradients. However, their spatial distribution is primarily controlled by meridional gas motions, which are typically much faster than the dust–gas drift. In particular, the grains that have settled near the mid-plane are carried rapidly inwards by a fast accretion stream to the inner edges of the gas rings, where they are lifted up by the gas flows diverted away from the mid-plane by a strong poloidal magnetic field. The flow pattern in our simulation provides an attractive explanation for the meridional flows recently inferred in HD 163296 and other discs, including both ‘collapsing’ regions where the gas near the disc surface converges towards the mid-plane and a disc wind. Our study highlights the prevalence of the potentially observable meridional flows associated with the gas substructure formation in non-ideal MHD discs and their crucial role in generating rings and gaps in dust.

    more » « less

    Gaps imaged in protoplanetary discs are suspected to be opened by planets. We compute the present-day mass accretion rates $\dot{M}_{\rm p}$ of seven hypothesized gap-embedded planets, plus the two confirmed planets in the PDS 70 disc. The accretion rates are based on disc gas surface densities Σgas from C18O observations, and planet masses Mp from simulations fitted to observed gaps. Assuming accretion is Bondi-like, we find in eight out of nine cases that $\dot{M}_{\rm p}$ is consistent with the time-averaged value given by the current planet mass and system age, Mp/tage. As system ages are comparable to circumstellar disc lifetimes, these gap-opening planets may be undergoing their last mass doublings, reaching final masses of $M_{\rm p} \sim 10\rm{\!-\!}10^2 \, M_\oplus$ for the non-PDS 70 planets, and $M_{\rm p} \sim 1\!-\!10 \, M_{\rm J}$ for the PDS 70 planets. For another 15 gaps without C18O data, we predict Σgas by assuming their planets are accreting at their time-averaged $\dot{M}_{\rm p}$. Bondi accretion rates for PDS 70b and c are orders of magnitude higher than accretion rates implied by measured U-band and H α fluxes, suggesting most of the accretion shock luminosity emerges in as yet unobserved wavebands, or that the planets are surrounded by dusty, highly extincting, quasi-spherical circumplanetary envelopes. Thermal emission from such envelopes or from circumplanetary discs, on Hill sphere scales, peaks at wavelengths in the mid-to-far-infrared and can reproduce observed mm-wave excesses.

    more » « less
  3. null (Ed.)
    ABSTRACT Recent ALMA molecular line observations have revealed 3D gas velocity structure in protoplanetary discs, shedding light on mechanisms of disc accretion and structure formation. (1) By carrying out viscous simulations, we confirm that the disc’s velocity structure differs dramatically using vertical stress profiles from different accretion mechanisms. Thus, kinematic observations tracing flows at different disc heights can potentially distinguish different accretion mechanisms. On the other hand, the disc surface density evolution is mostly determined by the vertically integrated stress. The sharp disc outer edge constrained by recent kinematic observations can be caused by a radially varying α in the disc. (2) We also study kinematic signatures of a young planet by carrying out 3D planet–disc simulations. The relationship between the planet mass and the ‘kink’ velocity is derived, showing a linear relationship with little dependence on disc viscosity, but some dependence on disc height when the planet is massive (e.g. 10MJ). We predict the ‘kink’ velocities for the potential planets in DSHARP discs. At the gap edge, the azimuthally averaged velocities at different disc heights deviate from the Keplerian velocity at similar amplitudes, and its relationship with the planet mass is consistent with that in 2D simulations. After removing the planet, the azimuthally averaged velocity barely changes within the viscous time-scale, and thus the azimuthally averaged velocity structure at the gap edge is due to the gap itself and not directly caused to the planet. Combining both axisymmetric kinematic observations and the residual ‘kink’ velocity is needed to probe young planets in protoplanetary discs. 
    more » « less

    Radial substructures have now been observed in a wide range of protoplanetary discs (PPDs), from young to old systems; however, their formation is still an area of vigorous debate. Recent magnetohydrodynamic (MHD) simulations have shown that rings and gaps can form naturally in PPDs when non-ideal MHD effects are included. However, these simulations employ ad hoc approximations to the magnitudes of the magnetic diffusivities in order to facilitate ring growth. We replace the parametrization of these terms with a simple chemical network and grain distribution model to calculate the non-ideal effects in a more self-consistent way. We use a range of grain distributions to simulate grain formation for different disc conditions. Including ambipolar diffusion, we find that large grain populations (>1 $\mu$m), and those including a population of very small polyaromatic hydrocarbons (PAHs) facilitate the growth of periodic, stable rings, while intermediate-sized grains suppress ring formation. Including Ohmic diffusion removes the positive influence of PAHs, with only large grain populations still producing periodic ring and gap structures. These results relate closely to the degree of coupling between the magnetic field and the neutral disc material, quantified by the non-dimensional Elsasser number Λ (the ratio of magnetic forces to Coriolis force). For both the ambipolar-only and ambipolar-ohmic cases, if the total Elsasser number is initially of the order of unity along the disc mid-plane, ring and gap structures may develop.

    more » « less
  5. Abstract Recent observations have revealed a gallery of substructures in the dust component of nearby protoplanetary discs, including rings, gaps, spiral arms, and lopsided concentrations. One interpretation of these substructures is the existence of embedded planets. Not until recently, however, most of the modelling effort to interpret these observations ignored the dust back reaction to the gas. In this work, we conduct local-shearing-sheet simulations for an isothermal, inviscid, non-self-gravitating, razor-thin dusty disc with a planet on a fixed circular orbit. We systematically examine the parameter space spanned by planet mass (0.1Mth ≤ Mp ≤ 1Mth, where Mth is the thermal mass), dimensionless stopping time (10−3 ≤ τs ≤ 1), and solid abundance (0 < Z ≤ 1). We find that when the dust particles are tightly coupled to the gas (τs < 0.1), the spiral arms are less open and the gap driven by the planet becomes deeper with increasing Z, consistent with a reduced speed of sound in the approximation of a single dust-gas mixture. By contrast, when the dust particles are marginally coupled (0.1 ≲ τs ≲ 1), the spiral structure is insensitive to Z and the gap structure in the gas can become significantly skewed and unidentifiable. When the latter occurs, the pressure maximum radially outside of the planet is weakened or even extinguished, and hence dust filtration by a low-mass (Mp < Mth) planet could be reduced or eliminated. Finally, we find that the gap edges where the dust particles are accumulated as well as the lopsided large-scale vortices driven by a massive planet, if any, are unstable, and they are broken into numerous small-scale dust-gas vortices. 
    more » « less