A Note on Second Species Solutions Generated from p–q Resonant Orbits

The first image under the flow of the restricted three-body problem of the p—q resonant strips — that appear in the study of the p—q resonant orbits — do not have, in general, intersection with the strip. In this paper we show some particular situations in which the above intersections exist for some very simple p—q resonant orbits which, at the same time, are periodic second species solutions.This revised version was published online in October 2005 with corrections to the Cover Date.     

Families of Periodic Orbits Emanating From Homoclinic Orbits in the Restricted Problem of Three Bodies

We describe and comment the results of a numerical exploration on the evolution of the families of periodic orbits associated with homoclinic orbits emanating from the equilateral equilibria of the restricted three body problem for values of the mass ratio larger than μ ₁. This exploration is, in some sense, a continuation of the work reported in Henrard [Celes. Mech. Dyn. Astr. 2002, 83, 291]. Indeed it shows how, for values of μ. larger than μ ₁, the Trojan web described there is transformed into families of periodic orbits associated with homoclinic orbits. Also we describe how families of periodic orbits associated with homoclinic orbits can attach (or detach) themselves to (or from) the best known families of symmetric periodic orbits. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Systematic Study of the Stability of Symmetric Periodic Orbits in the Planar, Circular, Restricted Three-Body Problem

We investigate symmetric periodic orbits in the framework of the planar, circular, restricted, three-body problem. Having fixed the mass of the primary equal to that of Jupiter, we determine the linear stability of a number of periodic orbits for different values of the eccentricity. A systematic study of internal resonances, with frequency p/q with 2p 9, 1 q 5 and 4/3 p/q 5, offers an overall picture of the stability character of inner orbits. For each resonance we compute the stability of the two possible periodic orbits. A similar analysis is performed for some external periodic orbits.Furthermore, we let the mass of the primary vary and we study the linear stability of the main resonances as a function of the eccentricity and of the mass of the primary. These results lead to interesting conclusions about the stability of exosolar planetary systems. In particular, we study the stability of Earth-like planets in the planetary systems HD168746, GI86, 47UMa,b and HD10697.     

Newtonian and relativistic periodic orbits around two fixed black holes

We compare families of simple periodic orbits of test particles in the Newtonian and relativistic problems of two fixed centers (black holes). The Newtonian problem is integrable, while the relativistic problem is highly non-integrable.The orbits are calculated on the meridian plane through the fixed centersM ₁ (atz=+1) andM ₂ (atz=–1) for energies smaller than the escape energyE=1. We use prolate spheroidal coordinates (, , =const) and also the variables =cosh and =–cos . The orbits are inside a curve of zero velocity (CZV). The Newtonian orbits are also limited by an ellipse and a hyperbola, or by two eillipses. There are 3 main types of periodic orbits (1) elliptic type (around both centers), (2) hyperbolic-type, and (3) resonant-type.The elliptic type orbits are stable in the Newtonian case and both stable and unstable in the relativistic case. From the stable orbits bifurcate double period orbits both symmetric and asymmetric with respect to thez-axis. There are also higher order bifurcations. The hyperbolic-type orbits are unstable. The Newtonian resonant orbits are defined by the ratiot _µ/t =n/m of oscillations along and during one period, and they are all marginally unstable. The corresponding relativistic orbits are stable, or unstable. The main families are figure eight orbits aroundM ₁, or aroundM ₂ (3/1 orbits); gamma, or inverse gamma orbits (4/2); higher resonant families 5/1,7/1,...,8/2,12/2,...;, more complicated orbits, like 5/3, and bifurcations from the above orbits. Satellite orbits aroundM ₁, orM ₂, and their bifurcations (e.g. double period) exist in the relativistic case but not in the Newtonian case. The characteristics of the various families are quite different in the Newtonian and the relativistic cases. The sizes of the orbits and their stabilities are also quite different in general. In the Appendix we study the various types of straight line orbits and prove that some subcases introduced by Charlier (1902) are impossible.     

