Numerical simulation of the final stages of terrestrial planet formation

Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are presented. The mass and relative-velocity distributions of the bodies in these simulations are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies, therefore, arise in a natural way and the assumption of expressions for the relative velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and the Earth is possible, at least for the case of a two-dimensional system. Simulations assuming an initial uniform distribution of orbital eccentricities on the interval from 0 to e_max are found to produce final states containing too many bodies with masses which are too small when e_max < 0.10, while simulations with e_max > 0.20 result in too many catastrophic collisions between bodies thus preventing rapid accretion of planetary-size bodies. The e_max = 0.15 simulation ends with a state surprisingly similar to that of the present terrestrial planets and, therefore, provides a rough estimate of the range of radial sampling to be expected for the terrestrial planets.     

A co-accretional model of satellite formation

Numerical calculation of a simple accretion model including the effects of tidal friction indicate that coformation is tenable only if the planet's Q is less than about 10³. The parameter which most strongly affects the final mass ratio of the pair is the time at which the secondary embryo is introduced. Our model yields the proper Moon-Earth mass ratio if the Moon embryo is introduced when the Earth is only about $\text{1}\text{10}$ of its final mass. The lunar orbit remains at about 10 Earth radii throughout most of the growth.This model of satellite formation overcomes two difficulties of the “circumterrestrial cloud” model of Ruskol (1960, 1963, 1972): (1) The difficulty of accumulating a mass as great as the entire Moon before gravitational instability reduces the cloud to a small number of moonlets is removed. (2) The differences between terrestrial and outer planet satellite systems is easily understood in terms of the differences in Q between these planets. The high Q of the outer planets does not allow a satellite embryo to survive a significant portion of the accretion process, thus only small bodies which formed very late in the accumulation of the planet remain as satellites. The low Q of the terrestrial planets allows satellite embryos of these planets to survive during accretion, thus massive satellites such as the Earth's Moon are expected. The present lack of such satellites of the other terrestrial planets may be the result of tidal evolution, either infall following primary despinning (Burns, 1973) or escape due to increase in orbit eccentricity.     

On the process of accretion in the formation of the planets and comets

The physically reasonable assumption that the seed bodies which initiated the accretion of the individual asteroids, planets, and comets (subsequently these objects are collectively called planetoids) formed by stochastic processes requires a radius distribution function which is unique except for two scaling parameters: the total number of planetoids and their most probable radius. The former depends on the ease of formation of the seed bodies while the second is uniquely determined by the average pre-encounter velocity, V, of the accretable material relative to an individual planetoid. This theoretical radius function can be fit to the initial asteroid radius distribution which Anders (1965) derived from the present-day distribution by allowing for fragmentation collisions among the asteroids since their formation. Normalizing the theoretical function to this empirical distribution reveals that there were about 10² precollision asteroids and that V = (2?4) × 10^?2 km/sec which was presumably the turbulent velocity in the Solar Nebula. Knowing V we can determine the scale height of the dust in the Solar Nebula and consequently its space density. The density of accretable material determines the rate of accretion of the planetoids. From this we find, for example, that the Earth formed in about 8 × 10⁶ yr and it attained a maximum temperature through accretion of about 3 × 10³°K. From the total mass of the terrestrial planets and the theoretical radius function we find that about 2 × 10³ planetoids formed in the vicinity of the terrestrial planets. Except for the asteroids the smaller planetoids have since been accreted by the terrestrial planets. About 15% of the present mass of the terrestrial planets was accumulated by the secondary accretion of these smaller primary planetoids. There are far fewer primary planetoids than craters on the Moon or Mars. The craters were likely produced by the collisional breakup of a few primary planetoids with masses between one-tenth and one lunar mass. This deduction comes from comparing the collision cross sections of the planetoids in this mass range to that of the terrestrial planets. This comparison shows that two to three collisions leading to the breakup of four to six objects likely occurred among these objects before their accretion by the terrestrial planets. The number of these fragments is quite adequate to explain the lunar and Martin craters. Furthermore the mass spectrum of such fragments is a power-law distribution which results in a power-law distribution of crater radii of just the type observed on the Moon and Mars. Applying the same analysis to the planetoids which formed in the vicinity of the giant planets reveals that it is unlikely that any fragmentation collisions took place among them before they were accreted by these planets due to the integrated collision cross section of the giant planets being about three orders of magnitude greater than that of the terrestrial planets. We can thus anticipate a marked scarcity of impact craters on the satellites of these outer planets. This prediction can be tested by future space probes. Our knowledge of the radius function of the comets is consistent with their being primary planetoids. The primary difference between the radius function of the planetoids which formed in the inner part of the solar system and that of the comets results from the fact that the seed bodies which grew into the comets formed far more easily than those which grew into the asteroids and the terrestrial planets. Thus in the outer part of the Solar Nebula the principal solid material (water and ammonia snow) accreted into a huge (～10¹²⁺) number of relatively small objects (comets) while in the inner part of the nebula the solid material (hard-to-stick refractory substances) accumulated into only a few (～10³) large objects (asteroids and terrestrial planets). Uranus and Neptune presumably formed by the secondary accretion of the comets.     

