The scattered disk population as a source of Oort cloud comets: evaluation of its current and past role in populating the Oort cloud

We have integrated the orbits of the 76 scattered disk objects (SDOs), discovered through the end of 2002, plus 399 clones for 5 Gyr to study their dynamical evolution and the probability of falling in one of the following end states: reaching Jupiter's influence zone, hyperbolic ejection, or transfer to the Oort cloud. We find that nearly 50% of the SDOs are transferred to the Oort cloud (i.e., they reach heliocentric distances greater than 20,000 AU in a barycentric elliptical orbit), from which about 60% have their perihelia beyond Neptune's orbit (31 AU<q<36 AU) at the moment of reaching the Oort cloud. This shows that Neptune acts as a dynamical barrier, scattering most of the bodies to near-parabolic orbits before they can approach or cross Neptune's orbit in non-resonant orbits (that may allow their transfer to the planetary region as Centaurs via close encounters with Neptune). Consequently, Neptune's dynamical barrier greatly favors insertion in the Oort cloud at the expense of the other end states mentioned above. We found that the current rate of SDOs with radii R>1 km incorporated into the Oort cloud is about 5 yr⁻¹, which might be a non-negligible fraction of comet losses from the Oort cloud (probably around or even above 10%). Therefore, we conclude that the Oort cloud may have experienced and may be even experiencing a significant renovation of its population, and that the trans-neptunian belt—via the scattered disk—may be the main feeding source.     

The origin and distribution of the Centaur population

We analyze the Centaur population as a group of objects with perihelion distances (q) of less than 30 AU and heliocentric distances outside the orbit of Jupiter, formed by objects entering this region from the Scattered Disk (SD). We perform a numerical integration of 95 real Scattered Disk Objects (SDOs) extracted from the Minor Planet Center database and of 905 synthetic SDOs compensating for observational biases. SDOs have in the Centaur zone a mean lifetime of 72 Myr, though this number falls with a decrease of q. After this incursion, 30% of them enter the zone interior to Jupiter's orbit. We find that the contribution to the Centaur population from the SD gives a total of ∼2.8×¹⁰⁸ Centaurs with a radius R>1 km. We also propose a model for the intrinsic distribution of orbital elements of Centaurs and their distance and apparent magnitude distribution.     

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 Dynamics of Objects in the Inner Edgeworth–Kuiper Belt

Objects in 3:2 mean motion resonance with Neptune are protected from close encounters with Neptune by the resonance. Bodies in orbits with semi-major axis between 39.5 and about 42 AU are not protected by the resonance; indeed due to overlapping secular resonances, the eccentricities of orbits in this region are driven up so that a close encounter with Neptune becomes inevitable. It is thus expected that such orbits are unstable. The list of known Trans-Neptunian objects shows a deficiency in the number of objects in this gap compared to the 43–50 AU region, but the gap is not empty. We numerically integrate models for the initial population in the gap, and also all known objects over the age of the Solar System to determine what fraction can survive. We find that this fraction is significantly less than the ratio of the population in the gap to that in the main belt, suggesting that some mechanism must exist to introduce new members into the gap. By looking at the evolution of the test body orbits, we also determine the manner in which they are lost. Though all have close encounters with Neptune, in most cases this does not lead to ejection from the Solar System, but rather to a reduced perihelion distance causing close encounters with some or all of the other giant planets before being eventually lost from the system, with Saturn appearing to be the cause of the ejection of most of the objects.     

Dynamical classification of trans-neptunian objects: Probing their origin, evolution, and interrelation

The orbital structure of trans-neptunian objects (TNOs) in the trans-neptunian belt (Edgeworth-Kuiper belt) and scattered disk provides important clues to understand the origin and evolution of the Solar System. To better characterize these populations, we performed computer simulations of currently observed objects using long-arc orbits and several thousands of clones. Our preliminary analysis identified 622 TNOs, and 65 non-resonant objects whose orbits penetrate that of at least one of the giant planets within 1 Myr (the centaurs). In addition, we identified 196 TNOs locked in resonances with Neptune, which, sorted by distance from the Sun, are 1:1 (Neptune trojans), 5:4, 4:3, 11:8, 3:2, 18:11, 5:3, 12:7, 19:11, 7:4, 9:5, 11:6, 2:1, 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4. Kozai resonant TNOs are found inside the 3:2, 5:3, 7:4, and 2:1 resonances. We present detailed general features for the resonant populations (i.e., libration amplitude angles, libration centers, Kozai libration amplitudes, etc.). Taking together the simulations of Lykawka and Mukai [Lykawka, P.S., Mukai, T., 2007. Icarus 186, 331-341], an improved classification scheme is presented revealing five main classes: centaurs, resonant, scattered, detached and classical TNOs. Scattered and detached TNOs (non-resonant) have q (perihelion distance) <37 AU and q>40 AU, respectively. TNOs with 37 AU<q<40 AU occupy an intermediate region where both classes coexist. Thus, there are no clear boundaries between the scattered and detached regions. We also securely identified a total of 9 detached TNOs by using 4-5 Gyr orbital integrations. Classical objects are non-resonant TNOs usually divided into cold and hot populations. Their boundaries are as follows: cold classical TNOs (i?5°) are located at 37 AU<a<40 AU (q>37 AU) and 42 AU<a<47.5 AU (q>38 AU), and hot classical TNOs (i>5°) occupy orbits with 37 AU<a<47.5 AU (q>37 AU). However, a more firm classification is found with i>10° for hot classical TNOs. Lastly, we discuss some implications of our classification scheme comparing all TNOs with our model and other past models.     

