Analytical development of the lunisolar disturbing function and the critical inclination secular resonance

We provide a detailed derivation of the analytical expansion of the lunar and solar disturbing functions. Although there exist several papers on this topic, many derivations contain mistakes in the final expansion or rather (just) in the proof, thereby necessitating a recasting and correction of the original derivation. In this work, we provide a self-consistent and definite form of the lunisolar expansion. We start with Kaula’s expansion of the disturbing function in terms of the equatorial elements of both the perturbed and perturbing bodies. Then we give a detailed proof of Lane’s expansion, in which the elements of the Moon are referred to the ecliptic plane. Using this approach the inclination of the Moon becomes nearly constant, while the argument of perihelion, the longitude of the ascending node, and the mean anomaly vary linearly with time. We make a comparison between the different expansions and we profit from such discussion to point out some mistakes in the existing literature, which might compromise the correctness of the results. As an application, we analyze the long-term motion of the highly elliptical and critically-inclined Molniya orbits subject to quadrupolar gravitational interactions. The analytical expansions presented herein are very powerful with respect to dynamical studies based on Cartesian equations, because they quickly allow for a more holistic and intuitively understandable picture of the dynamics.     

Chaotic motion of comets in near-parabolic orbit: Mapping approaches

There exist many comets with near-parabolic orbits in the Solar System. Among various theories proposed to explain their origin, the Oort cloud hypothesis seems to be the most reasonable (Oort, 1950). The theory assumes that there is a cometary cloud at a distance 10³ – 10⁵ AU from the Sun and that perturbing forces from planets or stars make orbits of some of these comets become of near-parabolic type. Concerning the evolution of these orbits under planetary perturbations, we can raise the question: Will they stay in the Solar System forever or will they escape from it? This is an attractive dynamical problem. If we go ahead by directly solving the dynamical differential equations, we may encounter the difficulty of long-time computation. For the orbits of these comets are near-parabolic and their periods are too long to study on their long-term evolution. With mapping approaches the difficulty will be overcome. In another aspect, the study of this model has special meaning for chaotic dynamics. We know that in the neighbourhood of any separatrix i.e. the trajectory with zero frequency of the unperturbed motion of an Hamiltonian system, some chaotic motions have to be expected. Actually, the simplest example of separatrix is the parabolic trajectory of the two body problem which separates the bounded and unbounded motion. From this point of view, the dynamical study on near-parabolic motion is very important. Petrosky's elegant but more abstract deduction gives a Kepler mapping which describes the dynamics of the cometary motion (Petrosky, 1988). In this paper we derive a similar mapping directly and discuss its dynamical characters.     

A symplectic mapping for the synchronous spin-orbit problem

We derive a symplectic mapping model based on Hadjidemetriou’s method for the synchronous spin-orbit problem with and without the additional precession of the nodes. The mapping is derived from the averaged potential of the spin-orbit dynamical model and includes the main spin-orbit interactions, i.e. the non-zero obliquity and wobble motion of the rotating body. In addition the orbit of the perturbing body allows non-zero inclination and eccentricity. To obtain the equilibrium configuration we calculate the position and stability of the fixed points in the 1:1 spin-orbit resonance and relate them to the equilibria of the continuous system. We use the mapping equations to investigate the long-term stability close to the fixed point solutions of the mapping. We also apply the mapping method to the case of the moon Titan and validate the mapping approach by means of numerical integrations. The mapping model reproduces all the characteristics of Deprit’s model of free rotation as well as the dynamical features of Henrard’s averaged model of spin-orbit interaction with great precision.     

Where are the Uranus Trojans?

