Transformation of Trojans into quasi-satellites during planetary migration and their subsequent close-encounters with the host planet

We use numerical integrations to investigate the dynamical evolution of resonant Trojan and quasi-satellite companions during the late stages of migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. Our migration simulations begin with Jupiter and Saturn on orbits already well separated from their mutual 2:1 mean-motion resonance. Neptune and Uranus are decoupled from each other and have orbital eccentricities damped to near their current values. From this point we adopt a planet migration model in which the migration speed decreases exponentially with a characteristic timescale τ (the e-folding time). We perform a series of numerical simulations, each involving the migrating giant planets plus test particle Trojans and quasi-satellites. We find that the libration frequencies of Trojans are similar to those of quasi-satellites. This similarity enables a dynamical exchange of objects back and forth between the Trojan and quasi-satellite resonances during planetary migration. This exchange is facilitated by secondary resonances that arise whenever there is more than one migrating planet. For example, secondary resonances may occur when the circulation frequencies, f, of critical arguments for the Uranus-Neptune 2:1 mean-motion near-resonance are commensurate with harmonics of the libration frequency of the critical argument for the Trojan and quasi-satellite 1:1 mean-motion resonance . Furthermore, under the influence of these secondary resonances quasi-satellites can have their libration amplitudes enlarged until they undergo a close-encounter with their host planet and escape from the resonance. High-resolution simulations of this escape process reveal that ≈80% of jovian quasi-satellites experience one or more close-encounters within Jupiter’s Hill radius (R_H) as they are forced out of the quasi-satellite resonance. As many as ≈20% come within R_H/4 and ≈2.5% come within R_H/10. Close-encounters of escaping quasi-satellites occur near or even below the 2-body escape velocity from the host planet. Finally, the exchange and escape of Trojans and quasi-satellites continues to as late as 6-9τ in some simulations. By this time the dynamical evolution of the planets is strongly dominated by distant gravitational perturbations between the planets rather than the migration force. This suggests that exchange and escape of Trojans and quasi-satellites may be a contemporary process associated with the present-day near-resonant configuration of some of the giant planets in our Solar System.     

How Long-Lived Are the Hypothetical Trojan Populations of Saturn, Uranus, and Neptune?

We investigate the possibility that fractions of the primordial populations at the triangular Lagrangian points of Saturn, Uranus, and Neptune have survived to the present and form (as yet unobserved) clusters of bodies coorbiting with these planets. Such leftovers would be analogs of the numerous objects (Trojans) leading and trailing the revolutions of Jupiter around the Sun. We focus on the dynamical stability of such populations over the age of the Solar System, assuming the current configuration of planets, and also discuss effects of the early radial migration of the outer planets. Our results suggest that, while Saturn's and Uranus' primordial Trojan populations should have been depleted by a factor of 100, Neptune may retain 50% of its original population of Trojans. A population of neptunian Trojans comparable to, or even larger than, Jupiter's Trojan population cannot be ruled out by existing observations. We compute the present-day sky densities of the hypothetical Trojans of the outer planets which can be used to guide observational surveys. Finally, we propose that the long-term instabilities that cause some jovian Trojans to escape the region of the Lagrange points at the present are due to three-body resonances.     

Dynamics of Mars Trojans

In this paper we explore the dynamical stability of the Mars Trojan region applying mainly Laskar's Frequency Map Analysis. This method yields the chaotic diffusion rate of orbits and allows to determine the most stable regions. It also gives the frequencies which are responsible for the instability of orbits. The most stable regions are found for inclinations between about 15° and 30°. For inclinations smaller than 15°, we confirm, by applying a synthetic secular theory, that the secular resonances ν₃, ν₄, ν₁₃, ν₁₄ rapidly excite asteroid orbits within a few Myrs, or even faster. The asteroids are removed from the Trojan region after a close encounter with Mars. For large inclinations, the secular resonance ν₅ clears a small region around 30° while the Kozai resonance rapidly removes bodies for inclinations larger than 35°. The dynamical lifetimes of the three L5 Trojans, (5261) Eureka, 1998 VF31, 2001 DH47, and the only L4 Trojan 1999 UJ7 are determined by numerically integrating clouds of corresponding clones over the age of the Solar System. All four Trojans reside in the most stable region with smallest diffusion coefficients. Their dynamical half-lifetime is of the order of the age of the Solar System. The Yarkovsky force has little effect on the known Trojans but for bodies smaller than about 1-5 m the drag is strong enough to destabilize Trojans on a timescale shorter than 4.5 Gyr.     

Irregular satellite capture during planetary resonance passage

