Stability of planetary orbits in binary systems

The possible existence of stable orbits is investigated in binary systems using Hill's method. Analytical stability conditions are established for satellites, for inner planets and for outer planets, allowing arbitrary values for the mass-ratio of the binary.Presented at the Symposium Star Catalogues, Positional Astronomy and Celestial Mechanics, held in honor of Paul Herget at the U.S. Naval Observatory, Washington, November 30, 1978.     

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.     

Stability of resonant planetary orbits in binary stars

This paper contains a numerical study of the stability of resonant orbits in a planetary system consisting of two planets, moving under the gravitational attraction of a binary star. Its results are expected to provide us with useful information about real planetary systems and, at the same time, about periodic motions in the general four-body problem (G4) because the above system is a special case of G4 where two bodies have much larger masses than the masses of the other two (planets). The numerical results show that the main mechanism which generates instability is the destruction of the Jacobi integrals of the massless planets when their masses become nonzero and that resonances in the motion of planets do not imply, in general, instability. Considerable intervals of stable resonant orbits have been found. The above quantitative results are in agreement with the existing qualitative predictions     

Synthetic secular theories of the planetary orbits: Regular and chaotic behavior

Modern computer technology allows dynamical astronomers to investigate the long term stability of real systems as thoroughly as ever. However, the process is not straightforward and new problems need to be solved. This work deals with only one such problem: the construction-from the numerical integration- of a secular perturbation theory that is able to describe the dynamical behavior of the system. The discussion refers to the outer planets and is based on the knowledge acquired by the author during her participation in project LONGSTOP. A digital filter is used in order to reduce the output and eliminate short periodic terms. Filtering uncovers long term variations in the semimajor axes. From the filtered output a secular perturbation theory is constructed in the assumption that the solution is regular, as secular perturbation theories can only be constructed for regular solutions. If we succeed, this means that the solution is indeed regular for the computed span of time; if not-and this can be established in a rigorous way-it has to be concluded a posteriori that the solution is not regular. The LONGSTOP 1A and 1B integrations show well that as the timespan of the integration increases it is possible to detect the non-regular behavior of the solution. This happens in the eccentricity of Saturn at the 10^–4 level.     

Planetary nebulae in planetary and binary systems

The problem of the survival of a low-mass secondary orbiting a primary that becomes a planetary nebula is studied. The values of the mass of the primary are 1.0, 1.5, and 2.0M ; the values of the mass of the secondary 0.001M , 0.01M and 0.1M . The orbital decay and mass of the secondary due to accretion and gravitational drag in the common envelope are presented. The possible application of the results to V471 Tau, UU Sge, WZ Sge and the Sun-Jupiter system are discussed.     

On periodic orbits and resonance in extrasolar planetary systems

We study the evolution of an extrasolar planetary system with two planets, for planar motion, starting from an exact resonant periodic motion and increasing the deviation from the equilibrium solution. We keep the semimajor axes and the eccentricities of the two planets fixed and we change the initial conditions by rotating the orbit of the outer planet by Δω. In this way the resonance is preserved, but we deviate from the exact periodicity and there is a transition from order to chaos as the deviation increases. There are three different routes to chaos, as far as the evolution of (ω ₂ ? ω ₁) is concerned: (a) Libration → rotation → chaos, with intermittent transition from libration to rotation in between, (b) libration → chaos and (c) libration → intermittent interchange between libration and rotation → chaos. This indicates that resonant planetary systems where the angle (ω ₂ ? ω ₁) librates or rotates are not different, but are closely connected to the exact periodic motion.     

The influence of shapes of the bodies on the dynamics of planetary and star systems

The dynamics of planetary and star systems including perturbing forces due to the flattering and quadrupole distortion of the bodies is studied. The analytical model describing the perturbations which influence on the orbital motion of extrasolar planetary systems is presented. The calculations of the secular evolution of the mean orbital elements have shown that the effects related to the shape of the body are more important than the ones due to the quadrupole distortion.     

