Co-orbital satellites of Saturn: congenital formation

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

On the stability of hypothetical satellites coorbital to Mimas or Enceladus

In the present work, we study the stability of hypothetical satellites that are coorbital with Enceladus and Mimas. We performed numerical simulations of 50 particles around the triangular Lagrangian equilibrium points of Enceladus and Mimas taking into account the perturbation of Mimas, Enceladus, Tethys, Dione, Titan and the oblateness of Saturn. All particles remain on tadpole orbits after 10 000 yr of integration. Since in the past the orbit of Enceladus and Mimas expanded due to the tidal perturbation, we also simulated the system with Enceladus and Mimas at several different values of semimajor axes. The results show that in general the particles remain on tadpole orbits. The exceptions occur when Enceladus is at semimajor axes that correspond to 6:7, 5:6 and 4:5 resonances with Mimas. Therefore, if Enceladus and Mimas had satellites librating around their Lagrangian triangular points in the past, they would have been removed if Enceladus crossed one of these first-order resonances with Mimas.

Numerical exploration of resonant dynamics in the system of Saturnian major satellites

We numerically investigate the long-term dynamics of the Saturnian system by analyzing the Fourier spectra of ensembles of orbits taken around the current orbits of Mimas, Enceladus, Tethys, Rhea and Hyperion. We construct dynamical maps around the current position of these satellites in their respective phase spaces. The maps are the result of a great deal of numerical simulations where we adopt dense sets of initial conditions and different satellite configurations. Several structures associated to the current two-body mean-motion resonances, unstable regions associated to close approaches between the satellites, and three-body mean-motion resonances in the system, are identified in the map. (C) 2010 Elsevier Ltd. All rights reserved.

Long-term dynamics of Methone, Anthe and Pallene

We numerically investigate the long-term dynamics of the Saturn's small satellites Methone (S/2004 S1), Anthe (S/2007 S4) and Pallene (S/2004 S2). In our numerical integrations, these satellites are disturbed by non-spherical shape of Saturn and the six nearest regular satellites. The stability of the small bodies is studied here by analyzing long-term evolution of their orbital elements.We show that long-term evolution of Pallene is dictated by a quasi secular resonance involving the ascending nodes (12) and longitudes of pericentric distances (pi) of Mimas (subscript 1) and Pallene (subscript 2), which critical argument is pi(2) - pi(1) - Omega(1) + Omega(2) Long-term orbital evolution of Methone and Anthe are probably chaotic since: i) their orbits randomly cross the orbit of Mimas in time scales of thousands years); ii) long-term numerical simulations involving both small satellites are strongly affected by small changes in the initial conditions.

Deformation and tidal evolution of close-in planets and satellites using a Maxwell viscoelastic rheology

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

Effects of the planetary migration on some primordial satellites of the outer planets

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

Analysis of 25 mutual eclipses and occultations between the Galilean satellites observed from Brazil in 2009

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

Orbital evolution of the mu and nu dust ring particles of Uranus

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

Time analysis for temporary gravitational capture: Satellites of Uranus

We study the problem of gravitational capture in the framework of the Sun-Uranus-particle system. Part of the space of initial conditions is systematically explored, and the duration of temporary gravitational capture is measured. The location and size of different capture-time regions are given in terms of diagrams of initial semimajor axis versus eccentricity. The other initial orbital elements - inclination (i), longitude of the node (Ω), argument of pericenter (ω), and time of pericenter passage (τ) - are first taken to be zero. Then we investigate the cases with ω = 90°, 180°, and 270°. We also present a sample of results for Ω = 90°, considering the cases i = 60°, 120°, 150°, and 180°. Special attention is given to the influence of the initial orbital inclination, taking orbits initially in opposition at pericenter. In this case, the initial inclination is varied from 0° to 180° in steps of 10°. The success of the final stage of the capture problem, which involves the transformation of temporary captures into permanent ones, is highly dependent on the initial conditions associated with the longest capture times. The largest regions of the initial-conditions space with the longest capture times occur at inclinations of 60°-70° and 160°. The regions of possible stability as a function of initial inclination are also delimited. These regions include not only a known set of retrograde orbits, but also a new sort of prograde orbit with inclinations greater than zero.

ORBITAL PERTURBATIONS OF THE GALILEAN SATELLITES DURING PLANETARY ENCOUNTERS

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

The role of Mab as a source for the mu ring of Uranus (Research Note)

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

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.

Nebular gas drag and planetary accretion with eccentric high-mass planets

Aims.We investigate the dynamics of pebbles immersed in a gas disk interacting with a planet on an eccentric orbit. The model has a prescribed gap in the disk around the location of the planetary orbit, as is expected for a giant planet with a mass in the range of 0.1-1 Jupiter masses. The pebbles with sizes in the range of 1 cm to 3 m are placed in a ring outside of the giant planet orbit at distances between 10 and 30 planetary Hill radii. The process of the accumulation of pebbles closer to the gap edge, its possible implication for the planetary accretion, and the importance of the mass and the eccentricity of the planet in this process are the motivations behind the present contribution. Methods. We used the Bulirsch-Stoer numerical algorithm, which is computationally consistent for close approaches, to integrate the Newtonian equations of the planar (2D), elliptical restricted three-body problem. The angular velocity of the gas disk was determined by the appropriate balance between the gravity, centrifugal, and pressure forces, such that it is sub-Keplerian in regions with a negative radial pressure gradient and super-Keplerian where the radial pressure gradient is positive. Results. The results show that there are no trappings in the 1:1 resonance around the L 4 and L5 Lagrangian points for very low planetary eccentricities (e2 < 0.07). The trappings in exterior resonances, in the majority of cases, are because the angular velocity of the disk is super-Keplerian in the gap disk outside of the planetary orbit and because the inward drift is stopped. Furthermore, the semi-major axis location of such trappings depends on the gas pressure profile of the gap (depth) and is a = 1.2 for a planet of 1 MJ. A planet on an eccentric orbit interacts with the pebble layer formed by these resonances. Collisions occur and become important for planetary eccentricity near the present value of Jupiter (e 2 = 0.05). The maximum rate of the collisions onto a planet of 0.1 MJ occurs when the pebble size is 37.5 cm ≤ s < 75 cm; for a planet with the mass of Jupiter, it is15 cm ≤ s < 30 cm. The accretion stops when the pebble size is less than 2 cm and the gas drag dominates the motion. © 2013 ESO.

Irregular satellites of Jupiter: capture configurations of binary-asteroids

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

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.