Instability in the 3rd Dimension of the Extended Störmer Problem

The problem of motions of charged particles or charged grains of matter in the field of two magnetic dipoles in rotation about a common center are studied by means of computing families of periodic solutions and their stability. Ten such families are found and their stability study shows that only two out of the ten have short arcs of stable members, all with low inclination to the plane of the two dipoles, while the rest of the families consist exclusively of unstable solutions. Hence the impact of adding one dipole to the configuration of the Störmer problem, leads to the conclusion that in the modified system the three-dimensional motions are basically unstable unlike those of the one-dipole case. The model magnetic field used in this study is a first approximation of the Z3model of the magnetic field of Saturn as computed from the Pioneer 11 and Voyager 1 and 2 measurements.     

Search for and study of nearly periodic orbits in the plane problem of three equal-mass bodies

We analyze nearly periodic solutions in the plane problem of three equal-mass bodies by numerically simulating the dynamics of triple systems. We identify families of orbits in which all three points are on one straight line (syzygy) at the initial time. In this case, at fixed total energy of a triple system, the set of initial conditions is a bounded region in four-dimensional parameter space. We scan this region and identify sets of trajectories in which the coordinates and velocities of all bodies are close to their initial values at certain times (which are approximately multiples of the period). We classify the nearly periodic orbits by the structure of trajectory loops over one period. We have found the families of orbits generated by von Schubart’s stable periodic orbit revealed in the rectilinear three-body problem. We have also found families of hierarchical, nearly periodic trajectories with prograde and retrograde motions. In the orbits with prograde motions, the trajectory loops of two close bodies form looplike structures. The trajectories with retrograde motions are characterized by leafed structures. Orbits with central and axial symmetries are identified among the families found.     

Mass ratios of three bodies for stable low-inclination three-dimensional periodic motion

The intervals of possible stability, on the -axis, of the basic families of three-dimensional periodie motions of the restricted three-body problem (determined in an earlier paper) are extended into regions of the -m ₃ parameter space of the general three-body problem. Sample three-dimensional periodic motions corresponding to these regions are computed and tested for stability. Six regions, corresponding to the vertical-critical orbitsl1v, m1v,m2v, andilv, survive this preliminary stability test-therefore, emerging as the mass parameters regions allowing the simplest types of stable low inclination three-dimensional motion of three massive bodies.     

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.     

Numerical continuation of families of frozen orbits in the zonal problem of artificial satellite theory

In the zonal problem of a satellite around the Earth, we continue numerically natural families of periodic orbits with the polar component of the angular momentum as the parameter. We found three families; two of them are made of orbits with linear stability while the third one is made of unstable orbits. Except in a neighborhood of the critical inclination, the stable periodic (or frozen) orbits have very small eccentricities even for large inclinations.     

Restriction of motions in wide pairs in the Galactic field

The motions of the components of wide binary stars in the solar neighborhood in the regular Galactic gravitational field on time scales ～10¹⁰ yr have been studied numerically. The regions of restricted motions of the components in wide pairs have been found depending on the initial conditions: the magnitude of the relative velocity of the components, their mutual distance, and the inclination of the relative velocity vector to the Galactic plane. The size of the main part of the region of restricted motions is approximately equal to the tidal radius. Profound changes in the eccentricity of the binary orbit occur at inclinations close to 90°, which can lead to close approaches of the stars with a pericenter distance less than 1 AU. In the case of retrograde motions (the binary rotates in a direction opposite to the Galactic rotation), there is a region of restricted motions extending at least to 10 pc. Examples of the trajectories of relative motion of the stars and the change in osculating orbital elements are given for systems with restricted motions.     

Frozen orbits for satellites close to an Earth-like planet

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.     

Simple three-dimensional periodic orbits in a galactic-type potential

We study the families of simple periodic orbits in a three-dimensional system that represents the inner parts of a perturbed triaxial galaxy. The perturbations depend on two control parameters. We find the regions where each family is stable, simply unstable, doubly unstable, or complex unstable. the stable and simply unstable families produce other families by bifurcation. Several families reach a maximum (or minimum) perturbation and then are continued by other families. The bifurcations are direct or inverse. The transition from one type of bifurcation to the other is theoretically explained. Another important phenomenon is the splitting of one family into two, or the joining of two families into one. We do not have any complex instability in the limiting cases of two-dimensional motions (when one control parameter is zero).The two main families of periodic orbits are in most cases stable when the energy is smaller than the escape energy. Most high energy orbits are unstable. However, we found stable orbits even for energies about four times larger than the escape energy.     

Periodic attitudes and bifurcations of a rigid spacecraft in the second degree and order gravity field of a uniformly rotating asteroid

