Dynamic habitability for Earth-like planets in 86 extrasolar planetary systems

In this paper we estimate the likelihood to find habitable Earth-like planets on stable orbits for 86 selected extrasolar planetary systems, where luminosity, effective temperature and stellar age are known. For determining the habitable zone (HZ) an integrated system approach is used taking into account a variety of climatological, biogeochemical, and geodynamical processes. Habitability is linked to the photosynthetic activity on the planetary surface. We find that habitability strongly depends on the age of the stellar system and the characteristics of a virtual Earth-like planet. In particular, the portion of land/ocean coverages plays an important role. We approximated the conditions for orbital stability using a method based on the Hill radius. Almost 60% of the investigated systems could harbour habitable Earth-like planets on stable orbits. In 18 extrasolar systems we find even better prerequisites for dynamic habitability than in our own solar system. In general our results are comparable to those with an HZ determination based only on climatic constraints. However, there are remarkable differences for land worlds and for systems older than about 7 Gyr.     

Mass-radius curve for extrasolar Earth-like planets and ocean planets

By comparison with the Earth-like planets and the large icy satellites of the Solar System, one can model the internal structure of extrasolar planets. The input parameters are the composition of the star (Fe/Si and Mg/Si), the Mg content of the mantle (Mg# = Mg/[Mg + Fe]), the amount of H₂O and the total mass of the planet. Equation of State (EoS) of the different materials that are likely to be present within such planets have been obtained thanks to recent progress in high-pressure experiments. They are used to compute the planetary radius as a function of the total mass. Based on accretion models and data on planetary differentiation, the internal structure is likely to consist of an iron-rich core, a silicate mantle and an outer silicate crust resulting from magma formation in the mantle. The amount of H₂O and the surface temperature control the possibility for these planets to harbor an ocean. In preparation to the interpretation of the forthcoming data from the CNES led CoRoT (Convection Rotation and Transit) mission and from ground-based observations, this paper investigates the relationship between radius and mass. If H₂O is not an important component (less than 0.1%) of the total mass of the planet, then a relation (R/REarth)=ab(M/MEarth) is calculated with (a,b)=(1,0.306) and (a,b)=(1,0.274) for ¹⁰⁻²MEarth<M<MEarth and MEarth<M<10MEarth, respectively. Calculations for a planet that contains 50% H₂O suggest that the radius would be more than 25% larger than that based on the Earth-like model, with (a,b)=(1.258,0.302) for ¹⁰⁻²MEarth<M<MEarth and (a,b)=(1.262,0.275) for MEarth<M<10MEarth, respectively. For a surface temperature of 300 K, the thickness of the ocean varies from 150 to 50 km for planets 1 to 10 times the Earth's mass, respectively. Application of this algorithm to bodies of the Solar System provides not only a good fit to most terrestrial planets and large icy satellites, but also insights for discussing future observations of exoplanets.     

Astrometric detection of Earth-like planets with OSI

The Orbiting Stellar Interferometer (OSI) is a space-based astrometric interferometer designed primarily for wide-angle astrometry. OSI is potentially capable of achieving astrometric accuracies of 1as in narrow angle (1°) astrometry. This paper discusses the implications for astrometric planet detection, specifically the detection of Earth-like planets around nearby stars. OSI has the potential to detect a limited number of planetary systems with Earths, if a number of technical problems are solved.     

The maximal runaway temperature of Earth-like planets

In Simpson’s (Simpson, G.C. [1927]. Mem. R. Meteorol. Soc. II (16), 69–95) classical derivation of the temperature of the Earth in the semi-gray model, the surface temperature diverges as the fourth root of the thermal radiation’s optical depth. No resolution to this apparent paradox was yet obtained under the strict semi-gray approximation. Using this approximation and a simplified approach, we study the saturation of the runaway greenhouse effect.First we generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption.Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and that any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200–1300 K irrespective of the nature of the greenhouse gas.We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of λ_cut, the wavelength above which the atmosphere becomes optically thick. Unless λ_cut is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has an interesting observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.     

Tidal interactions in multi-planet systems

We study systems of close orbiting planets evolving under the influence of tidal circularization. It is supposed that a commensurability forms through the action of disk induced migration and orbital circularization. After the system enters an inner cavity or the disk disperses the evolution continues under the influence of tides due to the central star which induce orbital circularization. We derive approximate analytic models that describe the evolution away from a general first order resonance that results from tidal circularization in a two planet system and which can be shown to be a direct consequence of the conservation of energy and angular momentum. We consider the situation when the system is initially very close to resonance and also when the system is between resonances. We also perform numerical simulations which confirm these models and then apply them to two and four planet systems chosen to have parameters related to the GJ 581 and HD 10180 systems. We also estimate the tidal dissipation rates through effective quality factors that could result in evolution to observed period ratios within the lifetimes of the systems. Thus the survival of, or degree of departure from, close commensurabilities in observed systems may be indicative of the effectiveness of tidal disipation, a feature which in turn may be related to the internal structure of the planets involved.     

The internal activity and thermal evolution of Earth-like planets

We model the internal thermal evolution of planets with Earth-like composition and masses ranging from 0.1 to 10 Earth masses over a period of 10 billion years. We also characterize the internal activity of the planets by the velocity of putative tectonic plates, the rate at which mantle material is processed through melting zones, and the time taken to process one mantle mass. The more massive the planet the larger its processing rate (?), which scales approximately as ?∝M0.8-1.0. The processing times for all the planets increase with time as they cool and become less active. As would be expected, the surface heat flow scales with planet mass. All planets have similar declines in mantle temperature except for the largest, in which pressure effects cause a larger decline. The larger planets have higher mantle temperatures over all times. The less massive the planet, the larger the decrease in core temperature with time. The core heat flow is also found to decrease more rapidly for smaller planet masses. Finally, rough predictions are made for the time required to generate an atmosphere from estimates of the time to degas water and carbon dioxide in mantle melting zones. The degassing times depend strongly on the initial temperature of the planet, but for the temperatures used in our model all the planets degas within ∼32 Ma after their formation.     

Searching for exoplanets with HEPS:I.detection probability of Earth-like planets in multiple systems

The astrometry method has great advantages in searching for exoplanets in the habitable zone around solar-like stars. However, the presence of multiple planets may cause a problem with degeneracy when trying to compute accurate planet parameters from observation data and reduce detectability. The degeneracy problem is extremely critical, especially in a space mission which has limited observation time and cadence. In this series of papers, we study the detectability of habitable Earth-mass planets in different types of multi-planet systems, aiming to find the most favorable targets for the potential space mission–Habitable ExoPlanet Survey(HEPS). In the first paper, we present an algorithm to find planets in the habitable zone around solar-like stars using astrometry. We find the detectability can be well described by planets' signal-to-noise ratio(SNR) and a defined parameter S = M2/(T1-T2)2, where M2 and T2are the mass and period of the second planet, respectively. T1 is the period of the planet in the habitable zone. The parameter S represents the influence of planetary architectures. We fit the detectability as a function of both the SNR of the planet in the habitable zone and the parameter S. An Earth-like planet in a habitable zone is harder to detect(with detectability PHP 80%) in a system with a hot Jupiter or warm Jupiter(within2 AU), in which the parameter S is large. These results can be used in target selections and to determine the priority of target stars for HEPS, especially when we select and rank nearby planet hosts with a single planet.     

Orbital stability of systems of closely-spaced planets

An investigation of the stability of systems of 1 M_⊕ (Earth-mass) bodies orbiting a Sun-like star has been conducted for virtual times reaching 10 billion years. For the majority of the tests, a symplectic integrator with a fixed timestep of between 1 and 10 days was employed; however, smaller timesteps and a Bulirsch-Stoer integrator were also selectively utilized to increase confidence in the results. In most cases, the planets were started on initially coplanar, circular orbits, and the longitudinal initial positions of neighboring planets were widely separated. The ratio of the semimajor axes of consecutive planets in each system was approximately uniform (so the spacing between consecutive planets increased slowly in terms of distance from the star). The stability time for a system was taken to be the time at which the orbits of two or more planets crossed. Our results show that, for a given class of system (e.g., three 1 M_⊕ planets), orbit crossing times vary with planetary spacing approximately as a power law over a wide range of separation in semimajor axis. Chaos tests indicate that deviations from this power law persist for changed initial longitudes and also for small but non-trivial changes in orbital spacing. We find that the stability time increases more rapidly at large initial orbital separations than the power-law dependence predicted from moderate initial orbital separations. Systems of five planets are less stable than systems of three planets for a specified semimajor axis spacing. Furthermore, systems of less massive planets can be packed more closely, being about as stable as 1 M_⊕ planets when the radial separation between planets is scaled using the mutual Hill radius. Finally, systems with retrograde planets can be packed substantially more closely than prograde systems with equal numbers of planets.     

On the protection of extrasolar Earth-like planets around K/M stars against galactic cosmic rays

Previous studies have shown that extrasolar Earth-like planets in close-in habitable zones around M-stars are weakly protected against galactic cosmic rays (GCRs), leading to a strongly increased particle flux to the top of the planetary atmosphere. Two main effects were held responsible for the weak shielding of such an exoplanet: (a) For a close-in planet, the planetary magnetic moment is strongly reduced by tidal locking. Therefore, such a close-in extrasolar planet is not protected by an extended magnetosphere. (b) The small orbital distance of the planet exposes it to a much denser stellar wind than that prevailing at larger orbital distances. This dense stellar wind leads to additional compression of the magnetosphere, which can further reduce the shielding efficiency against GCRs. In this work, we analyse and compare the effect of (a) and (b), showing that the stellar wind variation with orbital distance has little influence on the cosmic ray shielding. Instead, the weak shielding of M star planets can be attributed to their small magnetic moment. We further analyse how the planetary mass and composition influence the planetary magnetic moment, and thus modify the cosmic ray shielding efficiency. We show that more massive planets are not necessarily better protected against galactic cosmic rays, but that the planetary bulk composition can play an important role.     

About the Hill stability of extrasolar planets in stellar binary systems

We use a three dimensional generalization of Szebehely’s invariant relation obtained by us (Makó and Szenkovits, Celest. Mech. Dyn. Astron. 90, 51, 2004) in the elliptic restricted three-body problem, to establish more accurate criterion of the Hill stability. By using this criterion, the Hill stability of four extrasolar planets (γ Cephei Ab, Gliese 86 Ab, HD 41004 Ab and HD 41004 Bb) is investigated.     

Orbital stability of systems of closely-spaced planets,II: configurations with coorbital planets

We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 10⁸, 5 × 10⁵, and 2 × 10⁵ years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t _c /t ₀) = b β + c, where t _c is the crossing time, t ₀ = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.     

