The Hill stability of inclined small mass binary systems in three-body systems with special application to triple star systems, extrasolar planetary systems and Binary Kuiper Belt systems

The dynamical stability of a bound triple system composed of a small binary or minor planetary system moving on a orbit inclined to a central third body is discussed in terms of Hill stability for the full three-body problem. The situation arises in the determination of stability of triple star systems against disruption and component exchange and the determination of stability of extrasolar planetary systems and minor planetary systems against disruption, component exchange or capture. The Hill stability criterion is applied to triple star systems and extrasolar planetary systems, the Sun-Earth-Moon system and Kuiper Belt binary systems to determine the critical distances for stable orbits. It is found that increasing the inclination of the third body decreases the Hill regions of stability. Increasing the eccentricity of the binary also produces similar effects.These type of changes make exchange or disruption of the component masses more likely. Increasing the eccentricity of the binary orbit relative to the third body substantially decreases stability regions as the eccentricity reaches higher values. The Kuiper Belt binaries were found to be stable if they move on circular orbits. Taking into account the eccentricity, it is less clear that all the systems are stable.     

The Hill stability of a binary or planetary system during encounters with a third inclined body

The dynamical interaction of a binary or planetary system and a third body moving on a parabolic orbit inclined to the system is discussed in terms of Hill stability for the full three-body problem. The situation arises in binary star disruption and exchange, in extrasolar planetary system disruption, exchange and capture. It is found that increasing the inclination of the third body decreases the Hill regions of stability. This makes exchange or disruption of the component masses more likely as does increasing the eccentricity of the binary.
The stability criteria are applied to determine possible disruption and capture distances for currently known extrasolar planetary systems.     

The Hill stability of a binary or planetary system during encounters with a third inclined body moving on a hyperbolic orbit

The dynamical interaction of a binary or planetary system and a third body moving on a hyperbolic orbit inclined to the system is discussed in terms of Hill stability for the full three-body problem. The situation arises in binary star disruption and exchange and planetary system exchange or capture. It is found that increasing the inclination of the third body decreases the Hill regions of stability. Increasing the eccentricity of the third body also produces similar effects. These type of changes make exchange or disruption of the component masses more likely as also does increasing the eccentricity of the binary.

The critical distances and Hill stability ranges associated with the possible formation of roughly equal mass trans-Neptunian binaries from three-body interactions are determined for a range of secondary component masses.     

Limits on the Orbits of Possible Eccentric and Inclined Moons of Extrasolar Planets Orbiting Single Stars

Limits are placed on the range of orbits and masses of possible moons orbiting extrasolar planets which orbit single central stars. The Roche limiting radius determines how close the moon can approach the planet before tidal disruption occurs; while the Hill stability of the star–planet–moon system determines stable orbits of the moon around the planet. Here the full three-body Hill stability is derived for a system with the binary composed of the planet and moon moving on an inclined, elliptical orbit relative the central star. The approximation derived here in Eq. (17) assumes the binary mass is very small compared with the mass of the star and has not previously been applied to this problem and gives the criterion against disruption and component exchange in a closed form. This criterion was applied to transiting extrasolar planetary systems discovered since the last estimation of the critical separations (Donnison in Mon Not R Astron Soc 406:1918, 2010a) for a variety of planet/moon ratios including binary planets, with the moon moving on a circular orbit. The effects of eccentricity and inclination of the binary on the stability of the orbit of a moon is discussed and applied to the transiting extrasolar planets, assuming the same planet/moon ratios but with the moon moving with a variety of eccentricities and inclinations. For the non-zero values of the eccentricity of the moon, the critical separation distance decreased as the eccentricity increased in value. Similarly the critical separation decreased as the inclination increased. In both cases the changes though very small were significant.     

The Hill stability of low mass binaries in hierarchical triple systems

The Hill stability of the low mass binary system in the presence of a massive third body moving on a wider inclined orbit is investigated analytically. It is found that, in the case of the third body being on a nearly circular orbit, the region of Hill stability expands as the binary/third body mass ratio increases and the inclination (i) decreases. This i-dependence decreases very quickly with increasing eccentricity (e ₂) of the third body relative to the binary barycentre. In fact, if e ₂ is not extremely small, the Hill stable region can be approximately expressed in a closed form by setting i = 90°, and it contracts with increasing e ₂ as ${e_2^2}$ for sufficiently low mass binary. Our analytic results are then applied to the observed triple star systems and the Kuiper belt binaries.     

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.     