Periodic orbits and bifurcations in the Sitnikov four-body problem

We study the existence, linear stability and bifurcations of what we call the Sitnikov family of straight line periodic orbits in the case of the restricted four-body problem, where the three equal mass primary bodies are rotating on a circle and the fourth (small body) is moving in the direction vertical to the center mass of the other three. In contrast to the restricted three-body Sitnikov problem, where the Sitnikov family has infinitely many stability intervals (hence infinitely many Sitnikov critical orbits), as the “family parameter” ż₀ varies within a finite interval (while z ₀ tends to infinity), in the four-body problem this family has only one stability interval and only twelve 3-dimensional (3D) families of symmetric periodic orbits exist which bifurcate from twelve corresponding critical Sitnikov periodic orbits. We also calculate the evolution of the characteristic curves of these 3D branch-families and determine their stability. More importantly, we study the phase space dynamics in the vicinity of these orbits in two ways: First, we use the SALI index to investigate the extent of bounded motion of the small particle off the z-axis along its interval of stable Sitnikov orbits, and secondly, through suitably chosen Poincaré maps, we chart the motion near one of the 3D families of plane-symmetric periodic orbits. Our study reveals in both cases a fascinating structure of ordered motion surrounded by “sticky” and chaotic orbits as well as orbits which rapidly escape to infinity.     

The elliptic restricted problem at the 3 : 1 resonance

Four 3 : 1 resonant families of periodic orbits of the planar elliptic restricted three-body problem, in the Sun-Jupiter-asteroid system, have been computed. These families bifurcate from known families of the circular problem, which are also presented. Two of them, I _c , II _c bifurcate from the unstable region of the family of periodic orbits of the first kind (circular orbits of the asteroid) and are unstable and the other two, I _e , II _e , from the stable resonant 3 : 1 family of periodic orbits of the second kind (elliptic orbits of the asteroid). One of them is stable and the other is unstable. All the families of periodic orbits of the circular and the elliptic problem are compared with the corresponding fixed points of the averaged model used by several authors. The coincidence is good for the fixed points of the circular averaged model and the two families of the fixed points of the elliptic model corresponding to the families I _c , II _c , but is poor for the families I _e , II _e . A simple correction term to the averaged Hamiltonian of the elliptic model is proposed in this latter case, which makes the coincidence good. This, in fact, is equivalent to the construction of a new dynamical system, very close to the original one, which is simple and whose phase space has all the basic features of the elliptic restricted three-body problem.     

Homoclinic and heteroclinic orbits in the photogravitational restricted three-body problem

We study numerically the asymptotic homoclinic and heteroclinic orbits around the hyperbolic Lyapunov periodic orbits which emanate from Euler's critical points L ₁ and L ₂, in the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these Lyapunov orbits, are also presented. Poincaré surface of sections of these manifolds on appropriate planes and several homoclinic and heteroclinic orbits for the gravitational case as well as for varying radiation factor q ₁, are displayed. Homoclinic-homoclinic and homoclinic-heteroclinic-homoclinic chains which link the interior with the exterior Hill's regions, are illustrated. We adopt the Sun-Jupiter system and assume that only the larger primary radiates. It is found that for small deviations of its value from the gravitational case (q ₁ = 1), the radiation pressure exerts a significant impact on the Hill's regions and on these asymptotic orbits.     

Periodic orbits in the general three-body problem and the relationship between them

We study the regions of finite motions in the vicinity of three simple stable periodic orbits in the general problem of three equal-mass bodies with a zero angular momentum. Their distinctive feature is that one of the moving bodies periodically passes through the center of mass of the triple system. We consider the dynamical evolution of plane nonrotating triple systems for which the initial conditions are specified in such a way that one of the bodies is located at the center of mass of the triple system. The initial conditions can then be specified by three parameters: the virial coefficient k and the two angles, φ₁ and φ₂, that characterize the orientation of the velocity vectors for the bodies. We scanned the region of variation in these parameters k∈(0, 1); φ₁, φ₂∈(0, π) at steps of δk=0.01; δφ₁=δφ₂=1° and identified the regions of finite motions surrounding the periodic orbits. These regions are isolated from one another in the space of parameters (k, φ₁, φ₂). There are bridges that correspond to unstable orbits with long lifetimes between the regions. During the evolution of these metastable systems, the phase trajectory can “stick” to the vicinity of one of the periodic orbits or move from one vicinity to another. The evolution of metastable systems ends with their breakup.     

