Confined chaotic motion in three-body resonances: trapping of trans-Neptunian material induced by gas-drag

Jiang & Yeh proposed gas-drag-induced resonant capture as a mechanism able to explain the dominant 3:2 resonance observed in the trans-Neptunian belt. Using a model of a disc–star–planet system they concluded that gaseous drag in a protoplanetary disc can trap trans-Neptunian object (TNO) embryos into the 3:2 resonance rather easily although it could not trap objects into the 2:1 resonance. Here we further investigate this scenario using numerical simulations within the context of the planar restricted four-body problem by including both present-day Uranus and Neptune. Our results show that mean motion and corotation resonances are possible and trapping into both the 3:2 and 2:1 resonances as well as other resonances is observed. The associated corotation centres may easily form larger planetesimals from smaller ones. Corotation resonances evolve into pure Lindblad resonances in a time-scale of 0.5 Myr. The non-linear corotation and mean motion resonances produced are very size selective. The 3:2 resonance is dominant for submetric particles but for larger particles the 2:1 resonance is stronger. In summary, our calculations show that confined chaotic motion around the resonances not only increases trapping efficiency but also the orbital eccentricities of the trapped material, modifying the relative abundance of trapped particles in different resonances. If we assume a more compact planetary system, instead of using the present-day values of the orbital elements of Uranus and Neptune, our results remain largely unchanged.     

Resonant Periodic Orbits of Trans-Neptunian Objects

We study two and three-dimensional resonant periodic orbits, usingthe model of the restricted three-body problem with the Sun andNeptune as primaries. The position and the stability character ofthe periodic orbits determine the structure of the phase space andthis will provide useful information on the stability and longterm evolution of trans-Neptunian objects. The circular planarmodel is used as the starting point. Families of periodic orbitsare computed at the exterior resonances 1/2, 2/3 and 3/4 withNeptune and these are used as a guide to select the energy levelsfor the computation of the Poincaré maps, so that all basicresonances are included in the study. Using the circular planarmodel as the basic model, we extend our study to more realisticmodels by considering an elliptic orbit of Neptune and introducingthe inclination of the orbit. Families of symmetric periodicorbits of the planar elliptic restricted three-body problem andthe three-dimensional problem are found. All these orbitsbifurcate from the families of periodic orbits of the planarcircular problem. The stability of all orbits is studied. Althoughthe resonant structure in the circular problem is similar for allresonances, the situation changes if the eccentricity of Neptuneor the inclination of the orbit is taken into account. All theseresults are combined to explain why in some resonances there aremany bodies and other resonances are empty.     

Drag-induced resonant capture in a multiplanet scenario: An application to 55 Cancri A

The binary star 55 Cancri harbors at least five giant planets. The discovery of a possible dust disk around the primary of the system was announced in 1998 although it was later dismissed as an observational artifact induced by the presence of three background galaxies. However, the possible existence of an asteroid belt beyond 6 AU from the primary could not be excluded. The actual properties of a hypothetical outer asteroid belt around 55 Cancri A are likely to be the result of the gravitational influence of the outermost planets and gas drag during the late stages of the formation of the planetary system. Gaseous drag within a protoplanetary disk in a multiplanet context can induce resonant capture of solid material rather easily, triggering the formation of asteroid belts similar to our own transneptunian belt. In this paper we investigate this scenario using numerical simulations within the framework of the planar restricted four-body problem to find possible stable debris locations or additional planets around 55 Cancri A. Our results indicate that, beyond 55 Cancri d, only the 1:2 and 1:3 mean-motion resonances may be possible although the details of gas drag-induced resonant trapping depend strongly on the size of the particles, with larger particles being preferentially trapped farther away from the host star. For a primary of mass 0.94 M and i=53° we find solid particles submitted to confined chaos at 10.14 and 13.22 AU with average eccentricities 0.14 and 0.20, respectively. The argument of pericentre of the trapped particles is found to librate around π/2 with the pericentres of the two resonant structures nearly aligned. The resonance responsible for the inner ring effectively halts the migration of the particles. The outer structure resembles a rosette or separatrix orbit. Our calculations suggest that hypothetical planets or asteroid belts in these locations would be, in principle, dynamically stable.     

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.     

Model orbits of asteroid 3040 Kozai

