Slow and Fast Diffusion in Asteroid-Belt Resonances: A Review

This paper reviews recent advances in several topics of resonant asteroidal dynamics as the role of resonances in the transportation of asteroids and asteroidal debris to the inner and outer solar system; the explanation of the contrast of a depleted 2/1 resonance (Hecuba gap) and a high-populated 3/2 resonance (Hilda group); the overall stochasticity created in the asteroid belt by the short-period perturbations of Jupiter's orbit, with emphasis in the formation of significant three-period resonances, the chaotic behaviour of the outer asteroid belt, and the depletion of the Hecuba gap. This revised version was published online in July 2006 with corrections to the Cover Date.     

Origin and sustainability of the population of asteroids captured in the exterior resonance 1:2 with Mars

At present, approximately 1500 asteroids are known to evolve inside or sticked to the exterior 1:2 resonance with Mars at a ? 2.418 AU, being (142) Polana the largest member of this group. The effect of the forced secular modes superposed to the resonance gives rise to a complex dynamical evolution. Chaotic diffusion, collisions, close encounters with massive asteroids and mainly orbital migration due to the Yarkovsky effect generate continuous captures to and losses from the resonance, with a fraction of asteroids remaining captured over long time scales and generating a concentration in the semimajor axis distribution that exceeds by 20% the population of background asteroids. The Yarkovsky effect induces different dynamics according to the asteroid size, producing an excess of small asteroids inside the resonance. The evolution in the resonance generates a signature on the orbits, mainly in eccentricity, that depends on the time the asteroid remains captured inside the resonance and on the magnitude of the Yarkovsky effect. The greater the asteroids, the larger the time they remain captured in the resonance, allowing greater diffusion in eccentricity and inclination. The resonance generates a discontinuity and mixing in the space of proper elements producing misidentification of dynamical family members, mainly for Vesta and Nysa-Polana families. The half-life of resonant asteroids large enough for not being affected by the Yarkovsky effect is about 1 Gyr. From the point of view of taxonomic classes, the resonant population does not differ from the background population and the excess of small asteroids is confirmed.     

Surfaces of section in the Miranda-Umbriel 3:1 inclination problem

The recent numerical simulations of Tittemore and Wisdom (1988, 1989, 1990) and Dermottet al. (1988), Malhotra and Dermott (1990) concerning the tidal evolution through resonances of some pairs of Uranian satellites have revealed interesting dynamical phenomena related to the interactions between close-by resonances. These interactions produce chaotic layers and strong secondary resonances. The slow evolution of the satellite orbits in this dynamical landscape is responsible for temporary capture into resonance, enhancement of eccentricity or inclination and subsequent escape from resonance. The present contribution aims at developing analytical tools for predicting the location and size of chaotic layers and secondary resonances. The problem of the 3:1 inclination resonance between Miranda and Umbriel is analysed.     

The dynamics of Neptune Trojan – I. The inclined orbits

The stability of Trojan type orbits around Neptune is studied. As the first part of our investigation, we present in this paper a global view of the stability of Trojans on inclined orbits. Using the frequency analysis method based on the fast Fourier transform technique, we construct high-resolution dynamical maps on the plane of initial semimajor axis a ₀ versus inclination i ₀. These maps show three most stable regions, with i ₀ in the range of  (0°, 12°), (22°, 36°)  and  (51°, 59°),  respectively, where the Trojans are most probably expected to be found. The similarity between the maps for the leading and trailing triangular Lagrange points L ₄ and L ₅ confirms the dynamical symmetry between these two points. By computing the power spectrum and the proper frequencies of the Trojan motion, we figure out the mechanisms that trigger chaos in the motion. The Kozai resonance found at high inclination varies the eccentricity and inclination of orbits, while the  ν₈  secular resonance around   i ₀∼ 44°  pumps up the eccentricity. Both mechanisms lead to eccentric orbits and encounters with Uranus that introduce strong perturbation and drive the objects away from the Trojan like orbits. This explains the clearance of Trojan at high inclination  (>60°)  and an unstable gap around  44°  on the dynamical map. An empirical theory is derived from the numerical results, with which the main secular resonances are located on the initial plane of  ( a ₀, i ₀)  . The fine structures in the dynamical maps can be explained by these secular resonances.     

Review of the dynamics in the Kirkwood gaps

