Measures of basins of attraction in spin-orbit dynamics

We consider a dissipative spin-orbit model where it is assumed that the orbit of the satellite is Keplerian, the obliquity is zero, and the dissipative effects depend linearly on the relative angular velocity. The measure of the basins of attraction associated to periodic and quasi-periodic attractors is numerically investigated. The results depend on the interaction among the physically relevant parameters, namely, the orbital eccentricity, the equatorial oblateness and the dissipative constant. In particular, it appears that, for astronomically relevant parameter values, for low eccentricities (as in the Moon’s case) about 96% of the initial data belong to the basin of attraction of the 1/1 spin-orbit resonance; for larger values of the eccentricities higher order spin-orbit resonances and quasi-periodic attractors become dominant providing a mechanism for explaining the observed state of Mercury into the 3/2 resonance.     

Construction of invariant tori for the spin-orbit problem in the mercury-sun system

The stability of spin-orbit resonances, namely commensurabilities between the periods of rotation and revolution of an oblate satellite orbiting around a primary body, is investigated using perturbation theory. We reduce the system to a model described by a one-dimensional, time-dependent Hamiltonian function. By means of KAM theory we rigorously construct bidimensional invariant surfaces, which separate the three dimensional phase space. In particular with a suitable choice of the rotation numbers of the invariant tori we are able to trap the periodic orbit associated with a given resonance in a finite region of the phase space. This technique is applied to the Mercury-Sun system. A connection with the probability of capture in a resonance is also provided.     

Nebular drag and capture into spin-orbit resonance

The majority of planetary satellites whose spin period is known are observed to be in synchronous spin-orbit resonance. The commonly accepted explanation for this observation is that it is due to the effects of tidal evolution. However, cosmogonic theories state that the formation of planetary and satellite systems occurs within a primordial solar nebula and circumplanetary nebulae, respectively. In this paper the influence of nebular drag on the capture into spin-orbit resonance is analysed. The results show that the torques generated are important for these resonances in a wide range of cases. Using the protojovian nebula model by Lunine and Stevenson (1982), conservative estimates of the despinning time scales for the Galilean satellites are computed. In comparison the despinning time scale from tidal effects are several orders of magnitude larger.on leave of absence from Departamento de Matemática, Faculdade de Engenharia de Guaratinguetá, UNESP, CP 205, 12500-000, Guaratinguetá, SP, Brazil     

Lunisolar Resonances Revisited

Lunisolar resonances arise in the artificial satellite problem without short-periodic terms. The basic model including the Earth's J ₂and a Hill-type model for the Sun or the Moon admits 20 different periodic terms which may lead to a resonance involving the satellite's perigee, node and the longitude of the perturbing body. Some of the resonances have been studied separately since 1960s. The present paper reviews all single resonances, attaching an appropriate fundamental model to each case. Only a part of resonances match known fundamental models. An extended fundamental model is proposed to account for some complicated phenomena. Most of the double resonance cases still remain unexplored.     

Hamiltonian Stability of Spin–Orbit Resonances in Celestial Mechanics

The behaviour of ‘resonances’ in the spin-orbit coupling in celestial mechanics is investigated in a conservative setting. We consider a Hamiltonian nearly-integrable model describing an approximation of the spin-orbit interaction. The continuous system is reduced to a mapping by integrating the equations of motion through a symplectic algorithm. We study numerically the stability of periodic orbits associated to the above mapping by looking at the eigenvalues of the matrix of the linearized map over the full cycle of the periodic orbit. In particular, the value of the trace of the matrix is related to the stability character of the periodic orbit. We denote by ε_* (p/q) the value of the perturbing parameter at which a given elliptic periodic orbit with frequency p/q becomes unstable. A plot of the critical function ε_* (p/q) versus the frequency at different orbital eccentricities shows significant peaks at the synchronous resonance (for low eccentricities) and at the synchronous and 3:2 resonances (at higher eccentricities) in good agreement with astronomical observations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.     

Planar periodic orbits in exterior resonances with Neptune

In the framework of the planar restricted three-body problem we study a considerable number of resonances associated to the basic dynamical features of Kuiper belt and located between 30 and 48 a.u. Our study is based on the computation of resonant periodic orbits and their stability. Stable periodic orbits are surrounded by regular librations in phase space and in such domains the capture of trans-Neptunian object is possible. All the periodic orbits found are symmetric and there is an indication of the existence of asymmetric ones only in a few cases. In the present work first, second and third order resonances are under consideration. In the planar circular case we found that most of the periodic orbits are stable. The families of periodic orbits are temporarily interrupted by collisions but they continue up to relatively large values of the Jacobi constant and highly eccentric regular motion exists for all cases. In the elliptic problem and for a particular eccentricity value of the primary bodies, the periodic orbits are isolated. The corresponding families, where they belong to, bifurcate from specific periodic orbits of the circular problem and seem to continue up to the rectilinear problem. Both stable and unstable orbits are obtained for each case. In the elliptic problem, the unstable orbits found are associated with narrow chaotic domains in phase space. The evolution of the orbits, which are located in such chaotic domains, seems to be practically regular and bounded for long time intervals.     

