Exploring the 7:4 mean motion resonance—II: Scattering evolutionary paths and resonance sticking

In our preliminary study, we have investigated basic properties and dynamical evolution of classical TNOs around the 7:4 mean motion resonance with Neptune (a∼43.7 AU), motivated by observational evidences that apparently present irregular features near this resonance (see [Lykawka and Mukai, 2005a. Exploring the 7:4 mean motion resonance—I. Dynamical evolution of classical trans-Neptunian objects (TNOs). Space Planet. Sci. 53, 1175-1187]; hereafter “Paper I”). In this paper, we aim to explore the dynamical long-term evolution in the scattered disk (but not its early formation) based on the computer simulations performed in Paper I together with extra computations. Specifically, we integrated the orbital motion of test particles (totalizing a bit more than 10,000) placed around the 7:4 mean motion resonance under the effect of the four giant planets for the age of the Solar System. In order to investigate chaotic diffusion, we also conducted a special simulation with on-line computation of proper elements following tracks in phase space over 4-5 Gyr. We found that: (1) A few percent (1-2%) of the test particles survived in the scattered disk with direct influence of other Neptunian mean motion resonances, indicating that resonance sticking is an extremely common phenomenon and that it helps to enhance scattered objects longevity. (2) In the same region, the so-called extended scattered TNOs are able to form via very long resonance trapping under certain conditions. Namely, if the body spends more than about 80% of its dynamical lifetime trapped in mean motion resonance(s) and there is the action of a k+1 or (k+2)/2 mean motion resonance (e.g., external mean motion resonances with Neptune described as (j+k)/j with j=1 and 2, respectively). According to this hypothetical mechanism, 5-15% of current scattered TNOs would possess thus probably constituting a significant part of the extended scattered disk. (3) Moreover, considering hot orbital initial conditions, it is likely that the trans-Neptunian belt (or Edgeworth-Kuiper belt) has been providing members to the scattered disk, so that scattered TNOs observed today would consist of primordial scattered bodies mixed with TNOs that came from unstable regions of the trans-Neptunian belt in the past.Considering the three points together, our results demonstrated that the scattered disk has been evolving continuously since early times until present.     

两个可能存在大质量类冥王星的区域

1992年以来，在海王星外的太阳系发现了近千个小天体，称为Kuiper带天体(KBO)或Edgeworth—Kuiper带天体，其中有一部分偏心率和倾角较大的小天体与海王星之间存在3：2平运动共振，轨道特征类似冥王星，命名为类冥王星，自KBO发现以来，天文学家们进行了多次小天区的搜索，发现了几个质量较大的KBO，通过数值计算，在轨道参数空间发现了两个和冥王星一样同时具有3种共振的区域，在这两个区域里的小天体既避免了海王星的强摄动又不会与冥王星密切交会，轨道非常稳定，因此有可能在其中发现质量较大的类冥王星。     

Resonant motion of trans-Neptunian objects in high-eccentricity orbits

A substantial fraction of the Edgeworth-Kuiper belt objects are presently known to move in resonance with Neptune (the principal commensurabilities are 1/2, 3/5, 2/3, and 3/4). We have found that many of the distant (with orbital semimajor axes a > 50 AU) trans-Neptunian objects (TNOs) also execute resonant motions. Our investigation is based on symplectic integrations of the equations of motion for all multiple-opposition TNOs with a > 50 AU with allowance made for the uncertainties in their initial orbits. Librations near such commensurabilities with Neptune as 4/9, 3/7, 5/12, 2/5, 3/8, 4/27, and others have been found. The largest number of distant TNOs move near the 2/5 resonance with Neptune: 12 objects librate with a probability higher than 0.75. The multiplicity of objects moving in 2/5 resonance and the longterm stability of their librations suggest that this group of resonant objects was formed at early formation stages of the Solar system. For most of the other resonant objects, the librations are temporary. We also show the importance of asymmetric resonances in the large changes in TNO perihelion distances.     

The Kuiper Belt luminosity function from mR= 22 to 25

