Constructing the secular architecture of the solar system I: The giant planets

Constructing the secular architecture of the solar system I: The giant planets

Key Words.:
Solar System: formation

Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations. Moreover, the amplitudes of the eigenmodes characterising the current secular evolution of the eccentricities of Jupiter and Saturn would not be reproduced correctly; only one eigenmode is excited by resonance-crossing. We show that, at the very least, encounters between Saturn and one of the ice giants (Uranus or Neptune) need to have occurred, in order to reproduce the current secular properties of the giant planets, in particular the amplitude of the two strongest eigenmodes in the eccentricities of Jupiter and Saturn.

1 Introduction

The formation and evolution of the solar system is a longstanding open problem. Of particular importance is the issue of the origin of the orbital eccentricities of the giant planets. Even though these are small compared to those of most extra-solar planets discovered so far, they are nevertheless large compared to what is expected from formation and evolution models.

Giant planets are expected to be born on quasi-circular orbits because low relative velocities with respect to the planetesimals in the disk are a necessary condition to allow the rapid formation of their cores (Kokubo & Ida, 1996, 1998; Goldreich et al., 2004). Once the giant planets have formed, their eccentricities evolve under the effects of their interactions with the disc of gas. These interactions can in principle enhance the eccentricities of very massive planets (Goldreich & Sari, 2003), but for moderate-mass planets they have a damping effect. In fact, numerical hydro-dynamical simulations (Kley & Dirksen, 2006; D’Angelo et al., 2006) show that only planets of masses larger than 2–3 Jupiter masses that are initially on circular orbits are able to excite an eccentricity in the disk and, in response, to become eccentric themselves. Planets of Jupiter-mass or less have their eccentricities damped. Accounting for turbulence should not change the result significantly: the eccentricity excitation due to turbulence is only of the order of 0.01 for a 10 Earth mass planet and decreases rapidly with increasing mass of said planet (Nelson, 2005). By comparison, the mean eccentricities of Jupiter and Saturn are 0.045 and 0.05 respectively.

In addition, the interactions between Jupiter and Saturn, as they evolve and migrate in the disk of gas, should not lead to a significant enhancement of their eccentricities. Figure 1 shows a typical evolution of the Jupiter-Saturn pair, from Masset & Snellgrove (2001). The top panel shows the evolution of the semi major axes, where Saturn’ semi-major axis is depicted by the upper curve and that of Jupiter is the lower trajectory. Initially far away, Saturn swiftly approaches Jupiter, possibly passing across their mutual 2:1 resonance (at approximately 9 000 yr in the figure), and is eventually trapped in the 3:2 resonance. At this point, the migration of both planets slows down slightly and then reverses. Morbidelli & Crida (2007) argued that this dynamical evolution explains why Jupiter did not migrate all the way to the Sun in our System. Pierens & Nelson (2008) convincingly demonstrated that the trapping in the 3:2 resonance is the only possible outcome for the Jupiter-Saturn pair. The lower panel of Figure 1 shows the evolution of the eccentricities of both planets, where Saturn’s eccentricity is depicted by crosses and that of Jupiter by bullets. Both eccentricities remain low all the time. The burst of the eccentricities associated to the passage through the 2:1 resonance at approximately 9 000 yr is rapidly damped. Once trapped in the 3:2 resonance, the equilibrium eccentricities are approximately 0.003 for Jupiter and 0.01 for Saturn, i.e. five to ten times smaller than their current values.

Figure 1: The evolution of Jupiter (bullets) and Saturn (crosses) in the gas disk. Taken from Morbidelli & Crida (2007), but reproducing the evolution shown in Masset & Snellgrove (2001).
Figure 2: The evolution of the eccentricities of Jupiter (bullet) and Saturn (crosses) as they migrate through a 50 Earth masses planetesimal disk. From Gomes et al. (2004).

Once the gas has dispersed from the system, the giant planets are still expected to migrate, due to their interaction with a planetesimal disk (Fernandez & Ip, 1984; Malhotra, 1993, 1995; Hahn & Malhotra, 1999; Gomes et al., 2004). While migration in the gas disk causes the planets to approach each other (Morbidelli et al., 2007), migration in the planetesimal disk causes the planets to diverge i.e. it increases the ratio between the orbital periods (Fernandez & Ip, 1984). In this process, the orbital eccentricities are damped by a mechanism known as “dynamical friction” (e.g. Stewart & Wetherill, 1988). Figure 2 provides a example of the eccentricity evolution of Jupiter (bullets) and Saturn (crosses) who are initially at 5.4 and 8.7 AU with their current eccentricities (initial conditions typical of Malhotra, 1993, 1995; Hahn & Malhotra 1999). They migrate, together with Uranus and Neptune, through a planetesimal disk carrying in total 50. The disk is simulated using 10,000 tracers (see Gomes et al., 2004, for details). The figure shows the evolution of the eccentricities of Jupiter and Saturn: both are rapidly damped below 0.01. Thus, a smooth radial migration through the planetesimal disk, as originally envisioned by Malhotra (1995) cannot explain the current eccentricities (nor the inclinations) of the orbits of the giant planets.

Then, how did Jupiter and Saturn acquire their current eccentricities? In Tsiganis et al. (2005), the foundation paper for a comprehensive model of the evolution of the outer Solar System – often called the Nice model – it is argued that the current eccentricities were achieved when Jupiter and Saturn passed across their mutual 2:1 resonance, while migrating in divergent directions under the interactions with a planetesimal disk. They indeed showed that the mean eccentricities of Jupiter and Saturn are adequately reproduced during the resonance crossing (see electronic supplement of Tsiganis et al., 2005, or figure 5 below), as well as their orbital separations and mutual inclinations. However, the mean values of the eccentricities do not properly describe the secular dynamical architecture of a planetary system: the eccentricities of the planets oscillate with long periods, because of the mutual secular interactions among the planets. A system of planets has fundamental frequencies in the secular evolution of the eccentricities, and the amplitude of each mode – or, at least that of the dominant ones – should be reproduced in a successful model. We remind that Tsiganis et al. (2005) never checked if the Nice model reproduces the secular architecture of the giant planets (we will show below that it does) nor if this is achieved via the 2:1 resonance crossing (we will show here that it is not).

Actually, in this paper we make an abstraction of the Nice model, and investigate which events in the evolution of the giant planets are needed to achieve the current secular architecture of the giant planet system. We start in section 2 by reviewing what this secular architecture is and how it evolves during migration, in the case where no mean motion resonances are crossed. In section 3 we investigate the effect of the passage through the 2:1 resonance on the secular architecture of the Jupiter-Saturn pair. As we will see, this resonance crossing alone, although reproducing the mean eccentricities of both planets, does not reproduce the frequency decomposition of the secular system. In section 4 we discuss the effect of multiple mean motion resonance crossings between Jupiter and Saturn, showing that this is still not enough to achieve the good secular solution. In section 5 we examine the role of a third planet, with a mass comparable to that of Uranus or Neptune. We first consider the migration of this third planet on a circular orbit, then on an eccentric orbit and finally we discuss the consequences of encounters between this planet and Saturn. We show that encounters of Saturn with the ice giant lead to the correct secular evolution for the eccentricities of Jupiter and Saturn. In section 6 we return to the Nice model, verify its ability to reproduce the current secular architecture of the planetary system and discuss other models that could in principle be equally successful in this respect. Although this paper is mostly focused on the Jupiter-Saturn pair and the evolution of their eccentricities, in section 7 we briefly discuss the fate of Uranus and Neptune and the excitation of inclinations. The case of the terrestrial planets will be discussed in a second paper. The results are then summarised in section 8.

