On possible circumbinary configurations of the planetary systems of α Centauri and EZ Aquarii

Possible configurations of the planetary systems of the binary stars α Cen A–BandEZAqr A–C are analyzed. The P-type orbits—circumbinary ones, i.e., the orbits around both stars of the binary, are studied. The choice of these systems is dictated by the fact that α Cen is closest to us in the Galaxy, while EZ Aqr is the closest system whose circumbinary planets, as it turns out, may reside in the “habitability zone.” The analysis has been performed within the framework of the planar restricted three-body problem. The stability diagrams of circumbinary motion have been constructed: on representative sets of initial data (in the pericentric distance–eccentricity plane), we have computed the Lyapunov spectra of planetary motion and identified the domains of regular and chaotic motion through their statistical analysis. Based on present views of the dynamics and architecture of circumbinary planetary systems, we have determined the most probable planetary orbits to be at the centers of the main resonance cells, at the boundary of the dynamical chaos domain around the parent binary star, which allows the semimajor axes of the orbits to be predicted. In the case of EZ Aqr, the orbit of the circumbinary planet is near the habitability zone and, given that the boundary of this zone is uncertain, may belong to it.     

Stochasticity of Planetary Orbits in Double Star Systems

Cosmogonical theories as well as recent observations allow us to expect the existence of planets around many stars other than the Sun. On an other hand, double and multiple star systems are established to be more numerous than single stars (such as the Sun), at least in the solar neighborhood. We are then faced to the following dynamical problem: assuming that planets can form in a binary early environment (I do not deal here with), does long-term stability for planetary orbits exist in double star systems.Although preliminary studies were rather pessimistic about the possibility of existence of stable planetary orbits in double or multiple star systems, modern computation have shown that many such stable orbits do exist (but possible chaotic behavior), either around the binary as a whole (P-type) or around one component of the binary (S-type), this latter being explored here.The dynamical model is the elliptic plane restricted three-body problem; the phase space of initial conditions is systematically explored, and limits for stability have been established. Stable S-type planetary orbits are found up to distance of their "sun" of the order of half the periastron distance of the binary; moreover, among these stable orbits, nearly-circular ones exist up to distance of their "sun" of the order of one quarter the periastron distance of the binary; finally, among the nearly-circular stable orbits, several stay inside the "habitable zone", at least for two nearby binaries which components are nearly of solar type.Nevertheless, we know that chaos may destroy this stability after a long time (sometimes several millions years). It is therefore important to compute indicators of chaos for these stable planetary orbits to investigate their actual very long-term stability. Here we give an example of such a computation for more than a billion years.     

Ordered and chaotic Bohmian trajectories

We discuss the issue of ordered and chaotic trajectories in the Bohmian approach of Quantum Mechanics from points of view relevant to the methods of Celestial Mechanics. The Bohmian approach gives the same results as the orthodox (Copenhagen) approach, but it considers also underlying trajectories guided by the wave. The Bohmian trajectories are rather different from the corresponding classical trajectories. We give examples of a classical chaotic system that is ordered quantum-mechanically and of a classically ordered system that is mostly chaotic quantum mechanically. Then we consider quantum periodic orbits and ordered orbits, that can be represented by formal series of the “third integral” type, and we study their asymptotic properties leading to estimates of exponential stability. Such orbits do not approach the “nodal points” where the wavefunction ψ vanishes. On the other hand, when an orbit comes close to a nodal point, chaos is generated in the neighborhood of a hyperbolic point (called X-point). The generation of chaos is maximum when the X-point is close to the nodal point. Finally we remark that high order periodic orbits may behave as “effectively ordered” or “effectively chaotic” for long times before reaching the period.     

Precession of the Orbit of a Planet around Stars Revolving along: Low-Eccentricity Orbits in a Binary System: Analytical Solution

In our previous paper (hereafter, paper I) we presented analytical results on the non-planar motion of a planet around a binary star for the cases of the circular orbits of the components of the binary. We found that the orbital plane of the planet (the plane containing the “unperturbed” elliptical orbit of the planet), in addition to precessing about the angular momentum of the binary, undergoes simultaneously the precession within the orbital plane. We demonstrated that the analytically calculated frequency of this additional precession is not the same as the frequency of the precession of the orbital plane about the angular momentum of the binary, though the frequencies of both precessions are of the same order of magnitude. In the present paper we extend the analytical results from paper I by relaxing the assumption that the binary is circular – by allowing for a relatively small eccentricity ε of the stars orbits in the binary. We obtain an additional, ε-dependent term in the effective potential for the motion of the planet. By analytical calculations we demonstrate that in the particular case of the planar geometry (where the planetary orbit is in the plane of the stars orbits), it leads to an additional contribution to the frequency of the precession of the planetary orbit. We show that this additional, ε-dependent contribution to the precession frequency of the planetary orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency. Besides, we also obtain analytical results for another type of the non-planar configuration corresponding to the linear oscillatory motion of the planet along the axis of the symmetry of the circular orbits of the stars. We show that as the absolute value of the energy increases, the period of the oscillations decreases.     

Discrete Virial Theorem

