Dynamical evolution of a cometary swarm in the outer planetary region

The dynamical evolution of bodies under the gravitational influence of the accreting proto-Uranus and proto-Neptune is investigated. The main aim of this study is to analyze the interrelations between the accretion of Uranus and Neptune with other processes of cosmological importance as, for example, the formation of a cometary reservoir from bodies placed into near-parabolic orbits by planetary perturbations and the scattering of bodies to the region of the terrestrial planets. Starting with a mass ratio (initial mass/present mass) of 0.1, Uranus and Neptune acquire masses close to their present ones in a time scale of 10⁸ years. Neptune is found to be the most important contributor of comets to the cometary reservoir. The time scale of bodies scattered by Neptune to reach near-parabolic orbits (semimajor axes a > 10⁴ AU)is about 10⁹ years. The contribution of Uranus was partially inhibited because a large part of the residual bodies of its accretion zone fell under the strong gravitational influence of Jupiter and Saturn. A significant fraction of the bodies dispersed by Uranus and Neptune reached the region of the terrestrial planets in a time scale of some 10⁸ years.     

Evolution of comet orbits under the perturbing influence of the giant planets and nearby stars

The orbital evolution of 500 hypothetical comets during 10⁹ years is studied numerically. It is assumed that the birthplace of such comets was the region of Uranus and Neptune from where they were deflected into very elongated orbits by perturbations of these planets. Then, we adopted the following initial orbital elements: perihelion distances between 20 and 30 AU, inclinations to the ecliptic plane smaller than 20°, and semimajor axes from 5 × 10³ to 5 × 10⁴ AU. Gravitational perturbations by the four giant planets and by hypothetical stars passing at distances from the Sun smaller than 5 × 10⁵ AU are considered. During the simulation, somewhat more than 50% of the comets were lost from the solar system due to planetary or stellar perturbations. The survivors were removed from the planetary region and left as members of what is generally known as the cometary cloud. At the end of the studied period, the semimajor axes of the surviving comets tend to be concentrated in the interval 2 × 10⁴ < a < 3 × 10⁴ AU. The orbital planes of the comets with initial $\text{a ≧ 3 × 10}{}^4\text{AU}$ acquired a complete randomization while the others still maintain a slight predominance of direct orbits. In addition, comet orbits with final a < 6 × 10⁴AU preserve high eccentricities with an average value greater than 0.8 Most “new” comets from the sample entering the region interior to Jupiter's orbit had already registered earlier passages through the planetary region. By scaling up the rate of paritions of hypothetical new comets with the observed one, the number of members of the cometary cloud is estimated to be about 7 × 10¹⁰ and the conclusion is drawn that Uranus and Neptune had to remove a number of comets ten times greater.     

Uranus and Neptune: Shape and rotation

Both Uranus and Neptune are thought to have strong zonal winds with velocities of several 100 m s⁻¹. These wind velocities, however, assume solid-body rotation periods based on Voyager 2 measurements of periodic variations in the planets’ radio signals and of fits to the planets’ magnetic fields; 17.24 h and 16.11 h for Uranus and Neptune, respectively. The realization that the radio period of Saturn does not represent the planet’s deep interior rotation and the complexity of the magnetic fields of Uranus and Neptune raise the possibility that the Voyager 2 radio and magnetic periods might not represent the deep interior rotation periods of the ice giants. Moreover, if there is deep differential rotation within Uranus and Neptune no single solid-body rotation period could characterize the bulk rotation of the planets. We use wind and shape data to investigate the rotation of Uranus and Neptune. The shapes (flattening) of the ice giants are not measured, but only inferred from atmospheric wind speeds and radio occultation measurements at a single latitude. The inferred oblateness values of Uranus and Neptune do not correspond to bodies rotating with the Voyager rotation periods. Minimization of wind velocities or dynamic heights of the 1 bar isosurfaces, constrained by the single occultation radii and gravitational coefficients of the planets, leads to solid-body rotation periods of ∼16.58 h for Uranus and ∼17.46 h for Neptune. Uranus might be rotating faster and Neptune slower than Voyager rotation speeds. We derive shapes for the planets based on these rotation rates. Wind velocities with respect to these rotation periods are essentially identical on Uranus and Neptune and wind speeds are slower than previously thought. Alternatively, if we interpret wind measurements in terms of differential rotation on cylinders there are essentially no residual atmospheric winds.     

Origin of the structure of the Kuiper belt during a dynamical instability in the orbits of Uranus and Neptune