The area of stable motion for fictitious Trojan asteroids around Uranus’ equilateral equilibrium points is investigated with respect to the inclination of the asteroid’s orbit to determine the size of the regions and their shape. For this task we used the results of extensive numerical integrations of orbits for a grid of initial conditions around the points L ₄ and L ₅, and analyzed the stability of the individual orbits. Our basic dynamical model was the Outer Solar System (Jupiter, Saturn, Uranus and Neptune). We integrated the equations of motion of fictitious Trojans in the vicinity of the stable equilibrium points for selected orbits up to the age of the Solar system of 5 × 10⁹ years. One experiment has been undertaken for cuts through the Lagrange points for fixed values of the inclinations, while the semimajor axes were varied. The extension of the stable region with respect to the initial semimajor axis lies between 19.05 ≤ a ≤ 19.3 AU but depends on the initial inclination. In another run the inclination of the asteroids’ orbit was varied in the range 0° < i < 60° and the semimajor axes were fixed. It turned out that only four ‘windows’ of stable orbits survive: these are the orbits for the initial inclinations 0° < i < 7°, 9° < i < 13°, 31° < i < 36° and 38° < i < 50°. We postulate the existence of at least some Trojans around the Uranus Lagrange points for the stability window at small and also high inclinations.     

Finding how many isolating integrals of motion an orbit obeys

The correlation dimension, that is the dimension obtained by computing the correlation function of pairs of points of a trajectory in phase space, is a numerical technique introduced in the field of non-linear dynamics in order to compute the dimension of the manifold in which an orbit moves, without the need of knowing the actual equations of motion that give rise to the trajectory. This technique has been proposed in the past as a method to measure the dimension of stellar orbits in astronomical potentials, that is the number of isolating integrals of motion the orbits obey. Although the algorithm can in principle yield that number, some care has to be taken in order to obtain good results. We studied the relevant parameters of the technique, found their optimal values, and tested the validity of the method on a number of potentials previously studied in the literature, using the Smaller Alignment Index (SALI), Lyapunov exponents and spectral dynamics as gauges.     

Osculating and Superosculating Intermediate Orbits and their Applications

A method of construction of intermediate orbits for approximating the real motion of celestial bodies in the initial part of trajectory is proposed. The method is based on introducing a fictitious attracting centre with a time-variable gravitational parameter. The variation of thisparameter is assumed to obey the Eddington–Jeans mass-variationlaw. New classes of orbits having first-, second-, and third-order tangency to the perturbed trajectory at the initial instant of time are constructed. For planar motion, the tangency increases by one or two orders. The constructed intermediate orbits approximate the perturbed motion better than the osculating Keplerian orbit and analogous orbits of otherauthors. The applications of the orbits constructed in Encke's methodfor special perturbations and in the procedure for predicting themotion in which the perturbed trajectory is represented by a sequenceof short arcs of the intermediate orbits are suggested.The use of the constructed orbits is especially advantageous in the investigation of motion under the action of large perturbations.     

Sun-synchronous orbits near critical inclination

The long period dynamics of Sun-synchronous orbits near the critical inclination 116.6° are investigated. It is known that, at the critical inclination, the average perigee location is unchanged by Earth oblateness. For certain values of semimajor axis and eccentricity, orbit plane precession caused by Earth oblateness is synchronous with the mean orbital motion of the apparent Sun (a Sun-synchronism). Sun-synchronous orbits have been used extensively in meteorological and remote sensing satellite missions. Gravitational perturbations arising from an aspherical Earth, the Moon, and the Sun cause long period fluctuations in the mean argument of perigee, eccentricity, inclination, and ascending node. Double resonance occurs because slow oscillations in the perigee and Sun-referenced ascending node are coupled through the solar gravity gradient. It is shown that the total number and infinitesimal stability of equilibrium solutions can change abruptly over the Sun-synchronous range of semimajor axis values (1.54 to 1.70 Earth radii). The effect of direct solar radiation pressure upon certain stable equilibria is investigated.     

A numerical-analytical method for studying the orbital evolution of distant planetary satellites

We describe an approximate numerical-analytical method for calculating the perturbations of the elements of distant satellite orbits. The model for the motion of a distant satellite includes the solar attraction and the eccentricity and ecliptic inclination of the orbit of the central planet. In addition, we take into account the variations in planetary orbital elements with time due to secular perturbations. Our work is based on Zeipel’s method for constructing the canonical transformations that relate osculating satellite orbital elements to the mean ones. The corresponding transformation of the Hamiltonian is used to construct an evolution system of equations for mean elements. The numerical solution of this system free from rapidly oscillating functions and the inverse transformation from the mean to osculating elements allows the evolution of distant satellite orbits to be studied on long time scales on the order of several hundred or thousand satellite orbital periods.     