In this work, periodic attitudes and bifurcations of periodic families are investigated for a rigid spacecraft moving on a stationary orbit around a uniformly rotating asteroid. Under the second degree and order gravity field of an asteroid, the dynamical model of attitude motion is formulated by truncating the integrals of inertia of the spacecraft at the second order. In this dynamical system, the equilibrium attitude has zero Euler angles. The linearised equations of attitude motion are utilised to study the stability of equilibrium attitude. It is found that there are three fundamental types of periodic attitude motions around a stable equilibrium attitude point. We explicitly present the linear solutions around a stable equilibrium attitude, which can be used to provide the initial guesses for computing the true periodic attitudes in the complete model. By means of a numerical approach, three fundamental families of periodic attitudes are studied, and their characteristic curves, distribution of eigenvalues, stability curves and stability distributions are determined. Interestingly, along the characteristic curves of the fundamental families, some critical points are found to exist, and these points correspond to tangent and period-doubling bifurcations. By means of a numerical approach, the bifurcated families of periodic attitudes are identified. The natural and bifurcated families constitute networks of periodic attitude families.     

Motions in wide pairs at various galactocentric distances

The motions of the components of wide binary stars in the regular Galactic gravitational field on time scales ~10¹⁰ yr at various Galactocentric distances R ₀ have been studied numerically. Near the Galactic center, the influence of the bar has been taken into account. The regions of restricted motions of the components in wide pairs have been found depending on the initial conditions: the magnitude of the relative velocity of the components, their mutual distance, and the inclination of the relative velocity vector to the Galactic plane. The shape and sizes of these regions are shown to depend significantly on R ₀: the sizes of the region of initial conditions corresponding to restricted motions increase with R ₀. Profound changes in the eccentricity of the binary orbit occur at inclinations close to 90?, which can lead to close approaches of the stars with a pericenter distance less than 1 AU. For retrograde motions (the binary rotates in a direction opposite to the Galactic rotation), there are elongated branches of the region of restricted motions in all cases, which extend at least to 10 pc in some cases.     

A symplectic mapping model for the 3/4 exterior

We present a symplectic mapping model to study the evolution of a small body at the 3/4 exterior resonance with Neptune, for planar and for three dimensional motion. The mapping is based on the averaged Hamiltonian close to this resonance and is constructed in such a way that the topology of its phase space is similar to that of the Poincaré map of the elliptic restricted three-body problem. Using this model we study the evolution of a small object near the 3/4 resonance. Both chaotic and regular motions are found, and it is shown that the initial phase of the object plays an important role on the appearance of chaos. In the planar case, objects that are phase-protected from close encounters with Neptune have regular orbits even at eccentricities up to 0.44. On the other hand objects that are not phase protected show chaotic behaviour even at low eccentricities. The introduction of the inclination to our model affects the stable areas around the 3/4 mean motion resonance, which now become thinner and thinner and finally at is=10° the whole resonant region becomes chaotic. This may justify the absence of a large population of objects at this resonance.     

The problem of the critical inclination revisited

The behaviour of the argument of the pericentre is investigated for the orbit of an artificial satellite which is moving under the potential when the inclination of the orbit is close to thecritical value tan^?1 2. The theory is developed to first order and it is applicable to all values of the eccentricity, with the exception of those in the neighbourhood of zero and unity. Four principal types of behaviour are noted and these are illustrated in appropriate phase-plane diagrams. It is shown that the two types which exhibit double libration in the argument of the pericentre are restricted to a relatively small domain in the (a, e)-plane of possible motions. Moreover, it is demonstrated that for double libration to occur it is necessary, but not sufficient, that \(e > \sqrt 6/13\) . The ranges of values of the inclination for which libration of the pericentre is a possibility are given for the more important cases. The general results are applied to the specific case of artificial Earth satellites whose orbits are inclined to the equator at angles close to the value of the critical inclination.     

Resonant periodic orbits in the exoplanetary systems

The planetary dynamics of 4/3, 3/2, 5/2, 3/1 and 4/1 mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to 60^° of mutual planetary inclination, but in most families, the stability does not exceed 20^°–30^°, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the 4/3 and 5/2 resonance, respectively.     

The Sitnikov family and the associated families of 3D periodic orbits in the photogravitational RTBP with oblateness

We consider the photogravitational restricted three-body problem with oblateness and study the Sitnikov motions. The family of straight line oscillations exists only in the case where the primaries are of equal masses as in the classical Sitnikov problem and have the same oblateness coefficients and radiation factors. A perturbation method based on Floquet theory is applied in order to study the stability of the motion and critical orbits are determined numerically at which families of three-dimensional periodic orbits of the same or double period bifurcate. Many of these families are computed.     

Influence of periodic orbits on the formation of giant planetary systems