A model for accretion of the terrestrial planets

One possible origin of the terrestrial planets involves their formation by gravitational accretion of particles originally in Keplerian orbits about the sun. Some implications of this theory are considered. A formal expression for the rate of mass accretion by a planet is developed. The formal singularity of the gravitational collision cross-section for low relative velocities is shown to be without physical significance when the accreting bodies are in heliocentric orbits. The distribution of particle velocities relative to an accreting planet is considered; the mean velocity increases with time. The internal temperature of an accreting planet is shown to depend simply on the accretion rate. A simple and physically reasonable approximate expression for a planetary accretion rate is proposed.     

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

Accretion disks around Jupiter and Saturn at the stage of regular satellite formation

Modern models of the formation of the regular satellites of giant planets, constructed with consideration for their structure and composition suggest that this process lasted for a considerable period of time (0.1–1 Myr) and developed in gas-dust circumplanetary disks at the final stage of giant planet formation. The parameters of protosatellite disks (e.g., the radial distribution of surface density and temperature) serve as important initial conditions for such models. Therefore, the development of protosatellite disk models that take into account currently known cosmochemical and physical restrictions remains a pressing problem. It is this problem that is solved in the paper. New models of the accretion disks of Jupiter and Saturn were constructed with consideration for the disk heating by viscous dissipation of turbulent motions, by accretion of material from the surrounding region of the solar nebula, and by radiation from the central planets. The influence of a set of input model parameters (the total rate of mass infall onto the disk, the turbulent viscosity and opacity of disk material, and the centrifugal radius of the disk) on thermal conditions in the accretion disks was studied. The dependence of opacity on temperature and the abundance and size of solid particles present in the disk was taken into account. Those constructed models that satisfy the existing constraints limit the probable values of input parameters (primarily rates of mass infall onto the disks of Jupiter and Saturn at the stage of regular satellite formation and, to a lesser extent, the disk opacities). Constraints on the location of the regions of formation of the major satellites of Jupiter and Saturn are suggested based on the constructed models and simple analytical estimates concerning the formation of satellites in the accretion disks. It is shown that Callisto and Titan could hardly be formed at significantly greater distances from their planets.     

Origin and evolution of the solar system,II

(7)Formation of celestial bodies. The basic concepts of the accretional process are discussed, and the inadequacy of the contractional model is pointed out. A comparison is made between the general pre-planetary state on the one hand and the present state in the asteroidal region on the other. A model for accretion of resonance-captured grains leading to the formation of resonance-captured planets and satellites is suggested.(8)Spin and accretion. The relation between the accretional process and the spin of planets is analyzed.(9)Accretion of planets and satellites. It is shown that jet streams are a necessary intermediate stage in the formation of celestial bodies. The time sequence of planet formation is analyzed, and it is shown that the newly accreted bodies have a characteristic internal heat structure; the cases of the Earth and the Moon are considered in detail. A region of high initial temperature is found at 0.4 of the present Earth radius, whereas the culminating temperature of the Moon is near its present surface. An accretional heat wave is found to proceed outwards, and may produce the observed differentiation features.     