The passage of Jupiter and Saturn through mutual 1:2 mean-motion resonance has recently been put forward as explanation for their relatively high eccentricities [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] and the origin of Jupiter's Trojans [Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. Additional constraints on this event based on other small-body populations would be highly desirable. Since some outer satellite orbits are known to be strongly affected by the near-resonance of Jupiter and Saturn (“the Great Inequality”; ?uk, M., Burns, J.A., 2004b. Astron. J. 128, 2518-2541), the irregular satellites are natural candidates for such a connection. In order to explore this scenario, we have integrated 9200 test particles around both Jupiter and Saturn while they went through a resonance-crossing event similar to that described by Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]. The test particles were positioned on a grid in semimajor axes and inclinations, while their initial pericenters were put at just 0.01 AU from their parent planets. The goal of the experiment was to find out if short-lived bodies, spiraling into the planet due to gas drag (or alternatively on orbits crossing those of the regular satellites), could have their pericenters raised by the resonant perturbations. We found that about 3% of the particles had their pericenters raised above 0.03 AU (i.e. beyond Iapetus) at Saturn, but the same happened for only 0.1% of the particles at Jupiter. The distribution of surviving particles at Saturn has strong similarities to that of the known irregular satellites. If saturnian irregular satellites had their origin during the 1:2 resonance crossing, they present an excellent probe into the early Solar System's evolution. We also explore the applicability of this mechanism for Uranus, and find that only some of the uranian irregular satellites have orbits consistent with resonant pericenter lifting. In particular, the more distant and eccentric satellites like Sycorax could be stabilized by this process, while closer-in moons with lower eccentricity orbits like Caliban probably did not evolve by this process alone.     

Trojan orbits in secular resonances

A near equality between the nodal rates of suitably defined Trojan orbits and Jupiter represents an important type of a secular resonance. This case is realized by the model Sun-Jupiter-Saturn-Trojan, referred to the invariable plane. A second theoretical example is based on the elliptic three-body problem Sun-Jupiter-Trojan, where the vanishing nodal rate of a special Trojan orbit and the vanishing rate of Jupiter's longitude of perihelion define a secular resonance.We investigate the perturbations in the asteroidal inclinations and the nodes and consider the possibility of a libration.     

Computer studies of the evolution of planetary and satellite systems

This paper describes two computer experiments carried out with a CDC-Cyber 74 program (Barricelliet al., 1979) for computer simulation of a large number of objects in orbit about a central body or primary. The first experiment was started with 125 planets of which the two largest ones had coplanar orbits and masses comparable to those of Jupiter and Saturn, respectively. Their semimajor axes and eccentricities were, however, much larger. The smaller planets had a distribution promoting the formation of an axial meeting area. The experiment gives information relevant to the question of focusing of planetary orbits into a common plane and to the question of the formation and stability of an axial meeting area. Together with the next experiment, it also gives information about the development of commensurabilities (or resonances) with the largest planets.The second experiment started with 55 planets none of them with a mass greater than about 20% of Jupiter's but several of them with orbits close to a common plane. The aim of the experiment was to investigate whether successive captures followed by planetary fusion (Barricelli, 1972a) could lead to the formation of major planets comparable to Jupiter and Saturn, and in similar orbits. Also this experiment gives information relevant to the commensurability problem.     

Physical properties and orbital stability of the Trojan asteroids

All the Trojan asteroids orbit about the Sun at roughly the same heliocentric distance as Jupiter. Differences in the observed visible reflection spectra range from neutral to red, with no ultra-red objects found so far. Given that the Trojan asteroids are collisionally evolved, a certain degree of variability is expected. Additionally, cosmic radiation and sublimation are important factors in modifying icy surfaces even at those large heliocentric distances. We search for correlations between physical and dynamical properties, we explore relationships between the following four quantities; the normalised visible reflectivity indexes (S^′), the absolute magnitudes, the observed albedos and the orbital stability of the Trojans. We present here visible spectroscopic spectra of 25 Trojans. This new data increase by a factor of about 5 the size of the sample of visible spectra of Jupiter Trojans on unstable orbits. The observations were carried out at the ESO-NTT telescope (3.5 m) at La Silla, Chile, the ING-WHT (4.2 m) and NOT (2.5 m) at Roque de los Muchachos observatory, La Palma, Spain. We have found a correlation between the size distribution and the orbital stability. The absolute-magnitude distribution of the Trojans in stable orbits is found to be bimodal, while the one of the unstable orbits is unimodal, with a slope similar to that of the small stable Trojans. This supports the hypothesis that the unstable objects are mainly byproducts of physical collisions. The values of S^′ of both the stable and the unstable Trojans are uniformly distributed over a wide range, from 0%/1000 Å to about 15%/1000 Å. The values for the stable Trojans tend to be slightly redder than the unstable ones, but no significant statistical difference is found.     