High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

The Leonard Award Address Presented 1998 July 27, Dublin,Ireland: On the difficulties of making Earth-like planets

Abstract— Here I discuss the series of events that led to the formation and evolution of our planet to examine why the Earth is unique in the solar system. A multitude of factors are involved: These begin with the initial size and angular momentum of the fragment that separated from a molecular cloud; such random factors are crucial in determining whether a planetary system or a double star develops from the resulting nebula. Another requirement is that there must be an adequate concentration of heavy elements to provide the 2% “rock” and “ice” components of the original nebula. An essential step in forming rocky planets in the inner nebula is the loss of gas and depletion of volatile elements, due to early solar activity that is linked to the mass of the central star. The lifetime of the gaseous nebula controls the formation of gas giants. In our system, fine timing was needed to form the gas giant, Jupiter, before the gas in the nebula was depleted. Although Uranus and Neptune eventually formed cores large enough to capture gas, they missed out and ended as ice giants. The early formation of Jupiter is responsible for the existence of the asteroid belt (and our supply of meteorites) and the small size of Mars, whereas the gas giant now acts as a gravitational shield for the terrestrial planets. The Earth and the other inner planets accreted long after the giant planets, from volatile-depleted planetesimals that were probably already differentiated into metallic cores and silicate mantles in a gas-free, inner nebula. The accumulation of the Earth from such planetesimals was essentially a stochastic process, accounting for the differences among the four rocky inner planets—including the startling contrast between those two apparent twins, Earth and Venus. Impact history and accretion of a few more or less planetesimals were apparently crucial. The origin of the Moon by a single massive impact with a body larger than Mars accounts for the obliquity (and its stability) and spin of the Earth, in addition to explaining the angular momentum, orbital characteristics, and unique composition of the Moon. Plate tectonics (unique among the terrestrial planets) led to the development of the continental crust on the Earth, an essential platform for the evolution of Homo sapiens. Random major impacts have punctuated the geological record, accentuating the directionless course of evolution. Thus a massive asteroidal impact terminated the Cretaceous Period, resulted in the extinction of at least 70% of species living at that time, and led to the rise of mammals. This sequence of events that resulted in the formation and evolution of our planet were thus unique within our system. The individual nature of the eight planets is repeated among the 60-odd satellites—no two appear identical. This survey of our solar system raises the question whether the random sequence of events that led to the formation of the Earth are likely to be repeated in detail elsewhere. Preliminary evidence from the “new planets” is not reassuring. The discovery of other planetary systems has removed the previous belief that they would consist of a central star surrounded by an inner zone of rocky planets and an outer zone of giant planets beyond a few astronomical units (AU). Jupiter-sized bodies in close orbits around other stars probably formed in a similar manner to our giant planets at several astronomical units from their parent star and, subsequently, migrated inwards becoming stranded in close but stable orbits as “hot Jupiters”, when the nebula gas was depleted. Such events would prevent the formation of terrestrial-type planets in such systems.     

