Dynamics of two satellites in the 2/1 Mean-Motion resonance: application to the case of Enceladus and Dione

The dynamics of a pair of satellites similar to Enceladus-Dione is investigated with a two-degrees-of-freedom model written in the domain of the planar general three-body problem. Using surfaces of section and spectral analysis methods, we study the phase space of the system in terms of several parameters, including the most recent data. A detailed study of the main possible regimes of motion is presented, and in particular we show that, besides the two separated resonances, the phase space is replete of secondary resonances.

Effect of Jupiter's mass growth on satellite capture

We study the effects of Jupiter mass growth in order to permanently capture prograde satellites. Adopting the restricted three-body problem, Sun-Jupiter-Particle, we performed numerical simulations backward in time while considering the decrease in Jupiter's mass. We considered the particle's initial conditions to be prograde, at pericenter, in the region 100R(4) <= a <= 400R(4) and 0 <= e <= 0.5. The results give Jupiter's mass at the moment when the particle escapes from the planet. Such values give an indication of the conditions that are necessary for capture. An analysis of these results shows that prograde satellite capture is more complex than a retrograde one. It occurs in a two-step process. First, when the particles get inside about 0.85R(Hill) (Hills' radius), they become weakly bound to Jupiter. Then, they keep migrating toward the planet with a strong decrease in eccentricity, while the planet is growing. The radial oscillation of the particles reduces significantly when they reach a radial distance that is less than about 0.45R(Hill) from the planet. Three-dimensional simulations for the known prograde satellites of Jupiter were performed. The results indicate that Leda, Himalia, Lysithea, and Elara could have been permanently captured when Jupiter had between 50% and 60% of its present mass.

The effect of Jupiter's mass growth on satellite capture - Retrograde case

Gravitational capture can be used to explain the existence of the irregular satellites of giants planets. However, it is only the first step since the gravitational capture is temporary. Therefore, some kind of non-conservative effect is necessary to to turn the temporary capture into a permanent one. In the present work we study the effects of Jupiter mass growth for the permanent capture of retrograde satellites. An analysis of the zero velocity curves at the Lagrangian point L-1 indicates that mass accretion provides an increase of the confinement region ( delimited by the zero velocity curve, where particles cannot escape from the planet) favoring permanent captures. Adopting the restricted three-body problem, Sun-Jupiter-Particle, we performed numerical simulations backward in time considering the decrease of M-4. We considered initial conditions of the particles to be retrograde, at pericenter, in the region 100 R-4 less than or equal to a less than or equal to 400 R-4 and 0 less than or equal to e < 0.5. The results give Jupiter's mass at the moment when the particle escapes from the planet. Such values are an indication of the necessary conditions that could provide capture. An analysis of these results shows that retrograde satellites would be captured as soon as they get inside the Hills' radius and after that they keep migrating toward the planet while it is growing. For the region where the orbits of the four old retrograde satellites of Jupiter ( Ananke, Carme, Pasiphae and Sinope) are located we found that such satellites could have been permanently captured when Jupiter had between 62% and 93% of its present mass.

Stable satellites around extrasolar giant planets

In this work, we study the stability of hypothetical satellites of extrasolar planets. Through numerical simulations of the restricted elliptic three-body problem we found the borders of the stable regions around the secondary body. From the empirical results, we derived analytical expressions of the critical semimajor axis beyond which the satellites would not remain stable. The expressions are given as a function of the eccentricities of the planet, e(P), and of the satellite, e(sat). In the case of prograde, satellites, the critical semimajor axis, in the units of Hill's radius, is given by a(E) approximate to 0.4895 (1.0000 - 1.0305e(P) - 0.2738e(sat)). In the case of retrograde satellites, it is given by a(E) approximate to 0.9309 (1.0000 - 1.0764e(P) - 0.9812e(sat)). We also computed the satellite stability region (a(E)) for a set of extrasolar planets. The results indicate that extrasolar planets in the habitable zone could harbour the Earth-like satellites.

Dynamics of a spacecraft and normalization around Lagrangian points in the Neptune-Triton system