2 Secular eccentricity evolution of the Jupiter-Saturn pair

One can study the secular dynamics of a pair of planets as described in Michtchenko & Malhotra (2004). In that case the planets are assumed to evolve on the same plane and are far from mutual mean motion resonances. The Hamiltonian describing their interaction is averaged over the mean longitudes of the planets. This averaged Hamiltonian describes a two-degrees-of-freedom system, whose angles are the longitudes of perihelia of the two planets: and . The D’Alembert rules (see Chapter 1 of Morbidelli, 2002), ensure that the Hamiltonian depends only on the combination . Thus the system is effectively reduced to one degree of freedom, which is integrable. This means that, in addition to the value of the averaged Hamiltonian itself, which will improperly be called “energy” hereafter, the system must have a second constant of motion. Simple algebra on canonical transformations of variables allows one to prove that this constant is


where , , is the mass of the star, is the gravitational constant and and are the mass, semi major axis an eccentricity of each planet. Note that , called angular momentum deficit, actually measures the deviation in angular momentum for an eccentric two-planet system, with respect to a system of two planets, both on circular orbits, with the same values of . Now the global secular dynamics of the system can be illustrated by plotting level curves of the energy over manifolds defined by the condition that is constant.

Figure 3: Global illustration of the secular dynamics of the Jupiter-Saturn system. The bullets represent the current values of , and .

Figure 3 shows the result for the Jupiter-Saturn system. The value of that we have chosen corresponds to the current masses, semi major axes and eccentricities of these planets. The left panel illustrates the dynamics in the coordinates and , while the right panel uses the coordinates and , where refers to the eccentricity of Jupiter and to the eccentricity of Saturn. The bullets represent the current configuration of the Jupiter-Saturn system. We stress that the two panels are just two representations of the same dynamics. The same level curves of the energy are plotted in both panels. Thus, the level curve counting from the triangle in the left panel corresponds to the level curve counting from the triangle in the right panel. Indeed, the dot representing the current Jupiter-Saturn configuration is on the 5 level curve away from the triangle on each panel. The secular evolution of the system has to follow the energy level curve that passes through the dot. The other energy curves show the secular evolution that Jupiter and Saturn would have had, if the system were modified relative to the current configuration, preserving the current value of . We warn the reader that the dynamics illustrated in this figure is not very accurate from a quantitative point of view, because we have neglected the effects of the nearby 5/2 mean motion resonance between Jupiter and Saturn. Nevertheless all the qualitative aspects of the real dynamics are correctly reproduced.

We remark that the global secular dynamics of the Jupiter-Saturn system is characterised by the presence of two stable equilibrium points, one at (marked by a triangle in figure 3) and one at (marked by a cross). Thus, there are three kinds of energy level curves along which the Jupiter-Saturn system could evolve: those along which librates around (type I), those along which circulates (i.e. assumes all values from 0 to ; type II) and those along which librates around (type III). Notice that while type II curves wrap around the stable equilibrium at in the left panel, the curves wrap around the stable equilibrium at in the right panel. This means that during the circulation of , the eccentricity of Jupiter has a maximum when while that of Saturn has a maximum when . The real Jupiter-Saturn system has this type of evolution.

We stress that there is no critical curve (separatrix) separating the evolutions of type I, II and III. By critical curve we mean a trajectory passing though (at least) one unstable equilibrium point, along which the travel time is infinite; an example is the curve separating the libration and circulation regimes in a pendulum. In this respect, speaking of “resonance” when librates, as it is sometimes done when discussing the secular dynamics of extra-solar planets, is misleading because the word “resonance”, in the classical dynamical systems and celestial mechanics terminology, implies the existence of such a critical curve.

Frequency Value (/yr) Phase ()
4.26 30.67
28.22 128.11
Table 1: Frequencies and phases for the secular evolution of Jupiter and Saturn on their current orbits.
5 6
5 0.0442 0.0157
6 0.0330 0.0482
Table 2: Coefficients of the Lagrange–Laplace solution for the Jupiter-Saturn system. The coefficients of the terms with frequencies other than and are omitted.

In addition to using phase portraits, the secular dynamics of the Jupiter-Saturn system, or any pair of planets with small eccentricities, can also be described using the classical Lagrange-Laplace theory (see Chapter 7 in Murray & Dermott, 1999). This theory, which is in fact the solution of the averaged problem described above, in the linear approximation, states that the eccentricities and longitudes of perihelia of the pair of planets evolve as:


where and . Here and are the eigenfrequencies of the system, while and are their phases at . In equation (5) all . Tables 1 and 2 report the values of all the coefficients, obtained from the Fourier analysis of the complete 8-planet numerical solution (Nobili et al., 1989).

There is a one-to-one correspondence between the relative amplitudes of the coefficients and the three types of secular evolution illustrated in figure 3. We detail this relationship below, in order to achieve a better understanding of the planetary evolutions illustrated in the next sections.

From equation (5) the evolution of and (the quantities plotted on the -axes of the panels in figure 3) are:






When one has


where sign() is equal to if the argument of the function is negative, if it is positive and 0 if it is zero. Instead, when one has


Now, suppose that the amplitude corresponding to the frequency is zero, i.e. . Then, the dependence of equation (9) on vanishes and, from equation (13) or equation (15), one sees that the system is located at the equilibrium point at , so that and (the point marked by a cross in figure 3).

Let us now gradually increase the amplitudes of the mode, relative to that of the one. This implies increasing and at the same rate, while keeping and fixed. Initially, when and are small compared to and , all quantities in equations (13) and (15) are negative, and therefore the evolution of the system follows an energy level curve of type I, along which librates around . The distance of this curve from the equilibrium point, which we call amplitude of oscillation hereafter, is directly proportional to or .

Since and , then when one has . This implies that increasing the amplitude of the mode eventually brings us to the situation where becomes larger than , but is still less than . Now the value of at i.e.  becomes positive. When additionally its value remains negative, i.e. . Thus, the system now evolves on an energy curve of type II, along which circulates. Notice that the value of at remains negative, while the value at jumps from to . Thus, the level curve in the right panel of figure 3 flips from one looping around the equilibrium point at to one looping around the equilibrium at .

Further increasing the amplitude of the mode relative to that of the mode, eventually results in also becoming larger than . Now, all quantities in equations (13) and (15) are positive, which means the the system follows an energy level curve of type III, along which librates around 0. Notice that the value of at jumps from to , which means that the the level curve in the left panel of figure 3 flips from one going around the equilibrium point at to one going around the equilibrium at .

Finally, when the amplitude of the mode is zero, the system is on the stable equilibrium at (the cross in figure 3).

