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.     

Exoplanetary systems: The role of an equilibrium at high mutual inclination in shaping the global behavior of the 3-D secular planetary three-body problem

On the basis of a high-order (order 12) expansion of the perturbative potential in powers of the eccentricities and the inclinations, we analyze the secular interactions of two non-coplanar planets which are not in mean-motion resonance. The model is based on the planetary three-body problem which can be reduced to two degrees of freedom by the well-known elimination of the nodes [Jacobi, C.G.J., 1842. Astron. Nachr. XX, 81-102]. We introduce non-singular canonical variables which bring forward the symmetries of the problem. The main dynamical features depend on the location and stability of the equilibria which are easily found with our analytical model. We find that there exists an equilibrium when both eccentricities are zero. When the mutual inclination is small, this equilibrium is stable, but for larger mutual inclination it becomes unstable, generating a large chaotic zone and, by bifurcation, two regular regions, the so-called Kozai resonances. This analytical study which depends on only two parameters (the ratio of the semi-major axes and the mass ratio of the planets) makes possible a large survey of the problem and enables us to identify and quantify its main dynamical features, periodic orbits, regular and chaotic zones, etc. The results of our analytical model are illustrated and confirmed by numerical integrations.     

Modeling the 3-D secular planetary three-body problem: Discussion on the outer υ Andromedae planetary system

The three-dimensional secular behavior of a system composed of a central star and two massive planets is modeled semi-analytically in the frame of the general three-body problem. The main dynamical features of the system are presented in geometrical pictures allowing us to investigate a large domain of the phase space of this problem without time-expensive numerical integrations of the equations of motion and without any restriction on the magnitude of the planetary eccentricities, inclinations and mutual distance. Several regimes of motion of the system are observed. With respect to the secular angle Δ?, possible motions are circulations, oscillations (around 0° and 180°), and high-eccentricity/inclination librations in secular resonances. With respect to the arguments of pericenter, ω₁ and ω₂, possible motions are direct circulation and high-inclination libration around ±90° in the Lidov-Kozai resonance. The regions of transition between domains of different regimes of motion are characterized by chaotic behavior. We apply the analysis to the case of the two outer planets of the υ Andromedae system, observed edge-on. The topology of the 3-D phase space of this system is investigated in detail by means of surfaces of section, periodic orbits and dynamical spectra, mapping techniques and numerical simulations. We obtain the general structure of the phase space, and the boundaries of the spatial secular stability. We find that this system is secularly stable in a large domain of eccentricities and inclinations.     

Secular apsidal configuration of non-resonant exoplanetary systems

Using a high-order (order 12) expansion of the perturbative potential in powers of eccentricities [Libert, A.-S., Henrard, J., 2005. Celest. Mech. Dynam. Astron. 93, 187-200], we study the secular effects of two coplanar planets which are not in mean motion resonances. The main results concern eccentricity variations, oscillation amplitude of the angular difference of the apsidal lines (Δ?) and frequency of such an oscillation. We show that this analytical approach describes correctly the behaviour of most of the exosystems and underlines the known limitations of the linear Laplace-Lagrange theory. Apsidal configuration of υ Andromedae, HD 168443, HD 169830, HD 38529, HD 74156 and HD 12661 are examined. We also point out the great sensitivity of the υ Andromedae system to the initial values (e₁(0), e₂(0) or Δ?(0)).     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.     

Secular dynamics of the three-body problem: application to the υ Andromedae planetary system

The discovery of extra-solar planetary systems with multiple planets in highly eccentric orbits (∼0.1-0.6), in contrast with our own Solar System, makes classical secular perturbation analysis very limited. In this paper, we use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in co-planar orbits. We show that the secular dynamics of this system can be described using only two parameters, the ratios of the semi-major axes and the planetary masses. The main dynamical features of the system are presented in geometrical pictures that allows us to investigate a large domain of the phase space of this three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. The topology of the phase space is also investigated in detail by means of spectral map techniques, which allow us to detect the separatrix of a non-linear secular apsidal resonance. Finally, the qualitative study is supplemented by direct numerical integrations. Three different regimes of secular motion with respect to the secular angle Δ? are possible: they are circulation, oscillation (around 0° and 180°), and high eccentricity libration in a non-linear secular resonance. The first two regimes are a continuous extension of the classical linear secular perturbation theory; the last is a new feature, hitherto unknown, in the secular dynamics of the three-body problem. We apply the analysis to the case of the two outer planets in the υ Andromedae system, and obtain its periodic and ordinary orbits, the general structure of its secular phase space, and the boundaries of its secular stability; we find that this system is secularly stable over a large domain of eccentricities. Applying this analysis to a wide range of planetary mass and semi-major axis ratios (centered about the υ Andromedae parameters), we find that apsidal oscillation dominates the secular phase space of the three-body coplanar system, and that the non-linear secular resonance is also a common feature.     

Tidal interactions and disruptions of giant planets on highly eccentric orbits

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.     

Insolation patterns on synchronous exoplanets with obliquity

