Existence of equilibrium points and their linear stability in the generalized photogravitational Chermnykh-like problem with power-law profile

We consider the modified restricted three body problem with power-law density profile of disk, which rotates around the center of mass of the system with perturbed mean motion. Using analytical and numerical methods, we have found equilibrium points and examined their linear stability. We have also found the zero velocity surface for the present model. In addition to five equilibrium points there exists a new equilibrium point on the line joining the two primaries. It is found that L ₁ and L ₃ are stable for some values of inner and outer radius of the disk while other collinear points are unstable, but L ₄ is conditionally stable for mass ratio less than that of Routh’s critical value. Lastly, we have studied the effects of radiation pressure, oblateness and mass of the disk on the motion and stability of equilibrium points.     

The equilibrium solutions of the restricted problem of three rigid bodies with variable mass

In this paper we consider the circular planar restricted problem of three rigid bodiesS _i(i=1, 2, 3), two of them are axisymmetric ellipsoids and a third bodyS ₃ is a spherical satellite with decreasing mass, under the gravitational forces. The effect of small perturbations in the Coriolis force and the centrifugal forces on the location of equilibrium points has been studied. It is found only in the case when the primaries have equal differences between their respective principal moments of inertial the pointsL ₄ andL ₅ form nearly equilateral tringles with the primaries. The equilibrium pointsL ₁,L ₂,L ₃ remain collinear an ies on the line joining the primaries.     

Asymptotic Orbits at the Triangular Equilibria in the Photogravitational Restricted Three-Body Problem

We study numerically the asymptotic homoclinic and heteroclinic orbits associated with the triangular equilibrium points L ₄ and L ₅, in the gravitational and the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these critical points, are also presented. Hundreds of asymptotic orbits for equal mass of the primaries and for various values of the radiation pressure are computed and the most interesting of them are illustrated. In the Copenhagen case, which the problem is symmetric with respect to the x- and y-axis, we found and present non-symmetric heteroclinic asymptotic orbits. So pairs of heteroclinic connections (from L ₄ to L ₅ and vice versa) form non-symmetric heteroclinic cycles. The termination orbits (a combination of two asymptotic orbits) of all the simple families of symmetric periodic orbits, in the Copenhagen case, are illustrated.     

Non-symmetric doubly-asymptotic orbits at the outer collinear equilibrium pointL 3

In the framework of the solar system case (with only the larger primary radiating) of the photogravitational restricted three-body problem we compute and present some non-symmetric asymptotic orbits connecting the outer collinear equilibrium pointL ₃ with the neighbourhood of one of the triangular equilibrium pointsL _{4, 5}. Such orbits have not been found previously in the restricted problem.     

Periodic motion around stable collinear equilibrium point in the photogravitational restricted problem of three bodies

In a binary system with both bodies being luminous, the inner collinear equilibrium pointL ₁ becomes stable for values of the mass ratio and radiation pressure parameters in a certain region. The kind of periodic motions aroundL ₁ is examined in this case. Second-order parametric expansions are given and the families of periodic orbits generated fromL ₁ are numerically determined for several sets of values of the parameters. Short- and long-period solutions are identified showing a similarity in the character of periodicity with that aroundL ₄. It is also found that the finite periodic solutions in the vicinity ofL ₁ are stable.     

Artificial halo orbits for low-thrust propulsion spacecraft

We consider periodic halo orbits about artificial equilibrium points (AEP) near to the Lagrange points L ₁ and L ₂ in the circular restricted three body problem, where the third body is a low-thrust propulsion spacecraft in the Sun–Earth system. Although such halo orbits about artificial equilibrium points can be generated using a solar sail, there are points inside L ₁ and beyond L ₂ where a solar sail cannot be placed, so low-thrust, such as solar electric propulsion, is the only option to generate artificial halo orbits around points inaccessible to a solar sail. Analytical and numerical halo orbits for such low-thrust propulsion systems are obtained by using the Lindstedt Poincaré and differential corrector method respectively. Both the period and minimum amplitude of halo orbits about artificial equilibrium points inside L ₁ decreases with an increase in low-thrust acceleration. The halo orbits about artificial equilibrium points beyond L ₂ in contrast show an increase in period with an increase in low-thrust acceleration. However, the minimum amplitude first increases and then decreases after the thrust acceleration exceeds 0.415 mm/s². Using a continuation method, we also find stable artificial halo orbits which can be sustained for long integration times and require a reasonably small low-thrust acceleration 0.0593 mm/s².     