Repeat Ground Track Orbits of the Earth Tesseral Problem as Bifurcations of the Equatorial Family of Periodic Orbits

Orbits repeating their ground track on the surface of the earth are found to be members of periodic-orbit families (in a synodic frame) of the tesseral problem of the Earth artificial satellite. Families of repeat ground track orbits appear as vertical bifurcations of the equatorial family of periodic orbits, and they evolve from retrograde to direct motion throughout the 180 degrees of inclination. These bifurcations are always close to the resonances of the Earth's rotation rate and the mean motion of the orbiter.     

Tidally Induced Volcanism

The dissipation of tidal energy causes the ongoing silicate volcanism on Jupiter's satellite, Io, and cryovolcanism almost certainly has resurfaced parts of Saturn's satellite, Enceladus, at various epochs distributed over the latter's history. The maintenance of tidal dissipation in Io and the occurrence of the same on Enceladus depends crucially on the maintenance of the respective orbital eccentricities by the existence of mean motion resonances with nearby satellites. A formation of the resonances among the Galilean satellites by differential expansion of the satellite orbits from tides raised on Jupiter by the satellites means the onset of the volcanism on Io could be relatively recent. If, on the other hand, the resonances formed by differential migration from resonant interactions of the satellites with the disk of gas and particles from which they formed, Io would have been at least intermittently volcanically active throughout its history. Either means of assembling the Galilean satellite resonances lead to the same constraint on the dissipation function of Jupiter Q _J 10⁶, where the currently high heat flux from Io seems to favor episodic heating as Io's eccentricity periodically increases and decreases. Either of the two models might account for sufficient tidal dissipation in the icy satellite Enceladus to cause at least occasional cryovolcanism over much of its history. However, both models are assumption-dependent and not secure, so uncertainty remains on how tidal dissipation resurfaced Enceladus.     

Theory of satellite orbit-orbit resonance

On the basis of the strong mathematical and physical parallels between orbit-orbit and spin-orbit resonances, the dynamics of mutual orbit perturbations between two satellites about a massive planet are examined, exploiting an approach previously adopted in the study of spin-orbit coupling. The satellites are assumed to have arbitrary mass ratio and to move in non-intersecting orbits of arbitrary size and eccentricity. Resonances are found to exist when the mean orbital periods are commensurable with respect to some rotating axis, which condition also involves the apsidal and nodal motions of both satellites. In any resonant state the satellites are effectively trapped in separate potential wells, and a single variable is found to describe the simultaneous librations of both satellites. The librations in longitude are 180° out-of-phase, with fixed amplitude ratio that depends only on their relative masses and semimajor axes. At the same time the stroboscopic longitude of conjunction also librates about the commensurate axis with the same period. The theory is applicable to Saturn's resonant pairs Titan-Hyperion and Mimas-Tethys, and in these cases our calculated libration periods are in reasonably good agreement with the observed periods.This research supported under a grant from the California Institute of Technology President's Fund and NASA Contract NAS 7-100.     

The motion of axisymmetric satellite with drag and radiation pressure

The axisymmetric satellite problem including radiation pressure and drag is treated. The equations of motion of the satellite are derived. The energy-like and Laplace-like invariants of motion have been derived for a general drag force function of the polar angle, and the Laplace-like invariant is used to find the orbit equation in the case of a spherical satellite. Then using the small parameter, the orbit of the satellite is determined for an axisymmetric satellite.     

Periodic and quasi-periodic Earth satellite orbits

This paper reports the existence of three types of satellite orbit periodic in coordinates rotating with the Earth. Results on linear orbital stability are presented. Also included is a survey of the author's results on quasi-periodic orbits.     

Symmetric and asymmetric librations in extrasolar planetary systems: a global view

We present a global view of the resonant structure of the phase space of a planetary system with two planets, moving in the same plane, as obtained from the set of the families of periodic orbits. An important tool to understand the topology of the phase space is to determine the position and the stability character of the families of periodic orbits. The region of the phase space close to a stable periodic orbit corresponds to stable, quasi periodic librations. In these regions it is possible for an extrasolar planetary system to exist, or to be trapped following a migration process due to dissipative forces. The mean motion resonances are associated with periodic orbits in a rotating frame, which means that the relative configuration is repeated in space. We start the study with the family of symmetric periodic orbits with nearly circular orbits of the two planets. Along this family the ratio of the periods of the two planets varies, and passes through rational values, which correspond to resonances. At these resonant points we have bifurcations of families of resonant elliptic periodic orbits. There are three topologically different resonances: (1) the resonances (n + 1):n, (2:1, 3:2, ...), (2) the resonances (2n + 1):(2n-1), (3:1, 5:3, ...) and (3) all other resonances. The topology at each one of the above three types of resonances is studied, for different values of the sum and of the ratio of the planetary masses. Both symmetric and asymmetric resonant elliptic periodic orbits exist. In general, the symmetric elliptic families bifurcate from the circular family, and the asymmetric elliptic families bifurcate from the symmetric elliptic families. The results are compared with the position of some observed extrasolar planetary systems. In some cases (e.g., Gliese 876) the observed system lies, with a very good accuracy, on the stable part of a family of resonant periodic orbits.     

