Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

On the periodic orbits and the integrability of the regularized Hill lunar problem

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

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.

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.

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.

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.

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.

Resonance and chaos .2. Exterior resonances and asymmetric libration

The motion of a test particle in the vicinity of exterior resonances is examined in the context of the planar, circular, restricted three-body problem. The existence of asymmetric periodic orbits associated with the 1 : n resonances (where n = 2, 3, 4, 5) is confirmed; there is also evidence of asymmetric resonances associated with larger values of n. A detailed examination of the evolution of the family of orbits associated with the 1:2 resonance shows the sequence that leads to asymmetric libration. on the basis of numerical studies of the phase space it is concluded that the existence of asymmetric libration means that the region exterior to the perturbing mass is more chaotic than the interior region. The apparent absence of 'particles' in 1 : n resonances in the solar system may reflect this inherent bias.

Third body perturbation using a single averaged model considering elliptic orbits for the disturbing body

In the present work, we expanded the study done by Solorzanol(1) including the eccentricity of the perturbing body. The assumptions used to develop the single-averaged analytical model are the same ones of the restricted elliptic three-body problem. The disturbing function was expanded in Legendre polynomials up to fourth-order. After that, the equations of motion are obtained from the planetary equations and we performed a set of numerical simulations. Different initial eccentricities for the perturbing and perturbed body are considered. The results obtained perform an analysis of the stability of a near-circular orbits and investigate under which conditions this orbit remain near-circular.

Resonance and chaos: II. Exterior resonances and asymmetric libration

The motion of a test particle in the vicinity of exterior resonances is examined in the context of the planar, circular, restricted three-body problem. The existence of asymmetric periodic orbits associated with the 1 : n resonances (where n = 2, 3, 4, 5) is confirmed; there is also evidence of asymmetric resonances associated with larger values of n. A detailed examination of the evolution of the family of orbits associated with the 1:2 resonance shows the sequence that leads to asymmetric libration. On the basis of numerical studies of the phase space it is concluded that the existence of asymmetric libration means that the region exterior to the perturbing mass is more chaotic than the interior region. The apparent absence of 'particles' in 1 : n resonances in the solar system may reflect this inherent bias.

The stability evolution of a family of simply periodic lunar orbits

The present work deals with a family of simply periodic orbits around the Moon in the rotating Earth Moon-particle system. Taking the framework of the planar, circular, restricted three-body problem, we follow the evolution of this family of periodic orbits using the numerical technique of Poincaré surface of section. The maximum amplitude of oscillation about the periodic orbits are determined and can be used as a parameter to measure the degree of stability in the phase space for such orbits. Despite the fact that the whole family of periodic orbits remain stable, there is a dichotomy in the quasi-periodic ones at the Jacobi constant Cj = 2.85. The quasi-periodic orbits with Cj < 2.85 oscillate around the periodic orbits in a different way from those with Cj > 2.85. © 1999 Elsevier Science Ltd. All rights reserved.

The use of the two-body energy to study problems of escape/capture Ernesto Vieira Neto1

The problem of escape/capture is encountered in many problems of the celestial mechanics -the capture of the giants planets irregular satellites, comets capture by Jupiter, and also orbital transfer between two celestial bodies as Earth and Moon. To study these problems we introduce an approach which is based on the numerical integration of a grid of initial conditions. The two-body energy of the particle relative to a celestial body defines the escape/capture. The trajectories are integrated into the past from initial conditions with negative two-body energy. The energy change from negative to positive is considered as an escape. By reversing the time, this escape turns into a capture. Using this technique we can understand many characteristics of the problem, as the maximum capture time, stable regions where the particles cannot escape from, and others. The advantage of this kind of approach is that it can be used out of plane (that is, for any inclination), and with perturbations in the dynamics of the n-body problem. © 2005 International Astronomical Union.

Quantum propagator for some classes of three-dimensional three-body systems

In this work we solve exactly a class of three-body propagators for the most general quadratic interactions in the coordinates, for arbitrary masses and couplings. This is done both for the constant as the time-dependent couplings and masses, by using the Feynman path integral formalism. Finally, the energy spectrum and the eigenfunctions are recovered from the propagators. © 2005 Elsevier Inc. All rights reserved.