In summer 1999, we performed a survey optimized for the discovery of irregular satellites of Uranus and Neptune. We imaged 11.85 deg² of sky and discovered 66 new outer Solar system objects (not counting the three new Uranian satellites). Given the very short orbital arcs of our observations, only the heliocentric distance can be reliably determined. We were able to model the radial distribution of trans-Neptunian objects (TNOs). Our data support the idea of a strong depletion in the surface density beyond 45 au.
After fully characterizing this survey's detection efficiency as a function of object magnitude and rate of motion, we find that the apparent luminosity function of the trans-Neptunian region in the range   m_R = 22–25  is steep with a best-fitting cumulative power-law index of  α≃ 0.76  with one object per deg² estimated at magnitude   R _o= 23.3  . This steep slope, corresponding to a differential size index of   q ≃ 4.8  , agrees with other older and more recent analyses for the luminosity function brighter than 25 mag. A double power-law fit to the new data set turns out to be statistically unwarrented; this large and homogeneous data set provides no evidence for a break in the power-law slope, which must eventually occur if the Bernstein et al. sky density measurements are correct.     

Discovery of a New Class of Predicted Resonance Objects beyond Jupiter

In the regions of mean diurnal motions between the orbits of Jupiter and Saturn, predicted earlier by the authors, five asteroids have been discovered that move in 1:2 and 2:3 Lindblad orbital resonances with Jupiter (external orbital commensurability) and in 2:1 resonance with Saturn (internal version of commensurability). In addition to this, in the precalculated stable resonance zones between the giant planets Saturn and Uranus, three objects have been found that possess third-order (2:5) orbital commensurability with Saturn; nine objects have been discovered between the orbits of Uranus and Neptune, whose mean motions are in 1:3 and 1:4 orbital resonances with Saturn, and more than 200 libration-stable objects, linked by lower-order orbital resonances with Neptune and Uranus have been found in the Kuiper belt.     

The scattered trans-Neptunian object 1998 XY95

On 1998 December 12 a new trans-Neptunian object, 1998 XY95, was discovered as part of a deep search. Recent observations of this object have placed it amongst the class of objects known as the scattered trans-Neptunian objects (TNOs). A total of 39 CCD images of 1998 XY95 were taken over two nights, and these were used to search for a light curve, but no significant periodicity was found. An examination of the possible orbital evolution gives no indication of how it may have arrived on its present orbit. The current best-fitting orbit is unstable, but remains within a band of semi-major axis approximately 2 au wide. The bottom of this band is due to 3:1 mean motion resonance with Neptune, while the reason for the top of the band remains unclear.     

Enrichment in volatiles in the giant planets of the Solar System

We propose an interpretation of the enrichments in volatiles observed in the four giant planets with respect to the solar abundance. It is based on the assumption that volatiles were trapped in the form of solid clathrate hydrates and incorporated in planetesimals embedded in the feeding zones of each of the four giant planets. The mass of trapped volatiles is then held constant with time. The mass of hydrogen and of not trapped gaseous species continuously decreased with time until the formation of the planet was completed, resulting in an increase in the ratio of the mass of trapped volatiles to the mass of hydrogen (Gautier et al., Astrophys. J. 550 (2001) L227). The efficiency of the clathration depends upon the amount of ice available in the early feeding zone. The quasi-uniform enrichment in Ar, Kr, Xe, C, N, and S observed in Jupiter is reproduced because all volatiles were trapped. The non-uniform enrichment observed in C, N and S in Saturn is due to the fact that CH₄, NH₃, and H₂S were trapped but not CO and N₂. The non-uniform enrichment in C, N and S in Uranus and Neptune results from the trapping of CH₄, CO, NH₃ and H₂S, while N₂ was not trapped. Our scenario permits us to interpret the strongly oversolar sulfur abundance inferred by various modelers to be present in Saturn, Uranus and Neptune for reproducing the microwave spectra of the three planets. Abundances of Ar, Kr and Xe in these three are also predicted. Only Xe is expected to be substantially oversolar. The large enrichment in oxygen in Neptune with respect to the solar abundance, calculated by Lodders and Fegley (Icarus 112 (1994) 368) from the detection of CO in the upper troposphere of the planet, is consistent with the trapping of volatiles by clathration. The upper limit of CO in Uranus does not exclude that this process also occurred in Uranus.     