A previous paper [Dobrovolskis, A.R., 2007. Icarus 192, 1-23] showed that eccentricity can have profound effects on the climate, habitability, and detectability of extrasolar planets. This complementary study shows that obliquity can have comparable effects.The known exoplanets exhibit a wide range of orbital eccentricities, but those within several million kilometers of their suns are generally in near-circular orbits. This fact is widely attributed to the dissipation of tides in the planets. Tides in a planet affect its spin even more than its orbit, and such tidally evolved planets often are assumed to be in synchronous rotation, so that their rotation periods are identical to their orbital periods. The canonical example of synchronous spin is the way that our Moon always keeps nearly the same hemisphere facing the Earth.Tides also tend to reduce the planet’s obliquity (the angle between its spin and orbital angular velocities). However, orbit precession can cause the rotation to become locked in a “Cassini state”, where it retains a nearly constant non-zero obliquity. For example, our Moon maintains an obliquity of about 6.7° with respect to its orbit about the Earth. In comparison, stable Cassini states can exist for practically any obliquity up to ∼90° or more for planets of binary stars, or in multi-planet systems with high mutual inclinations, such as are produced by scattering or by the Kozai mechanism.This work considers planets in synchronous rotation with circular orbits, but arbitrary obliquity β; this affects the distribution of insolation over the planet’s surface, particularly near its poles. For β=0, one hemisphere bakes in perpetual sunshine, while the opposite hemisphere experiences eternal darkness. As β increases, the region of permanent daylight and the antipodal realm of endless night both shrink, while a more temperate area of alternating day and night spreads in longitude, and especially in latitude. The regions of permanent day or night disappear at β=90°. The insolation regime passes through several more transitions as β continues to increase toward 180°, but the surface distribution of insolation remains non-uniform in both latitude and longitude.Thus obliquity, like eccentricity, can protect certain areas of the planet from the worst extremes of temperature and solar radiation, and can improve the planet’s habitability. These results also have implications for the direct detectability of extrasolar planets, and for the interpretation of their thermal emissions.     

Obliquity evolution of extrasolar terrestrial planets

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ∼1 AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio (m/Mp) and initial obliquity (γ₀), through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized or larger bodies determine the planetary spin and form satellites. Since the giant impacts would be isotropic, tilted spins (sinγ₀∼1) are more likely to be produced than straight ones (sinγ₀∼0). The ratio m/Mp is dependent on the impact parameters and impactors' mass. However, most of previous studies on tidal evolution of the planet-satellite systems have focused on a particular case of the Earth-Moon systems in which m/Mp?0.0125 and γ₀∼10° or the two-body planar problem in which γ₀=0° and stellar torque is neglected. We numerically integrated the evolution of planetary spin and a satellite orbit with various m/Mp (from 0.0025 to 0.05) and γ₀ (from 0° to 180°), taking into account the stellar torques and precessional motions of the spin and the orbit. We start with the spin axis that almost coincides with the satellite orbit normal, assuming that the spin and the satellite are formed by one dominant impact. With initially straight spins, the evolution is similar to that of the Earth-Moon system. The satellite monotonically recedes from the planet until synchronous state between the spin period and the satellite orbital period is realized. The obliquity gradually increases initially but it starts decreasing down to zero as approaching the synchronous state. However, we have found that the evolution with initially tiled spins is completely different. The satellite's orbit migrates outward with almost constant obliquity until the orbit reaches the critical radius ∼10-20 planetary radii, but then the migration is reversed to inward one. At the reversal, the obliquity starts oscillation with large amplitude. The oscillation gradually ceases and the obliquity is reduced to ∼0° during the inward migration. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. We found that the character change of precession about total angular momentum vector into that about the planetary orbit normal is responsible for the oscillation with large amplitude and the reversal of migration. With the results of numerical integration and analytical arguments, we divided the m/Mp-γ₀ space into the regions of the qualitatively different evolution. The peculiar tidal evolution with initially tiled spins give deep insights into dynamics of extrasolar planet-satellite systems and discussions of surface environments of the planets.     

Internal structure of massive terrestrial planets

Planetary formation models predict the existence of massive terrestrial planets and experiments are now being designed that should succeed in discovering them and measuring their masses and radii. We calculate internal structures of planets with one to ten times the mass of the Earth (Super-Earths) to obtain scaling laws for total radius, mantle thickness, core size and average density as a function of mass. We explore different compositions and obtain a scaling law of R∝M0.267-0.272 for Super-Earths. We also study a second family of planets, Super-Mercuries with masses ranging from one mercury-mass to ten mercury-masses with similar composition to the Earth's but with a larger core mass fraction. We explore the effect of surface temperature and core mass fraction on the scaling laws for these planets. The scaling law obtained for the Super-Mercuries is R∝M∼0.3.     

The formation and habitability of terrestrial planets in the presence of close-in giant planets