We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. It successfully explains, for the first time, the current orbital architecture of the giant planets [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461], the existence of the Trojans populations of Jupiter and Neptune [Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465], and the origin of the late heavy bombardment of the terrestrial planets [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469]. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ∼35 AU to the Kuiper belt without excessive inclination excitation. In particular, we show that during the phase when Neptune's eccentricity is large, the region interior to its 1:2 mean motion resonance becomes unstable and disk particles can diffuse into this area. In addition, we perform numerical simulations where the planets are forced to evolve using fictitious analytic forces, in a way consistent with the direct N-body simulations of the Nice model. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ∼27 AU and an eccentricity of ∼0.3, and that subsequently Neptune's eccentricity damped in ∼1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level. In particular, our results explain, at least qualitatively: (1) the co-existence of resonant and non-resonant populations, (2) the eccentricity-inclination distribution of the Plutinos, (3) the peculiar semi-major axis—eccentricity distribution in the classical belt, (4) the outer edge at the 1:2 mean motion resonance with Neptune, (5) the bi-modal inclination distribution of the classical population, (6) the correlations between inclination and physical properties in the classical Kuiper belt, and (7) the existence of the so-called extended scattered disk. Nevertheless, we observe in the simulations a deficit of nearly-circular objects in the classical Kuiper belt.     

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.     

Analysis of the Hungaria asteroid population

There are approximately 5000 known asteroids in the Hungaria orbital space, a region defined by orbits with high inclination (16° < i < 34°), low eccentricities (e < 0.18), and semi-major axes 1.78 < a < 2.0 AU. We argue that this region is populated by a large number of asteroids formed after a catastrophic collision involving (434) Hungaria, the presumptive largest fragment of the Hungaria collisional family. The remaining objects form a background population that share orbital characteristics with the family members. Due to the general dynamic stability of the region, it is likely that most asteroids in Hungaria space (the Hungaria “group”) have been in this region since the formation of the Solar System or at least since the planets assumed their current orbital configuration. Our examination of the Hungaria group included comparing rotation rates, taxonomic classification, and orbital dynamics to determine the characteristics of the family and background populations. We first found there is an excess of slow rotators among the group but, otherwise, the distribution of spin frequencies is essentially uniform, i.e., that a plot of the cumulative number of objects over the range of 1 d⁻¹ < f < 9 d⁻¹ is nearly a straight line or, put another way, if the distribution over the range is binned by equal intervals of f (1-2 d⁻¹, 2-3 d⁻¹, etc.), the number of objects in each bin is statistically the same.There is a distinct family within the Hungaria group, centered at a semi-major axis of 1.940 AU, with a dispersion range that increases with decreasing size of members, as expected of an evolved collisional family. The larger members with well-determined taxonomic class, including (434) Hungaria itself, have flat spectra, mostly likely type E or similar. The degree of spreading versus size of family members is consistent with that expected from Yarkovsky thermal drift in roughly 0.5 Gyr, suggesting that age for the family. The Asteroid (434) Hungaria is displaced in semi-major axis by 0.004 AU from the center of the Hungaria family. The collision event that produced the family should not have left the largest body displaced by more than 0.001 AU from the original orbit, thus we infer that the displacement of (434) Hungaria is mainly due to Yarkovsky drift, and is consistent with the expected drift for that size body in ∼0.5 Gyr. Below ∼1.93 AU heliocentric distance the Hungaria family is perturbed by at least two secular resonances, 2g − g₅ − g₆ and one of the family of 4th or 6th order secular resonances near s ∼ −22.25 ″/year. Their combined effect results in larger inclination dispersion of the family members.     

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.     

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 scattered disk as a source of Halley-type comets

We investigate a new dynamical mechanism for producing Halley-type comets from the scattered disk of comets. Levison and Duncan [Levison, H., Duncan, M., 1997. Icarus 127, 13-32] and Duncan and Levison [Duncan, M., Levison, H., 1997. Science 276, 1670-1672] showed that a significant number of objects leave the scattered disk by evolving to semi-major axes greater than 1000 AU. We find that once these objects reach semi-major axes on the order of ¹⁰⁴ AU, a significant fraction immediately have their perihelia driven inward by the galactic tides. Approximately 0.01% of the objects that reach ¹⁰⁴ AU then evolve onto orbits similar to the observed Halley-like comets due to gravitational interactions with the giant planets. The orbital element distribution resulting from this process is statistically consistent with observations. We discuss the implications of this model for the number of objects in the scattered disk in the text. The model predicts a temporal variation in the influx of HTCs with a period of ∼118 Myr. At the peak, the model predicts that there should be roughly 10 times as many HTCs as currently observed (i.e., there should be weak HTC showers). However, the model may inflate the importance of these showers because it does not include the effects of passing stars and giant molecular clouds.     

What do the orbital motions of the outer planets of the Solar System tell us about the Pioneer anomaly?