We reexamine the classical virial theorem for bounded orbits of arbitrary autonomous Hamiltonian systems possessing both regular and chaotic orbits. New and useful forms of the virial theorem are obtained for natural Hamiltonian flows of arbitrary dimension. A discrete virial theorem is derived for invariant circles and periodic orbits of natural symplectic maps. A weak and a strong form of the virial theorem are proven for both flows and maps. While the Birkhoff Ergodic Theorem guarantees the existence of the relevant time averages for both regular and chaotic orbits, the convergence is very rapid for the former and extremely slow for the latter. This circumstance leads to a simple and efficient measure of chaoticity. The results are applied to several problems of current physical interest, including the Hénon–Heiles system, weak chaos in the standard map, and a 4D Froeschlé map.     

Onset of secular chaos in planetary systems: period doubling and strange attractors

As a result of resonance overlap, planetary systems can exhibit chaotic motion. Planetary chaos has been studied extensively in the Hamiltonian framework, however, the presence of chaotic motion in systems where dissipative effects are important, has not been thoroughly investigated. Here, we study the onset of stochastic motion in presence of dissipation, in the context of classical perturbation theory, and show that planetary systems approach chaos via a period-doubling route as dissipation is gradually reduced. Furthermore, we demonstrate that chaotic strange attractors can exist in mildly damped systems. The results presented here are of interest for understanding the early dynamical evolution of chaotic planetary systems, as they may have transitioned to chaos from a quasi-periodic state, dominated by dissipative interactions with the birth nebula.     

2/1 resonant periodic orbits in three dimensional planetary systems

We consider the general spatial three body problem and study the dynamics of planetary systems consisting of a star and two planets which evolve into 2/1 mean motion resonance and into inclined orbits. Our study is focused on the periodic orbits of the system given in a suitable rotating frame. The stability of periodic orbits characterize the evolution of any planetary system with initial conditions in their vicinity. Stable periodic orbits are associated with long term regular evolution, while unstable periodic orbits are surrounded by regions of chaotic motion. We compute many families of symmetric periodic orbits by applying two schemes of analytical continuation. In the first scheme, we start from the 2/1 (or 1/2) resonant periodic orbits of the restricted problem and in the second scheme, we start from vertical critical periodic orbits of the general planar problem. Most of the periodic orbits are unstable, but many stable periodic orbits have been, also, found with mutual inclination up to 50^?–60^?, which may be related with the existence of real planetary systems.     

Stochasticity of planetary orbits in double star systems. II

Cosmogonical theories as well as recent observations allow us to expect the existence of numerous exo-planets, including in binaries. Then arises the dynamical problem of stability for planetary orbits in double star systems. Modern computations have shown that many such stable orbits do exist, among which we consider orbits around one component of the binary (called S-type orbits). Within the framework of the elliptic plane restricted three-body problem, the phase space of initial conditions for fictitious S-type planetary orbits is systematically explored, and limits for stability had been previously established for four nearby binaries which components are nearly of solar type. Among stable orbits, found up to distance of their sun of the order of half the binarys periastron distance, nearly-circular ones exist for the three binaries (among the four) having a not too high orbital eccentricity. In the first part of the present paper, we compare these previous results with orbits around a 16 Cyg B-like binarys component with varied eccentricities, and we confirm the existence of stable nearly-circular S-type planetary orbits but for very high binarys eccentricity. It is well-known that chaos may destroy this stability after a very long time (several millions years or more). In a first paper, we had shown that a stable planetary orbit, although chaotic, could keep its stability for more than a billion years (confined chaos). Then, in the second part of the present paper, we investigate the chaotic behaviour of two sets of planetary orbits among the stable ones found around 16 Cyg B-like components in the first part, one set of strongly stable orbits and the other near the limit of stability. Our results show that the stability of the first set is not destroyed when the binarys eccentricity increases even to very high values (0.95), but that the stability of the second set is destroyed as soon as the eccentricity reaches the value 0.8.     

Resonant periodic orbits in the exoplanetary systems

The planetary dynamics of 4/3, 3/2, 5/2, 3/1 and 4/1 mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to 60^° of mutual planetary inclination, but in most families, the stability does not exceed 20^°–30^°, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the 4/3 and 5/2 resonance, respectively.     

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.     

Extrasolar planetary systems