Below we discuss how the secular dynamics of the planets changes as they migrate away from each other and are also submitted to dynamical friction, exerted by the planetesimal population.

2.1 Migration and the evolution of the secular dynamics

Let’s imagine two planets migrating, without passing through any major mean motion resonance. A good example could be Jupiter and Saturn migrating from a configuration with orbital period ratio slightly larger than 2 to their current configuration, with slightly smaller than 2.5.

The migration causes the semi-major axis ratio between the planets to change. This affects the values of the coefficients , since they depend explicitly on the above ratio; in turn, this affects the global portrait of secular dynamics. However, if the migration is slow enough, the amplitude of oscillation around the equilibrium point is preserved as an adiabatic invariant (Neishtadt, 1984; Henrard, 1993). More precisely, it can be demonstrated that the conserved quantity is


which is the action conjugate to , and where denotes the integral over a closed energy curve (i.e. a bounded trajectory) that characterises the secular motion of the two planets (as in figure 3), if migration is frozen. Therefore during migration the planets would react to the slow changes in the global dynamical portrait, by passing from one energy curve to another in such a way as to preserve the quantity in equation (16), i.e. the oscillation amplitude around the stable equilibrium point remains constant. Since for Jupiter and Saturn this amplitude is related to , it turns out that any smooth migration should not have changed this coefficient significantly. In the next section the additional effect of dynamical friction is analysed.

2.2 Dynamical friction and the evolution of the secular dynamics

Dynamical friction is the mechanism by which gravitating objects of different masses exchange energy so as to evolve towards an equipartition of energy of relative motion (Saslaw, 1985). For a system of planets embedded in a massive population of small bodies, the eccentricities and inclinations of the former are damped, while those of the latter are excited (Stewart & Wetherill, 1988).

In principle, each planet might suffer dynamical friction with different intensity, owing to a different location inside the planetesimal disk. However, the planets are connected to each other through their secular dynamics, so that dynamical friction, even if acting in an unbalanced way between the planets, turns out to have a systematic net effect.

Figure 4: Effect of eccentricity damping on the evolution of Jupiter and Saturn. In this experiment the damping force is applied to Saturn only. The top panel shows the coefficients (filled circles), (open circles), (filled triangles) and (open triangles) as a function of time. The bottom panel shows the eccentricities of Saturn (open circles) and Jupiter (filled circles).

To illustrate this point, consider again the Jupiter-Saturn system of figure 3 and suppose that dynamical friction is applied only to Saturn. The eccentricity of Saturn is damped, so the value of is reduced. Consequently, the location of the two stable equilibrium points has to move towards . If the adiabatic invariance of equation (16) held, the amplitude of oscillation around the equilibrium point would be preserved, eventually turning a libration of into a circulation. However, the adiabatic invariance does not hold in this case. The reason is that dynamical friction damps the eccentricity of Saturn, and therefore damps both the and the coefficients. Since the and coefficients are related, is also damped. In other words, the amplitude of oscillation around the equilibrium point is damped, and so the value of given in equation (16) decays with time.

As a check, we have run a simple numerical experiment. We have considered a Jupiter-Saturn system with semi major axes 5.4 and 8.85, with relatively eccentric orbits and large amplitude () of apsidal libration around . We have integrated the orbits using the Wisdom-Holman (Wisdom & Holman, 1991) method, with the code Swift-WHM (Levison & Duncan,1994). We used a time-step of 0.1 y and modified the equations of motion so that a damping term is included for the eccentricity of Saturn only. Figure 4 shows the result. The bottom panel shows the evolutions of the eccentricities of the two planets, where is represented by open circles and by filled circles: both are damped and decay with time at the same rate. The top panel shows the amplitudes of the coefficients of equation (5), i.e.  (filled circles), (open circles), (filled triangles) and (open triangles), computed at six different points in time. The computations were performed, by applying Fourier analysis to the time series produced in respectively six short-time integrations, where no eccentricity damping was applied. As one can see, all coefficients decrease with time, at comparable rates.

Thus, we conclude that, whatever the initial configuration of the planets, smooth migration and dynamical friction cannot increase the amplitude of the term and cannot turn the libration of into a circulation. This result will be relevant in the next section.

3 The effect of the passage through the 2:1 resonance

We now consider the migration of Jupiter and Saturn, initially on quasi-circular orbits, through their mutual 2:1 mean motion resonance. Tsiganis et al. (2005) argued that this passage through the resonance is responsible for the acquisition of the current eccentricities of the two planets.

Figure 5 shows the effect of the passage through this resonance, starting from circular orbits. The simulation is again done using the Swift-WHM integrator, but in this case the equations of motion are modified so as to induce radial migration to the planets, with a rate decaying as . No eccentricity damping is imposed. In practise, at every timestep the velocity of each planet is multiplied by a quantity , with being proportional to . For Jupiter is negative and for Saturn it is positive, so that the two planets migrate inwards and outwards respectively, as observed in realistic -body simulations (Fernandez & Ip, 1984; Hahn & Malhotra, 1999; Gomes et al., 2004). We choose Myr.

Figure 5: The evolution of Jupiter and Saturn, as they pass across their mutual 2:1 resonance. In the bottom panel, Saturn’s eccentricity is the upper curve and Jupiter’s is the lower one.
Figure 6: Similar to figure 3 but after the 2:1 resonance crossing.

The top panel of Figure 5 shows the ratio of the orbital periods of Saturn () and Jupiter () as a function of time. We stop the simulation well before achieves the current value, to emphasise the effect of the 2:1 resonance crossing. The middle panel shows as a function of time. The bottom panel shows the evolutions of the eccentricities of Jupiter and Saturn, where the lower trajectory corresponds to Jupiter and the upper curve corresponds to Saturn. We notice that the orbital period ratio abruptly jumps across the value of 2. Correspondingly, the eccentricities of Saturn and Jupiter jump to and , which, as noticed by Tsiganis et al. (2005), are quite close to the current mean eccentricities of the two planets. During the subsequent migration, the eccentricity of Saturn increases somewhat and that of Jupiter decreases respectively. The two planets enter into apsidal anti-alignment (i.e.  librates around 180 degrees) shortly before the resonance passage and the libration amplitude shrinks down to as the eccentricities o the two planets grow. The crossing of the 2:1 does not seem to significantly affect the libration amplitude, which remains of the order of during the post-resonance-crossing migration (see also Čuk, 2007). As a result of this narrow libration amplitude, the eccentricities of Jupiter and Saturn do not show any sign of secular oscillation. In practise, Jupiter and Saturn are located at the stable equilibrium point of their secular dynamics, as shown in figure 6. This is very different from the current situation (compare with figure 3). Table 3 reports the values of the coefficients of equation (5) at the end of the simulation. The and coefficients, related to the amplitude of oscillation around the equilibrium point as explained in section 2, are very small; they are more than an order of magnitude smaller than their current values. As discussed in section 2.1, they would not increase during the subsequent migration of the planets, because they behave as adiabatic invariants. Even worse, they would decrease if dynamical friction were applied.

