Periodic orbits in a dynamical system with three degrees of freedom

We use the analytical method of Lindstedt to make an inventory of the regular families of periodic orbits and to obtain approximate analytical solutions in a three-dimensional harmonic oscillator with perturbing cubic terms. We compare these solutions to the results of numerical computations at a specific orbital resonance.     

On the families of periodic orbits which bifurcate from the circular Sitnikov motions

In this paper we deal with the circular Sitnikov problem as a subsystem of the three-dimensional circular restricted three-body problem. It has a first analytical part where by using elliptic functions we give the analytical expressions for the solutions of the circular Sitnikov problem and for the period function of its family of periodic orbits. We also analyze the qualitative and quantitative behavior of the period function. In the second numerical part, we study the linear stability of the family of periodic orbits of the Sitnikov problem, and of the families of periodic orbits of the three-dimensional circular restricted three-body problem which bifurcate from them; and we follow these bifurcated families until they end in families of periodic orbits of the planar circular restricted three-body problem. We compare our results with the previous ones of other authors on this problem. Finally, the characteristic curves of some bifurcated families obtained for the mass parameter close to 1/2 are also described.     

High-order solutions of invariant manifolds associated with libration point orbits in the elliptic restricted three-body system

High-order analytical solutions of invariant manifolds, associated with Lissajous and halo orbits in the elliptic restricted three-body problem (ERTBP), are constructed in this paper. The equations of motion of ERTBP in the pulsating synodic coordinate system have five equilibrium points, and the three collinear libration points as well as the associated center manifolds are unstable. In our calculation, the general solutions of the invariant manifolds associated with Lissajous and halo orbits around collinear libration points are expressed as power series of five parameters: the orbital eccentricity, two amplitudes corresponding to the hyperbolic manifolds, and two amplitudes corresponding to the center manifolds. The analytical solutions up to arbitrary order are constructed by means of Lindstedt–Poincaré method, and then the center and invariant manifolds, transit and non-transit trajectories in ERTBP are all parameterized. Since the circular restricted three-body problem (CRTBP) is a particular case of ERTBP when the eccentricity is zero, the general solutions constructed in this paper can be reduced to describe the dynamics around the collinear libration points in CRTBP naturally. In order to check the validity of the series expansions constructed, the practical convergence of the series expansions up to different orders is studied.     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

On quasi-periodic motions around the collinear libration points in the real Earth–Moon system

Due to various perturbations, the collinear libration points of the real Earth–Moon system are not equilibrium points anymore. Under the assumption that the Moon’s motion is quasi-periodic, special quasi-periodic orbits called dynamical substitutes exist. These dynamical substitutes replace the geometrical collinear libration points as time-varying equilibrium points. In the paper, the dynamical substitutes of the three collinear libration points in the real Earth–Moon system are computed. For the points L ₁ and L ₂, linearized motions around the dynamical substitutes are described, and the variational equations of the dynamical substitutes are reduced to a form with a near constant coefficient matrix. Then higher order analytical formulae of the central manifolds are constructed. Using these analytical solutions as initial seeds, Lissajous orbits and halo orbits are computed with numerical algorithms.     

Relative motion of satellites exploiting the super-integrability of Kepler’s problem

This paper builds upon the work of Palmer and Imre exploring the relative motion of satellites on neighbouring Keplerian orbits. We make use of a general geometrical setting from Hamiltonian systems theory to obtain analytical solutions of the variational Kepler equations in an Earth centred inertial coordinate frame in terms of the relevant conserved quantities: relative energy, relative angular momentum and the relative eccentricity vector. The paper extends the work on relative satellite motion by providing solutions about any elliptic, parabolic or hyperbolic reference trajectory, including the zero angular momentum case. The geometrical framework assists the design of complex formation flying trajectories. This is demonstrated by the construction of a tetrahedral formation, described through the relevant conserved quantities, for which the satellites are on highly eccentric orbits around the Sun to visit the Kuiper belt.     

On quasi-satellite periodic motion in asteroid and planetary dynamics