Dynamical aspects of multi-round horseshoe-shaped homoclinic orbits in the RTBP

We consider the planar restricted three-body problem and the collinear equilibrium point L ₃, as an example of a center × saddle equilibrium point in a Hamiltonian with two degrees of freedom. We explore numerically the existence of symmetric and non-symmetric homoclinic orbits to L ₃, when varying the mass parameter μ. Concerning the symmetric homoclinic orbits (SHO), we study the multi-round, m-round, SHO for m ≥ 2. More precisely, given a transversal value of μ for which there is a 1-round SHO, say μ ₁, we show that for any m ≥ 2, there are countable sets of values of μ, tending to μ ₁, corresponding to m-round SHO. Some comments on related analytical results are also made.     

Control of chaos in the vicinity of the Earth‐Moon L5 Lagrangian point to keep a spacecraft in orbit

The concept of Space Manifold Dynamics is a new method of space research. We have applied it along with the basic idea of the method of Ott, Grebogi, and York (OGY method) to stabilize the motion of a spacecraft around the triangular Lagrange point L₅ of the Earth‐Moon system. We have determined the escape rate of the trajectories in the general three‐ and four‐body problem and estimated the average lifetime of the particles. Integrating the two models we mapped in detail the phase space around the L₅ point of the Earth‐Moon system. Using the phase space portrait our next goal was to apply a modified OGY method to keep a spacecraft close to the vicinity of L₅. We modified the equation of motions with the addition of a time dependent force to the motion of the spacecraft. In our orbit‐keeping procedure there are three free parameters: (i) the magnitude of the thrust, (ii) the start time, and (iii) the length of the control. Based on our numerical experiments we were able to determine possible values for these parameters and successfully apply a control phase to a spacecraft to keep it on orbit around L₅. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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 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.     

Numerical Exploration of Chermnykh’s Problem

We study the equilibrium points and the zero-velocity curves of Chermnykh’s problem when the angular velocity ω varies continuously and the value of the mass parameter is fixed. The planar symmetric simple-periodic orbits are determined numerically and they are presented for three values of the parameter ω. The stability of the periodic orbits of all the families is computed. Particularly, we explore the network of the families when the angular velocity has the critical value ω = 2√2 at which the triangular equilibria disappear by coalescing with the collinear equilibrium point L₁. The analytic determination of the initial conditions of the family which emanate from the Lagrangian libration point L₁ in this case, is given. Non-periodic orbits, as points on a surface of section, providing an outlook of the stability regions, chaotic and escape motions as well as multiple-periodic orbits, are also computed. Non-linear stability zones of the triangular Lagrangian points are computed numerically for the Earth–Moon and Sun–Jupiter mass distribution when the angular velocity varies.     

The effect of eccentricity and inclination on the motion near the Lagrangian points L4 and L5

We study the effect of eccentricity and inclination on small amplitude librations around the equilibrium points L ₄ and L ₅ in the circular restricted three-body problem. We show that the effective libration centres can be displaced appreciably from the equilateral configuration. The secular extrema of the eccentricity as a function of the argument of pericentre are shifted by ∼25 ° This revised version was published online in July 2006 with corrections to the Cover Date.     

Non-linear stability in the generalized restricted three-body problem

The non-linear stability of the triangular equilibrium point L ₄ in the generalized restricted three-body problem has been examined. The problem is generalized in the sense that the infinitesimal body and one of the primaries have been taken as oblate spheroids. It is found that the triangular equilibrium point is stable in the range of linear stability except for three mass ratios.     

Radiation Induced Dynamics of the Galactic Halo

We show results of numerical simulations of a three component plasma consisting of electrons, ions and dust with external gravitation and radiation fields. We perform simulation runs, starting from an analytic halo equilibrium, balancing pressure, gravitational, and radiative forces. Within these the equilibrium is perturbed by the radiation of a typical OB-star association. The perturbation has a total energy input of 10⁷ L _⊙ and a duration of 30 Myrs. After switching off the perturbation, the simulations are continued to further investigate the dynamics induced. We start with a self consistent one-fluid MHD model without background magnetic field and show for an asymmetric case that the system approaches a new equilibrium after switching on the perturbation. Later it relaxes into the starting configuration again, when the additional radiation is turned off. We then show, first by including a disk-parallel magnetic field and then by redoing the simulations with a full three-fluid code, the influence of magnetic fields and species separation on the plasma dynamics. With our computations we demonstrate that these features can be important for the explanation of the structures of galactic halos and large scale mass flows. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bifurcation points and intersections of families of periodic orbits in the three-dimensional restricted three-body problem

