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)     

Dynamics of possible Trojan planets in binary systems

This paper is devoted to the dynamical stability of possible Trojan planets in binaries and in binary systems where one of the substellar companions is not larger than a brown dwarf. Using numerical integrations, we investigated how the size of the stable region around the Lagrangian point L₄ depends on the mass parameter and the eccentricity of the secondary star. An additional goal of this work was to create a catalogue of all possible candidates, which could be useful for future observations to detect such objects.     

Detection of Earth-mass and super-Earth Trojan planets using transit timing variation method

We have carried out an extensive study of the possibility of the detection of Earth-mass and super-Earth Trojan planets using transit timing variation method with the Kepler space telescope. We have considered a system consisting of a transiting Jovian-type planet in a short period orbit, and determined the induced variations in its transit timing due to an Earth-mass/super-Earth Trojan planet. We mapped a large section of the phase space around the 1:1 mean-motion resonance and identified regions corresponding to several other mean-motion resonances where the orbit of the planet would be stable. We calculated transit timing variations (TTVs) for different values of the mass and orbital elements of the transiting and perturbing bodies as well as the mass of central star, and identified orbital configurations of these objects (ranges of their orbital elements and masses) for which the resulted TTVs would be within the range of the variations of the transit timing of Kepler’s planetary candidates. Results of our study indicate that in general, the amplitudes of the TTVs fall within the detectable range of timing precision obtained from the Kepler’s long-cadence data, and depending on the parameters of the system, their magnitudes may become as large as a few hours. The probability of detection is higher for super-Earth Trojans with slightly eccentric orbits around short-period Jovian-type planets with masses slightly smaller than Jupiter. We present the details of our study and discuss the implications of its results.     

On the Stability Regions of the Trojan Asteroids

The orbits of fictitious bodies around Jupiter’s stable equilibrium points L ₄ and L ₅ were integrated for a fine grid of initial conditions up to 100 million years. We checked the validity of three different dynamical models, namely the spatial, restricted three body problem, a model with Sun, Jupiter and Saturn and also the dynamical model with the Outer Solar System (Jupiter to Neptune). We determined the chaoticity of an orbit with the aid of the Lyapunov Characteristic Exponents (=LCE) and used also a method where the maximum eccentricity of an orbit achieved during the dynamical evolution was examined. The goal of this investigation was to determine the size of the regions of motion around the equilibrium points of Jupiter and to find out the dependance on the inclination of the Trojan’s orbit. Whereas for small inclinations (up to i=20°) the stable regions are almost equally large, for moderate inclinations the size shrinks quite rapidly and disappears completely for i>60^°. Additionally, we found a difference in the dynamics of orbits around L ₄ which – according to the LCE – seem to be more stable than the ones around L ₅.     

High-inclination planets and asteroids in multistellar systems

The Kozai mechanism often destabilizes high-inclination orbits. It couples changes in the eccentricity and inclination, and drives high inclination, circular orbits to low inclination, eccentric orbits. In a recent study of the dynamics of planetesimals in the quadruple star system HD 98800, there were significant numbers of stable particles in circumbinary polar orbits about the inner binary pair which are apparently able to evade the Kozai instability.
Here, we isolate this feature and investigate the dynamics through numerical and analytical models. The results show that the Kozai mechanism of the outer star is disrupted by a nodal libration induced by the inner binary pair on a shorter time-scale. By empirically modelling the period of the libration, a criteria for determining the high-inclination stability limits in general triple systems is derived. The nodal libration feature is interesting and, although affecting inclination and node only, shows many parallels to the Kozai mechanism. This raises the possibility that high-inclination planets and asteroids may be able to survive in multistellar systems.     

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.     

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.     

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.     

On the stability of Earth-like planets in multi-planet systems

