Bifurcation of central configuration in the 2N+1 body problem

The bifurcation of central configuration in the Newtonian N-body problem for any odd number N ≥ 7 is shown. We study a special case where 2n particles of mass m on the vertices of two different coplanar and concentric regular n-gons (rosette configuration) and an additional particle of mass m₀ at the center are governed by the gravitational law he 2n+1 body problem. This system is of two degrees of freedom and permits only one mass parameter μ =m ₀/m. This parameter μ controls the bifurcation. If n≥ 3, namely any odd N ≥ 7, then the number of central configurations is three when μ ≥ μ _c , and one when μ ≥ μ _c . By combining the results of the preceding studies and our main theorem, explicit examples of bifurcating central configuration are obtained for N ≤ 13, for any odd N ∈ [15,943], and for any N ≥ 945.     

EJECTION–COLLISION ORBITS AND INVARIANT PUNCTURED TORI IN A RESTRICTED FOUR‐BODY PROBLEM

We study the motion of an infinitesimal mass point under the gravitational action of three mass points of masses μ, 1–2μ and μ moving under Newton's gravitational law in circular periodic orbits around their center of masses. The three point masses form at any time a collinear central configuration. The body of mass 1–2μ is located at the center of mass. The paper has two main goals. First, to prove the existence of four transversal ejection–collision orbits, and second to show the existence of an uncountable number of invariant punctured tori. Both results are for a given large value of the Jacobi constant and for an arbitrary value of the mass parameter 0<μ≤1/2. This revised version was published online in July 2006 with corrections to the Cover Date.     

Chaotic Scattering in the Restricted Three‐Body Problem II. Small mass parameters

We study the scattering motion of the planar restricted three‐body problem for small mass parameters μ. We consider the symmetric periodic orbits of this system with μ = 0 that collide with the singularity together with the circular and parabolic solutions of the Kepler problem. These divide the parameter space in a natural way and characterize the main features of the scattering problem for small non‐vanishing μ. Indeed, continuation of these orbits yields the primitive periodic orbits of the system for small μ. For different regions of the parameter space, we present scattering functions and discuss the structure of the chaotic saddle. We show that for μ < μc and any Jacobi integral there exist departures from hyperbolicity due to regions of stable motion in phase space. Numerical bounds for μc are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body

We have studied periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body. We have determined the periodic orbits for different values of μ, h and A (h is energy constant, μ is mass ratio of the two primaries and A is an oblateness factor). These orbits have been determined by giving displacements along the tangent and normal to the mobile coordinates as defined by Karimov and Sokolsky (Celest. Mech. 46:335, 1989). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of oblateness by taking some fixed values of μ, A and h. As starters for our method, we use some known periodic orbits in the classical restricted three body problem.     

On the Two-body Problem with Slowly Decreasing Mass

We intend to present two approximate analytic solutions of the two-body problem with slowly decreasing mass which are obtained through the integration of the Hamilton equations. The law of mass variation used was m = -αn on the first case, and mi = αimin(i); i = 1,2 on the second which gives rise to a perturbed Keplerian problem dependent on one and two small parameters respectively. Practical applications are included. This revised version was published online in August 2006 with corrections to the Cover Date.     

Instability of the 2+2 body problem

The equations of motion of the 2+2 body problem (two interacting particles in the gravitational field of two much more massive primaries m₁ and m₂ in circular keplerian orbit) have an integral analogous to the Jacobi integral of the circular 2+1 body problem. We show here that with 2+2 bodies this integral does not give rise to Hill stability, i.e. to confinement for all time in a portion of the configuration space not allowing for some close approaches to occur. This is because all the level manifolds are connected and all exchanges of bodies between the regions surroundingm ₁,m ₂ and infinity do not contradict the conservation of the integral. However, it is worth stressing that some of these exchanges are physically meaningless, because they involve either unlimited extraction of potential energy from the binary formed by the small bodies (without taking into account their physical size) or significant mutual perturbations between the small masses without close approach, a process requiring, for the Sun-Jupiter-two asteroids system, timescales longer than the age of the Solar System.     