Accretion among preplanetary bodies: The many faces of runaway growth

When preplanetary bodies reach proportions of ∼1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. The numerical model is novel as it attempts to preserve the individual particle nature of the bodies (like N-body codes); yet it is statistical in nature since it must incorporate the very large number of planetesimals. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R₀. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We associate the initial runaway growth phase with increasing GF-factors for the largest body. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high-mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two-component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. We identify the peak of GF with the transition between the runaway growth and oligarchy accretion stages. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius, especially at the outer disk. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion, which could lead to a very fast growth.     

The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

N-Body simulations of growth from 1 km planetesimals at 0.4 AU

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M_⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10⁵ particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.     

When and where were the satellites of uranus formed?

Uranus exhibits an unusually large obliquity compared to other planets of the solar system; its equator is inclined by 98° to the plane of its orbit. However its five satellites are remarkably regular, with eccentricities and inclinations very nearly zero, but of course with orbit planes that are tilted by ～98° to the plane of the ecliptic. This circumstance is used here to relate the formation of satellites to planet formation. Six different cases are discussed, of which two can be ruled out and two others are highly improbable. In the analysis, use is made of the fact that satellites in near-equatorial orbits could not follow a rapid (“non-adiabatic”) change of the planet's obliquity. We assume, also, that the observed obliquity is the result of the last stages of planet accumulation. We can therefore exclude contemporaneous formation of planet and satellites, and conclude instead that the satellites were formed or acquired after the planet's axis had been tilted. A plausible scenario involves the tidal capture of a body having 5% to 10% of the planet's mass—sufficient to account for the tilt—followed by its accretion. However, tidal forces break up the body into chunks, slow the accretion, and allow ～1% of the chunks to form the satellites through interaction with a temporary dense atmosphere. The same reasoning may apply also for Saturn and Jupiter. It should be noted that the synchronous orbit it well within the Roche limit for all three planets.     

The internal activity and thermal evolution of Earth-like planets

We model the internal thermal evolution of planets with Earth-like composition and masses ranging from 0.1 to 10 Earth masses over a period of 10 billion years. We also characterize the internal activity of the planets by the velocity of putative tectonic plates, the rate at which mantle material is processed through melting zones, and the time taken to process one mantle mass. The more massive the planet the larger its processing rate (?), which scales approximately as ?∝M0.8-1.0. The processing times for all the planets increase with time as they cool and become less active. As would be expected, the surface heat flow scales with planet mass. All planets have similar declines in mantle temperature except for the largest, in which pressure effects cause a larger decline. The larger planets have higher mantle temperatures over all times. The less massive the planet, the larger the decrease in core temperature with time. The core heat flow is also found to decrease more rapidly for smaller planet masses. Finally, rough predictions are made for the time required to generate an atmosphere from estimates of the time to degas water and carbon dioxide in mantle melting zones. The degassing times depend strongly on the initial temperature of the planet, but for the temperatures used in our model all the planets degas within ∼32 Ma after their formation.     

Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.     

Satellite formation,II

The satellite formation model of Harris and Kaula (Icarus24, 516–524, 1975) is extended to include evolution of planetary ring material and elliptic orbital motion. This model is more satisfactory than the previous one in that the formation of the moon begins at a later time in the growth of the earth, and that a significant fraction of the lunar material is processed through a circumterrestrial debris cloud where volatiles might have been lost. Thus the chemical differences between the earth and moon are more plausibly accounted for. Satellites of the outer planets probably formed in large numbers throughout the growth of those planets. Because of rapid inward evolution of the orbits of small satellites, the present satellite systems represent only satellites formed in the last few percent of the growths of their primaries. The rings of Saturn and Uranus are most plausibly explained as the debris of satellites disrupted within the Roche limit. Because such a ring would collapse onto the planet in the course of any significant further accretion by the planet, the rings must have formed very near or even after the conclusion of accretion.     