We present a continuation of our numerical study on planetary systems with similar characteristics to the Solar System. This time we examine the influence of three giant planets on the motion of terrestrial-like planets in the habitable zone (HZ). Using the Jupiter–Saturn–Uranus configuration we create similar fictitious systems by varying Saturn’s semi-major axis from 8 to 11 AU and increasing its mass by factors of 2–30. The analysis of the different systems shows the following interesting results: (i) Using the masses of the Solar System for the three giant planets, our study indicates a maximum eccentricity (max-e) of nearly 0.3 for a test-planet placed at the position of Venus. Such a high eccentricity was already found in our previous study of Jupiter–Saturn systems. Perturbations associated with the secular frequency g ₅ are again responsible for this high eccentricity. (ii) An increase of the Saturn-mass causes stronger perturbations around the position of the Earth and in the outer HZ. The latter is certainly due to gravitational interaction between Saturn and Uranus. (iii) The Saturn-mass increased by a factor 5 or higher indicates high eccentricities for a test-planet placed at the position of Mars. So that a crossing of the Earth’ orbit might occur in some cases. Furthermore, we present the maximum eccentricity of a test-planet placed in the Earth’ orbit for all positions (from 8 to 11 AU) and masses (increased up to a factor of 30) of Saturn. It can be seen that already a double-mass Saturn moving in its actual orbit causes an increase of the eccentricity up to 0.2 of a test-planet placed at Earth’s position. A more massive Saturn orbiting the Sun outside the 5:2 mean motion resonance (a _S  ≥9.7 AU) increases the eccentricity of a test-planet up to 0.4.     

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.     

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.     

A SuperWASP search for additional transiting planets in 24 known systems

We present results from a search for additional transiting planets in 24 systems already known to contain a transiting planet. We model the transits due to the known planet in each system and subtract these models from light curves obtained with the SuperWASP (Wide Angle Search for Planets) survey instruments. These residual light curves are then searched for evidence of additional periodic transit events. Although we do not find any evidence for additional planets in any of the planetary systems studied, we are able to characterize our ability to find such planets by means of Monte Carlo simulations. Artificially generated transit signals corresponding to planets with a range of sizes and orbital periods were injected into the SuperWASP photometry and the resulting light curves searched for planets. As a result, the detection efficiency as a function of both the radius and orbital period of any second planet is calculated. We determine that there is a good (>50 per cent) chance of detecting additional, Saturn-sized planets in   P ∼  10 d orbits around planet-hosting stars that have several seasons of SuperWASP photometry. Additionally, we confirm previous evidence of the rotational stellar variability of WASP-10, and refine the period of rotation. We find that the period of the rotation is  11.91 ± 0.05  d, and the false alarm probability for this period is extremely low  (∼10⁻¹³)  .     

Formal Integrals and Nekhoroshev Stability in a Mapping Model for the Trojan Asteroids

A symplectic mapping model for the co-orbital motion (Sándor et al., 2002, Cel. Mech. Dyn. Astr. 84, 355) in the circular restricted three body problem is used to derive Nekhoroshev stability estimates for the Sun–Jupiter Trojans. Following a brief review of the analytical part of Nekhoroshev theory, a direct method is developed to construct formal integrals of motion in symplectic mappings without use of a normal form. Precise estimates are given for the region of effective stability based on the optimization of the size of the remainder of the formal series. The stability region found for t=10¹⁰ yrs corresponds to a libration amplitude D_p=10.6°. About 30% of asteroids with accurately known proper elements (Milani, 1993, Cel. Mech. Dyn. Astron. 57, 59), at low eccentricities and inclinations, are included within this region. This represents an improvement with respect to previous estimates given in the literature. The improvement is due partly to the choice of better variables, but also to the use of a mapping model, which is a simplification of the circular restricted three body problem.     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

Formation of the regular satellite systems and rings of the major planets

In this paper, we apply the ideas presented by one of us (Prentice, 1978a, b) for the development of the proto-solar cloud into a system of Laplacian rings to the development of the protoplanetary clouds which ultimately led to Jupiter, Saturn and Uranus. We show that if one accepts this scenario — especially the idea of supersonic turbulence in the proto-planetary clouds — one can satisfactorily explain, on the basis of fixing a single adjustable parameter, both the geometric precession of the orbital radii of the regular satellite systems of these three planets and the chemical composition and mass distribution of these satellites. We suggest that thermal stirring in the proto-planetary cloud in the vicinity of the surface of the planet may be responsible for the smaller masses of some of the inner satellites as well as for the formation of the rocky rings of Uranus. The icy rings of Saturn are suggested to be the product of condensation processes in a continuous gaseous disc within the Roche limit of the planet.     