Tidal fluxes and evolution of binary star systems

Explicit analytical expressions for the tidal velocities in the stars which are components of the binary system are obtained, taking into account the viscosity of the stellar matter, the star's rotation, the eccentricity of the orbit and tilt of the axis of rotation to the orbital plane.Equations which determine the tidal evolution of the star's axial and orbital motions are obtained and analyzed. Numerical solutions of these equations are analyzed.     

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.     

Stability of outer planetary orbits around binary stars: A comparison of Hill's and Laplace's stability criteria

A comparison is made between the stability criteria of Hill and that of Laplace to determine the stability of outer planetary orbits encircling binary stars. The restricted, analytically determined results of Hill's method by Szebehely and co-workers and the general, numerically integrated results of Laplace's method by Graziani and Black are compared for varying values of the mass parameter =m ₂/(m ₁+m ₂). For 00.15, the closest orbit (lower limit of radius) an outer planet in a binary system can have and still remain stable is determined by Hill's stability criterion. For >0.15, the critical radius is determined by Laplace's stability criterion. It appears that the Graziani-Black stability criterion describes the critical orbit within a few percent for all values of .     

The Hill stability of inclined bound triple star and planetary systems

The dynamical stability of a triple system composed of a binary or planetary system and a bound third body moving on a orbit inclined to the system 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 planetary systems against disruption, component 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 binary also produces similar effects. These type of changes make exchange or disruption of the component masses more likely. Increasing the eccentricity of the third body initially increases the stability of the system then decreases stability as the eccentricity reaches higher values.The Hill stability criterion is applied to extrasolar planetary systems to determine the critical distances at which planets of the same mass as the observed extrasolar planet moving on a circular orbit could remain on a stable orbit. It was found that these distances were sufficiently short suggesting that the presence of further as yet unobserved stable extrasolar planets in observed systems was very likely.     

A search for planets in binary star systems: A new type of approach

A theoretical model, based on particular type of the restricted three-body problem, is here presented in order to demonstrate the existence of a possible planetary motion near the center of mass in a binary star system. The superposition principle is used, with the introduction of two fictitious negative masses in order to simulate the real two primary bodies system.     

Morphological and kinematic signatures of a binary central star in the planetary nebula Hu 2-1

We present H α , [N  ii ] and [O  iii ] ground-based and HST archive images, VLA–A 3.6-cm continuum and H92 α emission-line data and high-resolution long-slit [N  ii ] spectra of the planetary nebula Hu 2-1. A large number of structural components are identified in the nebula: an outer bipolar and an inner shell, two pairs of collimated bipolar structures at different directions, monopolar bow-shock-like structures, and an extended equatorial structure within a halo. The formation of Hu 2-1 appears to be dominated by anisotropic mass ejection during the late-AGB stage of the progenitor and by variable, 'precessing' collimated bipolar outflows during the protoplanetary nebula and/or early planetary nebula phases. Different observational results strongly support the existence of a binary central star in Hu 2-1, among them (1) the observed point-symmetry of the bipolar lobes and inner shell, and the departures from axial symmetry of the bipolar lobes, (2) the off-centre position of the central star, (3) the detection of mass ejection towards the equatorial plane, and (4) the presence of 'precessing' collimated outflows. In addition, (5) an analysis of the kinematics shows that the systemic velocity of the bipolar outflows does not coincide with the systemic velocity of the bipolar shell. We propose that this velocity difference is a direct evidence of orbital motion of the ejection source in a binary system. From a deduced orbital velocity of ∼10 km s⁻¹, a semimajor axis of ∼ 9–27 au and period of ∼ 25–80 yr are obtained, assuming a reasonable range of masses. These parameters are used to analyse the formation of Hu 2-1 within current scenarios of planetary nebulae with binary central stars.     

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

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

On the dynamics of extrasolar planetary systems under dissipation: Migration of planets

We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.