The web of three-planet resonances in the outer Solar System: II. A source of orbital instability for Uranus and Neptune

The motion of the giant planets from Jupiter to Neptune is chaotic with Lyapunov time of approximately 10 Myr. A recent theory explains the presence of this chaos with three-planet mean-motion resonances, i.e. resonances among the orbital periods of at least three planets. We find that the distribution of these resonances with respect to the semi-major axes of all the planets is compatible with orbital instability. In particular, they overlap in a region of 10⁻³ AU with respect to the variation of the semi-major axes of Uranus and Neptune. Fictitious planetary systems with initial conditions in this region can undergo systematic variations of semi-major axes. The true Solar System is marginally in this region, and Uranus and Neptune undergo very slow systematic variations of semi-major axes with speed of order 10⁻⁴ AU/Gyr.     

Long term dynamical evolution and classification of classical TNOs

Classical trans-Neptunian objects (TNOs) are believed to represent the most dynamically pristine population in the trans-Neptunian belt (TNB) offering unprecedented clues about the formation of our Solar System. The long term dynamical evolution of classical TNOs was investigated using extensive simulations. We followed the evolution of more than 17000 particles with a wide range of initial conditions taking into account the perturbations from the four giant planets for 4 Gyr. The evolution of objects in the classical region is dependent on both their inclination and semimajor axes, with the inner (a<45 AU) and outer regions (a>45 AU) evolving differently. The reason is the influence of overlapping secular resonances with Uranus and Neptune (40–42 AU) and the 5:3 (a∼ ∼42.3 AU), 7:4 (a∼ ∼43.7 AU), 9:5 (a∼ ∼44.5 AU) and 11:6 (a∼ ∼ 45.0 AU) mean motion resonances strongly sculpting the inner region, while in the outer region only the 2:1 mean motion resonance (a∼ ∼47.7 AU) causes important perturbations. In particular, we found: (a) A substantial erosion of low-i bodies (i<10°) in the inner region caused by the secular resonances, except those objects that remained protected inside mean motion resonances which survived for billion of years; (b) An optimal stable region located at 45 AU<a<47 AU, q>40 AU and i>5° free of major perturbations; (c) Better defined boundaries for the classical region: 42–47.5 AU (q>38 AU) for cold classical TNOs and 40–47.5 AU (q>35 AU) for hot ones, with i=4.5° as the best threshold to distinguish between both populations; (d) The high inclination TNOs seen in the 40–42 AU region reflect their initial conditions. Therefore they should be classified as hot classical TNOs. Lastly, we report a good match between our results and observations, indicating that the former can provide explanations and predictions for the orbital structure in the classical region.     

The origin of TNO 2004 XR190 as a primordial scattered object

Numerical integrations of the equations of motion of the giant planets and scattering particles show that there is a possible orbital itinerary that a particle may follow from a scattering mode up to a stable position near the orbit of 2004 XR₁₉₀. This orbital evolution requires that the particle gets trapped in a mean motion resonance with Neptune coupled with the Kozai resonance. Imposing migration on Neptune while a particle is experiencing both resonances can entail an escape from resonance at a low particle’s eccentricity. This eccentricity and the associated inclination are always similar to those of 2004 XR₁₉₀. I conclude that 2004 XR₁₉₀ was most likely a scattered object that went through those resonance processes and was eventually deposited at its current position. By the same argument, it is expected that there must exist several other objects with similar semimajor axis, eccentricity and inclination as those of 2004 XR₁₉₀.     

Reducing the probability of capture into resonance

