Families of Asymmetric Periodic Orbits in Hill’s Problem of Three Bodies

We present five families of periodic solutions of Hill’s problem which are asymmetric with respect to the horizontal ξ axis. In one of these families, the orbits are symmetric with respect to the vertical η axis; in the four others, the orbits are without any symmetry. Each family consists of two branches, which are mirror images of each other with respect to the ξ axis. These two branches are joined at a maximum of Γ, where the family of asymmetric periodic solutions intersects a family of symmetric (with respect to the ξ axis) periodic solutions. Both branches can be continued into second species families for Γ → − ∞.     

Families of Asymmetric Periodic Orbits in the Restricted Three-body Problem

This paper studies the asymmetric solutions of the restricted planar problem of three bodies, two of which are finite, moving in circular orbits around their center of masses, while the third is infinitesimal. We explore, numerically, the families of asymmetric simple-periodic orbits which bifurcate from the basic families of symmetric periodic solutions f, g, h, i, l and m, as well as the asymmetric ones associated with the families c, a and b which emanate from the collinear equilibrium points L ₁, L ₂ and L ₃ correspondingly. The evolution of these asymmetric families covering the entire range of the mass parameter of the problem is presented. We found that some symmetric families have only one bifurcating asymmetric family, others have infinity number of asymmetric families associated with them and others have not branching asymmetric families at all, as the mass parameter varies. The network of the symmetric families and the branching asymmetric families from them when the primaries are equal, when the left primary body is three times bigger than the right one and for the Earth–Moon case, is presented. Minimum and maximum values of the mass parameter of the series of critical symmetric periodic orbits are given. In order to avoid the singularity due to binary collisions between the third body and one of the primaries, we regularize the equations of motion of the problem using the Levi-Civita transformations.     

Families of symmetric relative periodic orbits originating from the circular Euler solution in the isosceles three-body problem

We study symmetric relative periodic orbits in the isosceles three-body problem using theoretical and numerical approaches. We first prove that another family of symmetric relative periodic orbits is born from the circular Euler solution besides the elliptic Euler solutions. Previous studies also showed that there exist infinitely many families of symmetric relative periodic orbits which are born from heteroclinic connections between triple collisions as well as planar periodic orbits with binary collisions. We carry out numerical continuation analyses of symmetric relative periodic orbits, and observe abundant families of symmetric relative periodic orbits bifurcating from the two families born from the circular Euler solution. As the angular momentum tends to zero, many of the numerically observed families converge to heteroclinic connections between triple collisions or planar periodic orbits with binary collisions described in the previous results, while some of them converge to “previously unknown” periodic orbits in the planar problem.     

Planar Symmetric Periodic Orbits in Four Dipole Problem

We have extend Stormer’s problem considering four magnetic dipoles in motion trying to justify the phenomena of extreme “orderlines” such as the ones observed in the rings of Saturn; the aim is to account the strength of the Lorentz forces estimating that the Lorentz field, co-acting with the gravity field of the planet, will limit the motion of all charged particles and small size grains with surface charges inside a layer of about 200 m thickness as that which is observed in the rings of Saturn. For this purpose our interest feast in the motion of charged particles with neglected mass where only electromagnetic forces accounted in comparison to the weakness of the Newtonian fields. This study is particularly difficult because in the regions we investigate these motions there is enormous three dimensional instability. Following the Poincare’s hypothesis that periodic solutions are ‘dense’ in the set of all solutions in Hamiltonian systems we try to calculate many families of periodic solutions and to study their stability. In this work we prove that in this environment charged particles can trace planar symmetric periodic orbits. We discuss these orbits in details and we give their symplectic relations using the Hamiltonian formulation which is related to the symplectic matrix. We apply numerical procedures to find families of these orbits and to study their stability. Moreover we give the bifurcations of these families with families of planar asymmetric periodic orbits and families of three dimensional symmetric periodic orbits.     

Asymmetric Periodic Solutions in the Restricted Problem of Three Bodies

A systematic numerical exploration of the families of asymmetric periodic orbits of the restricted three-body problem when a) the primary bodies are equal and b) for the Earth-Moon mass ratio, is presented. Decades families of asymmetric periodic solutions were found and three of the simplest ones, in the first case, and ten of the second one are illustrated. All of these families consist of periodic orbits which are asymmetric with respect to x-axis while are simple symmetric periodic orbits with respect to y-axis (i.e. the orbit has only one perpendicular intersection at half period with y-axis). Many asymmetric periodic orbits, members of these families, are calculated and plotted. We studied the stability of all the asymmetric periodic orbits we found. These families consist, mainly, of unstable periodic solutions but there exist very small, with respect to x, intervals where these families have stable periodic orbits. We also found, using appropriate Poincaré surface of sections, that a relatively large region of phase space extended around all these stable asymmetric periodic orbits shows chaotic motion.     

Families of asymmetric periodic solutions of the restricted problem of three bodies for the Sun-Jupiter mass ratio and their relationship with the symmetric families