Evidence for an Extended Scattered Disk

By telescopic tracking, we have established that the transneptunian object (TNO) 2000 CR₁₀₅ has a semimajor axis of 220±1 AU and perihelion distance of 44.14±0.02 AU, beyond the domain which has heretofore been associated with the “scattered disk” of Kuiper Belt objects interacting via gravitational encounters with Neptune. We have also firmly established that the TNO 1995 TL₈ has a high perihelion (of 40.08±0.02 AU). These objects, and two other recent discoveries which appear to have perihelia outside 40 AU, have probably been placed on these orbits by a gravitational interaction which is not strong gravitational scattering off of any of the giant planets on their current orbits. Their existence may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced these objects, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of TNOs in an “extended scattered disk” with perihelia above the previously suggested 38 AU boundary. The total population is difficult to estimate due to the ease with which such objects would have been lost. This illustrates the great value of frequent and well time-sampled recovery observations of trans-neptunian objects within their discovery opposition.     

Near-Sun asteroids

As follows from dynamical studies, in the course of evolution, most near-Earth objects reach orbits with small perihelion distances. Changes of the asteroids in the vicinity of the Sun should play a key role in forming the physical properties, size distribution, and dynamical features of the near-Earth objects. Only seven of the discovered asteroids are currently moving along orbits with perihelion distances q < 0.1 AU. However, due to the Kozai–Lidov secular perturbations, the asteroids, having recently passed near the Sun, could by now have moved to orbits farther from the Sun. In this study, we found asteroids that have been recently orbiting with perihelion distances q < 0.1 AU. Asteroids may be on such orbits for hundreds to tens of thousands of years. To carry out astrophysical observations of such objects is a high priority.     

The Dynamics of Low-Perihelion Meteoroid Streams

The Canadian Meteor Orbit Radar (CMOR) has collected information on a number of weak meteor showers that have not been well characterized in the literature. A subsample of these showers (1) do not show a strong orbital resemblance to any known comets or asteroids, (2) have highly inclined orbits, (3) are at low perihelion distances ( AU) and (4) are at small semimajor axes (<2 AU). Though one might conclude that the absence of a parent object could be the result of its disruption, it is unclear how this relatively inaccessible (dynamically speaking) region of phase space might have been populated by parents in the first place. It will be shown that the Kozai secular resonance and/or Poynting–Robertson drag can modify meteor stream orbits rapidly (on time scales comparable to a precession cycle) and may be responsible for placing some of these streams into their current locations. These same effects are also argued to act on these streams so as to contribute to the high-ecliptic latitude north and south toroidal sporadic meteor sources. There remain some differences between the simple model results presented here and observations, but there may be no need to invoke a substantial population of high-inclination parents for the observed high-inclination meteoroid streams with small perihelion distances.     

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).     

Secular Dynamics of Asteroids in the Inner Solar System

In the last three years we have carried out numerical and semi-analytical studies on the secular dynamical mechanisms in the region (semimajor axis a < 2 AU) where the NEA orbits evolve. Our numerical integrations (over a time span of a few Myr) have shown that: (i) the linear secular resonances with both the inner and the outer planets may play an important role in the dynamical evolution of NEAs; (ii) the apsidal secular resonance with Mars could provide an important dynamical transport mechanism by which asteroids in the Mars-crossing region eventually achieve Earth-crossing orbits; (iii) in this region, due to the interaction with the terrestrial planets, the Kozai resonance can occur at small inclinations, with the argument of perihelion ω librating around 0° or 180°, providing a temporary protection mechanism against close approaches to the planets. The location of the linear secular resonances in this zone has also been obtained by an automatic procedure using a semi-numerical method valid for all values of the inclinations and eccentricities of the small bodies, and also in the case of libration of the argument of perihelion. A map of the secular resonances in the (a, i) plane shows — in agreement with the numerical integrations — that all the resonances with the terrestrial and giant planets are present, and also that some of them overlap. Thus the way is now open to fully take into account secular resonances in modelling the dynamical evolution of NEAs. This revised version was published online in July 2006 with corrections to the Cover Date.     