Hamiltonian Stability of Spin–Orbit Resonances in Celestial Mechanics

The behaviour of ‘resonances’ in the spin-orbit coupling in celestial mechanics is investigated in a conservative setting. We consider a Hamiltonian nearly-integrable model describing an approximation of the spin-orbit interaction. The continuous system is reduced to a mapping by integrating the equations of motion through a symplectic algorithm. We study numerically the stability of periodic orbits associated to the above mapping by looking at the eigenvalues of the matrix of the linearized map over the full cycle of the periodic orbit. In particular, the value of the trace of the matrix is related to the stability character of the periodic orbit. We denote by ε_* (p/q) the value of the perturbing parameter at which a given elliptic periodic orbit with frequency p/q becomes unstable. A plot of the critical function ε_* (p/q) versus the frequency at different orbital eccentricities shows significant peaks at the synchronous resonance (for low eccentricities) and at the synchronous and 3:2 resonances (at higher eccentricities) in good agreement with astronomical observations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical exploration of the photogravitational restricted five-body problem

We study numerically the restricted five-body problem when some or all the primary bodies are sources of radiation. The allowed regions of motion as determined by the zero-velocity surface and corresponding equipotential curves, as well as the positions of the equilibrium points are given. We found that the number of the collinear equilibrium points of the problem depends on the mass parameter β and the radiation factors q _i , i=0,…,3. The stability of the equilibrium points are also studied. Critical masses associated with the number of the equilibrium points and their stability are given. The network of the families of simple symmetric periodic orbits, vertical critical periodic solutions and the corresponding bifurcation three-dimensional families when the mass parameter β and the radiation factors q _i vary are illustrated. Series, with respect to the mass (and to the radiation) parameter, of critical periodic orbits are calculated.     

The Web of Periodic Orbits at L 4

We describe and comment the results of a numerical exploration of the numerous natural families of periodic orbits associated with the L ₄ equilibrium point of the restricted problem of three bodies (and of course by symmetry those associated with the L ₅ equilibrium point). These families are organized in a very structured network or coweb and this structure evolves, when the mass ratio varies, in a very organized way.     

Bifurcation of Relative Equilibria in the Main Problem of Artificial Satellite Theory for a Prolate Body

In this paper, we study circular orbits of the J ₂ problem that are confined to constant-z planes. They correspond to fixed points of the dynamics in a meridian plane. It turns out that, in the case of a prolate body, such orbits can exist that are not equatorial and branch from the equatorial one through a saddle-center bifurcation. A closed-form parametrization of these branching solutions is given and the bifurcation is studied in detail. We show both theoretically and numerically that, close to the bifurcation point, quasi-periodic orbits are created, along with two families of reversible orbits that are homoclinic to each one of them.     

Motion close to the Hopf Bifurcation Of the vertical family of periodic orbits of L4

The paper deals with different kinds of invariant motions (periodic orbits, 2D and 3D invariant tori and invariant manifolds of periodic orbits) in order to analyze the Hamiltonian direct Hopf bifurcation that takes place close to the Lyapunov vertical family of periodic orbits of the triangular equilibrium point L₄ in the 3D restricted three-body problem (RTBP) for the mass parameter, μ greater than (and close to) μ_R (Routh’s mass parameter). Consequences of such bifurcation, concerning the confinement of the motion close to the hyperbolic orbits and the 3D nearby tori are also described.     

Vertical stability of periodic solutions around the triangular equilibrium points