The orbital evolutions of the asteroid 3040 Kozai and model asteroids with similar orbits have been investigated. Their osculating orbits for an epoch 1991 December 10 were numerically integrated forward within the interval of 20,000 years, using a dynamical model of the solar system consisting of all inner planets, Jupiter, and Saturn.The orbit of the asteroid Kozai is stable. Its motion is affected only by long-period perturbations of planets. With change of the argument of perihelion of the asteroid Kozai, the evolution of the model asteroid orbits changes essentially, too. The model orbits with the argument of perihelion changed by the order of 10% show that asteroids with such orbital parameters may approach the Earth orbit, while asteroids with larger changes may even cross it, at least after 10,000 years. Long-term orbital evolution of asteroids with these orbital parameters is very sensitive on their angular elements.     

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.     

Numerical Simulation of the Orbital Evolution of Near-Earth Asteroids Close to Mean Motion Resonances

The dynamics of near-Earth asteroids near mean motion resonances with the Earth or other planets is considered. The probability domains of the motion of some near-Earth asteroids close to low-order resonances are presented. The investigations have been carried out by means of a numerical integration of differential equations, taking into account the perturbations from the major planets and the Moon. For each investigated object an ensemble of 100 test particles with orbital elements nearby those of the nominal orbit has been constructed and its evolution has been retraced over the time interval (–3000, +3000 years). The initial set of orbits has been generated on the basis of probable variations of the initial orbital elements obtained from the least square analysis of observations.This revised version was published online in October 2005 with corrections to the Cover Date.     

The absolute magnitude distribution of trans-Neptunian objects

It is shown that the known trans-Neptunian objects (TNOs) have an absolute magnitude distribution index that increases as a function of orbital perihelion distance. In no perihelion range is the TNO index the same as that found for known short-period comets. However, the fact that the median diameters of the known members of these two populations (220 and 2.9 km respectively) differ by a factor of about 75 means that very small TNOs and short-period comets might still be related.     

Origin of orbits of secondaries in the discovered trans-Neptunian binaries

The dependences of inclinations of orbits of secondaries in the discovered trans-Neptunian binaries on the distance between the primary and the secondary, on the eccentricity of orbits of the secondary around the primary, on the ratio of diameters of the secondary and the primary, and on the elements of heliocentric orbits of these binaries are studied. These dependences are interpreted using the model of formation of a satellite system in a collision of two rarefied condensations composed of dust and/or objects less than 1 m in diameter. It is assumed in this model that a satellite system forms in the process of compression of a condensation produced in such a collision. The model of formation of a satellite system in a collision of two condensations agrees with the results of observations: according to observational data, approximately 40% of trans-Neptunian binaries have a negative angular momentum relative to their centers of mass.     

On the periodically evolving orbits in the singly averaged Hill problem

We continue to analyze the periodic solutions of the singly averaged Hill problem. We have numerically constructed the families of solutions that correspond to periodically evolving satellite orbits for arbitrary initial values of their eccentricities and inclinations to the plane of motion of the perturbing body. The solutions obtained are compared with the numerical solutions of the rigorous (nonaveraged) equations of the restricted circular three-body problem. In particular, we have constructed a periodically evolving orbit for which the well-known Lidov-Kozai mechanism manifests itself, just as in the doubly averaged problem.     

Orbital evolution of giant comet-like objects

We investigated by numerical integrations the long-term orbital evolution of four giant comets or comet-like objects. They are Chiron, P/Schwassmann-Wachmann 1 (SW1), Hidalgo, and 1992AD (5145), and their orbits were traced for 100–200 thousand years (kyr) toward both the past and the future. For each object, 13 orbits were calculated, one for the nominal orbital elements and other 12 with slightly modified elements based on the rms residual of the orbit determination and on the number of observations. As past studies indicate, their orbital evolution is found to be very chaotic, and thus can be described only in terms of probability. Plots of the semi-major axis (a) and perihelion distance (q) of the objects treated here seem to cross each other frequently, suggesting a possibility of their common evolutionary paths. About a half of all the calculated orbits showedq- ora-decreasing evolution. This indicates that, at least on the time scale in question, the giant comet-like objects are possibly on a dynamical track that can lead to capture from the outer solar system. We could hardly find the orbits with perihelia far outside the orbit of Saturn (q>15 AU). This is perhaps because the evolution of the orbits beyond Saturn is so slow that substantial orbital changes do not take place within 100–200 kyr.     

Theoretical meteor radiants for macroscopic Taurid Complex objects