5 6
5 0.00272 0.00275
6 0.0378 0.0854
Table 3: Coefficients of the Jupiter-Saturn secular system at the end of the simulation illustrated in figure 5.
1 3.55
2 4.41
5 7.49
10 7.02
20 5.05
Table 4: Values of in Jupiter after the 2:1 resonance crossing with Saturn as a function of the migration e-folding time, .

Simulations that we performed assuming larger values of (i.e. slower migration rates) lead to the same result. The eccentricities of Jupiter and Saturn after the 2:1 resonance crossing are about the same as in figure 5. The amplitude of libration of is always between 6 and 15 degrees, with no apparent correlation on . Table 4 recapitulates the results, for what concerns the values of the coefficient, which are always comparably small. Thus, we conclude that the 2:1 resonance crossing, although it explains the current mean values of the eccentricities of Jupiter and Saturn, cannot by itself explain the current secular dynamical structure of the system.

Let us now provide an interpretation of the behaviour observed in the simulation. The dynamics of two planets in the vicinity of a mean motion resonance can be studied following Michtchenko et al. (2008). For the 2:1 resonance, the fundamental angles of the problem are


The motion of the first angle is conjugated with the motion of the angular momentum deficit of the planets, defined in equation (1). The motion of the second angle is conjugated with the motion of the quantity . If the system is far from the resonance, the motion of can be averaged out. Then becomes a constant of motion and the secular dynamics described in the previous section is recovered. In particular, if the planets migrate slowly enough so that the adiabatic invariance of – see equation (16) – holds, the amplitude of oscillation around the stable equilibrium point of the secular dynamics has to remain constant. However, as migration continues, the approach to the mean motion resonance forces the location of the equilibrium point to shift to larger eccentricity. Thus, by virtue of a geometric effect, the apparent amplitude of libration of has to shrink as the eccentricities increase. This is visible in Fig. 5 in the phase before the resonance crossing.

As the planets are approaching the resonance, the angle can no longer be averaged out in a trivial way. However, an adiabatic invariant can still be introduced, as long as the timescale for the motion of is significantly shorter than that of and that migration changes the system on even longer time scales. This invariant is


where the integral is taken over a path describing the coupled evolution of and , which is closed if and are frozen (i.e. a trajectory of the so-called “frozen” system).

As shown in Michtchenko et al. (2008), for a pair of planets with the Jupiter-Saturn mass ratio and , the dynamics of the 2:1 resonance presents one critical curve, or separatrix, for the degree of freedom and no critical curve for the degree of freedom. Thus, when the resonance is reached during the migration, the invariance of is broken (Neishtadt, 1984), since the two time-scales of the motion are no more well separated. As the planets are migrating in divergent directions, they cannot be trapped in the resonance (Henrard & Lemaitre, 1983). The resonant angle has to switch from clockwise to anti-clockwise circulation. Correspondingly, the quantity has to jump to a larger value, that is the planets acquire an angular momentum deficit. How this new angular momentum deficit is partitioned between the two planets is difficult to compute a priori, because the dynamics are fully four-dimensional (i.e. two degrees of freedom) and, therefore, not integrable.

Figure 7: Left: the evolution of and in polar coordinates. Right: the evolution of and . Before the resonance is reached, both planets have nearly zero eccentricities (region I, black dots). When the resonance is reached, both planets jump to the corresponding Region III, along the -axis and through the “X”-point of the critical curve (grey dots). From the supplementary material of Tsiganis et al. (2005)

As shown in figure 7, there are clearly two regimes of motion in the portraits and (remember that and ). Before the resonance crossing (black dots), the dynamics are confined in a narrow, off-centred region close to the origin of the polar coordinates. After the resonance crossing (grey dots), the dynamics fill a wide annulus, also asymmetric relative to the axis . The curves in each panel are free hand illustrations of the dynamics near and inside a first-order mean motion resonance. The region filled by the black dots is bounded by the inner loop of the critical curve (labelled I in the plot). The annulus filled by the grey dots is adhesive, at its inner edge, to the outer loop of the critical curve. The jump in eccentricity observed at the resonance crossing corresponds to the passage from the region inside the inner loop to that outside the outer loop.

Thus, in practise, it is as if each planet saw its own resonance: the one with critical angle for Jupiter and the one with critical angle for Saturn. The two resonances are just two slices of the same resonance, because only one critical curve exists (Michtchenko et al., 2008). This is the reason why the eccentricities of both planets jump simultaneously. The 2:1 resonance is structured by the presence of a periodic orbit, along which and remain constant and are equal to 0 and respectively (Michtchenko et al., 2008). As , the phase portrait of the resonance and that of the resonance are rotated by 180 degrees, with respect to each other. Thus, reaches a maximum when and when . Consequently, when the planets reach their maximal eccentricities and the angles start to circulate anti-clockwise, has to be (see figure 7). It is evident that the result of this transition through the resonance depends just on the resonance topology and not on the migration rate, as long as the latter is slow compared to the motion of the angles (i.e. the adiabatic approximation holds; Neishtadt, 1984).

After the resonance crossing, one can again average over and reduce the system to a one-degree of freedom secular system. As , Jupiter has to be on the negative -axis of a diagram like that of the left panel of figure 6. Thus the secular evolution of Jupiter will be an oscillation around the stable equilibrium at . The amplitude of this oscillation depends on the value of the eccentricity of Jupiter acquired at resonance crossing, relative to the value of the stable equilibrium of the secular problem. It turns out that, for the masses of Jupiter and Saturn these two values are almost the same. Thus, the amplitude of oscillation is very small. We think that this is a coincidence and that, in principle, it does not have to be that way. In fact, we have verified numerically that the result depends on the individual masses of both planets, even for the same mass ratio. For instance, if the masses of Jupiter and Saturn are both reduced by a factor 100, the eccentricities of both planets jump to at resonance crossing, and this puts Jupiter on a secular trajectory that brings down to 0; that is, the secular motion is now at the boundary between type-I and type-II, as defined in section 2. This is caused by a different scaling of the jumps in and with respect to the planetary masses and to a different global secular dynamics in the vicinity of the resonance, which in turn is caused by a different relative importance of the quadratic terms in the masses.

Once the planets are placed relative to the portrait of their secular dynamics, their destiny is fixed. As they move away from the mean motion resonance, the secular portrait can change, in particular because the near-resonant perturbation terms that are quadratic in the masses rapidly decrease in amplitude. Hence the location of the equilibrium points can change, but the planets have to follow them adiabatically. This explains the slow monotonic growth of the eccentricity of Saturn and the decay of that of Jupiter, observed in the top panel of figure 5, while the amplitude of libration of does not change (middle panel).

4 Passage through multiple resonances

Since the passage of Jupiter and Saturn through the 2:1 resonance, starting from initially circular orbits, produces a secular system that is incompatible with the current one, we now explore the effects of the passage of these planets through a series of resonances. This is done to determine whether or not such evolution could increase the value of in both planets.

We set Jupiter and Saturn initially on quasi-circular orbits just outside their mutual 3:2 resonance. The choice of these initial conditions is motivated by the result that during the gas-disk phase, Saturn should have been trapped in the 3:2 resonance with Jupiter (Morbidelli et al., 2007; Pierens & Nelson, 2008). Once the gas disappeared from the system, the two planets should have been extracted from the resonance at low eccentricity, by the interaction with the planetesimals, and subsequently start to migrate.