Families of periodic orbits in the restricted four-body problem

In this paper, families of simple symmetric and non-symmetric periodic orbits in the restricted four-body problem are presented. Three bodies of masses m ₁, m ₂ and m ₃ (primaries) lie always at the apices of an equilateral triangle, while each moves in circle about the center of mass of the system fixed at the origin of the coordinate system. A massless fourth body is moving under the Newtonian gravitational attraction of the primaries. The fourth body does not affect the motion of the three bodies. We investigate the evolution of these families and we study their linear stability in three cases, i.e. when the three primary bodies are equal, when two primaries are equal and finally when we have three unequal masses. Series, with respect to the mass m ₃, of critical periodic orbits as well as horizontal and vertical-critical periodic orbits of each family and in any case of the mass parameters are also calculated.     

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.     

Periodic orbits in the generalized perturbed restricted three-body problem

This paper studies the existence and stability of equilibrium points under the influence of small perturbations in the Coriolis and the centrifugal forces, together with the non-sphericity of the primaries. The problem is generalized in the sense that the bigger and smaller primaries are respectively triaxial and oblate spheroidal bodies. It is found that the locations of equilibrium points are affected by the non-sphericity of the bodies and the change in the centrifugal force. It is also seen that the triangular points are stable for 0<μ<μ _c and unstable for m_c ￡ m < \frac12\mu_{c}\le\mu <\frac{1}{2}, where μ _c is the critical mass parameter depending on the above perturbations, triaxiality and oblateness. It is further observed that collinear points remain unstable.     

The manifolds of families of 3D periodic orbits associated to Sitnikov motions in the restricted three-body problem

This paper deals with the Sitnikov family of straight-line motions of the circular restricted three-body problem, viewed as generator of families of three-dimensional periodic orbits. We study the linear stability of the family, determine several new critical orbits at which families of three dimensional periodic orbits of the same or double period bifurcate and present an extensive numerical exploration of the bifurcating families. In the case of the same period bifurcations, 44 families are determined. All these families are computed for equal as well as for nearly equal primaries (μ = 0.5, μ = 0.4995). Some of the bifurcating families are determined for all values of the mass parameter μ for which they exist. Examples of families of three dimensional periodic orbits bifurcating from the Sitnikov family at double period bifurcations are also given. These are the only families of three-dimensional periodic orbits presented in the paper which do not terminate with coplanar orbits and some of them contain stable parts. By contrast, all families bifurcating at single-period bifurcations consist entirely of unstable orbits and terminate with coplanar orbits.     

2D calculations of the collapse of magnetized gas cloud

All the families of planar symmetric simple-periodic orbits of the photogravitational restricted plane circular three-body problem, are determined numerically in the case when the primaries are of equal mass and radiate with equal radiation factors (q ₁=q₂=q). We obtain a global view of the possible patterns of periodic three-body motion while the full range of values of the common radiation factor is explored, from the gravitational case (q=1) down to near the critical value at which the triangular equilibria disappear by coalescing with the inner equilibrium pointL ₁ on the rotating axis of the primaries. It is found that for large deviations of its value from the gravitational case the radiation factorq can have a strong effect on the structure of the families.     

Stable Manifolds and Homoclinic Points Near Resonances in the Restricted Three-Body Problem

The restricted three-body problem describes the motion of a massless particle under the influence of two primaries of masses 1− μ and μ that circle each other with period equal to 2π. For small μ, a resonant periodic motion of the massless particle in the rotating frame can be described by relatively prime integers p and q, if its period around the heavier primary is approximately 2π p/q, and by its approximate eccentricity e. We give a method for the formal development of the stable and unstable manifolds associated with these resonant motions. We prove the validity of this formal development and the existence of homoclinic points in the resonant region. In the study of the Kirkwood gaps in the asteroid belt, the separatrices of the averaged equations of the restricted three-body problem are commonly used to derive analytical approximations to the boundaries of the resonances. We use the unaveraged equations to find values of asteroid eccentricity below which these approximations will not hold for the Kirkwood gaps with q/p equal to 2/1, 7/3, 5/2, 3/1, and 4/1. Another application is to the existence of asymmetric librations in the exterior resonances. We give values of asteroid eccentricity below which asymmetric librations will not exist for the 1/7, 1/6, 1/5, 1/4, 1/3, and 1/2 resonances for any μ however small. But if the eccentricity exceeds these thresholds, asymmetric librations will exist for μ small enough in the unaveraged restricted three-body problem.     