Message derived a method to detect bifurcations of a family of asymmetric periodic solutions from a family of symmetric periodic solutions in the restricted problem of three bodies for the limiting case when the second body has zero mass. This is used to examine several small integer commensurabilities. A total of 21 exterior and 21 interior small integer commensurabilities are examined and bifurcations (two in number) are found to exist only for exterior commensurabilities (q+1):1,q=1, 2,, 7. On investigating other commensurabilities of this form for values ofq up to 50 two bifurcations are still found to exist for each. The eccentricities of the two bifurcation orbits are given for eachq up to 20. For a Sun-Jupiter mass ratio the complete family of asymmetric periodic solutions associated withq=1, 2,..., 5, and the initial segments of the asymmetric family withq=6, 7,..., 12, have been numerically determined. The family associated withq=5 contains some unstable orbits but all orbits in the other four complete families are stable. The five complete families each begin and end on the same symmetric family. The network of asymmetric and symmetric families close to the commensurabilities (q+1):1,q=1, 2,..., 5 is discussed.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

Symmetric and asymmetric 3:1 resonant periodic orbits with an application to the 55Cnc extra-solar system

We study the dynamics of 3:1 resonant motion for planetary systems with two planets, based on the model of the general planar three body problem. The exact mean motion resonance corresponds to periodic motion (in a rotating frame) and the basic families of symmetric and asymmetric periodic orbits are computed. Four symmetric families bifurcate from the family of circular orbits of the two planets. Asymmetric families bifurcate from the symmetric families, at the critical points, where the stability character changes. There exist also asymmetric families that are independent of the above mentioned families. Bounded librations exist close to the stable periodic orbits. Therefore, such periodic orbits (symmetric or asymmetric) determine the possible stable configurations of a 3:1 resonant planetary system, even if the orbits of the two planets intersect. For the masses of the system 55Cnc most of the periodic orbits are unstable and they are associated with chaotic motion. There exist however stable symmetric and asymmetric orbits, corresponding to regular trajectories along which the critical angles librate. The 55Cnc extra-solar system is located in a stable domain of the phase space, centered at an asymmetric periodic orbit.     

Asymptotic orbits in the (<Emphasis Type="Italic">N</Emphasis>+1)-body ring problem

In this paper we study the asymptotic solutions of the (N+1)-body ring planar problem, N of which are finite and ν=N−1 are moving in circular orbits around their center of masses, while the Nth+1 body is infinitesimal. ν of the primaries have equal masses m and the Nth most-massive primary, with m ₀=β m, is located at the origin of the system. We found the invariant unstable and stable manifolds around hyperbolic Lyapunov periodic orbits, which emanate from the collinear equilibrium points L ₁ and L ₂. We construct numerically, from the intersection points of the appropriate Poincaré cuts, homoclinic symmetric asymptotic orbits around these Lyapunov periodic orbits. There are families of symmetric simple-periodic orbits which contain as terminal points asymptotic orbits which intersect the x-axis perpendicularly and tend asymptotically to equilibrium points of the problem spiraling into (and out of) these points. All these families, for a fixed value of the mass parameter β=2, are found and presented. The eighteen (more geometrically simple) families and the corresponding eighteen terminating homo- and heteroclinic symmetric asymptotic orbits are illustrated. The stability of these families is computed and also presented.     

Motion Around The Triangular Equilibrium Points Of The Restricted Three-Body Problem Under Angular Velocity Variation

We study numerically the asymmetric periodic orbits which emanate from the triangular equilibrium points of the restricted three-body problem under the assumption that the angular velocity ω varies and for the Sun–Jupiter mass distribution. The symmetric periodic orbits emanating from the collinear Lagrangian point L ₃, which are related to them, are also examined. The analytic determination of the initial conditions of the long- and short-period Trojan families around the equilibrium points, is given. The corresponding families were examined, for a combination of the mass ratio and the angular velocity (case of equal eigenfrequencies), and also for the critical value ω = 2

Numerical extrapolation of periodic solutions

Numerical procedures are established for the continuation of families of periodic solutions of non-integrable dynamical systems. They are based on the use of the previous known members of a family for non-linear prediction of the next member to be determined. Both symmetric and asymmetric periodic solutions are considered. The procedures are applied and compared in the case of the restricted three-body problem. They are shown to lead to considerable saving of computer time.     

Classes of Families of Generalized Periodic Solutions to the Beletsky Equation

For the equation describing plane oscillations and rotations of a satellite, we consider families of symmetric generalized periodic solutions with integral rotation number p. We give new confirmations of the hypothesis: there are only four classes of these families with topologically different structures, namely, the classes of families of periodic solutions with p≥ 1, p= 0, p=−1, and p≤−2. Besides, we demonstrate that the vertices of cusps of these families are placed on some analytical curves, and the same is true for the multiple intersections of these families with other families.     

Evolution of the General Solution of the Restricted Problem Covering Symmetric and Escape Solutions