In the above setting, our planets are forced to migrate through the 5:3, 7:4, 2:1, 9:4 and 7:3 resonances, ending up close to their current location in semi major axis (i.e. slightly interior to the 5/2 resonance). In the migration equations we set so that it takes about 40 Myr to reach although, as we saw before, the migration timescale has little influence on the resonant effects. The result of this experiment is shown in figure 8, which is similar in format to figure 5. The horizontal lines in the top panel denote the positions of the resonances mentioned above. Notice a distinct jump in the eccentricities of both planets at each resonance crossing. In order to prevent the system from becoming unstable, we applied eccentricity damping to Saturn, so to mimic the effect of dynamical friction and reach final eccentricities that are similar to the current mean values of the two planets. The parameters for the simulation depicted in figure 8 are Myr and yr. The effect of damping is visible in the eccentricity evolution, after the 2:1 resonance crossing; we will discuss the effect on the motion of below.

Figure 8: Like figure 5, but for Jupiter and Saturn evolving through a sequence of mean motion resonances, from just outside their mutual 3:2 commensurability up to their current location.

The middle panel of figure 8 shows that the passage through the 7:4 resonance significantly increases the libration amplitude of . In fact, the amplitude of the term in (5), increases from before the resonance crossing, to after the crossing. The passage through the 2:1 resonance, however, shrinks the amplitude of libration of . This happens because the 2:1 resonance crossing, as we have seen in the previous section, does not enhance (it remains equal to in this simulation) but does enhance the overall eccentricities of the planets. As explained earlier, this causes the amplitude of libration of to decrease.

Notice from figure 8 that, after the 2:1 resonance crossing, the amplitude of libration of starts to increase, slowly and monotonically. This is caused not by an enhancement of the amplitude of oscillation around the equilibrium point of the secular dynamics, but by the damping of the eccentricities of the planets. It is the opposite of what was just described before: a geometrical effect. In reality, the value of is decreased to (from 0.019) during this evolution. Hence, at the end of the simulation, the amplitude of the mode is about a factor of 3 smaller than in the real secular dynamics of Jupiter and Saturn. Without eccentricity damping, the amplitude of the mode would have remained equal to , still much smaller than in the current Jupiter-Saturn secular dynamics.

Several other experiments, changing the initial conditions slightly or the migration speed , lead essentially to the same result. Thus, we conclude that the migration of Jupiter and Saturn through a sequence of mean motion resonances is not enough to achieve their current secular configuration. A richer dynamics is required, likely involving interactions with a third planet.

5 Three-planet dynamics

From the discussions and the examples reported above, it is quite clear that, to enhance the amplitudes of the mode, it is necessary that the eccentricity of Saturn receives a kick that is not counterbalanced by a corresponding increase in the eccentricity of Jupiter (or vice versa). This would indeed move the planets away from the stable equilibrium point of their secular dynamics, thus enhancing the amplitude of oscillation around this point and, consequently, . Given that Jupiter and Saturn are not alone in the outer solar system, in this section we investigate the effect that interactions with a third planet with a mass comparable to that of Uranus and Neptune, which we simply refer to as ’Uranus’, has on the Jupiter-Saturn pair. We first address the effects of the migration of Uranus on a quasi-circular orbit. Then we study the effects of its migration on an initially eccentric orbit and, finally, we address the problem of encounters among the planets.

5.1 Migration of Uranus on a quasi-circular orbit

The main mean motion resonance with Saturn that Uranus can go through is the 2:1. Thus, this is the resonance crossing that we focus on here. Given that, as we have seen in the previous sections, the effect of a passage through a mean motion resonance is quite insensitive to the migration rate, the initial location of the planets etc., the main issue that may potentially lead to different results is whether the crossing of the 2:1 Saturn-Uranus resonance happened before or after the putative crossing of the 2:1 Jupiter-Saturn resonance. Below we investigate each of these two cases.

To have the crossing of the Saturn-Uranus resonance happen first, we have performed a numerical experiment, with Saturn and Jupiter having initially a small orbital period ratio (), and Uranus and Saturn having an orbital period ratio . The exact initial locations are not important, as long as they do not change the order of the resonance crossings. All planets initially have circular orbits. The three planets are forced to migrate to their current positions with Myr. Eccentricity damping is imposed on Saturn and Uranus, to mimic dynamical friction, with forces tuned such that, at the end of the simulation, Uranus approximately reaches its current eccentricity.

Figure 9 shows the result. The top panel shows the pericentre and apocentre of the planets which, from top to bottom, are Uranus, Saturn and Jupiter respectively. The separation among these curves gives a visual measure of the eccentricity of the orbit of the respective planet. The middle panel shows for Jupiter and Saturn.

Figure 9: A three-planet migration simulation, in which both the 2:1 resonance between Uranus and Saturn and the 2:1 resonance between Saturn and Jupiter are crossed.

In this simulation, Uranus crosses the 2:1 resonance with Saturn at Myr. This gives a kick to the eccentricity of Uranus (its -curves abruptly separate from each other) and, to a lesser extent, to the eccentricity of Saturn. This sudden increase in the eccentricity of Saturn moves the stable equilibrium point of the Jupiter-Saturn secular dynamics away from . However, the eccentricity of Jupiter does not receive an equivalent kick by this resonance crossing, so it remains close to zero. Consequently, Jupiter must start evolving secularly along a trajectory close to the boundary between type-I and type-II curves (see section 2); in other words . This is the reason why the amplitude of libration of changes abruptly at the Uranus-Saturn resonance crossing, and reaches an amplitude of .

The interim between 0.9 and 3.3 Myr is characterised by large, long-periodic, oscillations of the eccentricity of Uranus, which correlate with the modulation of the amplitude of libration of . These oscillations have a frequency equal to , where is the new fundamental frequency that characterises the extension of (5) to a three-planet system. Soon after the Uranus-Saturn resonance crossing, is small and therefore the oscillations have large amplitude. As Uranus departs from the resonance with Saturn, decreases; at the same time, increases, since Jupiter approaches the 2:1 resonance with Saturn. Hence, the oscillation with frequency becomes faster and its amplitude deceases. This sequence of increasingly shorter oscillations reduces the overall amplitude of libration of to approximately 40 degrees.

At  Myr, Jupiter and Saturn cross their mutual 2:1 mean motion resonance, which has the effects that we discussed before. The amplitude of is roughly preserved in this resonance crossing, as already illustrated in section 4. Its value at  Myr is 0.010, much larger than in section 3 but still about a factor of four smaller than the current value. The coefficient receives an additional small enhancement at  Myr, when Jupiter and Saturn cross their 7:3 resonance, but this does not change the substance of the result.