Motion and stability of triangular equilibrium points in elliptical restricted three body problem under the radiating primaries

This paper studies the motion and orbital stability of the infinitesimal mass in the vicinity of the equilateral (triangular) Lagrangian points of the elliptic restricted three body problem, considering photo gravitational effects of both the primaries. The stability of the triangular points is studied under the effects of radiating primaries around the binary system (Achird, Luyten, αCen AB, Kruger-60, Xi-Bootis); using simulation technique by drawing different curves of zero velocity.     

Periodic orbits in the photogravitational restricted problem with the smaller primary an oblate body

In this paper, we have studied periodic orbits generated by Lagrangian solutions of the restricted three body problem when more massive body is a source of radiation and the smaller primary is an oblate body. We have determined periodic orbits for fixed values of μ, σ and different values of p and h (μ mass ratio of the two primaries, σ oblate parameter, p radiation parameter and h energy constant). These orbits have been determined by giving displacements along the tangent and normal to the mobile co-ordinates as defined by Karimov and Sokolsky (in Celest. Mech. 46:335, 1989). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of radiation pressure on the periodic orbits by taking some fixed values of μ and σ.     

Periodic elliptic motions in a planar restricted (N+1)-body problem

LetN2 mass points (primaries) move on a collinear solution of relative equilibrium of theN-body problem; i.e. suitably fixed on a uniformly rotating straight line. Consider the motion of a massless particle in the gravitational field of these primaries with arbitrarily given masses. An existence proof for periodic solutions (i.e. closed trajectories in a rotating coordinate system) will be given, in which the particle performs nearly keplerian elliptic motions about (and close to) any one of the primaries.     

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.     

Effect of oblateness, perturbations, radiation and varying masses on the stability of equilibrium points in the restricted three-body problem

This paper investigates the combined effect of small perturbations ε,ε′ in the Coriolis and centrifugal forces, radiation pressure q _i , and changing oblateness of the primaries A _i (t) (i=1,2) on the stability of equilibrium points in the restricted three body problem in which the primaries is a supergiant eclipsing binary system which consists of a pair of bright oblate stars having the appearance of a giant peanut in space and their masses assumed to vary with time in the absence of reactive forces. The equations of motion are derived and the equilibrium points are obtained. For the autonomized system, it is seen that there are more than a pair of the triangular points as κ→∞; κ being the arbitrary sum of the masses of the primaries. In the case of the collinear points, two additional equilibrium points exist on the line joining the primaries when simultaneously κ+ε′<0 and both primaries are oblate, i.e., 0<α _i ?1. So there are five collinear equilibrium points in this case. Two non-planar equilibrium points exist for κ>1. Hence, there are at least nine equilibrium points of the system. The stability of these points is explored analytically and numerically. It is seen that the collinear and triangular points are stable with respect to certain conditions controlled by κ while the non-planar equilibrium points are unstable.     

The stability of inner collinear equilibrium points in the photogravitational elliptic restricted problem

The linear stability of the inner collinear equilibrium point of 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. The case of equal radiation factors of the two primaries is considered and the full range of values of the common radiation factor is explored, from the caseq ₁ =q ₂ =q = 1/8 at which the triangular equilibria disappear by coalescing on the rotating axis of the primaries transferring their stability to the collinear point, down toq = 0 at which value the stability regions in theµ - e plane disappear by shrinking down to zero size. 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. They evolve as in the case of the triangular equilibrium point considered in a previous paper. There exist values of the common radiation factor, in the range considered, for which the collinear equilibrium point is stable for the entire range of mass distribution among the primaries and for large eccentricities of their orbits.     

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.