The study of mean motion resonance dynamics was motivated by the search for an explanation for the puzzling problem of the Kirkwood gaps. The most important contributions in this field within the last 32 years are reviewed here. At the beginning of that period, which coincides with the first long-term numerical investigations of resonant motion, different hypotheses (collisional, gravitational, statistical and cosmological) to explain the origin of the gaps were still competing with each other. At present, a general theory, based on gravitational mechanisms only, is capable of explaining in a uniform way all the Kirkwood gaps except the 2/1 one. Indeed, in the 4/1, 3/1, 5/2 and 7/3 mean motion commensurabilities, the overlap of secular resonances leads to almost overall chaos where asteroids undergo large and wild variations in their orbital elements. Such asteroids, if not thrown directly into the Sun, are sooner or later subject to strong close encounters with the largest inner planets, the typical time scale of the whole process being of the order of a few million years. Unfortunately, this mechanism is not capable of explaining the 2/1 gap where the strong chaos produced by the overlapping secular resonances does not attain orbits with moderate eccentricity, of low inclination and with low to moderate amplitude of libration. In the light of the most recent studies, it appears that the 2/1 gap is the global consequence of slow diffusive processes. At present, the origin of these processes remains under study.     

On the Asteroidal Population of the First-Order Jovian Resonances

The frequency map analysis was applied to the fairly realistic models of the 2/1, 3/2, and 4/3 jovian resonances and the results were compared with the asteroidal distribution at these commensurabilities. The presence of the Hecuba gap at the 2/1 and of the Hilda group in the 3/2 is explained on the basis of different rates of the chaotic transport (diffusion) in these resonances. The diffusion in the most stable 2/1-resonant region is almost two orders in magnitude faster than the diffusion in the region which accommodates the Hildas. In the 2/1 commensurability there are two possible locations for long-surviving asteroids: the one centered at an eccentricity of 0.3 near the libration stable centers with small libration amplitude and the other at a slightly lower eccentricity with a moderate libration amplitude (∼90°). Surprisingly, all asteroids observed in the 2/1 resonance (8 numbered and multi-opposition objects in Bowell's catalog from 1994) occupy the moderate-libration area and avoid the area in a close vicinity of the libration stable centers. Possible explanations of this fact were discussed. Concerning the 4/3 resonance, the only asteroid in the corresponding stable region is 279 Thule, in spite of the fact that this region is almost as regular (although not as extensive) as the one where the Hilda group in the 3/2, with 79 members, is found.     

Atlas of the mean motion resonances in the Solar System

The aim of this work is to present a systematic survey of the strength of the mean motion resonances (MMRs) in the Solar System. We know by applying simple formulas where the resonances with the planets are located but there is no indication of the strength that these resonances have. We propose a numerical method for the calculation of this strength and we present an atlas of the MMRs constructed with this method. We found there exist several resonances unexpectedly strong and we look and find in the small bodies population several bodies captured in these resonances. In particular in the inner Solar System we find one asteroid in the resonance 6:5 with Venus, five asteroids in resonance 1:2 with Venus, three asteroids in resonance 1:2 with Earth and six asteroids in resonance 2:5 with Earth. We find some new possible co-orbitals of Earth, Mars, Saturn, Uranus and Neptune. We also present a discussion about the behavior of the resonant disturbing function and where the stable equilibrium points can be found at low and high inclination resonant orbits.     

A mechanism of formation for the Kirkwood gaps

In this paper an analytical model describing the effect of a displacement of the Jovian resonances in the asteroid belt is analyzed. It is found that small displacement can transform a truncated uniform density distribution of asteroids into a gap. As a possible explanation for the displacement, the effect of the removal of an accretion disk in the early stage of the solar system is investigated. It is found that removal of a disk containing a few percent of the solar mass between the orbit of the asteroids and the orbit of Jupiter is sufficient to account for the observed Hecuba gap.     

Yarkovsky origin of the unstable asteroids in the 2/1 mean motion resonance with Jupiter

The 2/1 mean motion resonance with Jupiter, intersecting the main asteroid belt at ≈3.27  au, contains a small population of objects. Numerical investigations have classified three groups within this population: asteroids residing on stable orbits (i.e. Zhongguos), those on marginally stable orbits with dynamical lifetimes of the order of 100 Myr (i.e. Griquas), and those on unstable orbits. In this paper, we reexamine the origin, evolution and survivability of objects in the 2/1 population. Using recent asteroid survey data, we have identified 100 new members since the last search, which increases the resonant population to 153. The most interesting new asteroids are those located in the theoretically predicted stable island A, which until now had been thought to be empty. We also investigate whether the population of objects residing on the unstable orbits could be resupplied by material from the edges of the 2/1 resonance by the thermal drag force known as the Yarkovsky effect (and by the YORP effect, which is related to the rotational dynamics). Using N -body simulations, we show that test particles pushed into the 2/1 resonance by the Yarkovsky effect visit the regions occupied by the unstable asteroids. We also find that our test bodies have dynamical lifetimes consistent with the integrated orbits of the unstable population. Using a semi-analytical Monte Carlo model, we compute the steady-state size distribution of magnitude   H < 14  asteroids on unstable orbits within the resonance. Our results provide a good match with the available observational data. Finally, we discuss whether some 2/1 objects may be temporarily captured Jupiter-family comets or near-Earth asteroids.     

Stable Chaos versus Kirkwood Gaps in the Asteroid Belt: A Comparative Study of Mean Motion Resonances

We have shown, in previous publications, that stable chaos is associated with medium/high-order mean motion resonances with Jupiter, for which there exist no resonant periodic orbits in the framework of the elliptic restricted three-body problem. This topological “defect” results in the absence of the most efficient mechanism of eccentricity transport (i.e., large-amplitude modulation on a short time scale) in three-body models. Thus, chaotic diffusion of the orbital elements can be quite slow, while there can also exist a nonnegligible set of chaotic orbits which are semiconfined (stable chaos) by “quasi-barriers” in the phase space. In the present paper we extend our study to all mean motion resonances of order q≤9 in the inner main belt (1.9-3.3 AU) and q≤7 in the outer belt (3.3-3.9 AU). We find that, out of the 34 resonances studied, only 8 possess resonant periodic orbits that are continued from the circular to the elliptic three-body problem (regular families), namely, the 2/1, 3/1, 4/1, and 5/2 in the inner belt and the 7/4, 5/3, 11/7, and 3/2 in the outer belt. Numerical results indicate that the 7/3 resonance also carries periodic orbits but, unlike the aforementioned resonances, 7/3-periodic orbits belong to an irregular family. Note that the five inner-belt resonances that carry periodic orbits correspond to the location of the main Kirkwood gaps, while the three outer-belt resonances correspond to gaps in the distribution of outer-belt asteroids noted by Holman and Murray (1996, Astron. J.112, 1278-1293), except for the 3/2 case where the Hildas reside. Fast, intermittent eccentricity increase is found in resonances possessing periodic orbits. In the remaining resonances the time-averaged elements of chaotic orbits are, in general, quite stable, at least for times t∼250 Myr. This slow diffusion picture does not change qualitatively, even if more perturbing planets are included in the model.     

Ringlet–satellite interactions

We study the interaction of a satellite and a nearby ringlet on eccentric and inclined orbits. Secular torques originate from mean motion resonances and the secular interaction potential which represents the m  = 1 global modes of the ring. The torques act on the relative eccentricity and inclination. The resonances damp the relative eccentricity. The inclination instability owing to the resonances is turned off by a finite differential eccentricity of the order of 0.27 for nearly coplanar systems. The secular potential torque damps the eccentricity and inclination and does not affect the relative semi-major axis; also, it suppresses the inclination instability that persists at small differential eccentricities. The damping of the relative eccentricity and inclination forces an initially circular and planar small mass ringlet to reach the eccentricity and inclination of the satellite. When the planet is oblate, the interaction of the satellite damps the proper precession of a small mass ringlet so that it precesses at the satellite's rate independently of their relative distance. The oblateness of the primary modifies the long-term eccentricity and inclination magnitudes and introduces a constant shift in the apsidal and nodal lines of the ringlet with respect to those of the satellite. These results are applied to Saturn's F-ring, which orbits between the moons Prometheus and Pandora.     

Collective resonant phenomena on small bodies in the solar system

For both asteroids and meteor streams, and also for comets, resonances play a major role for their orbital evolutions but on different time scales. For asteroids both mean motion resonances and secular resonances not only structure the phase space of regular orbits but are mainly at the origin for the inherent chaos of planet crosser objects.For comets and their chaotic routes temporary trapping into orbital resonances is a well known phenomenon. In addition for slow diffusion through the Kuiper belt resonances are the only candidates for originating a slow chaos.Like for asteroids, resonances with Jupiter play a major role for the orbital evolution of meteor streams. Crossing of separatrix like zones appears to be crucial for the formation of arcs and for the dissolution of streams. In particular the orbital inclination of a meteor stream appears to be a critical parameter for arc formation. Numerical results obtained in an other context show that the competition between the Poynting-Robertson drag and the gravitational interaction of grains near the 2/1 resonance might be very important in the long run for the structure of meteor streams.     

A survey of near-mean-motion resonances between Venus and Earth

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively, to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore, we did numerical experiments in different dynamical models with different initial conditions—on one hand the couple Venus–Earth was set close to different mean motion resonances (MMR), and on the other hand Venus’ orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40°). The couple Venus–Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth’s semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in initial eccentricity and/or inclination of Venus, and even escapes for the innermost planet are possible which may happen quite rapidly.     

