No Arabic abstract
As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planets luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planets luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.
A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a new source that can generate outward migration of low-mass planets. While the Lindblad and corotation torques depend mostly on the properties of the protoplanetary disc and on the planet mass, the thermal torque depends also on the luminosity of the planet, arising mainly from the accretion of solids. Thus, the accretion of solids plays an important role not only in the formation of the planet but also in its migration process. In a previous work, we evaluated the thermal torque effects on planetary growth and migration mainly in the planetesimal accretion paradigm. In this new work, we study the role of the thermal torque within the pebble accretion paradigm. Computations are carried out consistently in the framework of a global model of planet formation that includes disc evolution, dust growth and evolution, and pebble formation. We also incorporate updated prescriptions of the thermal torque derived from high resolution hydrodynamical simulations. Our simulations show that the thermal torque generates extended regions of outward migration in low viscosity discs. This has a significant impact in the formation of the planets.
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.
Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii with zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.
We study torques on migrating low-mass planets in locally isothermal discs. Previous work on low-mass planets generally kept the planet on a fixed orbit, after which the torque on the planet was measured. In addition to these static torques, when the planet is allowed to migrate it experiences dynamical torques, which are proportional to the migration rate and whose sign depends on the background vortensity gradient. We show that in discs a few times more massive than the Minimum Mass Solar Nebula, these dynamical torques can have a profound impact on planet migration. Inward migration can be slowed down significantly, and if static torques lead to outward migration, dynamical torques can take over, taking the planet beyond zero-torque lines set by saturation of the corotation torque in a runaway fashion. This means the region in non-isothermal discs where outward migration is possible can be larger than what would be concluded from static torques alone.
Using linear perturbation theory, we investigate the torque exerted on a low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity. When the planet does not release energy into the ambient disc, the main effect of thermal diffusion is the softening of the enthalpy peak near the planet, which results in the appearance of two cold and dense lobes on either side of the orbit, of size smaller than the thickness of the disc. The lobes exert torques of opposite sign on the planet, each comparable in magnitude to the one-sided Lindblad torque. When the planet is offset from corotation, the lobes are asymmetric and the planet experiences a net torque, the `cold thermal torque, which has a magnitude that depends on the relative value of the distance to corotation to the size of the lobes $simsqrt{chi/Omega_p}$, $chi$ being the thermal diffusivity and $Omega_p$ the orbital frequency. We believe that this effect corresponds to the phenomenon named `cold finger recently reported in numerical simulations, and we argue that it constitutes the dominant mode of migration of sub-Earth-mass objects. When the planet is luminous, the heat released into the ambient disc results in an additional disturbance that takes the form of hot, low-density lobes. They give a torque, named heating torque in previous work, that has an expression similar, but of opposite sign, to the cold thermal torque.