Intersections of families of three-dimensional periodic orbits which define bifurcation points are studied. The existence conditions for bifurcation points are discussed and an algorithm for the numerical continuation of such points is developed. Two sequences of bifurcation points are given concerning the family of periodic orbits which starts and terminates at the triangular equilibrium pointsL ₄,L ₅. On these sequences two trifurcation points are identified forµ = 0.124214 andµ = 0.399335. The caseµ = 0.5 is studied in particular and it is found that the space families originating at the equilibrium pointsL ₂,L ₃,L ₄,L ₅ terminate on the same planar orbitm _1v of the familym.     

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.     

Chaotic Diffusion And Effective Stability of Jupiter Trojans

It has recently been shown that Jupiter Trojans may exhibit chaotic behavior, a fact that has put in question their presumed long term stability. Previous numerical results suggest a slow dispersion of the Trojan swarms, but the extent of the ‘effective’ stability region in orbital elements space is still an open problem. In this paper, we tackle this problem by means of extensive numerical integrations. First, a set of 3,200 fictitious objects and 667 numbered Trojans is integrated for 4 Myrs and their Lyapunov time, T_L, is estimated. The ones following chaotic orbits are then integrated for 1 Gyr, or until they escape from the Trojan region. The results of these experiments are presented in the form of maps of T_Land the escape time, T_E, in the space of proper elements. An effective stability region for 1 Gyr is defined on these maps, in which chaotic orbits also exist. The distribution of the numbered Trojans follows closely the T_E=1 Gyr level curve, with 86% of the bodies lying inside and 14% outside the stability region. This result is confirmed by a 4.5 Gyr integration of the 246 chaotic numbered Trojans, which showed that 17% of the numbered Trojans are unstable over the age of the solar system. We show that the size distributions of the stable and unstable populations are nearly identical. Thus, the existence of unstable bodies should not be the result of a size-dependent transport mechanism but, rather, the result of chaotic diffusion. Finally, in the large chaotic region that surrounds the stability zone, a statistical correlation between T_LandT_E is found.     

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

Effect of radiation on the stability of equilibrium points in the binary stellar systems: RW-Monocerotis,Krüger 60

We have studied the stability of location of various equilibrium points of a passive micron size particle in the field of radiating binary stellar system within the framework of circular restricted three body problem. Influence of radial radiation pressure and Poynting-Robertson drag (PR-drag) on the equilibrium points and their stability in the binary stellar systems RW-Monocerotis and Krüger-60 has been studied. It is shown that both collinear and off axis equilibrium points are linearly unstable for increasing value of β ₁ (ratio of radiation to gravitational force of the massive component) in presence of PR-drag for the binary systems. Further we find that out of plane equilibrium points (L _i , i=6,7) may exists for range of values of β ₁>1 for these binary systems in the presence of PR-drag. Our linear stability analysis shows that the motion near the equilibrium points L _6,7 of the binary systems is unstable both in the absence and presence of PR-drag.     

Spacecraft trajectories to the L<Subscript>3</Subscript> point of the Sun–Earth three-body problem

Of the three collinear libration points of the Sun–Earth Circular Restricted Three-Body Problem (CR3BP), L₃ is that located opposite to the Earth with respect to the Sun and approximately at the same heliocentric distance. Whereas several space missions have been launched to the other two collinear equilibrium points, i.e., L₁ and L₂, taking advantage of their dynamical and geometrical characteristics, the region around L₃ is so far unexploited. This is essentially due to the severe communication limitations caused by the distant and permanent opposition to the Earth, and by the gravitational perturbations mainly induced by Jupiter and the close passages of Venus, whose effects are more important than those due to the Earth. However, the adoption of a suitable periodic orbit around L₃ to ensure the necessary communication links with the Earth, or the connection with one or more relay satellites located at L₄ or L₅, and the simultaneous design of an appropriate station keeping-strategy, would make it possible to perform valuable fundamental physics and astrophysics investigations from this location. Such an opportunity leads to the need of studying the ways to transfer a spacecraft (s/c) from the Earth’s vicinity to L₃. In this contribution, we investigate several trajectory design methods to accomplish such a transfer, i.e., various types of two-burn impulsive trajectories in a Sun-s/c two-body model, a patched conics strategy exploiting the gravity assist of the nearby planets, an approach based on traveling on invariant manifolds of periodic orbits in the Sun–Earth CR3BP, and finally a low-thrust transfer. We examine advantages and drawbacks, and we estimate the propellant budget and time of flight requirements of each.