Stability criteria in many-body systems

Celestial Mechanics and Dynamical Astronomy - In previous papers of this series the stability of hierarchical many-body dynamical systems has been considered in terms of parameters which are a...     

Stability criteria in many-body systems

The analytical stability criterion applicable to coplanar hierarchical three-body systems described in the first paper of this series, Walkeret al. (1980), is modified to give an exact representation ofHill-type stability in all such cases. The dependence of the stability on all orbital parameters (in the coplanar case) is taken into account. The criterion for stability is now dependant upon the participating masses, the elements of the initial osculating Keplerian orbits of the system (viz. the orbits ofm ₂ aboutm ₁ andm ₃ about the mass-centre of the (m ₁,m ₂) system) and the positions within these orbits.The behaviour of the stability of such systems is demonstrated (both analytically and numerically) with respect to certain of the parameters involved to consider effects not dealt with in the above-mentioned paper. In particular two interesting real cases of triple systems in the Solar System are discussed, namely Sun-Jupiter-Saturn and Earth-Moon-Sun. The results of the present paper are compared with those of past authors who considered the same systems.Finally some general features arising out of our analysis are discussed.     

Stability criteria in many-body systems

An expansion of the force function ofn-body dynamical systems, where the equations of motion are expressed in the Jacobian coordinate system, is shown to give rise naturally to a set of (n–1) (n–2) dimensionless parameters ^ki _li {i = 2,...,n;k = 2,...,i – 1 (i 3);l =i + 1,...,n (i n – 1)}, representative of the size of the disturbances on the Keplerian orbits of the various bodies. The expansion is particularized to the casen=3 which involves the consideration of only two parameters ²³ and ₃₂. Further, the work of Szebehely and Zare (1977) is reviewed briefly with reference to a sufficient condition for the stability of corotational coplanar three-body systems, in which two of the bodies form a binary system. This condition is sufficient in the sense that it precludes any possibility of an exchange of bodies, i.e. Hill type stability, however, it is not a necessary condition. These two approaches are then combined to yield regions of stability or instability in terms of the parameters ²³ and ₃₂ for any system of given masses and orbital characteristics (neglecting eccentricities and inclinations) with the following result: that there is a readily applicable rule to assess the likelihood of stability or instability of any given triple system in terms of ²³ and ₃₂.Treating a system ofn bodies as a set of disturbed three-body systems we use existing data from the solar system, known triple systems and numerical experiments in the many-body problem to plot a large number of triple systems in the ²³, ₃₂ plane and show the results agree well with the ²³, ₃₂ analysis above (eccentricities and inclinations as appropriate to most real systems being negligible). We further deal briefly with the extension of the criteria to many-body systems wheren>4, and discuss several interesting cases of dynamical systems.     

Stability criteria in many-body systems

The stability parameters developed and discussed in the first paper of this series (Walkeret al., 1980) are used to determine empirically, by means of numerical integration experiment, regions of stability for corotational, coplanar, hierarchical three-body systems. The initially circular case of these systems is studied: the components of the close binary are taken to move initially in circular orbits with respect to their common mass-centre, the third mass initially moving in a circular orbit with respect to the same mass-centre such that its orbit lies wholly outside those of the former two masses. The stability of these systems is then studied by reference to the empirical stability parameters and the initial ratio of the semi-major axes of the orbit of the close binary to that of the third mass about the binary's mass-centre, which is less than unity. For given values of the stability parameters it is determined how the stability of a system is affected by changes in the ratio of the semi-major axes. It is found that an upper limit to this ratio exists which determines the region of stability for such systems. It is also found possible, in the region of instability, to predict how unstable a system will be i.e. crudely speaking, the number of orbits it may be expected to execute before some gross instability sets in. The effect commensurabilities in mean motion have on the stability of these systems is also considered. It is generally found that these commensurabilities enhance the stability of these systems. The predictive powers of the method are then tested: using many test cases it is seen how accurately the stability or instability of a system may be predicted.