1:1 Ground-track resonance in a uniformly rotating 4th degree and order gravitational field

Using a gravitational field truncated at the 4th degree and order, the 1:1 ground-track resonance is studied. To address the main properties of this resonance, a 1-degree of freedom (1-DOF) system is firstly studied. Equilibrium points (EPs), stability and resonance width are obtained. Different from previous studies, the inclusion of non-spherical terms higher than degree and order 2 introduces new phenomena. For a further study about this resonance, a 2-DOF model which includes a main resonance term (the 1-DOF system) and a perturbing resonance term is studied. With the aid of Poincaré sections, the generation of chaos in the phase space is studied in detail by addressing the overlap process of these two resonances with arbitrary combinations of eccentricity (e) and inclination (i). Retrograde orbits, near circular orbits and near polar orbits are found to have better stability against the perturbation of the second resonance. The situations of complete chaos are estimated in the \(e-i\) plane. By applying the maximum Lyapunov Characteristic Exponent (LCE), chaos is characterized quantitatively and similar conclusions can be achieved. This study is applied to three asteroids 1996 HW1, Vesta and Betulia, but the conclusions are not restricted to them.     

The strange case of the missing apocentric librators in the 3:2 resonance

From a comparison of the 2:1 and 3:2 resonances (in the asteroidal belt) two possible explanations to the absence of 3:2 apocentric librators are suggested. The first one is that such 3:2 resonant motion is dynamically unstable. The second interpretation requires the absence of nearcircular orbits originally at 4 AU. The latter view, if correct, is inconsistent with cosmogonic models which predict the original orbits of the asteroids to be nearly circular.     