The stability of masses during three-body encounters

The dynamical interactio of a binary system and a third body not moving on a closed orbit arises in a large number of physical situations. The C²H condition for determining Hill stability of coplanar bound three-body systems is extended to cover situations where the outer body moves on a parabolic or hyperbolic orbit. Regions where such a body is stable against exchange or collision with other components of the system are determined for a number of important cases where closed solutions are possible.     

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.     

Stability of the coplanar planetary four-body system

We consider the coplanar planetary four-body problem, where three planets orbit a large star without the cross of their orbits. The system is stable if there is no exchange or cross of orbits. Starting from the Sundman inequality, the equation of the kinematical boundaries is derived. We discuss a reasonable situation, where two planets with known orbits are more massive than the third one. The boundaries of possible motions are controlled by the parameter c~2E. If the actual value of c~2E is less than or equal to a critical value(c~2 E)cr, then the regions of possible motions are bounded and therefore the system is stable.The criteria obtained in special cases are applied to the Solar System and the currently known extrasolar planetary systems. Our results are checked using N-body integrator.     

Ground-based detection of terrestrial extrasolar planets by photometry: the case for CM Draconis

The detection of extrasolar planets by measuring a photometric drop in the stellar brightness due to a planetary transit can be statistically improved by observing eclipsing binary systems and photometrically improved by observing small component systems. In particular the system CM Draconis, with two dM4 components, would allow the detection of extrasolar planets in the size range of Earth-to-Neptune requiring a ground-based photometric precision of about 0.08% to 1.1% (photometric precision of about 0.3% is routinely achievable with 1-meter class telescopes at the magnitude of CM Draconis, 11.07 inR-filter). In addition, the transit of extrasolar planets in a binary star system provides a unique, quasi-periodic signal that can be cross-correlated with the observational data to detect sub-noise signals. We examine the importance of making such observations to an understanding of the formation and evolution of terrestrial-type planets in main-sequence star systems. Terrestrial planets could have formed with substancially shorter periods in this lower luminosity system, for example, and might be expected to have accreted essentially in the binary orbital plane (however, non-coplanar planets may also eventually be detectable due to precession). We also report on a network of medium-sized telescopes at varying longitudes that have been organized to provide such constraints on terrestrial-planet formation processes and discuss the extention of near-term observations to other possible binary systems, as well. Finally, we discuss a more speculative, future observation that could be performed on the CM Draconis system that would be of exobiological as well as astrophysical interest.     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.     

A matched filter method for ground-based sub-noise detection of terrestrial extrasolar planets in eclipsing binaries: application to CM Draconis

The photometric detection of extrasolar planets by transits in eclipsing binary systems can be significantly improved by cross-correlating the observational light curves with synthetic models of possible planetary transit features, essentially a matched filter approach. We demonstrate the utility and application of this transit detection algorithm for ground-based detections of terrestrial-sized (Earth-to-Neptune radii) extrasolar planets in the dwarf M-star eclipsing binary system CM Draconis. Preliminary photometric observational data of this system demonstrate that the observational noise is well characterized as white and Gaussian at the observational time steps required for precision photometric measurements. Depending on planet formation scenarios, terrestrial-sized planets may form quite close to this low-luminosity system. We demonstrate, for example, that planets as small as 1.4 Earth radii with periods on the order of a few months in the CM Draconis system could be detected at the 99.9% confidence level in less than a year using 1-m class telescopes from the ground. This result contradicts commonly held assumptions limiting present ground-based efforts to, at best, detections of gas giant planets after several years of observation. This method can be readily extended to a number of other larger star systems with the utilization of larger telescopes and longer observing times. Its extension to spacecraft observations should also allow the determination of the presence of terrestrial-sized planets in nearly 100 other known eclipsing binary systems.     

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.     

Close encounters in young stellar clusters: implications for planetary systems in the solar neighbourhood

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N -body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction is not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the Solar System, the planets' orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet–planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extrasolar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of Solar System-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the Solar System around single stars, dynamical encounters in young stellar clusters may produce at least some of the extrasolar planetary systems observed in the solar neighbourhood.     

