A parametric study of stability and resonances around L 4 in the elliptic restricted three-body problem

The size of the stable region around the Lagrangian point L ₄ in the elliptic restricted three-body problem is determined by numerical integration as a function of the mass parameter and eccentricity of the primaries. The size distribution of the stable regions in the mass parameter-eccentricity plane shows minima at certain places that are identified with resonances between the librational frequencies of motions around L ₄. These are computed from an approximate analytical equation of Rabe relating the frequency, mass parameter and eccentricity. Solutions of this equation are determined numerically and the global behaviour of the frequencies depending on the mass parameter and eccentricity is shown and discussed. The minimum sizes of the stable regions around L ₄ change along the resonances and the relative strength of the resonances is analysed. Applications to possible Trojan exoplanets are indicated. Escape from L ₄ is also investigated.     

The three-dimensional motion of Trojan asteroids

The problem is considered within the framework of the elliptic restricted three-body problem. The asymptotic solution is derived by a three-variable expansion procedure. The variables of the expansion represent three time-scales of the asteroids: the revolution around the Sun, the libration around the triangular Lagrangian pointsL ₄,L ₅, and the motion of the perihelion. The solution is obtained completely in the first order and partly in the second order. The results are given in explicit form for the coordinates as functions of the true anomaly of Jupiter. As an example for the perturbations of the orbital elements the main perturbations of the eccentricity, the perihelion longitude and the longitude of the ascending node are given. Conditions for the libration of the perihelion are also discussed.     

The effect of solar radiation pressure on the Lagrangian points in the elliptic restricted three-body problem

In this paper the effect of solar radiation pressure on the location and stability of the five Lagrangian points is studied, within the frame of elliptic restricted three-body problem, where the primaries are the Sun and Jupiter acting on a particle of negligible mass. We found that the radiation pressure plays the rule of slightly reducing the effective mass of the Sun and changes the location of the Lagrangian points. New formulas for the location of the collinear libration points were derived. For large values of the force ratio β, we found that at β=0.12, the collinear point L₃ is stable and some families of periodic orbits can be drawn around it.     

Analysis of periodic orbits in the Saturn-Titan system using the method of Poincare section surfaces

We explore the periodic orbits and the regions of quasi-periodic motion around both the primaries in the Saturn-Titan system in the framework of planar circular restricted three-body problem. The location, nature and size of periodic and quasi-periodic orbits are studied using the numerical technique of Poincare surface of sections. The maximum amplitude of oscillations about the periodic orbits is determined and is used as a parameter to measure the degree of stability in the phase space for such orbits. It is found that the orbits around Saturn remain around it and their stability increases with the increase in the value of Jacobi constant C. The orbits around Titan move towards it with the increase in C. At C=3.1, the pericenter and apocenter are 358.2 and 358.5 km, respectively. No periodic or quasi-periodic orbits could be found by the present method around the collinear Lagrangian point L ₁ (0.9569373834…).     

Secondary resonances of co-orbital motions

The size distribution of the stability region around the Lagrangian point L ₄ is investigated in the elliptic restricted three-body problem as the function of the mass parameter and the orbital eccentricity of the primaries. It is shown that there are minimum zones in the size distribution of the stability regions, and these zones are connected with the secondary resonances between the frequencies of librational motions around L ₄. The results can be applied to hypothetical Trojan planets for predicting values of the mass parameter and the eccentricity for which such objects can be expected or their existence is less probable.     

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)     

The elliptical three-body problem. Triangular solution

The restricted three-body problem with eccentric orbit is reviewed and the positions of the triangular Lagrangian points (L₄, L₅) are determined. It is put in evidence the fact the fact L₄ and L₅ are situated at the corners of an isoscales triangle: AB = BC = 1 − e²/)1 + e cos ν )4/3 and AC = 1 − e²/)1 + e cos ν )     

The dynamics of inclined Neptune Trojans