Study and application of the resonant secular dynamics beyond Neptune

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (Celest Mech Dyn Astron, doi: 10.1007/s10569-016-9700-5, 2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or \(2012\text {VP}_{113}\), which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at \(a\in [100;300]\) AU, \(q\in [50;70]\) AU and \(I\in [30;50]^{\circ }\). In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.     

Origin of scattered disk resonant TNOs: Evidence for an ancient excited Kuiper belt of 50 AU radius

Resonance occupation of trans-neptunian objects (TNOs) in the scattered disk (>48 AU) was investigated by integrating the orbits of 85 observed members for 4 Gyr. Twenty seven TNOs were locked in the 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4 resonances. We then explored mechanisms for the origin of the resonant structure in the scattered disk, in particular the long-term 9:4, 5:2, and 8:3 resonant TNOs (median 4 Gyr), by performing large scale simulations involving Neptune scattering and planetary migration over an initially excited planetesimals disk (wide range of eccentricities and inclinations). To explain the formation of Gyr-resident populations in such distant resonances, our results suggest the existence of a primordial planetesimal disk of at least 45-50 AU radius that suffered a dynamical perturbation leading to 0.1-0.3 or greater eccentricities and a range of inclinations up to ∼20° during early stages of the Solar System history, before planetary migration.     

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.     

A Model for the Common Origin of Jupiter Family and Halley Type Comets

A numerical simulation of the Oort cloud is used to explain the observed orbital distributions and numbers of Jupiter-family (JF) and Halley-type (HT) short-period (SP) comets. Comets are given initial orbits with perihelion distances between 5 and 36 au, and evolve under planetary, stellar and Galactic perturbations for 4.5 Gyr. This process leads to the formation of an Oort cloud (which we define as the region of semimajor axes a > 1,000 au), and to a flux of cometary bodies from the Oort cloud returning to the planetary region at the present epoch. The results are consistent with the dynamical characteristics of SP comets and other observed cometary populations: the near-parabolic flux, Centaurs, and high-eccentricity trans-Neptunian objects. To achieve this consistency with observations, the model requires that the number of comets versus initial perihelion distance is concentrated towards the outer planetary region. Moreover, the mean physical lifetime of observable comets in the inner planetary region (q < 2.5 au) at the present epoch should be an increasing function of the comets’ initial perihelion distances. Virtually all observed HT comets and nearly half of observed JF comets come from the Oort cloud, and initially (4.5 Gyr ago) from orbits concentrated near the outer planetary region. Comets that have been in the Oort cloud also return to the Centaur (5 < q < 28 au, a < 1,000 au) and near-Neptune high-eccentricity regions. Such objects with perihelia near Neptune are hard to discover, but Centaurs with characteristics predicted by the model (e.g. large semimajor axes, above 60 au, or high inclinations, above 40°) are increasingly being found by observers. The model provides a unified picture for the origin of JF and HT comets. It predicts that the mean physical lifetime of all comets in the region q < 1.5 au is less than ～200 revolutions.     

Analytical Solution of the Kozai Resonance and its Application

When Kozai (1962) studied the secular resonance of asteroids, he found the so-called Kozai resonance and expressed the analytical solution with the use of Weierstrass ℘. Here we discuss the case where the disturber is outside a disturbed body and give the analytical solution of the eccentricity, the inclination and the argument of pericenter with the use of the Jacobi elliptic functions, which are more familiar than the Weierstrass ℘. Then we derive the Fourier expansion of the longitude of node and the mean anomaly. The analytical expressions obtained here can be used for any value of the eccentricity and the inclination. Finally we applied these analytical expressions to several dynamical systems – Nereid, that is a highly eccentric satellite of Neptune, and newly discovered retrograde satellites of Uranus. This revised version was published online in July 2006 with corrections to the Cover Date.     

The evidence of an early stellar encounter in Edgeworth-Kuiper belt

We investigate the influence of a stellar fly-by encounter on the Edgeworth-Kuiper belt objects through numerical orbital calculations, in order to explain both mass depletion and high orbital inclinations of the classical Edgeworth-Kuiper belt (CEKB) objects, which have semimajor axis of 42-48 AU and perihelia beyond 35 AU. The observationally inferred total mass of the CEKB is ∼1/10 Earth masses, which is only ∼0.02 of that extrapolated from the minimum-mass solar nebula model. The CEKB consists of bimodal population: “hot population” with inclinations i?0.2-0.6 radians and “cold population” with i?0.1. The observationally suggested difference in size and color of objects between the two populations may imply different origins of the two populations. We find that both the depletion of solid materials in the CEKB and the formation of the hot population are accounted for by a single close stellar encounter with pericenter distance of 80-100 AU and inclination relative to the initial protoplanetary disk ?50°-70°. Such a stellar encounter highly pumps up eccentricities of most objects in the CEKB and then their perihelia migrate within 35 AU. These objects would be removed by Neptune's perturbations after Neptune is formed at or migrates to the current position (30 AU). Less than 10% of the original objects remain in stable orbits with small eccentricities and perihelion distances larger than 35 AU, in the CEKB, which is consistent with the observation. We find that i of the remaining objects are as large as that of the observed hot population. The only problem is how to stop Neptune's migration at ∼30 AU, which is addressed in a separate paper. The depletion by the stellar encounter extends deeply into ∼30-35 AU, which provides the basis of the formation model for the cold population through Neptune's outward migration by Levison and Morbidelli (2003, Nature, 426, 419-421). The combination of our model with Levison and Morbidelli's model could consistently explain the mass depletion, truncation at 50 AU, bimodal distribution in i, and differences in size and color between the hot and the cold populations in the CEKB.     

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.     

Discovery Of Mg-Rich Pyroxenes In Comet C/1995 O1 (Hale-Bopp): Pristine Grains Revealed At Perihelion

The NASA Ames HIFOGS spectrometer observed comet C/1995 O1 (Hale-Bopp) at epochs including 96 Oct 7–14 UT (2.8 AU), 97 Feb 14–15 UT (1.2 AU), 97 Apr 11 UT (0.93 AU), and 97 Jun 22, 25 UT (1.7 AU). The HIFOGS 7.5–13.5 μm spectrophotometry (R = 360 - 180) of the silicate feature at 2.8 AU is identical in shape to the ISO SWS spectra of comet Hale-Bopp (Crovisier et al., 1997); the strong 11.2 μm peak in the structured silicate feature is identified as olivine. Upon close passage to the sun, the HIFOGS spectra at 1.2 AU and 0.93 AU reveals strong peaks at 9.3 μm and 10.0 μm. The post-perihelion 10 μm silicate feature at 1.7 AU is weaker but has nearly the same shape as the pre-perihelion spectra at 1.2 AU, reverting to its pre-perihelion shape: there is no change in the dust chemistry by close passage to the sun. The appearance of the strong peaks at 9.3 μm and 10.0 μm at r_h ≲ 1.7 AU is attributed to the rise in the contribution of pryoxenes (clino-pyroxene and orthopyroxene crystals) to the shape of the feature, and leads to the hypothesis that the pyroxenes are significantly cooler than the olivines. The pyroxenes are radiating on the Wien side of the blackbody at 2.8 AU and transition to the Rayleigh-Jeans tail of the blackbody upon closer approach to the Sun. Composite fits to the observed 10 μm silicate features using IDPs and laboratory minerals shows that a good empirical fit to the spectra is obtained when the pryoxenes are about 150 K cooler than the olivines. The pyroxenes, because they are cooler and contribute signficantly at perihelion, are more abundant than the olivines. The perihelion temperature of the pyroxenes implies that the pyroxenes are more Mg-rich than the other minerals including the olivines, amorphous olivines, and amorphous pyroxenes. The PUMA-1 flyby measurements of comet P/Halley also indicated an overabundance of Mg-rich pryoxenes compared to olivines. Comet Hale-Bopp's pyroxenes are similar to pyroxere IDPs from the ’Spray‘ class, known for their D-richness and their unaltered morphologies: Hale-Bopp's Mg-rich pyroxenes may be pristine relic ISM grains. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Simple Method for Numerical Calculations of the Evolution of Orbits of Near-Earth Asteroids

A simple method for numerical integration of the equations of motion of small bodies of the Solar System is proposed, which is especially efficient in studying the orbits with small perihelion distances. The evolution of orbits of 121 numbered asteroids with perihelion distances q < 1.2 AU is investigated over the time interval of years 2000–2100 with allowance made for the gravitational influence of nine planets and three largest asteroids. The circumstances of close encounters of asteroids with the Earth and other terrestrial planets are presented.