In this paper, we investigate the effects that an anomalous acceleration as that experienced by the Pioneer spacecraft after they passed the 20 AU threshold would induce on the orbital motions of the Solar System planets placed at heliocentric distances of 20 AU or larger as Uranus, Neptune and Pluto. It turns out that such an acceleration, with a magnitude of 8.74 × 10⁻¹⁰ m s⁻², would affect their orbits with secular and short-period signals large enough to be detected according to the latest published results by E.V. Pitjeva, even by considering errors up to 30 times larger than those released. The absence of such anomalous signatures in the latest data rules out the possibility that in the region 20–40 AU of the Solar System an anomalous force field inducing a constant and radial acceleration with those characteristics affects the motion of the major planets.     

The web of three-planet resonances in the outer Solar System

In this paper we numerically detect the web of three-planet resonances (i.e., resonances among mean anomalies, nodes and perihelia of three planets) with respect to the variation of the semi-major axis of Saturn and Jupiter, in a model including the planets from Jupiter to Neptune. The measure confirms the relevance of these resonances in the long-term evolution of the outer Solar System and provides a technique to identify some of the related coefficients.     

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

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.     

Uranus and Neptune interior models

This paper is concerned with the interior structure of Uranus and Neptune. Our approach is three-fold. First, a set of three-layer models for both Uranus and Neptune are constructed using a method similar to that used in the study of the terrestrial planets. The variations of the mass density (s) and flattening e(s) with fractional mean radius s for two representative models of Uranus and Neptune are calculated. The results are tabulated. A comparison of these models shows that these two planets are probably very similar to each other in their basic dynamical features. Such similarity is very seldom seen in our solar system. Secondly, we check the conformance between the theoretical results and observational data for the two planets. And thirdly, the 6th degree Stokes zonal parameters for Uranus and for Neptune are predicted, based on the interior models put forward in this paper.     

The origin of the Kuiper Belt high-inclination population

I simulate the orbital evolution of the four major planets and a massive primordial planetesimal disk composed of 10⁴ objects, which perturb the planets but not themselves. As Neptune migrates by energy and angular momentum exchange with the planetesimals, a large number of primordial Neptune-scattered objects are formed. These objects may experience secular, Kozai, and mean motion resonances that induce temporary decrease of their eccentricities. Because planets are migrating, some planetesimals can escape those resonances while in a low-eccentricity incursion, thus avoiding the return path to Neptune close encounter dynamics. In the end, this mechanism produces stable orbits with high inclination and moderate eccentricities. The population so formed together with the objects coming from the classical resonance sweeping process, originates a bimodal distribution for the Kuiper Belt orbits. The inclinations obtained by the simulations can attain values above 30° and their distribution resembles a debiased distribution for the high-inclination population coming from the real classical Kuiper Belt.     

Influence of fast interstellar gas flow on the dynamics of dust grains

The orbital evolution of a dust particle under the action of a fast interstellar gas flow is investigated. The secular time derivatives of Keplerian orbital elements and the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle’s orbit are derived. The secular time derivatives of the semi-major axis, eccentricity, and of the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle’s orbit constitute a system of equations that determines the evolution of the particle’s orbit in space with respect to the gas flow velocity vector. This system of differential equations can be easily solved analytically. From the solution of the system we found the evolution of the Keplerian orbital elements in the special case when the orbital elements are determined with respect to a plane perpendicular to the gas flow velocity vector. Transformation of the Keplerian orbital elements determined for this special case into orbital elements determined with respect to an arbitrary oriented plane is presented. The orbital elements of the dust particle change periodically with a constant oscillation period or remain constant. Planar, perpendicular and stationary solutions are discussed. The applicability of this solution in the Solar System is also investigated. We consider icy particles with radii from 1 to 10 μm. The presented solution is valid for these particles in orbits with semi-major axes from 200 to 3000 AU and eccentricities smaller than 0.8, approximately. The oscillation periods for these orbits range from 10⁵ to 2 × 10⁶ years, approximately.     

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

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

The dynamics of inclined Neptune Trojans

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

Accretion of terrestrial planets from oligarchs in a turbulent disk

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al. [Laughlin, G., Steinacker, A., Adams, F.C., 2004. Astrophys. J. 608, 489-496]. The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas is also included. Without any effect of disk gas, Earth-mass planets are formed in terrestrial planet regions in a disk comparable to MMSN but with too large orbital eccentricities to be consistent with the present eccentricities of Earth and Venus in our Solar System. With the eccentricity damping caused by the tidal interaction with a remnant gas disk, Earth-mass planets with eccentricities consistent with those of Earth and Venus are formed in a limited range of disk gas surface density (∼10⁻⁴ times MMSN). However, in this case, on average, too many (?6) planets remain in terrestrial planet regions, because the damping leads to isolation between the planets. We have carried out a series of N-body simulations including the random torques with different disk surface density and strength of turbulence. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2 AU, which is more consistent with Solar System, for relatively wide range of disk gas surface density (∼¹⁰⁻⁴-¹⁰⁻² times MMSN).