A migrating planet can capture planetesimals into mean motion resonances. However, resonant trapping can be prevented when the drift or migration rate is sufficiently high. Using a simple Hamiltonian system for first- and second-order resonances, we explore how the capture probability depends on the order of the resonance, drift rate and initial particle eccentricity. We present scaling factors as a function of the planet mass and resonance strength to estimate the planetary migration rate above which the capture probability drops to less than half. Applying our framework to multiple extrasolar planetary systems that have two planets locked in resonance, we estimate lower limits for the outer planet's migration rate, allowing resonance capture of the inner planet.
Mean motion resonances are comprised of multiple resonant subterms. We find that the corotation subterm can reduce the probability of capture when the planet eccentricity is above a critical value. We present factors that can be used to estimate this critical planet eccentricity. Applying our framework to the migration of Neptune, we find that Neptune's eccentricity is near the critical value that would make its 2 : 1 resonance fail to capture twotinos. The capture probability is affected by the separation between resonant subterms and so is also a function of the precession rates of the longitudes of periapse of both planet and particle near resonance.     

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 Dynamical Evolution of the Atens

In this investigation the orbits of 21 Atens (semimajor axes smaller than the Earth) are studied with the aid of numerical integrations over the time interval of one million years. The dynamical model was a 6-body Solar System, where the perturbations of Uranus and Neptune were ignored, and where Mercury's mass was added to the Sun's mass. Thus mean motion resonances, secular resonances and the Kozai resonance were fully taken into account. The evolution of the semimajor axes shows the typical step function like pattern which we know also from comets although some Atens have a very fuzzy development of the orbital elements, and some of them stay in a mean motion resonance for very long time. The evolution from Atens to Apollos (with semimajor axes larger than the Earth) and vice versa is also a phenomenon which we could observe. The main goal was the study of encounters of the Atens with the Earth and Venus. We found out that Venus encounters occur somewhat more often than Earth encounters (approximately one within the distance Earth - Moon every 40000 years with Venus, one every 50000 years with the Earth). This revised version was published online in July 2006 with corrections to the Cover Date.     

Resonant and Secular Families of the Kuiper Belt

We review ongoing efforts to identify occupants of mean-motion resonances(MMRs) and collisional families in the Edgeworth–Kuiper belt. Directintegrations of trajectories of Kuiper belt objects (KBOs) reveal the 1:1(Trojan), 5:4, 4:3, 3:2 (Plutino), 5:3, 7:4, 9:5, 2:1 (Twotino), and 5:2 MMRsto be inhabited. Apart from the Trojan, resonant KBOs typically have largeorbital eccentricities and inclinations. The observed pattern of resonanceoccupation is consistent with resonant capture and adiabatic excitation bya migratory Neptune; however, the dynamically cold initial conditions priorto resonance sweeping that are typically assumed by migration simulationsare probably inadequate. Given the dynamically hot residents of the 5:2 MMRand the substantial inclinations observed in all exterior MMRs, a fraction ofthe primordial belt was likely dynamically pre-heated prior to resonancesweeping. A pre-heated population may have arisen as Neptune gravitationallyscattered objects into trans-Neptunian space. The spatial distribution of Twotinosoffers a unique diagnostic of Neptune's migration history. The Neptunian Trojanpopulation may rival the Jovian Trojan population, and the former's existence isargued to rule out violent orbital histories for Neptune. Finally, lowest-order seculartheory is applied to several hundred non-resonant KBOs with well-measured orbitsto update proposals of collisional families. No convincing family is detected.     

Mass of the Kuiper belt

The Kuiper belt includes tens of thousand of large bodies and millions of smaller objects. The main part of the belt objects is located in the annular zone between 39.4 and 47.8 au from the Sun; the boundaries correspond to the average distances for orbital resonances 3:2 and 2:1 with the motion of Neptune. One-dimensional, two-dimensional, and discrete rings to model the total gravitational attraction of numerous belt objects are considered. The discrete rotating model most correctly reflects the real interaction of bodies in the Solar system. The masses of the model rings were determined within EPM2017—the new version of ephemerides of planets and the Moon at IAA RAS—by fitting spacecraft ranging observations. The total mass of the Kuiper belt was calculated as the sum of the masses of the 31 largest trans-Neptunian objects directly included in the simultaneous integration and the estimated mass of the model of the discrete ring of TNO. The total mass is \((1.97 \pm 0.35)\times 10^{-2} \ m_{\oplus }\). The gravitational influence of the Kuiper belt on Jupiter, Saturn, Uranus, and Neptune exceeds at times the attraction of the hypothetical 9th planet with a mass of \(\sim 10 \ m_{\oplus }\) at the distances assumed for it. It is necessary to take into account the gravitational influence of the Kuiper belt when processing observations and only then to investigate residual discrepancies to discover a possible influence of a distant large planet.     