To reverse the order of the resonance crossings, we have run a second experiment, in which we placed Jupiter and Saturn just beyond their 2:1 resonance (), on orbits typical of those achieved by the 2:1 resonance crossing (see section 3). This means that Jupiter and Saturn are in apsidal libration around 180, with an amplitude of the mode that is small in both planets, compared to the current value (see table 3). Uranus was placed on a circular orbit at . Again the planets were forced to migrate to their current locations, with Myr. The 2:1 resonance crossing between Uranus and Saturn again kicked the eccentricity of Saturn, which in turn enhanced the amplitude of the term. In this case, the final value of the coefficient was 0.014, i.e. of the same order as in the previous experiment.

Given the above results, we conclude that, no matter when the Uranus-Saturn 2:1 resonance crossing occurs, there is an enhancement of the amplitude of the term, as expected, but it is too small (by a factor of ) to explain the current Jupiter-Saturn secular system. It appears that the mass of Uranus is too small to provide enough eccentricity excitation on Saturn, when passing through a mean motion resonance with it.

5.2 Migration of Uranus on an eccentric orbit

In the current solar system, the proper frequency of perihelion of Uranus, , is smaller than (3.1 and 4.3 /yr, respectively). If Uranus was much closer to Saturn, however, had to be much larger too. For instance, if Uranus were just outside the 2:1 resonance with Saturn (say  AU and  AU), was /yr. One might then wonder if, during the migration of Uranus the secular resonance could have occurred. On the other hand, if Jupiter was closer to Saturn, would have been higher as well; /yr if . Thus, the occurrence of the resonance depends on the locations of both Jupiter and Uranus, relative to Saturn.

Figure 10: A three-planet simulation, with Uranus initially on an eccentric orbit and migrating to its current location. No migration is imposed on Jupiter and Saturn.

To see the effect of this secular resonance, we performed an idealised experiment, in which we placed Jupiter at  AU, Saturn at  AU and Uranus at  AU with an eccentricity of 0.25. The initial eccentricities and for Jupiter and Saturn were taken from a run, in which these planets passed through their mutual 2:1 resonance and migrated up to the locations reported above. Hence would librate with very small amplitude, in the absence of Uranus. In this experiment the latter was forced to migrate towards its current location, with Myr. No migration was imposed on Jupiter and Saturn. Since initially /yr and , the resonance crossing had to occur during this simulation. No eccentricity damping was applied to any of the planets in this run.

The result is shown in figure 10. In the first part of the simulation ( Myr) the amplitude of libration of , which is initially very small, suffers a large modulation, correlated with the oscillations of the eccentricity of Uranus. The dynamics here are in analogous to what we described before, for the interim between the two mean motion resonance crossings in figure 9. At  Myr, the resonance is crossed. This leads to an exchange of angular momentum between Uranus and Jupiter. The eccentricity of Uranus decreases a bit, while the value of the coefficient is enhanced. As a response, starts to circulate. The final value of is 0.04, essentially matching the current value.

Although Jupiter has to be far from Saturn to have a genuine secular resonance crossing, we have found that more realistic simulations, in which Jupiter and Saturn are initially much closer to each other – so to be able to migrate in the correct proportion with respect to Uranus – can lead to interesting results as well. The reason is that, although from the beginning , the two frequencies can become quasi-resonant; such interactions also allow for a significant transfer of eccentricity from Uranus to Jupiter and can excite the value of up to the current figure.

However, the reader should be aware that, while the effect of mean motion resonances is quite insensitive on parameters and initial conditions, in the case of a secular resonance, the outcome depends critically on a variety of issues. More precisely, the effects of the resonance, or quasi-resonance, must depend on the eccentricity of Uranus, the migration timescale and the position of relative to , immediately before this resonant interaction. The reason for the first dependence is that sets the strength of the secular resonance. The dependence on and has to do with the fact that the timescale associated with a secular resonance is very long, of the order of 1 Myr. Thus, migration through a secular resonance, unlike migration through a mean motion resonance, is not an adiabatic process, at least for values of up to 10 Myr that we are focusing on here. Thus, the time spent in the vicinity of the resonance and the values of the phases at which the planets enter the resonance have important impact on the resulting dynamics.

Given the above, we conclude that the secular interaction with an eccentric Uranus is a mechanism that is potentially capable of exciting the mode in the Jupiter-Saturn system to the observed level, but this mechanism is quite un-generic. Moreover, if we invoked an eccentric Uranus to explain the origin of the Jupiter-Saturn dynamical architecture, we would still need to explain how Uranus got so eccentric in first place. Finally, the secular resonance cannot alone explain the excitation of the planetary inclinations, which will be discussed in section 7. For all these reasons, we continue our search for a better mechanism and consider below the effect of planetary encounters.

5.3 Planetary encounters with Uranus

Close encounters between Uranus and Saturn could potentially be a very effective mechanism for kicking the eccentricity of Saturn and enhancing the amplitude of the mode.

To investigate this, we have run a series of twenty simulations, where the initial semi-major axes of Jupiter and Saturn were chosen such that these two planets are just outside their 2:1 mean motion resonance ( = 2.06), on orbits typical of those achieved during the 2:1 resonance crossing (in apsidal anti-alignment with negligible oscillation amplitude). Uranus was placed on an orbit with semi-major axis ranging from 11.8 AU to 13.4 AU at 0.2 AU intervals, with an initial eccentricity of 0.1. This value of the eccentricity is of the order of that achieved by Uranus, under the secular forcing induced by Jupiter and Saturn. The system was then allowed to evolve under the mutual gravitational forces and external migration forces. For all simulations, the e-folding time for the migration forces was set at 5 Myr. Eccentricity damping was applied to Uranus and Saturn, the values being /yr and /yr. The damping coefficient for Uranus was assumed to be six times larger than that of Saturn because, in principle, Uranus is more affected by the planetesimal disk than the gas-giants. The strength of the damping term was calibrated so that the post-encounter evolution of the eccentricity of Uranus follows the one seen in the full -body simulations of Tsiganis et al. (2005). These details are not very important, because we focus here on the final secular dynamics of Jupiter and not on the final orbit of Uranus. The latter is very sensitive to the prescription of damping, but not the former as we have seen in sect. 2.2. Uranus was typically found to be scattered by Saturn (and sometimes by Jupiter). The simulations were stopped once the phase of encounters among the planets ended, either because Uranus was decoupled from the giant planets, due to eccentricity damping, or because it was ejected from the system.

The simulations that yielded the best results in terms of final value are those with initial semi-major axes of Uranus AU. In these successful simulations, a total of four or 20%, the average final value of was approximately 0.04, in very good agreement with the current configuration. Figure 11 gives an example of evolution from one of these successful runs, and is similar to figure 9 and figure 10. As seen in our results, initially placing Uranus further away from Saturn results in encounters that are too weak to pump up the mode in Jupiter.

Figure 11: Example of an encounter between Uranus and Saturn. The plot is similar to figure 9.

In summary, we conclude that encounters between Saturn and Uranus constitute an effective and quite generic mechanism for achieving a final secular evolution of the Jupiter-Saturn system that is consistent with their current state. Compared to all other mechanisms investigated in this paper, which either do not work or work only for an ad-hoc set of conditions, planet-planet scattering is our favoured solution to the problem of the origin of the secular architecture of the giant planet system. In Section 7 we provide further arguments in favour of this conclusion.

6 The Nice model and its alternatives