The vertical stability character of the families of short and long period solutions around the triangular equilibrium points of the restricted three-body problem is examined. For three values of the mass parameter less than equal to the critical value of Routh (μ _R ) i.e. for μ = 0.000953875 (Sun-Jupiter), μ = 0.01215 (Earth-Moon) and μ = μ _R = 0.038521, it is found that all such solutions are vertically stable. For μ > (μ _R ) vertical stability is studied for a number of ‘limiting’ orbits extended to μ = 0.45. The last limiting orbit computed by Deprit for μ = 0.044 is continued to a family of periodic orbits into which the well known families of long and short period solutions merge. The stability characteristics of this family are also studied.     

Instabilities and bifurcations of the families of collision periodic orbits in the restricted three-body problem

By using Birkhoff's regularizing transformation, we study the evolution of some of the infinite j-k type families of collision periodic orbits with respect to the mass ratio μ as well as their stability and dynamical structure, in the planar restricted three-body problem. The μ-C characteristic curves of these families extend to the left of the μ-C diagram, to smaller values of μ and most of them go downwards, although some of them end by spiralling around the constant point S* (μ=0.47549, C=3) of the Bozis diagram (1970). Thus we know now the continuation of the families which go through collision periodic orbits of the Sun-Jupiter and Earth-Moon systems. We found new μ-C and x-C characteristic curves. Along each μ-C characteristic curve changes of stability to instability and vice versa and successive very small stable and very large unstable segments appear. Thus we found different types of bifurcations of families of collision periodic orbits. We found cases of infinite period doubling Feigenbaum bifurcations as well as bifurcations of new families of symmetric and non-symmetric collision periodic orbits of the same period. In general, all the families of collision periodic orbits are strongly unstable. Also, we found new x-C characteristic curves of j-type classes of symmetric periodic orbits generated from collision periodic orbits, for some given values of μ. As C varies along the μ-C or the x-C spiral characteristics, which approach their focal-terminating-point, infinite loops, one inside the other, surrounding the triangular points L₄ and L₅ are formed in their orbits. So, each terminating point corresponds to a collision asymptotic symmetric periodic orbit for the case of the μ-C curve or a non-collision asymptotic symmetric periodic orbit for the case of the x-C curve, that spiral into the points L₄ and L₅, with infinite period. All these are changes in the topology of the phase space and so in the dynamical properties of the restricted three-body problem.     

Symmetric and Asymmetric Librations in Planetary and Satellite Systems at the 2/1 Resonance

Families of asymmetric periodic orbits at the 2/1 resonance are computed for different mass ratios. The existence of the asymmetric families depends on the ratio of the planetary (or satellite) masses. As models we used the Io-Europa system of the satellites of Jupiter for the case m₁>m₂, the system HD82943 for the new masses, for the case m₁=m₂ and the same system HD82943 for the values of the masses m₁<m₂ given in previous work. In the case m₁≥ m₂ there is a family of asymmetric orbits that bifurcates from a family of symmetric periodic orbits, but there exist also an asymmetric family that is independent of the symmetric families. In the case m₁<m₂ all the asymmetric families are independent from the symmetric families. In many cases the asymmetry, as measured by and by the mean anomaly M of the outer planet when the inner planet is at perihelion, is very large. The stability of these asymmetric families has been studied and it is found that there exist large regions in phase space where we have stable asymmetric librations. It is also shown that the asymmetry is a stabilizing factor. A shift from asymmetry to symmetry, other elements being the same, may destabilize the system.     

Regular and irregular periodic orbits