We investigated the stable area for fictive Trojan asteroids around Neptune’s Lagrangean equilibrium points with respect to their semimajor axis and inclination. To get a first impression of the stability region we derived a symplectic mapping for the circular and the elliptic planar restricted three body problem. The dynamical model for the numerical integrations was the outer Solar system with the Sun and the planets Jupiter, Saturn, Uranus and Neptune. To understand the dynamics of the region around L ₄ and L ₅ for the Neptune Trojans we also used eight different dynamical models (from the elliptic problem to the full outer Solar system model with all giant planets) and compared the results with respect to the largeness and shape of the stable region. Their dependence on the initial inclinations (0° < i < 70°) of the Trojans’ orbits could be established for all the eight models and showed the primary influence of Uranus. In addition we could show that an asymmetry of the regions around L ₄ and L ₅ is just an artifact of the different initial conditions.     

Connection between Hill stability and weak stability in the elliptic restricted three-body problem

This paper provides a study on the connection between Hill stability and weak stability in the framework of the spatial elliptic restricted three-body problem. We determine a necessary condition for weak stability by giving an upper and a lower bound of qualitative measure of the Hill stability. The sufficient condition for weak stability and the symmetry of weak stable regions around the planets of the Solar System is also investigated.     

On the effect of the eccentricity of a planetary orbit on the stability of satellite orbits

The effect of the eccentricity of a planet’s orbit on the stability of the orbits of its satellites is studied. The model used is the elliptic Hill case of the planar restricted three-body problem. The linear stability of all the known families of periodic orbits of the problem is computed. No stable orbits are found, the majority of them possessing one or two pairs of real eigenvalues of the monodromy matrix, while a part of a family with complex instability is found. Two families of periodic orbits, bifurcating from the Lagrangian points L₁, L₂ of the corresponding circular case are found analytically. These orbits are very unstable and the determination of their stability coefficients is not accurate, so we compute the largest Liapunov exponent in their vicinity. In all cases these exponents are positive, indicating the existence of chaotic motions     

Infinite Feigenbaum sequences and spirals in the vicinity of the Lagrangian periodic solutions

We studied systematically cases of the families of non-symmetric periodic orbits in the planar restricted three-body problem. We took interesting information about the evolution, stability and termination of bifurcating families of various multiplicities. We found that the main families of simple non-symmetric periodic orbits present a similar dynamical structure and bifurcation pattern. As the Jacobi constant changes each branch of the characteristic of a main family spirals around a focal point-terminating point in x- at which the Jacobi constant is C _∞ = 3 and their periodic orbits terminate at the corotation (at the Lagrangian point L₄ or L₅). As the family approaches asymptotically its termination point infinite changes of stability to instability and vice versa occur along its characteristic. Thus, infinite bifurcation points appear and each one of them produces infinite inverse Feigenbaum sequences. That is, every bifurcating family of a Feigenbaum sequence produces the same phenomenon and so on. Therefore, infinite spiral characteristics appear and each one of them generates infinite new inner spirals and so on. Each member of these infinite sets of the spirals reproduces a basic bifurcation pattern. Therefore, we have in general large unstable regions that generate large chaotic regions near the corotation points L₄, L₅, which are unstable. As C varies along the spiral characteristic of every bifurcating family, which approaches its focal point, infinite loops, one inside the other, surrounding the unstable triangular points L₄ or L₅ are formed on their orbits. So, each terminating point corresponds to an asymptotic non-symmetric periodic orbit that spirals into the corotation points L₄, L₅ with infinite period. This is a new mechanism that produces very large degree of stochasticity. These conclusions help us to comprehend better the motions around the points L₄ and L₅ of Lagrange.     

The Phase Space Structure Around L4 in the Restricted Three-Body Problem

The phase space structure around L ₄ in the restricted three-body problem is investigated. The connection between the long period family emanating from L ₄ and the very complex structure of the stability region is shown by using the method of Poincarés surface of section. The zero initial velocity stability region around L ₄ is determined by using a method based on the calculation of finite-time Lyapunov characteristic numbers. It is shown that the boundary of the stability region in the configuration space is formed by orbits suffering slow chaotic diffusion.     

Stability of Trojan planets in multi-planetary systems