The work presented above shows that a combination of the effects provided by the 2:1 resonance crossing between Jupiter and Saturn and by encounters and/or secular interactions with an eccentric Uranus, can produce a Jupiter-Saturn system that behaves secularly like the real planets.

The 2:1 resonance crossing, the encounters among the planets and the high-eccentricity phases of Uranus and Neptune are essential ingredients of the Nice model (Tsiganis et al., 2005; Gomes et al., 2005; Morbidelli et al., 2007). Thus, we expect that this model not only reproduces the mean orbital eccentricities of the planets, as shown in Tsiganis et al. (2005), but also the correct architecture of secular modes. Curiously, this has never been properly checked before, and we do so in the following.

Figure 12: Sample evolution of the giant planets in the Nice model (from Gomes et al., 2005). Each planet is represented by a pair of curves, showing the time evolution of their perihelion and aphelion distance.

In Figure 12, the pericentre and apocentre distance of the four giant planets are plotted as a function of time, in a simulation taken from Gomes et al. (2005) that adequately reproduces the current positions of the giant planets. The curve starting around 5 AU represents Jupiter, the one around 9 AU is Saturn, the trajectory at 12 AU is Uranus (who ends up switching positions with Neptune) and the uppermost curve at 16 AU is Neptune (who ends up closer to the Sun than Uranus). This plot is a magnification around the time when Jupiter and Saturn cross their 2:1 resonance and the system becomes unstable. The plot shows the phase until all encounters had stopped, which in this case happened when the period ratio between Jupiter and Saturn was . The final semi-major axes of the four planets are , so that the largest “error” is in Saturn’s orbit. A Fourier spectrum of Jupiter’s eccentricity at the end of the simulation gives , and . The amplitude of the term is a bit small, but well within a factor of 2 from the real value. Another Nice-model run from the same Gomes et al. (2005) study gave an essentially identical result. However, a third simulation gave , which is higher than the real value. In a fourth simulation, in which not only Saturn but also Jupiter encountered an ice giant, a value close to the real one was again recoverd, namely . Given the chaotic nature of planetary encounters and that the resulting values are close to the real one (0.044) or even larger (e.g. 0.059), we conclude that the Nice model is able to reproduce the secular architecture of the Jupiter-Saturn system.

In principle, depending on the initial separations between Saturn and the innermost ice giant, the Nice model can also give planetary evolutions in which encounters between Saturn and an ice giant do not occur (only the ice giants encountering each other; see Tsiganis et al., 2005). These evolutions have been rejected already in Tsiganis et al. (2005), because they lead to final mean eccentricities (and inclinations) that are too small and an orbital separation between Saturn and Uranus that is too narrow et the end. They also lead to a value of that is much too small compared to the current value, because the 2:1 resonance crossing alone is not capable of pumping the excitation of the mode, as we have seen earlier.

At this point, one might wonder whether the Nice model, in the version with Saturn-Uranus encounters, is the only model capable of this result. From the study reported in this paper, it seems likely that encounters among planets might be sufficient to excite the modes of the final secular system, without any need for Jupiter and Saturn crossing their mutual 2:1 resonance. In other words, one might envision a model where Jupiter and Saturn formed on circular orbits, well separated from each other in the beginning, so that was always larger than 2. These planets then had close encounters with other planetary embryos, which at the end left them on eccentric orbits with both the and modes excited.

A single encounter of an embryo with one planet would not work because, by kicking the eccentricity of one planet and not of the other, it would produce a secular system with (if the embryo encountered Saturn) or with (if the embryo encountered Jupiter). The real system is different from these two extremes. However, multiple encounters with one planet or with both of them should do the job. To achieve an estimate of the mass of the planetary embryo that could excite the secular modes of Jupiter and Saturn up to the observed values, we have run four sets of four simulations each. In each run we considered Jupiter, Saturn and one embryo, initially on circular orbits. The mass of the embryo was 1, 5, 10 and 15 Earth masses respectively, for the four sets of simulations. The initial location of the embryo was  AU, 8.0 AU, 10.1 AU, 10.7 AU, for the four simulations in each set, whereas Jupiter and Saturn were initially at AU and AU in all cases. In most simulations, the embryo was eventually ejected from the Solar System: in two runs the embryo collided with Jupiter. The values of the and coefficients for each set of simulations are reported in table 5. It turns out that the putative embryo had to be massive, of the order of . We stress that multiple embryos with the same total mass would not do an equal job, because the geometries of the encounters would be randomized, rather leading to dynamical friction instead of excitation. In fact, we did the same experiment with 100 Mars mass objects instead of a unique 10 Earth-mass object; it resulted in and being smaller than 0.001, demonstrating that an ensemble of small objects could not have excited the relevant modes to their current states. In conclusion, for the excitation of the mode to reach its current value, Jupiter and Saturn should have encountered Uranus or Neptune or a putative third ice giant of comparable mass. Therefore, for what concerns the excitation of the secular Fourier modes of the planetary orbits, a generic scenario of global instability and mutual scattering of the four giant planets, as originally proposed by Thommes et al. (1999) would work; the passage of Jupiter and Saturn through their mutual 2:1 resonance, which is specific to the Nice model relative to Thommes et al. (1999) (or Thommes et al., 2007, in which Jupiter and Saturn are initially locked in the 2:1 resonance) is not necessary.

Pierens and Nelson (2008), however, showed that the only possible final configuration achieved by Jupiter and Saturn in the gas disk is in their mutual 3:2 resonance. Unless alternative evolutions have been missed in that work, this result invalidates the initial planetary configurations considered in Thommes et al. (1999; 2007) and supports the Nice model, in particular in its newest version described in Morbidelli et al. (2007), where the four giant planets start locked in a quadruple resonance with .

[] [AU]
1 7.2 4.41 5.92
1 8 1.34 3.17
1 10.1 0.475 8.91
1 10.7 8.97 8.56
5 7.2 5.23 85.3
5 8.0 15.0 57.3
5 10.1 7.78 12.0
5 10.7 7.49 20.2
10 7.2 38.0 88.5
10 8.0 22.5 15.1
10 10.1 11.0 32.8
10 10.7 66.0 28.7
15 7.2 32.3 41.0
15 8.0 20.3 28.9
15 10.1 20.4 80.9
15 10.7 14.8 152.1
Table 5: Values of in Jupiter and in Saturn after ejecting planetary embryos of various masses from various original locations.

7 Ice giants and inclination constraints

Up to now we have focused our discussion on the secular evolution of the eccentricities of Jupiter and Saturn. However, the outer Solar System has two additional planets: Uranus and Neptune, which introduce the additional frequencies and in the secular evolution of the eccentricities. Thus, one should also be concerned about the correct excitation of the and modes in all planets. In addition, the planets have a rich secular dynamics in inclination, associated with the precession of their nodes. The excitation of the correct modes in the inclinations is a problem as crucial as that of the eccentricities.

The reason that we did not discuss these issues so far is because considerations based solely on the secular evolution of the eccentricities of Jupiter and Saturn have proved enough to guide us towards a solution, which is also valid in the more general problem. That is, these planetary encounters that are necessary to explain the Jupiter-Saturn secular architecture, can also explain the excitation of the and modes i.e.  and , and of the inclinations of the planets.