We distinguish between regular orbits, that bifurcate from the main families of periodic orbits (those that exist also in the unperturbed case) and irregular periodic orbits, that are independent of the above. The genuine irregular families cannot be made to join the regular families by changing some parameters. We present evidence that all irregular families appear inside lobes formed by the asymptotic curves of the unstable periodic orbits. We study in particular a dynamical system of two degrees of freedom, that is symmetric with respect to the x-axis, and has also a triple resonance in its unperturbed form. The distribution of the periodic orbits (points on a Poincaré surface of section) shows some conspicuous lines composed of points of different multiplicities. The regular periodic orbits along these lines belong to Farey trees. But there are also lines composed mainly of irregular orbits. These are images of the x-axis in the map defined on the Poincaré surface of section. Higher order iterations of this map , close to the unstable triple periodic orbit, produce lines that are close to the asymptotic curves of this unstable orbit. The homoclinic tangle, formed by these asymptotic curves, contains many regular orbits, that were generated by bifurcation from the central orbit, but were trapped inside the tangle as the perturbation increased. We found some stable periodic orbits inside the homoclinic tangle, both regular and irregular. This proves that the homoclinic tangle is not completely chaotic, but contains gaps (islands of stability) filled with KAM curves.     

Periodic orbits and their bifurcations in a 3-D system

We study some simple periodic orbits and their bifurcations in the Hamiltonian . We give the forms of the orbits, the characteristics of the main families, and some existence diagrams and stability diagrams. The existence diagram of the family 1a contains regions that are stable (S), simply unstable (U), doubly unstable (DU) and complex unstable (). In the regionsS andU there are lines of equal rotation numberm/n. Along these lines we have bifurcations of families of periodic orbits of multiplicityn. When these lines reach the boundary of the complex unstable region, they are tangent to it. Inside the region there are linesm/n, along which the orbits 1a, describedn-times, are doubly unstable; however, along these lines there are no bifurcations ofn-ple periodic orbits. The families bifurcating from 1a exist only in certain regions of the parameter space (, ). The limiting lines of these regions join at particular points representing collisions of bifurcations. These collisions of bifurcations produce a nonuniqueness of the various families of periodic orbits. The complicated structure of the various bifurcations can be understood by constructing appropriate stability diagrams.     

Study of the transfer from the Earth to a halo orbit around the equilibrium pointL 1

The purpose of this paper is to study a transfer strategy from the vicinity of the Earth to a halo orbit around the equilibrium pointL ₁ of the Earth-Sun system. The study is done in the real solar system (we use the DE-118 JPL ephemeris in the simulations of motion) although some simplified models, such as the restricted three body problem (RTBP) and the bicircular problem, have been also used in order to clarify the geometrical aspects of the problem. The approach used in the paper makes use of the hyperbolic character of the halo orbits under consideration. The invariant stable manifold of these orbits enables the transfer to be achieved with, theoretically, only one manoeuvre: the one of insertion into the stable manifold. For the total v required, the figures obtained are similar to the ones given by the standard procedures of optimization.     

The structure of chaos in a potential without escapes

We study the structure of chaos in a simple Hamiltonian system that does no have an escape energy. This system has 5 main periodic orbits that are represented on the surface of section by the points (1)O(0,0), (2)C ₁,C ₂(±y _c, 0), (3)B ₁,B ₂(O,±1) and (4) the boundary . The periodic orbits (1) and (4) have infinite transitions from stability (S) to instability (U) and vice-versa; the transition values of are given by simple approximate formulae. At every transitionS U a set of 4 asymptotic curves is formed atO. For larger the size and the oscillations of these curves grow until they destroy the closed invariant curves that surroundO, and they intersect the asymptotic curves of the orbitsC ₁,C ₂ at infinite heteroclinic points. At every transitionU S these asymptotic curves are duplicated and they start at two unstable invariant points bifurcating fromO. At the transition itself the asymptotic curves fromO are tangent to each other. The areas of the lobes fromO increase with ; these lobes increase even afterO becomes stable again. The asymptotic curves of the unstable periodic orbits follow certain rules. Whenever there are heteroclinic points the asymptotic curves of one unstable orbit approach the asymptotic curves of another unstable orbit in a definite way. Finally we study the tangencies and the spirals formed by the asymptotic curves of the orbitsB ₁,B ₂. We find indications that the number of spiral rotations tends to infinity as . Therefore new tangencies between the asymptotic curves appear for arbitrarily large . As a consequence there are infinite new families of stable periodic orbits that appear for arbitrarily large .