‘Hot jupiters,’ giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form [P.J. Armitage, A reduced efficiency of terrestrial planet formation following giant planet migration, Astrophys. J. 582 (2003) L47-L50]. We present results of simulations of terrestrial planet formation in the presence of hot/warm jupiters, broadly defined as having orbital radii ?0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot/warm jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g., 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with close-in giant planets.     

On the possibility of earth-type habitable planets around 47 UMa

We investigate whether Earth-type habitable planets can in principle exist in the planetary system of 47 UMa. The system of 47 UMa consists of two Jupiter-size planets beyond the outer edge of the stellar habitable zone, and thus resembles our own Solar System most closely compared to all exosolar planetary systems discovered so far. Our study of habitability deliberately follows an Earth-based view according to the concept of Franck and colleagues, which assumes the long-term possibility of photosynthetic biomass production under geodynamic conditions. Consequently, a broad variety of climatological, biogeochemical, and geodynamical processes involved in the generation of photosynthesis-driven life conditions is taken into account. The stellar luminosity and the age of the star/planet system are of fundamental importance for planetary habitability. Our study considers different types of planetary continental growth models and takes into account a careful assessment of the stellar parameters. In the event of successful formation and orbital stability, two subjects of intense research, we find that Earth-type habitable planets around 47 UMa are in principle possible! The likelihood of those planets is increased if assumed that 47 UMa is relatively young (?6 Gyr) and has a relatively small stellar luminosity as permitted by the observational range of those parameters.     

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.     

Spin states and climates of eccentric exoplanets

The known extrasolar planets exhibit a wide range of orbital eccentricities e. This has a profound influence on their rotations and climates. Because of tides in their interiors, mostly solid exoplanets are expected eventually to despin to a state of spin-orbit resonance, where the orbital period is some integer or half-integer times the rotation period. The most important of these resonances is the synchronous state, where the planet's spin period exactly equals its orbital period (like Earth's Moon, and indeed most of the regular satellites in the Solar System). Such planets seem doomed to roast on one side and freeze on the other. However, synchronous planets rock back and forth by an angle of ∼2Arcsine with respect to the sub-stellar point. For e=0.055 (as for the Moon), this optical libration amounts to only ∼6°; but for a synchronous planet with e=0.50, for example, it would rise to ∼59°. This greatly expands the temperate “twilight zone” near the terminator and considerably improves the planet's prospects for habitability. For e?0.72389, the optical libration exceeds 90°; for such planets, the sector of permanent night vanishes, while the sunniest region splits in two. Furthermore, the synchronous state is not the only possible spin resonance. For example, Mercury (with e≈0.206) has an orbital period exactly 1.5 times its rotation period. A terrestrial exoplanet with e=0.40, say, is liable to have an orbital period of 2.0, 2.5, or 3.0 times its spin period. The corresponding insolation patterns are generally complicated, and all different from the synchronous state. Yet these non-synchronous resonances also protect certain longitudes from the worst extremes of temperature and solar radiation, and improve the planet's habitability, compared to non-resonant rotation. These results also have implications for the direct detectability of extrasolar planets, and the interpretation of their thermal emissions.     

The web of three-planet resonances in the outer Solar System

In this paper we numerically detect the web of three-planet resonances (i.e., resonances among mean anomalies, nodes and perihelia of three planets) with respect to the variation of the semi-major axis of Saturn and Jupiter, in a model including the planets from Jupiter to Neptune. The measure confirms the relevance of these resonances in the long-term evolution of the outer Solar System and provides a technique to identify some of the related coefficients.     

Oligarchic growth of giant planets

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as “oligarchic growth.” Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk (∼10 × minimum-mass) is required to produce giant planet core-sized bodies (∼10 M_⊕) within the lifetime of the nebular gas (?10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of noninteracting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System’s two outer giant planets at their current locations on any timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.     

Giant planet migration through the action of disk torques and planet-planet scattering

This paper presents a parametric study of giant planet migration through the combined action of disk torques and planet-planet scattering. The torques exerted on planets during Type II migration in circumstellar disks readily decrease the semi-major axes a, whereas scattering between planets increases the orbital eccentricities ?. This paper presents a parametric exploration of the possible parameter space for this migration scenario using two (initial) planetary mass distributions and a range of values for the time scale of eccentricity damping (due to the disk). For each class of systems, many realizations of the simulations are performed in order to determine the distributions of the resulting orbital elements of the surviving planets; this paper presents the results of ∼8500 numerical experiments. Our goal is to study the physics of this particular migration mechanism and to test it against observations of extrasolar planets. The action of disk torques and planet-planet scattering results in a distribution of final orbital elements that fills the a-? plane, in rough agreement with the orbital elements of observed extrasolar planets. In addition to specifying the orbital elements, we characterize this migration mechanism by finding the percentages of ejected and accreted planets, the number of collisions, the dependence of outcomes on planetary masses, the time spent in 2:1 and 3:1 resonances, and the effects of the planetary IMF. We also determine the distribution of inclination angles of surviving planets and the distribution of ejection speeds for exiled planets.     

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.