Resonant periodic orbits in the exoplanetary systems

The planetary dynamics of 4/3, 3/2, 5/2, 3/1 and 4/1 mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to 60^° of mutual planetary inclination, but in most families, the stability does not exceed 20^°–30^°, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the 4/3 and 5/2 resonance, respectively.     

Chaotic zone boundary for low free eccentricity particles near an eccentric planet

We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Through collisionless particle integration, we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semimajor axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7th law has been previously explained by estimating the semimajor at which the first-order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first-order corotation resonances that have strength proportional to the planet's eccentricity. However, we find that the corotation resonance width at low free eccentricity is small; also the first-order resonance width at zero free eccentricity is the same as that for a zero-eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.     

The occurrence of high-order resonances and Kozai mechanism in the scattered disk

By means of numerical methods we explore the relevance of the high-order exterior mean motion resonances (MMR) with Neptune that a scattered disk object (SDO) can experience in its diffusion to the Oort cloud. Using a numerical method for estimate the strength of these resonances we show that high-eccentricity or high-inclination resonant orbits should have evident dynamical effects. We investigate the properties of the Kozai mechanism (KM) for non-resonant SDO's and the conditions that generate the KM inside a MMR associated with substantial changes in eccentricity and inclination. We found that the KM inside a MMR is typical for SDO's with Pluto-like or greater inclinations and is generated by the oscillation of ω inside the mixed (e,i) resonant terms of the disturbing function. A SDO diffusing to the Oort cloud should experience temporary captures in MMR, preferably of the type 1:N, and when evolving inside a MMR and experiencing the KM it can reach regions where the strength of the resonance drops and consequently there is a possibility of being decoupled from the resonance generating by this way a long-lived high-perihelion scattered disk object (HPSDO).     

On the orbital dependence of the asteroidal collision process

Collision of asteroids with the main-belt asteroid population is considered with the effect of the impact kinetic energy taken into account. It is found that objects in eccentric orbits have a larger probability of destructive collision as compared to objects in orbits with mean values of eccentricity (e = 0.15) and inclination (i = 10°); also orbits with small semimajor axes (a ≈ 2.3 AU) are found to have peak values of the probability of destructive collision.     

Numerical investigation of motions of resonant asteroids in the three-dimensional space

Motions of asteroids in mean motion resonances with Jupiter are studied in three-dimensional space. Orbital changes of fictitious asteroids in the Kirkwood gaps are calculated by numerical integrations for 10⁵ – 10⁶ years. The main results are as follows: (1) There are various motions of resonant asteroids, and some of them are very complicated and chaotic and others are regular. (2) The eccentricity of some asteroids becomes very large, and the variation of the inclination is large while the eccentricity is large. (3) In the 3:1 resonance, there is a long periodic change in the variation of the inclination, when (7 : ) is a simple ratio (7: longitude of perihelion, : longitude of node). (4) In the 7:3 resonance, the variation of the inclination of some resonant asteroids is so large that prograde motion becomes retrograde. Some asteroids in the 7:3 resonance can collide with the Sun as well as with the inner planets.