On the stability of the inner planets,satellites and close binaries

I consider the range of Hill stability in the restricted circular problem of three bodies when the larger one of the two principal bodies has a finite oblateness. I show that the range r satisfies the equation

$\text{r =}\text{1? μ}\text{C}{}_{\mathrm{cr}}\text{? 3μ + μr}{}^2\text{? v (1? μ)r}{}^{\mathrm{?3}}\text{±(2 + 3v)}\text{(1 ? μ}\text{1 +}\text{3v}\text{r}{}^2\text{r}\text{,}$

where μ is the mass parameter and v is an oblateness parameter. This result is applied to the solar system, the Earth-Moon system and binary star systems. It is then shown that, all the inner planets of the solar system, the great majority of asteroids and some short-period comets are Hill stable, that direct artificial satellites of the Earth are more stable than retrograde ones, and that contact binaries possess cores between which no mass exchange takes place.     

On the dynamics of extrasolar planetary systems under dissipation: Migration of planets

We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.     

Stability of Trojan planets in multi-planetary systems

Today there are more than 340 extra-solar planets in about 270 extra-solar systems confirmed. Besides the observed planets there exists also the possibility of a Trojan planet moving in the same orbit as the Jupiter-like planet. In our investigation we take also into account the habitability of a Trojan planet and whether such a terrestrial planet stays in the habitable zone. Its stability was investigated for multi-planetary systems, where one of the detected giant planets moves partly or completely in the habitable zone. By using numerical computations, we studied the orbital behaviour up to 10⁷ years and determined the size of the stable regions around the Lagrangian equilibrium points for different dynamical models for fictitious Trojans. We also examined the interaction of the Trojan planets with a second or third giant planet, by varying its semimajor axis and eccentricity. We have found two systems (HD 155358 and HD 69830) that can host habitable Trojan planets. Another aim of this work was to determine the size of the stable region around the Lagrangian equilibrium points in the restricted three body problem for small mass ratios μ of the primaries μ ≤ 0.001 (e.g. Neptune mass of the secondary and smaller masses). We established a simple relation for the size depending on μ and the eccentricity.     

On the stability of hierarchical four-body systems

Generalized Jacobian coordinates can be used to decompose anN-body dynamical system intoN-1 2-body systems coupled by perturbations. Hierarchical stability is defined as the property of preserving the hierarchical arrangement of these 2-body subsystems in such a way that orbit crossing is avoided. ForN=3 hierarchical stability can be ensured for an arbitrary span of time depending on the integralz=c ² h (angular momentum squared times energy): if it is smaller than a critical value, defined by theL ₂ collinear equilibrium configuration, then the three possible hierarchical arrangements correspond to three disconnected subsets of the invariant manifold in the phase space (and in the configuration space as well; see Milani and Nobili, 1983a). The same definitions can be extended, with the Jacobian formalism, to an arbitrary hierarchical arrangement ofN≥4 bodies, and the main confinement condition, the Easton inequality, can also be extended but it no longer provides separate regions of trapped motion, whatever is the value ofz for the wholeN-body system,N≥4. However, thez criterion of hierarchical stability applies to every 3-body subsystem, whosez ‘integral’ will of course vary in time because of the perturbations from the other bodies. In theN=4 case we decompose the system into two 3-body subsystems whosec ² h ‘integrals’,z ₂₃ andz ₃₄, att=0 are assumed to be smaller than the corresponding critical values \(\tilde z_{23} \) and \(\tilde z_{34} \) , so that both the subsystems are initially hierarchically stable. Then the hierarchical arrangement of the 4 bodies cannot be broken until eitherz ₂₃ orz ₃₄ is changed by an amount \(\tilde z_{ij} - z_{ij} \left( 0 \right)\) ; that is the whole system is hierarchically stable for a time spain not shorter than the minimum between \(\Delta t_{23} = {{\left( {\tilde z_{23} - z_{23} \left( 0 \right)} \right)} \mathord{\left/ {\vphantom {{\left( {\tilde z_{23} - z_{23} \left( 0 \right)} \right)} {\dot z_{23} }}} \right. \kern-0em} {\dot z_{23} }}\) and \(\Delta t_{34} = {{\left( {\tilde z_{34} - z_{34} \left( 0 \right)} \right)} \mathord{\left/ {\vphantom {{\left( {\tilde z_{34} - z_{34} \left( 0 \right)} \right)} {\dot z_{34} }}} \right. \kern-0em} {\dot z_{34} }}\) . To estimate how long is this stability time, two main steps are required. First the perturbing potentials have to be developed in series; the relevant small parameters are some combinations of mass ratios and length ratios, the? _ij of Roy and Walker. When an appropriate perturbation theory is based on the? _ij , the asymptotic expansions are much more rapidly decreasing than the usual expansions in powers of the mass ratios (as in the classical Lagrange perturbation theory) and can be extended also to cases such as lunar theory or double binaries. The second step is the computation of the time derivatives \(\dot z_{ij} \) (we limit ourselves to the planar case). To assess the long term behaviour of the system, we can neglect the short-periodic perturbations and discuss only the long-periodic and the secular perturbations. By using a Poisson bracket formalism, a generalization of Lagrange theorem for semimajor axes and a generalization of the classical first order theories for eccentricities and pericenters, we prove that thez _ij do not undergo any secular perturbation, because of the interaction with the other subsystem, at the first order in the? _ik . After the long-periodic perturbations have been accounted for, and apart from the small divisors problems that could arise both from ordinary and secular resonances, only the second order terms have to be considered in the computation of Δt ₂₃, Δt ₃₄. A full second order perturbative theory is beyond the scope of this paper; however an order-of-magnitude lower estimate of the Δt _ij can be obtained with the very pessimistic assumption that essentially all the second order terms affect in a secular way thez _ij . The same method could be applied also toN≥5 body systems. Since almost everyN-body system existing in nature is strongly hierarchical, the product of two? _ij is very small for almost all the real astronomical problems. As an example, the hierarchical stability of the 4-body system Sun, Mercury, Venus, and Jupiter is investigated; this system turns out to be stable for at least 110 million years. Although this hierarchical stability time is ～10 times less than the real age of the Solar System, taking into account that many pessimistic assumptions have been done we can conclude that the stability of the Solar System is no more a forbidden problem for Celestial Mechanics.     

On the stability of circumbinary planetary systems

The dynamics of circumbinary planetary systems (the systems in which the planets orbit a central binary) with a small binary mass ratio discovered to date is considered. The domains of chaotic motion have been revealed in the “pericentric distance–eccentricity” plane of initial conditions for the planetary orbits through numerical experiments. Based on an analytical criterion for the chaoticity of planetary orbits in binary star systems, we have constructed theoretical curves that describe the global boundary of the chaotic zone around the central binary for each of the systems. In addition, based on Mardling’s theory describing the separate resonance “teeth” (corresponding to integer resonances between the orbital periods of a planet and the binary), we have constructed the local boundaries of chaos. Both theoretical models are shown to describe adequately the boundaries of chaos on the numerically constructed stability diagrams, suggesting that these theories are efficient in providing analytical criteria for the chaoticity of planetary orbits.