No Arabic abstract
We examine how the late divergent migration of Jupiter and Saturn may have perturbed the terrestrial planets. We identify six secular resonances between the nu_5 apsidal eigenfrequency of Jupiter and Saturn and the four eigenfrequencies of the terrestrial planets (g_{1-4}). We derive analytic upper limits on the eccentricity and orbital migration timescale of Jupiter and Saturn when these resonances were encountered to avoid perturbing the eccentricities of the terrestrial planets to values larger than the observed ones. If Jupiter and Saturn migrated with eccentricities comparable to their present day values, smooth migration with exponential timescales characteristic of planetesimal-driven migration (tau~5-10 Myr) would have perturbed the eccentricities of the terrestrial planets to values greatly exceeding the observed ones. This excitation may be mitigated if the eccentricity of Jupiter was small during the migration epoch, migration was very rapid (e.g. tau<~ 0.5 Myr perhaps via planet-planet scattering or instability-driven migration) or the observed small eccentricity amplitudes of the j=2,3 terrestrial modes result from low probability cancellation of several large amplitude contributions. Further, results of orbital integrations show that very short migration timescales (tau<0.5 Myr), characteristic of instability-driven migration, may also perturb the terrestrial planets eccentricities by amounts comparable to their observed values. We discuss the implications of these constraints for the relative timing of terrestrial planet formation, giant planet migration, and the origin of the so-called Late Heavy Bombardment of the Moon 3.9+/-0.1 Ga ago. We suggest that the simplest way to satisfy these dynamical constraints may be for the bulk of any giant planet migration to be complete in the first 30-100 Myr of solar system history.
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.
Many exoplanets in close-in orbits are observed to have relatively high eccentricities and large stellar obliquities. We explore the possibility that these result from planet-planet scattering by studying the dynamical outcomes from a large number of orbit integrations in systems with two and three gas-giant planets in close-in orbits (0.05 AU < a < 0.15 AU). We find that at these orbital separations, unstable systems starting with low eccentricities and mutual inclinations ($elesssim0.1$, $ilesssim0.1$) generally lead to planet-planet collisions in which the collision product is a planet on a low-eccentricity, low-inclination orbit. This result is inconsistent with the observations. We conclude that eccentricity and inclination excitation from planet-planet scattering must precede migration of planets into short-period orbits. This result constrains theories of planet migration: the semi-major axis must shrink by 1-2 orders of magnitude without damping the eccentricity and inclination.
The internal structure of gas giant planets may be more complex than the commonly assumed core-envelope structure with an adiabatic temperature profile. Different primordial internal structures as well as various physical processes can lead to non-homogenous compositional distributions. A non-homogenous internal structure has a significant impact on the thermal evolution and final structure of the planets. In this paper, we present alternative structure and evolution models for Jupiter and Saturn allowing for non-adiabatic primordial structures and the mixing of heavy elements by convection as these planets evolve. We present the evolution of the planets accounting for various initial composition gradients, and in the case of Saturn, include the formation of a helium-rich region as a result of helium rain. We investigate the stability of regions with composition gradients against convection, and find that the helium shell in Saturn remains stable and does not mix with the rest of the envelope. In other cases, convection mixes the planetary interior despite the existence of compositional gradients, leading to the enrichment of the envelope with heavy elements. We show that non-adiabatic structures (and cooling histories) for both Jupiter and Saturn are feasible. The interior temperatures in that case are much higher that for standard adiabatic models. We conclude that the internal structure is directly linked to the formation and evolution history of the planet. These alternative internal structures of Jupiter and Saturn should be considered when interpreting the upcoming Juno and Cassini data.
We study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar-nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets orbits, Delta a_SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inward or outward). We also provide an assessment of the planets orbital dynamics at different epochs of Saturns growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Delta a_SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the timescales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.
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 isothermal 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.