Superosculating intermediate orbits: Fourth-and fifth-order tangencies to trajectories of perturbed motion

The theory of superosculating intermediate orbits previously suggested by the author is developed. A new class of orbits with a fourth-order tangency to the actual trajectory of a celestial body at the initial time is constructed. Orbits with a fifth-order tangency have been constructed for the first time. The motion in the constructed orbits is represented as a combination of two motions: the motion of a fictitious attracting center with a variable mass and the motion relative to this center. The first motion is generally parabolic, while the second motion is described by the equations of the Gylden—Mestschersky problem. The variation in the mass of the fictitious center obeys Mestschersky’s first and combined laws. The new orbits represent more accurately the actual motion in the initial segment of the trajectory than an osculating Keplerian orbit and other existing analogues. Encke’s generalized methods of special perturbations in which the constructed intermediate orbits are used as reference orbits are presented. Numerical simulations using the approximations of the motions of Asteroid Toutatis and Comet P/Honda—Mrkos—Pajdu?áková as examples confirm that the constructed orbits are highly efficient. Their application is particularly beneficial in investigating strongly perturbed motion.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Quasi-satellite orbits in the general context of dynamics in the 1:1 mean motion resonance: perturbative treatment

Our investigation is motivated by the recent discovery of asteroids orbiting the Sun and simultaneously staying near one of the Solar System planets for a long time. This regime of motion is usually called the quasi-satellite regime, since even at the times of the closest approaches the distance between the asteroid and the planet is significantly larger than the region of space (the Hill’s sphere) in which the planet can hold its satellites. We explore the properties of the quasi-satellite regimes in the context of the spatial restricted circular three-body problem “Sun–planet–asteroid”. Via double numerical averaging, we construct evolutionary equations which describe the long-term behaviour of the orbital elements of an asteroid. Special attention is paid to possible transitions between the motion in a quasi-satellite orbit and the one in another type of orbits available in the 1:1 resonance. A rough classification of the corresponding evolutionary paths is given for an asteroid’s motion with a sufficiently small eccentricity and inclination.     

Dynamical Effects of Mars on Asteroidal Dust Particles

Asteroidal dust particles resulting from family-forming events migrate from their source locations in the asteroid belt inwards towards the Sun under the effect of Poynting-Robertson (PR) drag. Understanding the distribution of these dust particle orbits in the inner solar system is of great importance to determining the asteroidal contribution to the zodiacal cloud, the accretion rate by the Earth, and the threat that these particles pose to spacecraft and satellites in near-Earth space. In order to correctly describe this distribution of orbits in the inner solar system, we must track the dynamical perturbations that the dust particle orbits experience as they migrate inwards. In a seminal paper Öpik (1951) determines that very few of the μm-cm sized dust particles suffer a collision with the planet face as they decay inwards past Mars. Here we re-analyze this problem, considering additionally the likelihood that the dust particle orbits pass through the Hill sphere of Mars (to various depths) and experience potentially significant perturbations to their orbits. We find that a considerable fraction of dust particle orbits will enter the Hill sphere of Mars. Furthermore, we find that there is a bias with inclination, particle size, and eccentricity of the particle orbits that enter the Martian Hill sphere. In particular the bias with inclination may create a bias towards higher-inclination sources in the proportions of asteroid family particles that reach near-Earth space.     

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.     

A numerical investigation of coorbital stability and libration in three dimensions