Giant planet migration through the action of disk torques and planet-planet scattering

This paper presents a parametric study of giant planet migration through the combined action of disk torques and planet-planet scattering. The torques exerted on planets during Type II migration in circumstellar disks readily decrease the semi-major axes a, whereas scattering between planets increases the orbital eccentricities ?. This paper presents a parametric exploration of the possible parameter space for this migration scenario using two (initial) planetary mass distributions and a range of values for the time scale of eccentricity damping (due to the disk). For each class of systems, many realizations of the simulations are performed in order to determine the distributions of the resulting orbital elements of the surviving planets; this paper presents the results of ∼8500 numerical experiments. Our goal is to study the physics of this particular migration mechanism and to test it against observations of extrasolar planets. The action of disk torques and planet-planet scattering results in a distribution of final orbital elements that fills the a-? plane, in rough agreement with the orbital elements of observed extrasolar planets. In addition to specifying the orbital elements, we characterize this migration mechanism by finding the percentages of ejected and accreted planets, the number of collisions, the dependence of outcomes on planetary masses, the time spent in 2:1 and 3:1 resonances, and the effects of the planetary IMF. We also determine the distribution of inclination angles of surviving planets and the distribution of ejection speeds for exiled planets.     

Structure and evolutionary history of the solar system,IV

In this fourth and last part of our analysis, the first section (14) contains a study of the chemical composition of the planets and satellites. A sharp distinction is made between the large quantity of speculations about the interiors of the bodies and the rather meagerfacts known with a reasonable degree of certainty. It is shown, however, that the latter are sufficient todisprove the old concept of a Laplacian disc of homogeneous chemical composition. There is asystematic variations in the chemical composition of planets (and probably also of satellites) so that heavy elements are more abundant in the outermost and in the innermost regions of the systems. Section 15 containsa study of meteorites. These have earlier been interpreted in terms of ‘exploded planets’ and condensation processes in thermodynamic equilibrium. It is shown that such models are irreconcilable with the laws of physics and also with the meteoritic observations. These instead are found toprovide abundant information on the processes in jet streams and on early fractionation and condensation. Further work along these lines supplemented with other solar system materials studies may lead to a detailed reconstruction of important events in the evolution of the solar system. Section 16 demonstrates that the location of the different groups of secondary bodies is a result of a plasma phenomenon occurring at the critical velocity limit. These have recently been studied in detail in the laboratory but have not yet been fully applied to astrophysics.Groups of bodies in the planetary and the satellite systems related by the critical velocity shouldhave the same gravitational potential. There are large chemical differences between groups of different gravitational potential. This is reconcilable with the chemical differentiation found in Section 14. Finally, Section 17 deals with thestructure of the different groups of bodies and shows that the mass distributionis a function of the spin of the central body. Summarizing the properties and distribution of bodies in the solar system against this background, it is shown that there isno need for ‘missing planets’ or to explode hypothetical large bodies. Nor is there any justification for involvingdrastic ad hoc changes in the orbits of existing bodies. The scheme is complete in the sense that in all places where groups of bodies are expected, such bodies are actually found. All of the existing bodies are accounted for (with the exception of the small Martian satellites!). The general conclusion is that already with the empirical material now availableit is possible to suggest a series of basic processes leading to the present structure of planet and satellite systems in an internally consistent way. With the expected flow of data from space research the evolution of the solar system may eventually be described with about the same confidence and accuracy as the geological evolution of the Earth.     

Numerical experiment applicable to the latest stage of planet growth