The work presented in paper I (Papadakis, K.E., Goudas, C.L.: Astrophys. Space Sci. (2006)) is expanded here to cover the evolution of the approximate general solution of the restricted problem covering symmetric and escape solutions for values of μ in the interval [0, 0.5]. The work is purely numerical, although the available rich theoretical background permits the assertions that most of the theoretical issues related to the numerical treatment of the problem are known. The prime objective of this work is to apply the ‘Last Geometric Theorem of Poincaré’ (Birkhoff, G.D.: Trans. Amer. Math. Soc. 14, 14 (1913); Poincaré, H.: Rend. Cir. Mat. Palermo 33, 375 (1912)) and compute dense sets of axisymmetric periodic family curves covering the initial conditions space of bounded motions for a discrete set of values of the basic parameter μ spread along the entire interval of permissible values. The results obtained for each value of μ, tested for completeness, constitute an approximation of the general solution of the problem related to symmetric motions. The approximate general solution of the same problem related to asymmetric solutions, also computable by application of the same theorem (Poincaré-Birkhoff) is left for a future paper. A secondary objective is identification-computation of the compact space of escape motions of the problem also for selected values of the mass parameter μ. We first present the approximate general solution for the integrable case μ = 0 and then the approximate solution for the nonintegrable case μ = 10⁻³. We then proceed to presenting the approximate general solutions for the cases μ = 0.1, 0.2, 0.3, 0.4, and 0.5, in all cases building them in four phases, namely, presenting for each value of μ, first all family curves of symmetric periodic solutions that re-enter after 1 oscillation, then adding to it successively, the family curves that re-enter after 2 to 10 oscillations, after 11 to 30 oscillations, after 31 to 50 oscillations and, finally, after 51 to 100 oscillations. We identify in these solutions, considered as functions of the mass parameter μ, and at μ = 0 two failures of continuity, namely: 1. Integrals of motion, exempting the energy one, cease to exist for any infinitesimal positive value of μ. 2. Appearance of a split into two separate sub-domains in the originally (for μ = 0) unique space of bounded motions. The computed approximations of the general solution for all values of μ appear to fulfill the ‘completeness’ criterion inside properly selected sub-domains of the domain of bounded motions in the (x, C) plane, which means that these sub-domains are filled countably densely by periodic family curves, which form a laminar flow-line pattern. The family curves in this pattern may, or may not, be intersected by a ‘basic’ family curve segment of order from 1 up to 3. The isolated points generating asymptotic solutions resemble ‘sink’ points toward which dense sets of periodic family curves spiral. The points in the compact domain in the (x, C) plane resting outside the domain of bounded motions (μ = 0), including the gap between the two large sub-domains (μ > 0) created by the aforementioned split, generate escape motions. The gap between the two large sub-domains of bounded motions grows wider for growing μ. Also, a number of compact gaps that generate escape motions exist within the body of the two sub-domains of bounded motions. The approximate general solutions computed include symmetric, heteroclinic, asymptotic, collision and escape solutions, thus constituting one component of the full approximate general solution of the problem, the second and final component being that of asymmetric solutions.     

Numerical investigation of the manifold I qp in the Sun-Jupiter case of the restricted three-body problem

Message and Taylor (1978) have given values of the mean eccentricities and commen-surabilities which correspond to bifurcation orbits of families of symmetric periodic orbits with families of asymmetric periodic orbits in the limit as the mass ratio tends to zero. These bifurcations have been given in a way that they seem to be isolated and unrelated from the whole structure of the periodic orbits of the system.In this paper a numerical investigation of the horizontal stability of the family I and its branches reveals the above bifurcations orbits in the Sun-Jupiter case of the restricted three-body problem and associates these orbits with the whole structure of the system, giving extensive information on them.     

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.     

Periodic solutions of the restricted problem that are analytic continuations of periodic solutions of Hill's problem for small μ>0

This paper establishes the existence and first order perturbation approximation of an infinite number of one-parameter families of symmetric periodic solutions of the restricted three body problem that are analytic continuations of symmetric periodic solutions of Hill's problem for small values of the mass ratio μ>0.     

Symmetric and Nonsymmetric Periodic Orbits in the Exterior Mean Motion Resonances with Neptune

We present families of periodic orbits and their stability for the exterior mean motion resonances 1:2, 1:3 and 1:4 with Neptune in the framework of the planar circular restricted three-body problem. We found that in each resonance there exist two branches of symmetric elliptic periodic orbits with stable and unstable segments. Asymmetric periodic orbits bifurcate from the corresponding symmetric ones. Asymmetric periodic orbits are stable and the motion in their neighbourhood is a libration with respect to the resonant angle variable. In all the families of asymmetric periodic orbits the eccentricity extends to high values. Poincaré sections reveal the changes of the topology in phase space.     

Periodic orbits near infinity in the restrictedN-body problem

This paper shows that there exist two families of periodic solutions of the restrictedN-body problem which are close to large circular orbits of the Kepler problem. These solutions are shown to be of general elliptic type and hence are stable. If the restricted problem admits a symmetry, then there are symmetric periodic solutions which are close to large elliptic orbits of the Kepler problem.     

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.