The article deals with the occurrence of planetary systems in the Universe. In Section I, the terms “planet” and “planet-like objects” are defined. Two definitions proposed for the term “planetary system” are examined from the point of view (1) of the relation between planetary systems and binary and multiple star systems and (2) of planetary systems as abodes of intelligent beings. In Section II, the observational search for extrasolar planetary systems is described, as performable by earthbound optical telescopes, by space probes, by long baseline radio interferometry, and finally by inference from the reception of signals sent by intelligent beings in other worlds.In Section III we show that any planetary system must be preceded by a rotating disk of gas and dust around a central mass. Both observational evidence and theoretical reasons indicate the ease of formation of such disk structures in the cosmos. The time scale of collapse of a gaseous medium into a disk and that of the latter's dissipation are examined. This provides us with a new empirical approach and leads us to consider the problem of the frequency of occurrence of planetary systems to be ripe for scientific study. In Section IV, a brief review of theories of the formation of the solar system is given along with a proposed scheme for classification of these theories. In Section V, the evidence for magnetic activity in the early stages of stellar evolution is presented, as developed from six independent clues: the nuclear abundance of light elements, the behavior of flare stars, the intensities of H and K emission in stars, the nonthermal radiation of premain sequence stars, the properties of meteorites, and finally the existence of contact binaries. The magnetic braking theories of solar and stellar rotation are discussed in Section VI, thereby introducing the idea of formation of a rotating disk of gas and dust around stars in Section VII. From this disk a planetary system emerges.Section VIII gives an estimate for the frequency of occurrence of planetary systems in the Universe. It is based on the rotational behavior of main-sequence stars, and concludes that planetary systems have a far greater chance to appear around single main-sequence stars of spectral types later than F5 than around any other kind of star. The combined probability distribution of sizes and masses could be obtained. From physical considerations, it appears that sizes of planetary systems around stars of any given spectral type may not vary greatly from one to another.     

Some properties of a two-body system under the influence of the Galactic tidal field

In this paper the effect of the Galactic tidal field on a Sun–comet pair will be considered when the comet is situated in the Oort cloud and planetary perturbations can be neglected. First, two averaged models were created, one of which can be solved analytically in terms of Jacobi elliptic functions. In the latter system, switching between libration and circulation of the argument of perihelion is prohibited. The non-averaged equations of motion are integrated numerically in order to determine the regions of the ( e ,  i ) phase space in which chaotic orbits can be found, and an effort is made to explain why the chaotic orbits manifest in these regions only. It is evident that for moderate values of semimajor axis, a ∼50 000 au , chaotic orbits are found for ( e <0.15 , 40°≤ i ≤140°) as determined by integrating the evolution of the comet over a period of 10⁴ orbits. These regions of chaos increase in size with increasing semimajor axis. The typical e-folding times for these orbits range from around 600 Myr to 1 Gyr and thus are of little practical interest, as the time-scales for chaos arising from passing stars are much shorter. As a result of Galactic rotation, the chaotic regions in ( e ,  i ) phase space are not symmetric for prograde and retrograde orbits.     

Exploiting Unstable Periodic Orbits of a Chaotic Invariant Set for Spacecraft Control

The purpose of this work is to show that chaos control techniques (OGY, in special) can be used to efficiently keep a spacecraft around another body performing elaborate orbits. We consider a satellite and a spacecraft moving initially in coplanar and circular orbits, with slightly different radii, around a heavy central planet. The spacecraft, which is the inner body, has a slightly larger angular velocity than the satellite so that, after some time, they eventually go to a situation in which the distance between them becomes sufficiently small, so that they start to interact with one another. This situation is called as an encounter. In previous work we have shown that this scenario is a typical situation of a chaotic scattering for some well-defined range of parameters. Considering this scenario, we first show how it is possible to find the unstable periodic orbits that are located in the chaotic invariant set. From the set of unstable periodic orbits, we select the ones that can be combined to provide the desired elaborate orbit. Then, chaos control technique based on the OGY method is used to keep the spacecraft in the desired orbit. Finally, we analyze the results and make considerations regarding a realistic scenario of space exploration.     

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.     

The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincaré variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator Mean Exponential Growth factor of Nearby Orbits (MEGNO), a measure of the maximal Lyapunov exponent, that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best-fitting solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best-fitting solutions are found in dynamically active region of the phase space. The long-term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myr of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behaviour.     

FAST LYAPUNOV INDICATORS. APPLICATION TO ASTEROIDAL MOTION

We present a very simple and fast method to separate chaotic from regular orbits for non-integrable Hamiltonian systems. We use the standard map and the Hénon and Heiles potential as model problems and show that this method appears to be at least as sensitive as the frequency-analysis method. We also study the chaoticity of asteroidal motion. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

GLISSE: A GPU-optimized planetary system integrator with application to orbital stability calculations.

We present a GPU-accelerated numerical integrator specifically optimized for stability calculations of small bodies in planetary systems. Specifically, the integrator is designed for cases when large numbers of test particles (tens or hundreds of thousands) need to be followed for long durations (millions of orbits) to assess the orbital stability of their initially “close-encounter free” orbits. The GLISSE (Gpu Long-term Integrator for Solar System Evolution) code implements several optimizations to achieve a roughly factor of 100 speed increase over running the same code on a CPU. We explain how various hardware speed bottlenecks can be avoided by the careful code design, although some of the choices restrict the usage to specific types of application.As a first application, we study the long-term stability of small bodies initially on orbits between Uranus and Neptune. We map out in detail the small portion of the phase space in which small bodies can survive for 4.5 billion years of evolution; the ability to integrate large numbers of particles allow us to identify for the first time how instability-inducing mean-motion resonance overlaps sharply define the stable regions. As a second application, we map the boundaries of 4 Gyr stability for transneptunian objects in the 5:2 and 3:1 mean-motion resonances, demonstrating that long-term perturbations remove the initially stable Neptune-crossing members.