The resonant structure of Jupiter's Trojan asteroids – II. What happens for different configurations of the planetary system

In a previous paper, we have found that the resonance structure of the present Jupiter Trojan swarms could be split up into four different families of resonances. Here, in a first step, we generalize these families in order to describe the resonances occurring in Trojan swarms embedded in a generic planetary system. The location of these families changes under a modification of the fundamental frequencies of the planets and we show how the resonant structure would evolve during a planetary migration. We present a general method, based on the knowledge of the fundamental frequencies of the planets and on those that can be reached by the Trojans, which makes it possible to predict and localize the main events arising in the swarms during migration. In particular, we show how the size and stability of the Trojan swarms are affected by the modification of the frequencies of the planets. Finally, we use this method to study the global dynamics of the Jovian Trojan swarms when Saturn migrates outwards. Besides the two resonances found by Morbidelli et al. which could have led to the capture of the current population just after the crossing of the 2:1 orbital resonance, we also point out several sequences of chaotic events that can influence the Trojan population.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation when the two planets are in migration

In a previous paper (Hou et al. in Celest Mech Dyn Astron 119:119–142, 2014a), the problem of dynamical symmetry between two Jupiter triangular libration points (TLPs) with Saturn’s perturbation in the present configuration of the two planets was studied. A small short-time scale spatial asymmetry exists but gradually disappears with the time going, so the planar stable regions around the two Jupiter TLPs should be dynamically symmetric from a longtime perspective. In this paper, the symmetry problem is studied when the two planets are in migration. Several mechanisms that can cause asymmetries are discussed. Studies show that three important ones are the large short-time scale spatial asymmetry when Jupiter and Saturn are in resonance, the changing orbits of Jupiter and Saturn in the planet migration process, and the chaotic nature of Trojan orbits during the planet migration process. Their joint effects can cause an observable difference to the two Jupiter Trojan swarms. The thermal Yarkovsky effect is also found to be able to cause dynamical differences to the two TLPs, but generally they are too small to be practically observed.     

Survival of Trojan-type companions of Neptune during primordial planet migration

We investigate the survivability of Trojan-type companions of Neptune during primordial radial migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. We adopt the usual planet migration model in which the migration speed decreases exponentially with a characteristic time scale τ (the e-folding time). We perform a series of numerical simulations, each involving the migrating giant planets plus ∼1000 test particle Neptune Trojans with initial distributions of orbital eccentricity, inclination, and libration amplitude similar to those of the known jovian Trojans asteroids. We analyze these simulations to measure the survivability of Neptune's Trojans as a function of migration rate. We find that orbital migration with the characteristic time scale τ=10⁶ years allows about 35% of preexisting Neptune Trojans to survive to 5τ, by which time the giant planets have essentially reached their final orbits. In contrast, slower migration with τ=10⁷ years yields only a ∼5% probability of Neptune Trojans surviving to a time of 5τ. Interestingly, we find that the loss of Neptune Trojans during planetary migration is not a random diffusion process. Rather, losses occur almost exclusively during discrete prolonged episodes when Trojan particles are swept by secondary resonances associated with mean-motion commensurabilities of Uranus with Neptune. These secondary resonances arise when the circulation frequencies, f, of critical arguments for Uranus-Neptune mean-motion near-resonances (e.g., f^UN_1:2, f^UN_4:7) are commensurate with harmonics of the libration frequency of the critical argument for the Neptune-Trojan 1:1 mean-motion resonance (f^NT_1:1). Trojans trapped in the secondary resonances typically have their libration amplitudes amplified until they escape the 1:1 resonance with Neptune. Trojans with large libration amplitudes are susceptible to loss during sweeping by numerous high-order secondary resonances (e.g., f^UN_1:2≈11f^NT_1:1). However, for the slower migration, with τ=10⁷ years, even tightly bound Neptune Trojans with libration amplitudes below 10° can be lost when they become trapped in 1:3 or 1:2 secondary resonances between f^UN_1:2 and f^NT_1:1. With τ=10⁷ years the 1:2 secondary resonance was responsible for the single greatest episode of loss, ejecting nearly 75% of existing Neptune Trojans. This episode occurred during the late stages of planetary migration when the remnant planetesimal disk would have been largely dissipated. We speculate that if the number of bodies liberated during this event was sufficiently high they could have caused a spike in the impact rate throughout the Solar System.     