The excitation of the and modes does not appear to be very consraining. During the encounters between Saturn-Uranus and Uranus-Neptune encounters, the eccentricities of the ice giants become typically much larger than the current values. Thus, the combination of encounters and dynamical friction can produce a wide range of amplitudes of the and modes, including the current amplitudes. But we cannot exclude that other mechanisms, such as a sequence of mutual resonance crossing, could have led to the correct amplitudes, as well.

Conversely, the inclination excitation is particularly interesting, because it provides a strong, additional argument in favour of a violent evolution of the planets that involves mutual close encounters. In fact, the planets should form on essentially co-planar orbits, for the same reasons for which they should form on circular orbits: low relative velocities with respect to the planetesimals in the disk are necessary for the rapid formation of their cores. Once the planets are formed, tidal interactions with the gas disk damp the residual inclinations of the planets (Lubow & Ogilvie, 2001). After the disappearance of the gas, dynamical friction exerted by the remnant planetesimal disk would also damp the planetary inclinations. Thus, similar to the eccentricities, a relatively-late excitation mechanism is required to explain the inclinations of the planets. However, unlike the eccentricities, the passage across mean motion resonances does not significantly excite the inclinations because the resonant terms depending on the longitude of the nodes are at least quadratic in the inclinations. Secular resonances are also ineffective, if all planetary inclinations are initially small. The only mechanism by which inclinations can be efficiently increased is by close encounters. In fact, in the Nice model, the inclinations of Saturn, Uranus and Neptune, relative to the orbit of Jupiter, are well reproduced. A similar result holds also for the Thommes et al. (1999) model.

It is interesting to note that the eccentricities of the giant planets are about twice as large as the inclinations (respectively and radians or ). This is what one would expect, if both the eccentricities and the inclinations had been acquired by a combination of gravitational scattering and dynamical friction. Conversely, if encounters among the planets had never happened and the eccentricities had been acquired through specific resonance crossings, we would expect the planetary eccentricities to be much larger than the inclinations.

8 Conclusion

In this work we have demonstrated (see Fig. 2 and sects. 2.1 and 2.2) that the secular architecture of the giant planets, which have non negligible eccentricities and inclinations and specific amplitudes of the modes in the eccentricity evolution of Jupiter and Saturn, could not have been achieved if the planets migrated smoothly through a planetesimal disk, as originally envisioned in Malhotra (1993; see also Malhotra 1995; Hahn & Malhotra, 1999; Gomes et al., 2004). Thus, we believe that the community should no longer consider the smooth migration model as a valid template for the evolution of the solar system. Even repeated passages through mutual resonances, that would have happened if the planetary system was originally more compact than envisioned in Malhotra’s model, could not account for the orbital architecture of the giant planets, as we observe them today.

Instead, the outer planetary system had to evolve in a more violent way, in which the gas giants encountered the ice giants, or other rogue planets of equivalent mass. In this respect, the correct templates for the planetary evolution are the Nice model or the Thommes et al. model (or variants of these two). As we noted above, the Nice model seems to be, so far, more coherent and consistent in all its facets. We remark that encounters among the planets have been shown to also be an effective mechanism for the capture of the systems of irregular satellites (Nesvorny et al., 2007). These arguments, which are different and independent of those reported in this paper, also support a violent evolution scenarios of the outer solar system.

In a forthcoming paper we will investigate the orbital dynamics of the terrestrial planets in the context of the evolution of the giant planets that we have outlined in this work.

This work is part of the Helmholtz Alliance’s ’Planetary evolution and life’, which RB and AM thank for financial support. Exchanges between Nice and Thessaloniki have been funded by a PICS programme of France’s CNRS, and RB thanks the host KT for his hospitality during a recent visit. RG thanks Brasil’s CNPq and FAPERJ for financial support. HFL thanks NASA’s OSS and OPR programmes for support. Most of the simulations in this work were performed on the CRIMSON Beowulf cluster at OCA.

9 Bilbiography

Čuk, M. 2007, DPS 39, 60
D’Angelo, G., Lubow, S. & Bate, M., 2006 ApJ 652, 1698
Fernandez, J. A. & Ip, W.-H. 1984 Icarus 58, 109
Goldreich, P. & Sari, R. 2003, ApJ 585, 1024
Goldreich, P., Lithwick, Y. & Sari, R. 2004, ApJ 614, 497
Gomes, R., Morbidelli, A. & Levison, H. 2004 Icarus 170, 492
Gomes, R., Levison, H., Tsiganis, K. & Morbidelli, A. 2005 Nature 435,466
Hahn, J. & Malhotra, R. 1999 AJ 117,3041
Henrard, J. & Lemaître, A. 1983 Icarus55, 482
Henrard, J. 1993 in Dynamics Reported New series volume 2 117
Kley, W. & Dirksen, G. 2006, A&A 447, 369
Kokubo, E. & Ida, S. 1996, Icarus 123, 180
Kokubo, E. & Ida, S. 1998, Icarus 131, 171
Levison, H. & Duncan, M. 1994 Icarus 108, 18
Lubow, S. & Ogilvie, G. 2001 ApJ 560, 997
Malhotra, R. 1993 Nature 365, 819
Malhotra, R. 1995 AJ 110, 420
Masset F. & Snellgrove, M. 2001 MNRAS 320, 55
Michtchenko, T. & Malhotra, R. 2004 Icarus 168, 237
Michtchenko, T., Beaugé, C. & Ferraz-Mello S. 2008 MNRAS 391, 215
Morbidelli, A. 2002 Modern Celestial Mechanics - Aspects of Solar System Dynamics (Taylor & Francis, UK). Morbidelli A. & Crida, A. 2007 Icarus 191, 158
Morbidelli, A., Tsiganis, K., Crida, A., Levison, H. & Gomes, R. 2007 AJ 134,1790
Murray, C. & Dermott, S. 1999 Solar System Dynamics (Cambridge University Press, Cambridge, UK)
Neishtadt, A. 1984 J. App. Math. Mech 48, 133
Nelson, R. 2005 A&A 443 1067
Nesvorný, D., Vokrouhlický, D. & Morbidelli, A. 2007 AJ 133, 1962
Nobili, A., Milani, A. & Carpino, M. 1989 A&A 210, 313
Pierens, A. & Nelson, R. 2008 A&A 482, 333
Saslaw, W. 1985 ApJ 297, 49
Stewart, G. & Wetherill G. 1988 Icarus 74, 542
Tsiganis, K., Gomes, R., Morbidelli, A. & Levison, H. 2005 Nature 435, 459
Thommes, E., Duncan, M. & Levison, H. 1999 Nature 402, 635
Thommes, E., Nilsson, L & Murray, N. 2007 ApJ 656, 25
Wisdom, J. & Holman, M. 1991 AJ 102, 1528

Comments 0
Request comment
You are adding your first comment
How to quickly get a good reply:
  • Give credit where it’s Due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should ultimately help the author improve the paper.
  • Remember, the better we can be at sharing our knowledge with each other, the faster we can move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description