Motivated by the dynamics of resonance capture, we study numerically the coorbital resonance for inclination \(0\le I\le 180^\circ \) in the circular restricted three-body problem. We examine the similarities and differences between planar and three dimensional coorbital resonance capture and seek their origin in the stability of coorbital motion at arbitrary inclination. After we present stability maps of the planar prograde and retrograde coorbital resonances, we characterize the new coorbital modes in three dimensions. We see that retrograde mode I (R1) and mode II (R2) persist as we change the relative inclination, while retrograde mode III (R3) seems to exist only in the planar problem. A new coorbital mode (R4) appears in 3D which is a retrograde analogue to an horseshoe-orbit. The Kozai–Lidov resonance is active for retrograde orbits as well as prograde orbits and plays a key role in coorbital resonance capture. Stable coorbital modes exist at all inclinations, including retrograde and polar obits. This result confirms the robustness the coorbital resonance at large inclination and encourages the search for retrograde coorbital companions of the solar system’s planets.     

Periodic orbits in the generalized photogravitational Chermnykh-like problem with power-law profile

The orbits about Lagrangian equilibrium points are important for scientific investigations. Since, a number of space missions have been completed and some are being proposed by various space agencies. In light of this, we consider a more realistic model in which a disk, with power-law density profile, is rotating around the common center of mass of the system. Then, we analyze the periodic motion in the neighborhood of Lagrangian equilibrium points for the value of mass parameter $0<\mu\leq\frac{1}{2}$ . Periodic orbits of the infinitesimal mass in the vicinity of equilibrium are studied analytically and numerically. In spite of the periodic orbits, we have found some other kind of orbits like hyperbolic, asymptotic etc. The effects of radiation factor as well as oblateness coefficients on the motion of infinitesimal mass in the neighborhood of equilibrium points are also examined. The stability criteria of the orbits is examined with the help of Poincaré surfaces of section (PSS) and found that stability regions depend on the Jacobi constant as well as other parameters.     

Using FLI maps for preliminary spacecraft trajectory design in multi-body environments

Fast Lyapunov Indicator (FLI) maps are presented as a tool for solving spacecraft preliminary trajectory design problems in multi-body environments with long-term stability requirements. In particular, the FLI maps are shown to provide a global overview of the dynamics in the restricted three-body problem that can guide mission designers in selecting long-term stable regions of phase space which are inherently more robust to model parameter perturbations. The FLI is also shown to numerically detect the normally hyperbolic manifolds associated with unstable periodic orbits. These, in turn, provide a global map of the principal heteroclinic connections between the various resonance regions which form the basic backbone of dynamical transfers design. Examples of maps and transfers are provided in the restricted three-body problem modeling the Jupiter–Europa system.     

Energy and stability in the Full Two Body Problem

The conditions for relative equilibria and their stability in the Full Two Body Problem are derived for an ellipsoid–sphere system. Under constant angular momentum it is found that at most two solutions exist for the long-axis solutions with the closer solution being unstable while the other one is stable. As the non-equilibrium problem is more common in nature, we look at periodic orbits in the F2BP close to the relative equilibrium conditions. Families of periodic orbits can be computed where the minimum energy state of one family is the relative equilibrium state. We give results on the relative equilibria, periodic orbits and dynamics that may allow transition from the unstable configuration to a stable one via energy dissipation.      

Solutions and periodicity of satellite relative motion under even zonal harmonics perturbations

Finding relative satellite orbits that guarantee long-term bounded relative motion is important for cluster flight, wherein a group of satellites remain within bounded distances while applying very few formationkeeping maneuvers. However, most existing astrodynamical approaches utilize mean orbital elements for detecting bounded relative orbits, and therefore cannot guarantee long-term boundedness under realistic gravitational models. The main purpose of the present paper is to develop analytical methods for designing long-term bounded relative orbits under high-order gravitational perturbations. The key underlying observation is that in the presence of arbitrarily high-order even zonal harmonics perturbations, the dynamics are superintegrable for equatorial orbits. When only J ₂ is considered, the current paper offers a closed-form solution for the relative motion in the equatorial plane using elliptic integrals. Moreover, necessary and sufficient periodicity conditions for the relative motion are determined. The proposed methodology for the J ₂-perturbed relative motion is then extended to non-equatorial orbits and to the case of any high-order even zonal harmonics (J _2n , n ≥ 1). Numerical simulations show how the suggested methodology can be implemented for designing bounded relative quasiperiodic orbits in the presence of the complete zonal part of the gravitational potential.