The problem of a spacecraft orbiting the Neptune-Triton system is presented. The new ingredients in this restricted three body problem are the Neptune oblateness and the high inclined and retrograde motion of Triton. First we present some interesting simulations showing the role played by the oblateness on a Neptune's satellite, disturbed by Triton. We also give an extensive numerical exploration in the case when the spacecraft orbits Triton, considering Sun, Neptune and its planetary oblateness as disturbers. In the plane a x I (a = semi-major axis, I = inclination), we give a plot of the stable regions where the massless body can survive for thousand of years. Retrograde and direct orbits were considered and as usual, the region of stability is much more significant for the case of direct orbit of the spacecraft (Triton's orbit is retrograde). Next we explore the dynamics in a vicinity of the Lagrangian points. The Birkhoff normalization is constructed around L-2, followed by its reduction to the center manifold. In this reduced dynamics, a convenient Poincare section shows the interplay of the Lyapunov and halo periodic orbits, Lissajous and quasi-halo tori as well as the stable and unstable manifolds of the planar Lyapunov orbit. To show the effect of the oblateness, the planar Lyapunov family emanating from the Lagrangian points and three-dimensional halo orbits are obtained by the numerical continuation method. Published by Elsevier Ltd. on behalf of COSPAR.

Possibility of collision between co-orbital asteroids and the Earth

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Distant stable direct orbits around the Moon

In the present work we explore regions of distant direct stable orbits around the Moon. First, the location and size of apparently stable regions are searched for numerically, adopting the approach of temporary capture time presented in Vieira Neto & Winter (2001). The study is made in the framework of the planar, circular, restricted three-body problem, Earth-Moon-particle. Regions of the initial condition space whose trajectories are apparently stable are determined. The criterion adopted was that the trajectories do not escape from the Moon during an integration period of 10(4) days. Using Poincare surface of sections the reason for the existence of the two stable regions found is studied. The stability of such regions proved to be due to two families of simple periodic orbits, h1 and h2, and the associated quasi-periodic orbits that oscillate around them. The robustness of the stability of the larger region, h2, is tested with the inclusion of the solar perturbation. The size of the region decreases, but it is still significant in size and can be useful in spacecraft missions.

Time analysis for temporary gravitational capture - Stable orbits

In a previous work, Vieira Neto & Winter (2001) numerically explored the capture times of particles as temporary satellites of Uranus. The study was made in the framework of the spatial, circular, restricted three-body problem. Regions of the initial condition space whose trajectories are apparently stable were determined. The criterion adopted was that the trajectories do not escape from the planet during an integration of 10(5) years. These regions occur for a wide range of orbital initial inclinations (i). In the present work it is studied the reason for the existence of such stable regions. The stability of the planar retrograde trajectories is due to a family of simple periodic orbits and the associated quasi-periodic orbits that oscillate around them. These planar stable orbits had already been studied (Henon 1970; Huang & Innanen 1983). Their results are reviewed using Poincare surface of sections. The stable non-planar retrograde trajectories, 110 degrees less than or equal to i < 180, are found to be tridimensional quasi-periodic orbits around the same family of periodic orbits found for the planar case (i = 180 degrees). It was not found any periodic orbit out of the plane associated to such quasi-periodic orbits. The largest region of stable prograde trajectories occurs at i = 60 degrees. Trajectories in such region are found to behave as quasi-periodic orbits evolving similarly to the stable retrograde trajectories that occurs at i = 120 degrees.

The Liapunov exponent as a tool for exploring phase space

We have used the Liapunov exponent to explore the phase space of a dynamical system. Considering the planar, circular restricted three-body problem for a mass ratio mu = 10(-3) (close to the Jupiter/Sun case), we have integrated similar to 16,000 starting conditions for orbits started interior to that of the perturber and we have estimated the maximum Liapunov characteristic exponent for each starting condition. Despite the fact that the integrations, in general, are for only a few thousand orbital periods of the secondary, a comparative analysis of the Liapunov exponents for various values of the 'cut-off' gives a good overview of the structure of the phase space. It provides information about the diffusion rates of the various chaotic regions, the location of the regular regions associated with primary resonances and even details such as the location of secondary resonances that produce chaotic regions inside the regular regions of primary resonances.

On the scattering of comets by a planet