Today there are more than 340 extra-solar planets in about 270 extra-solar systems confirmed. Besides the observed planets there exists also the possibility of a Trojan planet moving in the same orbit as the Jupiter-like planet. In our investigation we take also into account the habitability of a Trojan planet and whether such a terrestrial planet stays in the habitable zone. Its stability was investigated for multi-planetary systems, where one of the detected giant planets moves partly or completely in the habitable zone. By using numerical computations, we studied the orbital behaviour up to 10⁷ years and determined the size of the stable regions around the Lagrangian equilibrium points for different dynamical models for fictitious Trojans. We also examined the interaction of the Trojan planets with a second or third giant planet, by varying its semimajor axis and eccentricity. We have found two systems (HD 155358 and HD 69830) that can host habitable Trojan planets. Another aim of this work was to determine the size of the stable region around the Lagrangian equilibrium points in the restricted three body problem for small mass ratios μ of the primaries μ ≤ 0.001 (e.g. Neptune mass of the secondary and smaller masses). We established a simple relation for the size depending on μ and the eccentricity.     

On the stability of triangular points in the elliptic R3BP under radiating and oblate primaries

This paper investigates the stability of triangular equilibrium points (L _4,5) in the elliptic restricted three-body problem (ER3BP), when both oblate primaries emit light energy simultaneously. The positions of the triangular points are seen to shift away from the line joining the primaries than in the classical case on account of the introduction of the eccentricity, semi-major axis, radiation and oblateness factors of both primaries. This is shown for the binary systems Achird, Luyten 726-8, Kruger 60, Alpha Centauri AB and Xi Bootis. We found that motion around these points is conditionally stable with respect to the parameters involved in the system dynamics. The region of stability increases and decreases with variability in eccentricity, oblateness and radiation pressures.     

On the non-linear stability of motions around L5 in the elliptic restricted problem of the three bodies

The non-linear stability of motions around L₅ in the elliptic restricted problem of the three bodies is investigated numerically with emphasis on the effect of the orbital eccentricity of the primaries on the shape of the established stability regions. It is shown that with increasing eccentricity, the width of these regions is decreasing.     

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.     

Comment on the paper “On the triangular libration points in photogravitational restricted three-body problem with variable mass” by Zhang, M.J., et al.

In a recent paper, published in Astrophys. Space Sci. (337:107, 2012) (hereafter paper ZZX) and entitled “On the triangular libration points in photogravitational restricted three-body problem with variable mass”, the authors study the location and stability of the generalized Lagrange libration points L ₄ and L ₅. However their study is flawed in two aspects. First they fail to write correctly the equations of motion of the variable mass problem. Second they attribute a variable mass to the third body of the restricted three-body model, a fact that is not compatible with the assumptions used in deriving the mathematical formulation of this model.     

Periodic orbits of the planetary type and their stability

A review is presented of periodic orbits of the planetary type in the general three-body problem and fourbody problem and the restricted circular and elliptic tnreebody problem. These correspond to planetary systems with one Sun and two or three planets (or a planet and its satellites), the motion of asteoids and also planetary systems with two Suns. The factors which affect the stability of the above configurations are studied in connection with resonance or additional perturbations. Finally, the correspondence of the periodic orbits in the restricted three-body problem with the fixed points obtained by the method of averaging or the method of surface of section is indicated.     

Mass ratios of three bodies for stable low-inclination three-dimensional periodic motion

The intervals of possible stability, on the -axis, of the basic families of three-dimensional periodie motions of the restricted three-body problem (determined in an earlier paper) are extended into regions of the -m ₃ parameter space of the general three-body problem. Sample three-dimensional periodic motions corresponding to these regions are computed and tested for stability. Six regions, corresponding to the vertical-critical orbitsl1v, m1v,m2v, andilv, survive this preliminary stability test-therefore, emerging as the mass parameters regions allowing the simplest types of stable low inclination three-dimensional motion of three massive bodies.     

Linear stability of the triangular equilibrium points in the photogravitational elliptic restricted problem,II

The linear stability of the triangular equilibrium points in the photogravitational elliptic restricted three-body problem is examined and the stability regions are determined in the space of the parameters of mass, eccentricity, and radiation pressure, in the case of equal radiation factors of the two primaries. The full range of values of the common radiation factor is explored, from the gravitational caseq ₁ =q ₂ =q = 1 down to the critical value ofq = 1/8 at which the triangular equilibria disappear by coalescing on the rotating axis of the primaries. It is found that radiation pressure exerts a significant influence on the stability regions. For certain intervals of radiation values these regions become qualitatively different from the gravitational case as well as the solar system case considered in Paper I. There exist values of the common radiation factor, in the range considered, for which the triangular equilibrium points are stable for the entire range of mass distribution among the primaries and for large eccentricities of their orbits.