The role of cluster evolution in disrupting planetary systems and discs: the Kozai mechanism

We examine the effects of dynamical evolution in clusters on planetary systems or protoplanetary discs orbiting the components of binary stars. In particular, we look for evidence that the companions of host stars of planetary systems or discs could have their inclination angles raised from zero to between the threshold angles (39.23° and 140.77°) that can induce the Kozai mechanism. We find that up to 20 per cent of binary systems have their inclination angles increased to within the threshold range. Given that half of all extrasolar planets could be in binary systems, we suggest that up to 10 per cent of extrasolar planets could be affected by this mechanism.     

天体运动的轨道偏心率研究概述---从恒星系统到行星系统

偏心率是描述天体运动轨道的重要参数之一, 能够为揭示天体的动力学演化提供重要线索, 进而帮助理解天体形成与演化的过程及背后的物理机制. 随着天文观测技术的不断发展, 人们对于天体运动轨道的研究已经走出太阳系, 包含的系统也从大质量端的恒星系统延伸到了低质量端的行星系统. 聚焦天体轨道偏心率研究, 回顾了目前在恒星系统(包括主序恒星、褐矮星以及致密星)和行星系统(包括太阳系外巨行星以及``超级地球''、``亚海王星''等小质量系外行星)方面取得的进展, 总结了不同尺度结构下偏心率研究的一些共同之处和待解决的问题. 并结合当下和未来的相关天文观测设备和项目, 对未来天体轨道偏心率方面的研究工作进行了展望.     

Dynamics Of 2/1 Resonant Extrasolar Systems Application to HD82943 and GLIESE876

A complete study is made of the resonant motion of two planets revolving around a star, in the model of the general planar three body problem. The resonant motion corresponds to periodic motion of the two planets, in a rotating frame, and the position and stability properties of the periodic orbits determine the topology of the phase space and consequently play an important role in the evolution of the system. Several families of symmetric periodic orbits are computed numerically, for the 2/1 resonance, and for the masses of some observed extrasolar planetary systems. In this way we obtain a global view of all the possible stable configurations of a system of two planets. These define the regions of the phase space where a resonant extrasolar system could be trapped, if it had followed in the past a migration process.The factors that affect the stability of a resonant system are studied. For the same resonance and the same planetary masses, a large value of the eccentricities may stabilize the system, even in the case where the two planetary orbits intersect. The phase of the two planets (position at perihelion or aphelion when the star and the two planets are aligned) plays an important role, and the change of the phase, other things being the same, may destabilize the system. Also, the ratio of the planetary masses, for the same total mass of the two planets, plays an important role and the system, at some resonances and some phases, is destabilized when this ratio changes.The above results are applied to the observed extrasolar planetary systems HD 82943, Gliese 876 and also to some preliminary results of HD 160691. It is shown that the observed configurations are close to stable periodic motion.     

Resonances in Multiple Planetary Systems

Since the first extrasolar planet was discovered about 10 years ago, a major point of dynamical investigations was the determination of stable regions in extrasolar planetary systems where additional planets may exist. Using numerical methods, we investigate the dynamical stability in known multiple planetary systems (HD74156, HD38529, HD168443, HD169830) with special interest on the region between the two known planets and on the mean motion resonances inside this region. As a dynamical model we take the restricted 4-body problem containing the host star, the two planets and massless test-planets. For our numerical integrations, we used the Lie-integrator and additionally the Fast Lyapunov Indicators as a tool for detecting chaotic motion. We also investigated the inner resonances with the outer planet and the outer resonances with the inner planet with test-planets located inside the resonances.     

Stability of S-type Orbits in Binaries

This numerical investigation is concerned with the stability of planets moving around one component of a double star system. Since the discovery of four extra solar planets moving in such orbits, there is a growing interest of stability studies thereto. We determined the stable regions in the elliptic restricted three body problem, for the whole range of mass-ratios from 0.1 to 0.9, by means of the Fast Lyapunov Indicators. The computations have been carried out for eccentricities of the binary and of the planet in the range 0–0.5. Therefore, we present for the first time the variation of the stable regions when the initial eccentricity of the planet is increased. We have found a correlation between the reduction of the stable zones if the eccentricity of the planet or of the binary is increased — of course the latter one is more effective.This revised version was published online in October 2005 with corrections to the Cover Date.