Embedded in the gaseous protoplanetary disk, Jupiter and Saturn naturally become trapped in 3:2 resonance and migrate outward. This serves as the basis of the Grand Tack model. However, previous hydrodynamical simulations were restricted to isotherma
l disks, with moderate aspect ratio and viscosity. Here we simulate the orbital evolution of the gas giants in disks with viscous heating and radiative cooling. We find that Jupiter and Saturn migrate outward in 3:2 resonance in modest-mass ($M_{disk} approx M_{MMSN}$, where MMSN is the minimum-mass solar nebula) disks with viscous stress parameter $alpha$ between $10^{-3}$ and $10^{-2} $. In disks with relatively low-mass ($M_{disk} lesssim M_{MMSN}$) , Jupiter and Saturn get captured in 2:1 resonance and can even migrate outward in low-viscosity disks ($alpha le 10^{-4}$). Such disks have a very small aspect ratio ($hsim 0.02-0.03$) that favors outward migration after capture in 2:1 resonance, as confirmed by isothermal runs which resulted in a similar outcome for $h sim 0.02$ and $alpha le 10^{-4}$. We also performed N-body runs of the outer Solar System starting from the results of our hydrodynamical simulations and including 2-3 ice giants. After dispersal of the gaseous disk, a Nice model instability starting with Jupiter and Saturn in 2:1 resonance results in good Solar Systems analogs. We conclude that in a cold Solar Nebula, the 2:1 resonance between Jupiter and Saturn can lead to outward migration of the system, and this may represent an alternative scenario for the evolution of the Solar System.
Planets close to their stars are thought to form farther out and migrate inward due to angular momentum exchange with gaseous protoplanetary disks. This process can produce systems of planets in co-orbital (Trojan or 1:1) resonance, in which two plan
ets share the same orbit, usually separated by 60 degrees. Co-orbital systems are detectable among the planetary systems found by the Kepler mission either directly or by transit timing variations. However, no co-orbital systems have been found within the thousands of Kepler planets and candidates. Here we study the orbital evolution of co-orbital planets embedded in a protoplanetary disk using a grid-based hydrodynamics code. We show that pairs of similar-mass planets in co-orbital resonance are disrupted during large-scale orbital migration. Destabilization occurs when one or both planets is near the critical mass needed to open a gap in the gaseous disk. A confined gap is opened that spans the 60 degree azimuthal separation between planets. This alters the torques imparted by the disk on each planet -- pushing the leading planet outward and the trailing planet inward -- and disrupts the resonance. The mechanism applies to systems in which the two planets masses differ by a factor of two or less. In a simple flared disk model the critical mass for gap opening varies from a few Earth masses at the inner edge of the disk to 1 Saturn-mass at 5 AU. A pair of co-orbital planets with masses in this range that migrates will enter a region where the planets are at the gap-opening limit. At that point the resonance is disrupted. We therefore predict an absence of planets on co-orbital configurations with masses in the super-Earth to Saturn mass range with similar masses.
Earth-mass bodies are expected to undergo Type I migration directed either inward or outward depending on the thermodynamical state of the protoplanetary disc. Zones of convergent migration exist where the Type I torque cancels out. We study the evol
ution of multiple protoplanets of a few Earth masses embedded in a non-isothermal protoplanetary disc. The protoplanets are located in the vicinity of a convergence zone located at the transition between two different opacity regimes. Inside the convergence zone, Type I migration is directed outward and outside the zone migration is directed inward. We used a grid-based hydrodynamical code that includes radiative effects. We performed simulations varying the initial number of embryos and tested the effect of including stochastic forces to mimic the effects resulting from turbulence. We also performed N-body runs calibrated on hydrodynamical calculations to follow the evolution on Myr timescales. For a small number of initial embryos (N = 5-7) and in the absence of stochastic forcing, the population of protoplanets migrates convergently toward the zero-torque radius and forms a stable resonant chain that protects embryos from close encounters. In systems with a larger initial number of embryos, or in which stochastic forces were included, these resonant configurations are disrupted. This in turn leads to the growth of larger cores via a phase of giant impacts, after which the system settles to a new stable resonant configuration. Giant planets cores with masses of 10 Earth masses formed in about half of the simulations with initial protoplanet masses of m_p = 3 Earth masses but in only 15% of simulations with m_p = 1 Earth mass. This suggests that if ~2-3 Earth mass protoplanets can form in less than ~1 Myr, convergent migration and giant collisions can grow giant planet cores at Type I migration convergence zones.
Earth-mass planets embedded in gaseous protoplanetary disks undergo Type I orbital migration. In radiative disks an additional component of the corotation torque scaling with the entropy gradient across the horseshoe region can counteract the general
inward migration, Type I migration can then be directed either inward or outward depending on the local disk properties. Thus, special locations exist in the disk toward which planets migrate in a convergent way. Here we present N-body simulations of the convergent migration of systems of low-mass (M=1-10 m_earth) planets. We show that planets do not actually converge in convergence zones. Rather, they become trapped in chains of mean motion resonances (MMRs). This causes the planets eccentricities to increase to high enough values to affect the structure of the horseshoe region and weaken the positive corotation torque. The zero-torque equilibrium point of the resonant chain of planets is determined by the sum of the attenuated corotation torques and unattenuated differential Lindblad torques acting on each planet. The effective convergence zone is shifted inward. Systems with several planets can experience stochastic migration as a whole due to continuous perturbations from planets entering and leaving resonances.