Resonance sticking in the scattered disk

We investigate the dynamical evolution of trans-neptunian objects (TNOs) in typical scattered disk orbits (scattered TNOs) by performing simulations using several thousand particles lying initially on Neptune-encountering orbits. We explore the role of resonance sticking in the scattered disk, a phenomenon characterized by multiple temporary resonance captures (‘resonances’ refers to external mean motion resonances with Neptune, which can be described in the form r:s, where the arguments r and s are integers). First, all scattered TNOs evolve through intermittent temporary resonance capture events and gravitational scattering by Neptune. Each scattered TNO experiences tens to hundreds of resonance captures over a period of 4 Gyr, which represents about 38% of the object's lifetime (mean value). Second, resonance sticking plays an important role at semimajor axes , where the great majority of such captures occurred. It is noteworthy that the stickiest (i.e., dominant) resonances in the scattered disk are located within this distance range and are those possessing the lowest argument s. This was evinced by r:1, r:2 and r:3 resonances, which played the greatest role during resonance sticking evolution, often leading to captures in several of their neighboring resonances. Finally, the timescales and likelihood of temporary resonance captures are roughly proportional to resonance strength. The dominance of low s resonances is also related to the latter. In sum, resonance sticking has an important impact on the evolution of scattered TNOs, contributing significantly to the longevity of these objects.     

Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals

The final stage of the accretion of Uranus and Neptune is numerically investigated. The four Jovian planets are considered with Jupiter and Saturn assumed to have reached their present sizes, whereas Uranus and Neptune are taken with initial masses 0.2 of their present ones. Allowance is made for the orbital variation of the Jovian planets due to the exchange of angular momentum with interacting bodies (“planetesimals”). Two possible effects that may have contributed to the accretion of Uranus and Neptune are incorporated in our model: (1) an enlarged cross section for accretion of incoming planetesimals due to the presence of extended gaseous envelopes and/or circumplanetary swarms of bodies; and (2) intermediate protoplanets in mid-range orbits between the Jovian planets. Significant radial displacements are found for Uranus and Neptune during their accretion and scattering of planetesimals. The orbital angular momentum budgets of Neptune, Uranus, and Saturn turn out to be positive; i.e., they on average gain orbital angular momentum in their interactions with planetesimals and hence they are displaced outwardly. Instead, Jupiter as the main ejector of bodies loses orbital angular momentum so it moves sunward. The gravitational stirring of planetesimals caused by the introduction of intermediate protoplanets has the effect that additional solid matter is injected into the accretion zones of Uranus and Neptune. For moderate enlargements of the radius of the accretion cross section (2–4 times), the accretion time scale of Uranus and Neptune are found to be a few 10⁸ years and the initial amount of solid material required to form them of a few times their present masses. Given the crucial role played by the size of the accretion cross section, questions as to when Uranus and Neptune acquired their gaseous envelopes, when the envelopes collapsed onto the solid cores, and how massive they were are essential in defining the efficiency and time scale of accretion of the two outer Jovian planets.     

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.     

Simulations for the Original Distribution of KBOs

Using the N-body dynamical model that includes the sun, the 8 planets, Pluto, UB313 and massless particles, we simulate the orbital evolution of 551 Kuiper Belt Objects (KBOs) with known parameters. The initial conditions of the simulations are the currently observed orbital parameters. The integration backtracks from now to -10×10⁸ yr. The results show that about 10×10⁸ years ago, more than 1/3 of the presently observed KBOs resided in the region of the present Kuiper main belt, a few were located inside the Neptune orbit, and the rest were beyond 50AU; and that about 4.5×10⁸ years ago, all the objects in the Kuiper main belt exhibited a rather good normal distribution, without so many objects concentrated in the Neptune's 3:2 resonance region, as at present time.     

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.