The dynamics of inclined Neptune Trojans

We investigated the stable area for fictive Trojan asteroids around Neptune’s Lagrangean equilibrium points with respect to their semimajor axis and inclination. To get a first impression of the stability region we derived a symplectic mapping for the circular and the elliptic planar restricted three body problem. The dynamical model for the numerical integrations was the outer Solar system with the Sun and the planets Jupiter, Saturn, Uranus and Neptune. To understand the dynamics of the region around L ₄ and L ₅ for the Neptune Trojans we also used eight different dynamical models (from the elliptic problem to the full outer Solar system model with all giant planets) and compared the results with respect to the largeness and shape of the stable region. Their dependence on the initial inclinations (0° < i < 70°) of the Trojans’ orbits could be established for all the eight models and showed the primary influence of Uranus. In addition we could show that an asymmetry of the regions around L ₄ and L ₅ is just an artifact of the different initial conditions.     

Jupiter's obliquity and a long-lived circumplanetary disk

It has been claimed [Canup, R.M., Ward, W.R., 2002. Astron. J. 124, 3404-3423; Ward, W.R., 2003. In: AGU, Fall Meeting 2003] that a long-lived minimum mass circumplanetary gas disk is inconsistent with Jupiter's low obliquity. Here we find that while Jupiter's obliquity may constrain its characteristics it does not rule out a long-lived massive (compared to the mass of the Galilean satellites) disk. This is because the argument assumes a Solar System much like that of the present day with the one exception of a circumjovian disk which is then allowed to dissipate on a long timescale (¹⁰⁶-¹⁰⁷ yr). Given that the sequence of events in Solar System history that fit known constraints is non-unique, we choose for the sake of clarity of exposition the orbital architecture framework of Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461], in which Jupiter and Saturn were once in compact, nearly coplanar orbits, and show that in this case Jupiter's low obliquity is consistent with the SEMM (solids-enhanced minimum mass) satellite formation model of Mosqueira and Estrada [Mosqueira, I., Estrada, P.R., 2003a. Icarus 163, 198-231; Mosqueira, I., Estrada, P.R., 2003b. Icarus 163, 232-255]. We suggest that a low inclination starting condition may apply, but stress that our SEMM satellite formation model could be compatible with Jupiter's obliquity even for mutually inclined giant planets.     

A planetary system with an escaping Mars

The chaotic behaviour of the motion of the planets in our Solar System is well established. In this work to model a hypothetical extrasolar planetary system our Solar System was modified in such a way that we replaced the Earth by a more massive planet and let the other planets and all the orbital elements unchanged. The major result of former numerical experiments with a modified Solar System was the appearance of a chaotic window at κ _E ∈ (4, 6), where the dynamical state of the system was highly chaotic and even the body with the smallest mass escaped in some cases. On the contrary for very large values of the mass of the Earth, even greater than that of Jupiter regular dynamical behaviour was observed. In this paper the investigations are extended to the complete Solar System and showed, that this chaotic window does still exist. Tests in different ‘Solar Systems’ clarified that including only Jupiter and Saturn with their actual masses together with a more ‘massive’ Earth (4 < κ _E < 6) perturbs the orbit of Mars so that it can even be ejected from the system. Using the results of the Laplace‐Lagrange secular theory we found secular resonances acting between the motions of the nodes of Mars, Jupiter and Saturn. These secular resonances give rise to strong chaos, which is the cause of the appearance of the instability window. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

An Upper Limit of Mass for a Stable Minor Planet in the Main Asteroid Belt

The stability of an imaginary planet located in the present main asteroid belt is studied with a 7-body model (Sun, Mars, Jupiter, Saturn, Uranus, Neptune and the imaginary planet). The fourth-order Hermite algorithm P(EC)³ is used, which has a very small secular energy error for the integration of periodic orbits with a constant time-step. The evolution of orbits is followed up to 10⁸ years. Our numerical results show that the low-order resonances with Jupiter can enhance the stability of the imaginary planet in some cases. The survival probability of the imaginary planet decreases with the planet mass. The upper limit of the imaginary planet's mass that can survive in the main belt is around 10²⁵ kg, i.e., about the Earth's mass.     

Implications of the Galilean satellites ice envelope explosions III. The origin of the Trojans and of some comets