Applying the method of analytical continuation of periodic orbits, we study quasi-satellite motion in the framework of the three-body problem. In the simplest, yet not trivial model, namely the planar circular restricted problem, it is known that quasi-satellite motion is associated with a family of periodic solutions, called family f, which consists of 1:1 resonant retrograde orbits. In our study, we determine the critical orbits of family f that are continued both in the elliptic and in the spatial models and compute the corresponding families that are generated and consist the backbone of the quasi-satellite regime in the restricted model. Then, we show the continuation of these families in the general three-body problem, we verify and explain previous computations and show the existence of a new family of spatial orbits. The linear stability of periodic orbits is also studied. Stable periodic orbits unravel regimes of regular motion in phase space where 1:1 resonant angles librate. Such regimes, which exist even for high eccentricities and inclinations, may consist dynamical regions where long-lived asteroids or co-orbital exoplanets can be found.     

Second-species solutions in the circular and elliptic restricted three-body problem

An analytical proof of the existence of some kinds of periodic orbits of second species of Poincaré, both in the Circular and Elliptic Restricted three-body problem, is given for small values of the mass parameter. The proof uses the asymptotic approximations for the solutions and the matching theory developed by Breakwell and Perko. In the paper their results are extended to the Elliptic problem and applied to prove the existence of second-species solutions generated by rectilinear ellipses in the Circular problem and nearly-rectilinear ones in the Elliptic case.     

Note on secular perturbations between a retrograde body and a prograde body

Hirayama (1927) studied secular perturbations between a retrograde body and a prograde body by considering that the mean motion of the retrograde body is negative. In this paper we discuss the same problem by measuring angle variables from the departure point and keeping the mean motions positive for both the retrograde body and the prograde body, and compare the analytical solutions with numerically integrated orbits.     

Analytical theories of satellite motion constructed with a new package of formula manipulation and compared with numerical integration

Specialized to the Lie series based perturbation method of Kirchgraber and Stiefel (1978) a new computer algebra package called ANALYTOS has been developed for constructing analytical orbital theories either in noncanonical or canonical form. We present results on the (extended) Main Problem of orbital theory of artificial earth satellites and related issues. The order of the solutions achieved is generally one order higher than those known from literature. Moreover, the analytical orbits have been checked succesfully against precise numerical ephemerides. This revised version was published online in July 2006 with corrections to the Cover Date.     

Refined model for the evolution of distant satellite orbits

We consider a model that describes the evolution of distant satellite orbits and that refines the solution of the doubly averaged Hill problem. Generally speaking, such a refinement was performed previously by J. Kovalevsky and A.A. Orlov in terms of Zeipel’s method by constructing a solution of the third order with respect to the small parameter m, the ratio of the mean motions of the planet and the satellite. The analytical solution suggested here differs from the solutions obtained by these authors and is closest in form to the general solution of the doubly averaged problem (∼m ²). We have performed a qualitative analysis of the evolutionary equations and conditions for the intersection of satellite orbits with the surface of a spherical planet with a finite radius. Using the suggested solution, we have obtained improved analytical time dependences of the elements of evolving orbits for a number of distant satellites of giant planets compared to the solution of the doubly averaged Hill problem and, thus, achieved their better agreement with the results of our numerical integration of the rigorous equations of perturbed motion for satellites.     

Über periodische Bahnen im eingeschränkten Dreikörperproblem und in der Nähe der Kommensurabilitäten vom Typus (k + 1)/k

The authors, collaborating during the winter 1966/67, performed numerical investigations on form and character of periodic solutions of the restricted problem of three bodies in the Sun–Jupiter-System and in the neighbourhood of the resonances of the (k + 1)/k type. The principal part of these computations which were achieved using a fast IBM-computer refers to the cases 2/1 (Hecuba-gap) and 3/2 (Hilda-group). The results gave new insight into the genealogy of periodic orbits in the regions which do not contain solutions of Poincaré's first genre. Moreover, L. Carpenter tried to develop an analytical theory (here not explicitely described) which leads to convergent trigonometric series for the coordinates of periodic orbits with restricted eccentricities, using the method of general perturbations in rectangular coordinates recently published by P. Musen and himself.     

Analytical solutions around planar equilibria of three dipole problem

Proceeding with our investigation into the motion of a particle influenced by the electromagnetic field of three celestial bodies of a magnetic-dipole nature we give here for the first time the analytical expressions of periodic solutions around a planar equilibrium point. These relations are expansions of the planar equations of motion in series of second order power of a parameter in the vicinity of equilibria. The above analytical expressions of periodic solutions give the first members of the family of periodic orbits which emanate from a stable equilibrium point. The whole family can then be calculated using a predictor-corrector algorithm.     