The late-stage formation of giant planetary systems is rich in interesting dynamical mechanisms. Previous simulations of three giant planets initially on quasi-circular and quasi-coplanar orbits in the gas disc have shown that highly mutually inclined configurations can be formed, despite the strong eccentricity and inclination damping exerted by the disc. Much attention has been directed to inclination-type resonance, asking for large eccentricities to be acquired during the migration of the planets. Here we show that inclination excitation is also present at small to moderate eccentricities in two-planet systems that have previously experienced an ejection or a merging and are close to resonant commensurabilities at the end of the gas phase. We perform a dynamical analysis of these planetary systems, guided by the computation of planar families of periodic orbits and the bifurcation of families of spatial periodic orbits. We show that inclination excitation at small to moderate eccentricities can be produced by (temporary) capture in inclination-type resonance and the possible proximity of the non-coplanar systems to spatial periodic orbits contributes to maintaining their mutual inclination over long periods of time.     

Analytical investigations of quasi-circular frozen orbits in the Martian gravity field

Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80 × 80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre is kept at 270° because J ₃ has the same sign as J ₂; while for a minority of ones with low altitude and low inclination, argument of pericentre can be kept at 90° because of the effect of the higher degree odd zonals; for the critical inclination cases, argument of pericentre can also be kept at 90°. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models.     

The Restricted Hill Four-Body Problem with Applications to the Earth–Moon–Sun System

Starting from the four-body problem a generalization of both the restricted three-body problem and the Hill three-body problem is derived. The model is time periodic and contains two parameters: the mass ratio ν of the restricted three-body problem and the period parameter m of the Hill Variation orbit. In the proper coordinate frames the restricted three-body problem is recovered as m → 0 and the classical Hill three-body problem is recovered as ν → 0. This model also predicts motions described by earlier researchers using specific models of the Earth–Moon–Sun system. An application of the current model to the motion of a spacecraft in the Sun perturbed Earth–Moon system is made using Hill's Variation orbit for the motion of the Earth–Moon system. The model is general enough to apply to the motion of an infinitesimal mass under the influence of any two primaries which orbit a larger mass. Using the model, numerical investigations of the structure of motions around the geometric position of the triangular Lagrange points are performed. Values of the parameter ν range in the neighborhood of the Earth–Moon value as the parameter m increases from 0 to 0.195 at which point the Hill Variation orbit becomes unstable. Two families of planar periodic orbits are studied in detail as the parameters m and ν vary. These families contain stable and unstable members in the plane and all have the out-of-plane stability. The stable and unstable manifolds of the unstable periodic orbits are computed and found to be trapped in a geometric area of phase space over long periods of time for ranges of the parameter values including the Earth–Moon–Sun system. This model is derived from the general four-body problem by rigorous application of the Hill and restricted approximations. The validity of the Hill approximation is discussed in light of the actual geometry of the Earth–Moon–Sun system. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the effect of the eccentricity of a planetary orbit on the stability of satellite orbits

The effect of the eccentricity of a planet’s orbit on the stability of the orbits of its satellites is studied. The model used is the elliptic Hill case of the planar restricted three-body problem. The linear stability of all the known families of periodic orbits of the problem is computed. No stable orbits are found, the majority of them possessing one or two pairs of real eigenvalues of the monodromy matrix, while a part of a family with complex instability is found. Two families of periodic orbits, bifurcating from the Lagrangian points L₁, L₂ of the corresponding circular case are found analytically. These orbits are very unstable and the determination of their stability coefficients is not accurate, so we compute the largest Liapunov exponent in their vicinity. In all cases these exponents are positive, indicating the existence of chaotic motions     

On Sitnikov-like motions generating new kinds of 3D periodic orbits in the R3BP with prolate primaries

The existence of new equilibrium points is established in the restricted three-body problem with equal prolate primaries. These are located on the Z-axis above and below the inner Eulerian equilibrium point L ₁ and give rise to a new type of straight-line periodic oscillations, different from the well known Sitnikov motions. Using the stability properties of these oscillations, bifurcation points are found at which new types of families of 3D periodic orbits branch out of the Z-axis consisting of orbits located entirely above or below the orbital plane of the primaries. Several of the bifurcating families are continued numerically and typical member orbits are illustrated.     

Long-Period Variations of the Eccentricity Vector Valid also for Near Circular Orbits around a Non-Spherical Body

This paper studies the long period variations of the eccentricity vector of the orbit of an artificial satellite, under the influence of the gravity field of a central body. We use modified orbital elements which are non-singular at zero eccentricity. We expand the long periodic part of the corresponding Lagrange equations as power series of the eccentricity. The coefficients characterizing the differential system depend on the zonal coefficients of the geopotential, and on initial semi-major axis, inclination, and eccentricity. The differential equations for the components of the eccentricity vector are then integrated analytically, with a definition of the period of the perigee based on the notion of “free eccentricity”, and which is also valid for circular orbits. The analytical solution is compared to a numerical integration. This study is a generalization of (Cook, Planet. Space Sci., 14, 1966): first, the coefficients involved in the differential equations depend on all zonal coefficients (and not only on the very first ones); second, our method applies to nearly circular orbits as well as to not too eccentric orbits. Except for the critical inclination, our solution is valid for all kinds of long period motions of the perigee, i.e., circulations or librations around an equilibrium point.