Effects of the Eccentricity of a Perturbing Third Body on the Orbital Correction Maneuvers of a Spacecraft

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Contact parameters in two dimensions for general three-body systems

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

On the evection resonance and its connection to the stability of outer satellites

In terms of stability around the primary, it is widely known that the semimajor axis of the retrograde satellites is much larger than the corresponding semimajor axis of the prograde satellites. Usually this conclusion is obtained numerically, since precise analytical derivation is far from being easy, especially, in the case of two or more disturbers. Following the seminal idea that what is unstable in the restricted three-body problem is also unstable in the general N-body problem, we present a simplified model which allows us to derive interesting resonant configurations. These configurations are responsible for cumulative perturbations which can give birth to strong instability that may cause the ejection of the satellite. Then we obtain, analytically, approximate bounds of the stability of prograde and retrograde satellites. Although we recover quite well previous results of other authors, we comment very briefly some weakness of these bounds. Copyright (c) 2008 Tadashi Yokoyama et al.

On the stability of the satellites of asteroid 87 Sylvia

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

Exploring S-type orbits in the Pluto-Charon binary system

This work generates, through a sample of numerical simulations of the restricted three-body problem, diagrams of semimajor axis and eccentricity which defines stable and unstable zones for particles in S-type orbits around Pluto and Charon. Since we consider initial conditions with 0 <= e <= 0.99, we found several new stable regions. We also identified the nature of each one of these newly found stable regions. They are all associated to families of periodic orbits derived from the planar circular restricted three-body problem. We have shown that a possible eccentricity of the Pluto-Charon system slightly reduces, but does not destroy, any of the stable regions.

Gravitational Capture of Asteroids by Gas Drag

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

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.

Stochastic variational approach applied to 2, 3 and 4 alpha particles

A stochastic variational method is applied to calculate the binding energies and root-mean-square radii of 2, 3 and 4 alpha particles using an S-wave Ali-Bodmer potential. The results agree with other calculations. We discuss the application of the present method to study the universality in weakly-bound three and four-body systems in the context of ultracold atomic traps.

Universality and Scaling Limit of Weakly-Bound Tetramers

The occurrence of a new limit cycle in few-body physics, expressing a universal scaling function relating the binding energies of two successive tetramer states, is revealed by considering a renormalized zero-range two-body interaction in bound state of four identical bosons. The tetramer energy spectrum is obtained by adding a boson to an Efimov bound state with energy B-3 in the unitary limit (for zero two-body binding energy or infinite two-body scattering length). Each excited N-th tetramer energy B-4((N)) is shown to slide along a scaling function as a short-range four-body scale is changed, emerging from the 3+1 threshold for a universal ratio B-4((N))/B-3 = 4.6, which does not depend on N. The new scale can also be revealed by a resonance in the atom-trimer recombination process.

Solutions of the bound-state Faddeev-Yakubovsky equations in three dimensions by using NN and 3N potential models

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

Binding and structure of tetramers in the scaling limit

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

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.

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.