In this paper, numerical simulations are made, using the three-dimensional restricted three-body problem as the mathematical model, to calculate the effects of a swing-by with the planet Saturn in the orbit of a comet. To show the results, the orbit of the comet is classified in four groups: elliptic direct, elliptic retrograde, hyperbolic direct and hyperbolic retrograde. Then, the modification in the orbit of the comet due to the close approach is shown in plots that specify from which group the comet's orbit is coming and to which group it is going. Several families of orbits are found and shown in detail. An analysis about the trends as parameters (position and velocity at the periapse) vary is performed and the influence of each of them is shown and explained. The result is a collection of maps that describe the evolution of the trajectory of the comet due to the close approach. Those maps can be used to estimate the probability of some events, like the capture or escape of a comet. An example of this technique is shown in the paper. (C) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Resonance and chaos .1. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted three-body problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincare surfaces of section with a mass ratio of 10(-3) (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. on theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Effects of the eccentricity on the primaries in the gravitational capture phenomenon

Gravitational capture is a characteristic of some dynamical systems in celestial mechanics, as in the elliptic restricted three-body problem that is considered in this paper. The basic idea is that a spacecraft (or any particle with negligible mass) can change a hyperbolic orbit with a small positive energy around a celestial body into an elliptic orbit with a small negative energy without the use of any propulsive system. The force responsible for this modification in the orbit of the spacecraft is the gravitational force of the third body involved in the dynamics. In this way, this force is used as a zero cost control, equivalent to a continuous thrust applied in the spacecraft. One of the most important applications of this property is the construction of trajectories to the Moon. The objective of the present paper is to study in some detail the effects of the eccentricity of the primaries in this maneuver.

Peculiar trajectories around the Lagrangian equilateral points

In the present work we analyse the behaviour of a particle under the gravitational influence of two massive bodies and a particular dissipative force. The circular restricted three body problem, which describes the motion of this particle, has five equilibrium points in the frame which rotates with the same angular velocity as the massive bodies: two equilateral stable points (L-4, L-5) and three colinear unstable points (L-1, L-2, L-3). A particular solution for this problem is a stable orbital libration, called a tadpole orbit, around the equilateral points. The inclusion of a particular dissipative force can alter this configuration. We investigated the orbital behaviour of a particle initially located near L4 or L5 under the perturbation of a satellite and the Poynting-Robertson drag. This is an example of breakdown of quasi-periodic motion about an elliptic point of an area-preserving map under the action of dissipation. Our results show that the effect of this dissipative force is more pronounced when the mass of the satellite and/or the size of the particle decrease, leading to chaotic, although confined, orbits. From the maximum Lyapunov Characteristic Exponent a final value of gamma was computed after a time span of 10(6) orbital periods of the satellite. This result enables us to obtain a critical value of log y beyond which the orbit of the particle will be unstable, leaving the tadpole behaviour. For particles initially located near L4, the critical value of log gamma is -4.07 and for those particles located near L-5 the critical value of log gamma is -3.96. (c) 2006 Elsevier B.V. All rights reserved.

Alternative paths to Earth-Moon transfer

The planar, circular, restricted three-body problem predicts the existence of periodic orbits around the Lagrangian equilibrium point L1. Considering the Earth-lunar-probe system, some of these orbits pass very close to the surfaces of the Earth and the Moon. These characteristics make it possible for these orbits, in spite of their instability, to be used in transfer maneuvers between Earth and lunar parking orbits. The main goal of this paper is to explore this scenario, adopting a more complex and realistic dynamical system, the four-body problem Sun-Earth-Moon-probe. We defined and investigated a set of paths, derived from the orbits around L1, which are capable of achieving transfer between low-altitude Earth (LEO) and lunar orbits, including high-inclination lunar orbits, at a low cost and with flight time between 13 and 15 days.

Irregular Satellites of Jupiter: a study of the capture direction

The irregular satellites of Jupiter are believed to be captured asteroids or planetesimals. In the present work is studied the direction of capture of these objects as a function of their orbital inclination. We performed numerical simulations of the restricted three-body problem, Sun-Jupiter-particle, taking into account the growth of Jupiter. The integration was made backward in time. Initially, the particles have orbits as satellites of Jupiter, which has its present mass. Then, the system evolved with Jupiter losing mass and the satellites escaping from the planet. The reverse of the escape direction corresponds to the capture direction. The results show that the Lagrangian points L1 and L2 mainly guide the direction of capture. Prograde satellites are captured through these two gates with very narrow amplitude angles. In the case of retrograde satellites, these two gates are wider. The capture region increases as the orbital inclination increases. In the case of planar retrograde satellites the directions of capture cover the whole 360 degrees around Jupiter. We also verified that prograde satellites are captured earlier in actual time than retrograde ones.