A three-dimensional numerical model was developed with the goal of studying limited dynamical problems relevant to the latest stage of planet growth in the accretion theory. A small number of large protoplanets (～ Moon size) of different masses, moving around the Sun, are considered. The dynamical evolution and growth of the population is studied under mutual gravitational perturbations, accretion, and collisional fragmentation processes. Gravitational encounters are treated exactly by numerical integration of the N-body problem. Outcomes of collisional fragmentation are modeled according to the results of R. Greenberg et al. (1978, Icarus, 35, 1–26). In the present work, we consider 25 protoplanets with uniform mass distribution in the range 2 × 10²⁵?4 × 10²⁶ g on heliocentric orbits in the Earth zone. These bodies are initially confined to a small volume of space to permit gravitational perturbations by close approaches and collisions within a finite length of integration time. The dynamical evolution of the swarm is followed for four different sets of initial ranges in semimajor axis, eccentricity, and inclination: Δa=0.01, 0.02, 0.04, 0.08 AU; Δe= 0.005, 0.01, 0.02, 0.04; Δi=0°3, 0°6, 1°2, 2°4. Among other results, it is found that average eccentricities and inclinations evolve toward a steady state such that $\text{i}\text{?}\text{1}\text{2}\text{,}\text{e}$; it is also found that, whatever the initial conditions, the population evolves toward a quasi-equilibrium relative velocity distribution corresponding to a Safronov parameter value $\text{θ?10}$. Moreover, the growth process of the growing planet presents very similar behavior in the four cases considered, except for the time scale of evolution, which increases with the initial range of orbital elements. Earlier works of this kind have been presented by L.P. Cox and J.S. Lewis (1980, Icarus, 44, 706–721) and by G.N. Wetherill (1980b, In Geol. Soc. Canad. Spec. Publ., p. 20), although a number of differences exist between the three approaches.     

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.     

Late stages of planetary accretion

We derive first-order differential equations for the late stages of planetary accretion (planetesimal mass >10¹³ g). The effect of gravitational encounters, energy exchange, collisions, and gas drag has been included. Two simple models are discussed, namely, (i) when all planetesimals have the same mass and (ii) when there is one large planetesimal and numerous small planetesmals. Gravitational two-body encounters are modeled according to Chandrasekhar's classical theory from stellar dynamics. It is shown that the velocity increase due to mutual encounters can be modeled according to the simple theory of random flights. We find analytical equations for the average velocity decrease due to collisions. Gas drag, if present, is modeled in averaged form up to the first order in the eccentricities and inclinations of the planetesimals. Characteristic time scales for the formation of the terrestrial planets are found for the most favorable models to be of order 10⁸ year. The calculated mass of rock and ice of the giant planets is too low as compared to the observed one. This difficulty of our model could be overcome by assuming a several times larger surface density, an enlarged accretion cross section, and gas accretion during the final stages of accretion of the solid cores of the giant planets. Analytical and numerical results are presebted, the evolutionary tracks showing satisfactory agreement with observations for some models.     

On the stability of a planet between Mars and the asteroid belt: Implications for the Planet V hypothesis

The stability of an additional planet between the orbit of Mars and the asteroid belt is examined in the context of the Planet V hypothesis. In this model, the Solar System initially contained a fifth terrestrial planet, “Planet V,” which was removed after ∼700 Myr, a possible trigger for the late heavy bombardment on the inner planets. The model is investigated using 96 N-body integrations of the 8 major planets with an additional body between Mars and the asteroid belt. In more than 1/4 of simulations, Planet V survives for 1000 Myr. In most other cases, Planet V collides with the Sun or hits another planet after several hundred Myr, leaving 4 surviving terrestrial planets. In 24/96 simulations, Planet V is lost by ejection or collision with the Sun while the other four terrestrial planets survive without undergoing a collision. In 18 cases, Planet V is removed at least 200 Myr after the beginning of the simulation. The endstate depends sensitively on the mass of Planet V. Collision with the Sun is likely when Planet V's mass is 0.25 Mars masses or less. When Planet V is more massive than this, collisions involving it and/or other terrestrial planets become commonplace. In unstable systems, the times of first encounter and first collision/ejection depend on the initial aphelion distance of Mars. Reducing Mars's aphelion distance increases these times and also increases the fraction of systems surviving for 1000 Myr. When Mars's current orbit is used, the stability of Planet V increases when these two planets are widely separated initially. Planet V's aphelion distance Q typically begins to cross the asteroid belt within a few tens to a few hundred Myr, and its orbit last leaves the belt several hundred Myr later in most cases. The total time spent with Q>2.1 AU is typically less than 200 Myr.