Stability of inclined orbits of terrestrial planets in habitable zones

In long-term stability studies of terrestrial planets moving in the habitable zone (HZ) of a sun-like star, we distinguish four different configurations: (i) planets moving in binary star systems, (ii) the inner type (where the gas giant moves outside the HZ), (iii) the outer type (where the gas giant is closer to the star, than the HZ) and (iv) the Trojan type (where the gas giant moves in the HZ). Since earlier calculations indicated, that the stability of the motion in the HZ also depends on the inclination of the terrestrial planet orbits, we present a detailed numerical investigation to show correlations between the eccentricity, the mass and the distance of the giant planet for various inclinations of the terrestrial planets. The orbital stability of the HZ was examined for all four configurations stated above. While we could find hardly any stable orbits for the first three types for inclinations higher than 40^°, the Trojan planets can be stable up to an inclination of 60^°. Additionally, we could also find some stabilizing effects of the inclination for the first three types. As dynamical model we used the elliptic restricted three-body problem, which consists of two massive and one mass-less body. This allows an application to all detected and future extrasolar single planet systems.     

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.     

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.     

On the habitability of the OGLE-2006-BLG-109L planetary system

We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analogue of the Solar system. Our work is inspired by the work of Malhotra & Minton. Using the linear Laplace–Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the subsystem of Jupiters may be found in a dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semimajor axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.     

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.     

The dynamical architecture and habitable zones of the quintuplet planetary system 55 Cancri

We perform numerical simulations to study the secular orbital evolution and dynamical structure of the quintuplet planetary system 55 Cancri with the self-consistent orbital solutions by Fischer and coworkers. In the simulations, we show that this sys-tem can be stable for at least 108 yr. In addition, we extensively investigate the planetary configuration of four outer companions with one terrestrial planet in the wide region of 0.790 AU ≤ a ≤ 5.900 AU to examine the existence of potential asteroid structure and Habitable Zones (HZs). We show that there are unstable regions for orbits about 4:1, 3:1 and 5:2 mean motion resonances (MMRs) of the outermost planet in the system, and sev-eral stable orbits can remain at 3:2 and 1:1 MMRs, which resembles the asteroid belt in the solar system. From a dynamical viewpoint, proper HZ candidates for the existence of more potential terrestrial planets reside in the wide area between 1.0 AU and 2.3 AU with relatively low eccentricities.     

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.     

Using long-term transit timing to detect terrestrial planets

We propose that the presence of additional planets in extrasolar planetary systems can be detected by long-term transit timing studies. If a transiting planet is on an eccentric orbit then the presence of another planet causes a secular advance of the transiting planet's pericentre over and above the effect of general relativity. Although this secular effect is impractical to detect over a small number of orbits, it causes long-term differences when future transits occur, much like the long-term decay observed in pulsars. Measuring this transit-timing delay would thus allow the detection of either one or more additional planets in the system or the first measurements of non-zero oblateness ( J ₂) of the central stars.     

Orbital Migration Of Giant Planets: Using Numerical Integration To Investigate Consequences For Other Bodies

A number of extrasolar planets have been detected in close orbits around nearby stars. It is probable that these planets did not form in these orbits but migrated from their formation locations beyond the ice line. Orbital migration mechanisms involving angular momentum transfer through tidal interactions between the planets and circumstellar gas-dust disks or by gravitational interaction with a residual planetesimal disk together with several means of halting inward migration have been identified. These offer plausible schemes to explain the orbits of observed extrasolar giant planets and giant planets within the Solar System. Recent advances in numerical integration methods and in the power of computer workstations have allowed these techniques to be applied to modelling directly the mechanisms and consequences of orbital migration in the Solar System. There is now potential for these techniques also to be applied to modelling the consequences of the orbital migration of planets in the observed exoplanetary systems. In particular the detailed investigation of the stability of terrestrial planets in the habitable zone of these systems and the formation of terrestrial planets after the dissipation of the gas disk is now possible. The stability of terrestrial planets in the habitable zone of selected exoplanetary systems has been established and the possibility of the accretion of terrestrial planets in these systems is being investigated by the author in collaboration with Barrie W. Jones (Open University), and with John Chambers (NASA-Ames) and Mark Bailey of Armagh Observatory, using numerical integration. The direct simulation of orbital migration by planetesimal scattering must probably await faster hardware and/or more efficient algorithms. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

系外类地行星的内部结构及大气成分研究进展

系外类地行星是目前搜寻地外生命的主要目标.随着观测仪器的发展,现在已经能探测到低于10个地球质量的系外行星.该文简要回顾了系外类地行星的形成与演化,介绍了当前研究它们内部结构的模型和方法,以及由此得出的类地行星质量-半径关系.同时,对应不同的行星初始物质成分,讨论了各种可能的大气结构.最后介绍了未来的空间任务在相关方面的工作.     

地外生命搜索和太阳系外的行星的发现

阐述地外生命搜索的意义，手段和现状，地外生命搜索的SETI计划及其发展，并着重介绍太阳系外的行星系统的发现，它的目的，成果以及未来。太阳系外的行星系统的发展是当代天文学最时髦的，也将是未来21世纪成果最丰富的研究领域之一，从1992年第一个确认了脉冲星PSR 1257＋12的行星系统以来，1995年确认主序星51Peg有一颗行星，至2001年5月，已经发现了60个太阳系外的行星系统，太阳系外的行星系统的发现与地外生命搜索研究是密不可分的，新的发现也提出了很多新的谜，这项研究近10年来发展很快，它的研究也促进了航天学，宇宙化学，天文生物学乃至哲学等其他学科的发展。     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