Constructive analytical solution of the evolution hill problem

A new analytical solution of the system of differential equations describing secular perturbations and long-period solar perturbations of mean orbits of outer satellites of giant planets was obtained. As distinct from other solutions, the solution constructed using von Zeipel’s method approximately takes into account, in the secular part of the perturbing function, the totality of fourth order with respect to the small parameter m of the ratio of the mean motions of the primary planet and the satellite. This enables us to describe more accurately the evolution of satellite orbits with large apocentric distances, which in the course of evolution may exceed the halved radius of the Hill sphere of the planet with respect to the Sun. Among these are the orbits of the two outermost Neptunian satellites N10 (Psamathe) and N13 (Neso). For these satellites, the parameter m amounts to 0.152 and 0.165, respectively. Different from a purely analytical solution, the proposed solution requires preliminary calculations for each satellite. More precisely, in doing so, we need to construct some simple functions to approximate more complex ones. This is why we use the phrase “constructive analytical.” To illustrate the solution, we compare it with the results of the numerical integration of the strict motion equations of the satellites N10 and N13 over time intervals 5–15 thousand years.     

Periodic Orbits in Analytical Planar Galactic Potentials

We study the regular families of periodic orbits in an analytical planar galactic potential, using the method of Lindstedt. We obtain analytical expressions describing these orbits, validity of which is not limited to small amplitudes. We can delimit, in the space of the parameters, the domain of existence of each family of orbits. This revised version was published online in July 2006 with corrections to the Cover Date.     

The long period behavior of the orbits of the Galilean satellites of Jupiter

Some of the results of an investigation into the long period behavior of the orbits of the Galilean satellites of Jupiter are presented. Special purpose computer programs were used to perform all the algebraic manipulations and series expansions that are necessary to describe the mutual interactions among the satellites.The disturbing function was expanded as a Poisson series in the modified Keplerian elements referred to a Jovicentric coordinate system. The differential equations for the modified Keplerian elements were then formed, and all short period perturbations were removed using Kamel's perturbation method. Approximate analytical solutions for these differential equations are derived, and the general form of the solutions are given.     

Intermediary LEO propagation including higher order zonal harmonics

Two new intermediary orbits of the artificial satellite problem are proposed. The analytical solutions include higher order effects of the geopotential, and are obtained by means of a torsion transformation applied to the quasi-Keplerian system resulting after the elimination of the parallax simplification, for the first intermediary, and after the elimination of the parallax and perigee simplifications, for the second one. The new intermediaries perform notably well for low Earth orbits propagation, are free from special functions, and result advantageous, both in accuracy and efficiency, when compared to the standard Cowell integration of the \(J_2\) problem, thus providing appealing alternatives for onboard, short-term, orbit propagation under limited computational resources.     

Qualitative study of the planar isosceles three-body problem

We consider the particular case of the planar three body problem obtained when the masses form an isosceles triangle for all time. Various authors [1, 2, 12, 8, 9, 13, 10] have contributed in the knowledge of the triple collision and of several families of periodic orbits in this problem. We study the flow on a fixed level of negative energy. First we obtain a topological representation of the energy manifold including the triple collision and infinity as boundaries of that manifold. The existence of orbits connecting the triple collision and infinity gives some homoclinic and heteroclinic orbits. Using these orbits and the homothetic solutions of the problem we can characterize orbits which pass near triple collision and near infinity by pairs of sequences. One of the sequences describes the regions visited by the orbit, the other refers to the behaviour of the orbit between two consecutive passages by a suitable surface of section. This symbolic dynamics which has a topological character is given in an abstract form and after it is applied to the isosceles problem. We try to keep globality as far as possible. This strongly relies on the fact that the intersection of some invariant manifolds with an equatorial plane (v=0) have nice spiraling properties. This can be proved by analytical means in some local cases. Numerical simulations given in Appendix A make clear that these properties hold globally.     

Stability of axial orbits in analytical galactic potentials

We study the stability of axial orbits in analytical galactic potentials as a function of the energy of the orbit and the ellipticity of the potential. The problem is solved by an analytical method, the validity of which is not limited to small amplitudes. The lines of neutral stability divide the parameter space in regions corresponding to different organizations of the main families of orbits in the symmetry planes.     

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange’s planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting–Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.