No Arabic abstract
We present global hydrodynamic and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (i) The accretion flow in the Hill sphere is intrinsically 3D for hydrodynamic and MHD models. Net inflow toward the planet is dominated by high latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (ii) Gap opening magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere and a CPD whose structure displays high levels of variability. (iii) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch a highly collimated, one-sided jet is observed. (iv) The CPDs surface density $sim30{rm,g,cm^{-2}}$, small enough for significant ionization and turbulence to develop. (v) The accretion rate onto the planet in the MHD simulation reaches a steady value $8 times 10^{-3} {rm M_oplus yr^{-1}}$, and is similar in the viscous hydrodynamic runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.
Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.
High contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their masses, from the observed gap properties, such as their depths and widths. Combing 2D and 3D hydrodynamics simulations with 3D radiative transfer simulations, we investigate the morphology of planet-opened gaps in near-infrared scattered light images. Quantitatively, we obtain correlations that directly link intrinsic gap depths and widths in the gas surface density to observed depths and widths in images of disks at modest inclinations under finite angular resolution. Subsequently, the properties of the surface density gaps enable us to derive the disk scale height at the location of the gap $h$, and to constrain the quantity $M_{rm p}^2/alpha$, where $M_{rm p}$ is the mass of the gap-opening planet and $alpha$ characterizes the viscosity in the gap. As examples, we examine the gaps recently imaged by VLT/SPHERE, Gemini/GPI, and Subaru/HiCIAO in HD 97048, TW Hya, HD 169142, LkCa 15, and RX J1615.3-3255. Scale heights of the disks and possible masses of the gap-opening planets are derived assuming each gap is opened by a single planet. Assuming $alpha=10^{-3}$, the derived planet mass in all cases are roughly between 0.1-1 $M_{rm J}$.
We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disks self-gravity affects the gap formation process and the migration of the giant planet. Two series of 1-D and 2-D hydrodynamic simulations are performed. We select several surface densities and focus on the gravitationally stable region. To obtain more reliable gravity torques exerted on the planet, a refined treatment of disks gravity is adopted in the vicinity of the planet. Our results indicate that the net effect of the disks self-gravity on the gap formation process depends on the surface density of the disk. We notice that there are two critical values, Sigma_I and Sigma_II. When the surface density of the disk is lower than the first one, Sigma_0 < Sigma_I, the effect of self-gravity suppresses the formation of a gap. When Sigma_0 > Sigma_I, the self-gravity of the gas tends to benefit the gap formation process and enlarge the width/depth of the gap. According to our 1-D and 2-D simulations, we estimate the first critical surface density Sigma_I approx 0.8MMSN. This effect increases until the surface density reaches the second critical value Sigma_II. When Sigma_0 > Sigma_II, the gravitational turbulence in the disk becomes dominant and the gap formation process is suppressed again. Our 2-D simulations show that this critical surface density is around 3.5MMSN. We also study the associated orbital evolution of a giant planet. Under the effect of the disks self-gravity, the migration rate of the giant planet increases when the disk is dominated by gravitational turbulence. We show that the migration timescale associates with the effective viscosity and can be up to 10^4 yr.
We carry out three-dimensional hydrodynamical simulations to study planet-disc interactions for inclined high mass planets, focusing on the discs secular evolution induced by the planet. We find that, when the planet is massive enough and the induced gap is deep enough, the disc inside the planets orbit breaks from the outer disc. The inner and outer discs precess around the systems total angular momentum vector independently at different precession rates, which causes significant disc misalignment. We derive the analytical formulae, which are also verified numerically, for: 1) the relationship between the planet mass and the depth/width of the induced gap, 2) the migration and inclination damping rates for massive inclined planets, and 3) the condition under which the inner and outer discs can break and undergo differential precession. Then, we carry out Monte-Carlo radiative transfer calculations for the simulated broken discs. Both disc shadowing in near-IR images and gas kinematics probed by molecular lines (e.g. from ALMA) can reveal the misaligned inner disc. The relationship between the rotation rate of the disc shadow and the precession rate of the inner disc is also provided. Using our disc breaking condition, we conclude that the disc shadowing due to misaligned discs should be accompanied by deep gaseous gaps (e.g. in Pre/Transitional discs). This scenario naturally explains both the disc shadowing and deep gaps in several systems (e.g. HD 100453, DoAr 44, AA Tau, HD 143006) and these systems should be the prime targets for searching young massive planets ($>M_J$) in discs.
Recent observations of protoplanetary disks have revealed ring-like structures that can be associated to pressure maxima. Pressure maxima are known to be dust collectors and planet migration traps. Most of planet formation works are based either on the pebble accretion model or on the planetesimal accretion model. However, recent studies proposed the possible formation of Jupiter by the hybrid accretion of pebbles and planetesimals. We aim to study the full process of planet formation consisting of dust evolution, planetesimal formation and planet growth at a pressure maximum in a protoplanetary disk. We compute, through numerical simulations, the gas and dust evolution, including dust growth, fragmentation, radial drift and particle accumulation at a pressure bump. We also consider the formation of planetesimals by streaming instability and the formation of a moon-size embryo that grows into a giant planet by the hybrid accretion of pebbles and planetesimals. We find that pressure maxima in protoplanetary disks are efficient collectors of dust drifting inwards. The condition of planetesimal formation by streaming instability is fulfilled due to the large amount of dust accumulated at the pressure bump. Then, a massive core is quickly formed (in $sim 10^4$ yr) by the accretion of pebbles. After the pebble isolation mass is reached, the growth of the core slowly continues by the accretion of planetesimals. The energy released by planetesimal accretion delays the onset of runaway gas accretion, allowing a gas giant to form after $sim$1 Myr of disk evolution. The pressure maximum also acts as a migration trap. Pressure maxima in protoplanetary disks are preferential locations for dust traps, planetesimal formation by streaming instability and planet migration traps. All these conditions allow the fast formation of a giant planet by the hybrid accretion of pebbles and planetesimals.