Tilting Styx and Nix but not Uranus with a Spin-Precession-Mean-motion resonance

A Hamiltonian model is constructed for the spin axis of a planet perturbed by a nearby planet with both planets in orbit about a star. We expand the planet–planet gravitational potential perturbation to first order in orbital inclinations and eccentricities, finding terms describing spin resonances involving the spin precession rate and the two planetary mean motions. Convergent planetary migration allows the spinning planet to be captured into spin resonance. With initial obliquity near zero, the spin resonance can lift the planet’s obliquity to near 90\(^\circ \) or 180\(^\circ \) depending upon whether the spin resonance is first or zeroth order in inclination. Past capture of Uranus into such a spin resonance could give an alternative non-collisional scenario accounting for Uranus’s high obliquity. However, we find that the time spent in spin resonance must be so long that this scenario cannot be responsible for Uranus’s high obliquity. Our model can be used to study spin resonance in satellite systems. Our Hamiltonian model explains how Styx and Nix can be tilted to high obliquity via outward migration of Charon, a phenomenon previously seen in numerical simulations.     

Making the Earth: Combining dynamics and chemistry in the Solar System

No terrestrial planet formation simulation completed to date has considered the detailed chemical composition of the planets produced. While many have considered possible water contents and late veneer compositions, none have examined the bulk elemental abundances of the planets produced as an important check of formation models. Here we report on the first study of this type. Bulk elemental abundances based on disk equilibrium studies have been determined for the simulated terrestrial planets of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58]. These abundances are in excellent agreement with observed planetary values, indicating that the models of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58] are successfully producing planets comparable to those of the Solar System in terms of both their dynamical and chemical properties. Significant amounts of water are accreted in the present simulations, implying that the terrestrial planets form “wet” and do not need significant water delivery from other sources. Under the assumption of equilibrium controlled chemistry, the biogenic species N and C still need to be delivered to the Earth as they are not accreted in significant proportions during the formation process. Negligible solar photospheric pollution is produced by the planetary formation process. Assuming similar levels of pollution in other planetary systems, this in turn implies that the high metallicity trend observed in extrasolar planetary systems is in fact primordial.     

Stability of fictitious Trojan planets in extrasolar systems

Our work deals with the dynamical possibility that in extrasolar planetary systems a terrestrial planet may have stable orbits in a 1:1 mean motion resonance with a Jovian like planet. We studied the motion of fictitious Trojans around the Lagrangian points L₄/L₅ and checked the stability and/or chaoticity of their motion with the aid of the Lyapunov Indicators and the maximum eccentricity. The computations were carried out using the dynamical model of the elliptic restricted three‐body problem that consists of a central star, a gas giant moving in the habitable zone, and a massless terrestrial planet. We found 3 new systems where the gas giant lies in the habitable zone, namely HD99109, HD101930, and HD33564. Additionally we investigated all known extrasolar planetary systems where the giant planet lies partly or fully in the habitable zone. The results show that the orbits around the Lagrangian points L₄/L₅ of all investigated systems are stable for long times (10⁷ revolutions). (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

How common are Earth-Moon planetary systems?

The Earth’s comparatively massive moon, formed via a giant impact on the proto-Earth, has played an important role in the development of life on our planet, both in the history and strength of the ocean tides and in stabilizing the chaotic spin of our planet. Here we show that massive moons orbiting terrestrial planets are not rare. A large set of simulations by Morishima et al. (Morishima, R., Stadel, J., Moore, B. [2010]. Icarus. 207, 517-535), where Earth-like planets in the habitable zone form, provides the raw simulation data for our study. We use limits on the collision parameters that may guarantee the formation of a circumplanetary disk after a protoplanet collision that could form a satellite and study the collision history and the long term evolution of the satellites qualitatively. In addition, we estimate and quantify the uncertainties in each step of our study. We find that giant impacts with the required energy and orbital parameters for producing a binary planetary system do occur with more than 1 in 12 terrestrial planets hosting a massive moon, with a low-end estimate of 1 in 45 and a high-end estimate of 1 in 4.     

The interaction of the stellar wind with an extrasolar planet—3D hybrid and drift-kinetic simulations

To be able to simulate the interaction of extrasolar planets with the stellar wind, a number of planetary parameters are required. Some of these (like planetary mass and radius) can be obtained directly from observational data. Other properties are not known very precisely. For example, up to now, there is no observation providing information on the strength of planetary magnetic moments. However, there is good reason to expect only very small magnetic moments for planets in very close orbits around their stars (like HD 209458 b and OGLE-TR-56 b). Thus, as a first step towards a more complete treatment, it seems reasonable to treat the interaction of the stellar wind with an unmagnetized planet. Calculations were performed for a nonconducting as well as for a weakly conducting planet. The interaction with the stellar wind and the resulting induced magnetosphere was simulated using a three dimensional hybrid code as well as in the drift-kinetic approximation. The effect of a interplanetary magnetic field oriented perpendicular to the incoming stellar wind was included. In the case of a weakly conducting body an asymmetrical Mach cone is formed, whereas for a nonconducting body no Mach cone is observed. These investigations will serve as the first step in the search for particular effects occurring at extrasolar planets, which could possibly lead to observable effects, e.g. radio emission. The results are also relevant for plasma structures near weakly conducting, unmagnetized bodies like the Earth's moon.