The calculation of theoretical meteor radiants is discussed for comets and asteroids whose orbits pass within, but at present do not necessarily intersect, that of the Earth, in particular from the perspective of developing a suitable method for application to Taurid Complex orbits. The main question addressed here is how to allow for dynamical evolution between epochs when an orbit isnot Earth-intersecting (as at present in most cases for macroscopic bodies) and those when itis (i.e., when meteors can actually be observed). This should be understood in terms of evolution in the past, such that meteoroids released some time ago have evolved differentially from the putative parents, allowing meteors to be detected now. Theoretical radiants for macroscopic Taurid objects are then presented and compared with observations of the nighttime and daytime Taurid meteor showers. These are found to be broadly similar in form, given the sparsity of some of the data, adding weight to the hypothesis that this sub-jovian complex contains kilometre-plus asteroids. A similar conclusion results for the group of objects in similar orbits to (2212) Hephaistos.     

Exploring the 7:4 mean motion resonance—I: Dynamical evolution of classical transneptunian objects

In the transneptunian classical region (), an unexpected orbital excitation in eccentricity and inclination, dynamically distinct populations and the presence of chaotic regions are observed. For instance, the 7:4 mean motion resonance () appears to have been causing unique dynamical excitation according to observational evidences, namely, an apparent shallow gap in number density and anomalies in the colour distribution, both features enhanced near the 7:4 mean motion resonance location. In order to investigate the resonance dynamics, we present extensive computer simulation results totalizing almost 10,000 test particles under the effect of the four giant planets for the age of the solar system. A chaotic diffusion experiment was also performed to follow tracks in phase space over 4-5 Gyr. The 7:4 mean motion resonance is weakly chaotic causing irregular eccentricity and inclination evolution for billions of years. Most 7:4 resonant particles suffered significant eccentricities and/or inclinations excitation, an outcome shared even by those located in the vicinity of the resonance. Particles in stable resonance locking are rare and usually had 0.25<e<0.3. For other regions, 7:4 resonants had quite large mobility in phase space typically leaving the resonance (and being scattered) after reaching a critical e∼0.2. The escape happened in 10⁸-10⁹ yr time scales. Concerning the inclination dependence for 7:4 resonants, we found strong instability islands for approximately i>10°. Taking into account those particles still locked in the resonance at the end of the simulations, we determined a retainability of 12-15% for real 7:4 resonant transneptunian objects (TNOs). Lastly, our results demonstrate that classical TNOs associated with the 7:4 mean motion resonance have been evolving continuously until present with non-negligible mixing of populations.     

Space weathering on near-Earth objects investigated by neutral-particle detection

The ion-sputtering (IS) process is active in many planetary environments in the solar system where plasma precipitates directly on the surface (for instance, Mercury, Moon and Europa). In particular, solar wind sputtering is one of the most important agents for the surface erosion of a near-Earth object (NEO), acting together with other surface release processes, such as photon stimulated desorption (PSD), thermal desorption (TD) and micrometeoroid impact vaporization (MIV). The energy distribution of the IS-released neutrals peaks at a few eVs and extends up to hundreds of eVs. Since all other release processes produce particles of lower energies, the presence of neutral atoms in the energy range above 10 eV and below a few keVs (sputtered high-energy atoms (SHEA)) identifies the IS process. SHEA easily escape from the NEO, due to NEO's extremely weak gravity. Detection and analysis of SHEA will give important information on surface-loss processes as well as on surface elemental composition. The investigation of the active release processes, as a function of the external conditions and the NEO surface properties, is crucial for obtaining a clear view of the body's present loss rate as well as for getting clues on its evolution, which depends significantly on space weather.In this work, an attempt to analyze processes that take place on the surface of these small airless bodies, as a result of their exposure to the space environment, has been realized. For this reason, a new space weathering model (space weathering on NEO-SPAWN) is presented. Moreover, an instrument concept of a neutral-particle analyzer specifically designed for the measurement of neutral density and the detection of SHEA from a NEO is proposed.     

On the stability of compact supermassive objects

Proceeded from the gravitation equations proposed by one of authors it was argued in a previous paper that there can exist supermassive compact configurations of degenerated Fermi‐gas without events horizon. In the present paper we consider the stability of these objects by a method like the one used in the theory of stellar structure. It is shown that the configurations are stable.     

The evolution of bodies orbiting in the trans-Neptunian region on to intermediate-type orbits

We investigate the dynamical evolution of 210 hypothetical massless bodies initially situated between 10 and 30 au from the Sun in order to determine the general characteristics of the evolved system. This is of particular relevance to the understanding of the origin of Edgeworth–Kuiper belt objects on scattered intermediate orbits, such as 1996 TL ₆₆, which have high eccentricity and semimajor axis but nevertheless have perihelion in the region between 30 and 50 au from the Sun.