No Arabic abstract
We carry out three-dimensional SPH simulations to study whether planets can survive in self-gravitating protoplanetary discs. The discs modelled here use a cooling prescription that mimics a real disc which is only gravitationally unstable in the outer regions. We do this by modelling the cooling using a simplified method such that the cooling time in the outer parts of the disc is shorter than in the inner regions, as expected in real discs. We find that both giant (> M_Sat) and low mass (< M_Nep) planets initially migrate inwards very rapidly, but are able to slow down in the inner gravitationally stable regions of the disc without needing to open up a gap. This is in contrast to previous studies where the cooling was modelled in a more simplified manner where regardless of mass, the planets were unable to slow down their inward migration. This shows the important effect the thermodynamics has on planet migration. In a broader context, these results show that planets that form in the early stages of the discs evolution, when they are still quite massive and self-gravitating, can survive.
Large-scale, dust-trapping vortices may account for observations of asymmetric protoplanetary discs. Disc vortices are also potential sites for accelerated planetesimal formation by concentrating dust grains. However, in 3D discs vortices are subject to destructive `elliptic instabilities, which reduces their viability as dust traps. The survival of vortices in 3D accretion discs is thus an important issue to address. In this work, we perform shearing box simulations to show that disc self-gravity enhances the survival of 3D vortices, even when self-gravity is weak in the classic sense (e.g. with a Toomre $Qsimeq5$). We find a 3D, self-gravitating vortex can grow on secular timescales in spite of the elliptic instability. The vortex aspect-ratio decreases as it strengthens, which feeds the elliptic instability. The result is a 3D vortex with a turbulent core that persists for $sim 10^{3}$ orbits. We find when gravitational and hydrodynamic stresses become comparable, the vortex may undergo episodic bursts, which we interpret as interaction between elliptic and gravitational instabilities. We estimate the distribution of dust particles in self-gravitating, turbulent vortices. Our results suggest large-scale vortices in protoplanetary discs are more easily observed at large radii.
Type-II migration of giant planets has a speed proportional to the discs viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instability at the edges of the gap opened by the planet. Migration is then determined by vortex-planet interactions. Our aim is to study migration in low viscosity 3D discs. We performed numerical simulations using 2D (including self-gravity) and 3D codes. After selecting disc masses for which self-gravity is not important, 3D simulations without self-gravity can be safely used. In our nominal simulation only numerical viscosity is present. We then performed simulations with prescribed viscosity to assess the threshold below which the new migration processes appear. We show that for alpha viscosity <= 1.e-5 two migration modes are possible which differ from classical Type-II migration, in the sense that they are not proportional to the discs viscosity. The first occurs when the gap opened by the planet is not very deep. This occurs in 3D simulations and/or when a big vortex forms at the outer edge of the planetary gap, diffusing material into the gap. We call this type of migration vortex-driven migration. This migration is very slow and cannot continue indefinitely, because eventually the vortex dissolves. The second migration mode occurs when the gap is deep so that the planets eccentricity grows to a value ~0.2 due to inefficient eccentricity damping by corotation resonances. This second, faster migration mode appears to be typical of 2D models in discs with slower damping of temperatures perturbations.
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.
Planets form in the discs of gas and dust that surround young stars. It is not known whether gas giant planets on wide orbits form the same way as Jupiter or by fragmentation of gravitationally unstable discs. Here we show that a giant planet, which has formed in the outer regions of a protostellar disc, initially migrates fast towards the central star (migration timescale ~10,000 yr) while accreting gas from the disc. However, in contrast with previous studies, we find that the planet eventually opens up a gap in the disc and the migration is essentially halted. At the same time, accretion-powered radiative feedback from the planet, significantly limits its mass growth, keeping it within the planetary mass regime (i.e. below the deuterium burning limit) at least for the initial stages of disc evolution. Giant planets may therefore be able to survive on wide orbits despite their initial fast inward migration, shaping the environment in which terrestrial planets that may harbour life form.
We present the results of hydrodynamical simulations of the orbital evolution of planets undergoing runaway gas accretion in radiative discs. We consider accreting disc models with constant mass flux through the disc, and where radiative cooling balances the effect of viscous heating and stellar irradiation. We assume that 20-30 $M_oplus$ giant planet cores are formed in the region where viscous heating dominates and migrate outward under the action of a strong corotation torque. In the case where gas accretion is neglected, we find evidence for strong dynamical torques in accreting discs with accretion rates ${dot M}gtrsim 7times 10^{-8} ;M_odot/yr$. Their main effect is to increase outward migration rates by a factor of $sim 2$ typically. In the presence of gas accretion, however, runaway outward migration is observed with the planet passing through the zero-torque radius and the transition between the viscous heating and stellar heating dominated regimes. The ability for an accreting planet to enter a fast migration regime is found to depend strongly on the planet growth rate, but can occur for values of the mass flux through the disc of ${dot M}gtrsim 5times 10^{-8} ;M_odot/yr$. We find that an episode of runaway outward migration can cause an accreting planet formed in the 5-10 AU region to temporarily orbit at star-planet separations as large as $sim$60-70 AU. However, increase in the amplitude of the Lindblad torque associated with planet growth plus change in the streamline topology near the planet systematically cause the direction of migration to be reversed. Our results indicate that a planet can reach large orbital distances under the combined effect of dynamical torques and gas accretion, but an alternative mechanism is required to explain the presence of massive planets on wide orbits.