Resonance transition periodic orbits in the circular restricted three-body problem

This work studies a special type of cislunar periodic orbits in the circular restricted three-body problem called resonance transition periodic orbits, which switch between different resonances and revolve about the secondary with multiple loops during one period. In the practical computation, families of multiple periodic orbits are identified first, and then the invariant manifolds emanating from the unstable multiple periodic orbits are taken to generate resonant homoclinic connections, which are used to determine the initial guesses for computing the desired periodic orbits by means of multiple-shooting scheme. The obtained periodic orbits have potential applications for the missions requiring long-term continuous observation of the secondary and tour missions in a multi-body environment.     

Some properties of the dumbbell satellite attitude dynamics

The dumbbell satellite is a simple structure consisting of two point masses connected by a massless rod. We assume that it moves around the planet whose gravity field is approximated by the field of the attracting center. The distance between the point masses is assumed to be much smaller than the distance between the satellite’s center of mass and the attracting center, so that we can neglect the influence of the attitude dynamics on the motion of the center of mass and treat it as an unperturbed Keplerian one. Our aim is to study the satellite’s attitude dynamics. When the center of mass moves on a circular orbit, one can find a stable relative equilibrium in which the satellite is permanently elongated along the line joining the center of mass with the attracting center (the so called local vertical). In case of elliptic orbits, there are no stable equilibrium positions even for small values of the eccentricity. However, planar periodic motions are determined, where the satellite oscillates around the local vertical in such a way that the point masses do not leave the orbital plane. We prove analytically that these planar periodic motions are unstable with respect to out-of-plane perturbations (a result known from numerical investigations cf. Beletsky and Levin Adv Astronaut Sci 83, 1993). We provide also both analytical and numerical evidences of the existence of stable spatial periodic motions.     

On the periodic solutions of a rigid dumbbell satellite in a circular orbit

The aim of the present paper is to provide sufficient conditions for the existence of periodic solutions of the perturbed attitude dynamics of a rigid dumbbell satellite in a circular orbit.     

Planetary perturbations on Mercury’s libration in longitude

Two space missions dedicated to Mercury (MESSENGER and BepiColombo) aim at understanding its rotation and confirming the existence of a liquid core. This double challenge requires much more accurate models for the spin-orbit resonant rotation of Mercury. The purpose of this paper is to introduce planetary perturbations on Mercury’s rotation using an analytical method and to analyse the influence of the perturbations on the libration in longitude. Applying a perturbation theory based on the Lie triangle, we were able to re-introduce short periodic terms into the averaged Hamiltonian and to compute the evolution of the rotational variables. The perturbations on Mercury’s forced libration in longitude mainly come from the orbital motion of Mercury (with an amplitude around 41 arcsec that depends on the momenta of inertia). It is completed by various effects from Jupiter (11.86 and 5.93 year-periods), Venus (with a 5.66 year-period), Saturn (14.73 year-period), and the Earth (6.58 year-period). The amplitudes of the oscillations due to Jupiter and Venus are approximately 33% and 10% of those from the orbital motion of Mercury and the amplitudes of the oscillations due to Saturn and the Earth are approximately 3% and 2%. We compare the analytical results with the solution obtained from the spin-orbit numerical model SONYR.     

Origin and continuation of 3/2, 5/2, 3/1, 4/1 and 5/1 resonant periodic orbits in the circular and elliptic restricted three-body problem

We consider a planetary system consisting of two primaries, namely a star and a giant planet, and a massless secondary, say a terrestrial planet or an asteroid, which moves under their gravitational attraction. We study the dynamics of this system in the framework of the circular and elliptic restricted three-body problem, when the motion of the giant planet describes circular and elliptic orbits, respectively. Originating from the circular family, families of symmetric periodic orbits in the 3/2, 5/2, 3/1, 4/1 and 5/1 mean-motion resonances are continued in the circular and the elliptic problems. New bifurcation points from the circular to the elliptic problem are found for each of the above resonances, and thus, new families continued from these points are herein presented. Stable segments of periodic orbits were found at high eccentricity values of the already known families considered as whole unstable previously. Moreover, new isolated (not continued from bifurcation points) families are computed in the elliptic restricted problem. The majority of the new families mainly consists of stable periodic orbits at high eccentricities. The families of the 5/1 resonance are investigated for the first time in the restricted three-body problems. We highlight the effect of stable periodic orbits on the formation of stable regions in their vicinity and unveil the boundaries of such domains in phase space by computing maps of dynamical stability. The long-term stable evolution of the terrestrial planets or asteroids is dependent on the existence of regular domains in their dynamical neighbourhood in phase space, which could host them for long-time spans. This study, besides other celestial architectures that can be efficiently modelled by the circular and elliptic restricted problems, is particularly appropriate for the discovery of terrestrial companions among the single-giant planet systems discovered so far.