Secondary explosions of the primary ice fragments ejected in the explosion of the electrolyzed massive ice envelopes of the Galilean satellites are capable of imparting velocities of up to ~5km s^–1 to the secondary fragments. As a result, the secondary fragments can enter the orbits of the irregular satellites (Agafonova and Drobyshevski, 1984b) and the Trojan libration orbits. In the latter case a perturbation velocity of V 0.3–2 km s^–1 is sufficient.The primary fragments ejected by the gravitational perturbations due to the Galilean satellites sunward from Jupiter's sphere of action move faster relative to the Sun than Jupiter does and therefore reach their first aphelion ahead of Jupiter in the neighborhood of L ₄. At the same time the fragments propelled from Jupiter's sphere of action beyond the planet's orbit approach it again in their perihelia behind Jupiter in the region of L ₅. The concentration of the fragments and, hence, the frequency of their collisions and explosions at L ₄ turn out to be much greater than those at L ₅. As a result, the number of the secondary fragments of diameter 15 km captured into libration orbits ahead of Jupiter can be as high as many hundreds and should exceed by more than a factor 3.5 that captured behind Jupiter.Since the icy mix of the fragments contains hydrocarbons and particulate material (silicates and the like), after ice sublimation from the surface layers the Trojans should reveal type C and RD spectra typical for Jupiter's irregular satellites, comet nuclei and other distant ice bodies of similar origin. Among the Trojans there cannot be rocky or metallic objects which are known to exist in the main asteroid belt.It is shown that a velocity perturbation of 150–200 m s^–1 resulting from a purely mechanical impact of two bodies may be sufficient to move collision fragments from the orbits of the Trojans to horseshoe-shaped trajectories with a subsequent transfer to the cometary orbits of Jupiter's family.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region     

The Trojan problem

There has been a renewed interest in the Trojan problem in recent years. Significant progress has been made in discovering and understanding dynamical features of motion of Jupiter's Trojan asteroids. The dynamics of hypothetical Trojan-type asteroids of other major planets has also been the subject of several recent investigations. This paper offers an overview on the current status of researches on real and hypothetical Trojan asteroids of the major planets. Results of analytical and numerical works are surveyed. Questions of dynamical properties, long-term evolution of orbits, stability regions around the triangular Lagrangian points are discussed among other problems of the Trojans.     

The primordial excitation and clearing of the asteroid belt—Revisited

We have performed new simulations of two different scenarios for the excitation and depletion of the primordial asteroid belt, assuming Jupiter and Saturn on initially circular orbits as predicted by the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. First, we study the effects of sweeping secular resonances driven by the depletion of the solar nebula. We find that these sweeping secular resonances are incapable of giving sufficient dynamical excitation to the asteroids for nebula depletion timescales consistent with estimates for solar-type stars, and in addition cannot cause significant mass depletion in the asteroid belt or produce the observed radial mixing of different asteroid taxonomic types. Second, we study the effects of planetary embryos embedded in the primordial asteroid belt. These embedded planetary embryos, combined with the action of jovian and saturnian resonances, can lead to dynamical excitation and radial mixing comparable to the current asteroid belt. The mass depletion driven by embedded planetary embryos alone, even in the case of an eccentric Jupiter and Saturn, is roughly 10-20× less than necessary to explain the current mass of the main belt, and thus a secondary depletion event, such as that which occurs naturally in the Nice Model, is required. We discuss the implications of our new simulations for the dynamical and collisional evolution of the main belt.     

The dynamics of Jupiter and Saturn in the gaseous protoplanetary disk

We study the possibility that the mutual interactions between Jupiter and Saturn prevented Type II migration from driving these planets much closer to the Sun. Our work extends previous results by Masset and Snellgrove [Masset, F., Snellgrove, M., 2001. Mon. Not. R. Astron. Soc. 320, L55-L59], by exploring a wider set of initial conditions and disk parameters, and by using a new hydrodynamical code that properly describes for the global viscous evolution of the disk. Initially both planets migrate towards the Sun, and Saturn's migration tends to be faster. As a consequence, they eventually end up locked in a mean motion resonance. If this happens in the 2:3 resonance, the resonant motion is particularly stable, and the gaps opened by the planets in the disk may overlap. This causes a drastic change in the torque balance for the two planets, which substantially slows down the planets' inward migration. If the gap overlap is substantial, planet migration may even be stopped or reversed. As the widths of the gaps depend on disk viscosity and scale height, this mechanism is particularly efficient in low viscosity, cool disks. The initial locking of the planets in the 2:3 resonance is a likely outcome if Saturn formed at the edge of Jupiter's gap, but also if Saturn initially migrated rapidly from further away. We also explore the possibility of trapping in other resonances, and the subsequent evolutions. We discuss the compatibility of our results with the initial conditions adopted in Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] and Gomes et al. [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469] to explain the current orbital architecture of the giant planets and the origin of the Late Heavy Bombardment of the Moon.