High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets - Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model”. We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size.     

Populating the asteroid belt from two parent source regions due to the migration of giant planets—“The Grand Tack”

Abstract– The asteroid belt is found today in a dramatically different state than that immediately following its formation. It is estimated that it has been depleted in total mass by a factor of at least 1000 since its formation, and that the asteroids’ orbits evolved from having near‐zero eccentricity and inclination to the complex distributions we find today. The asteroid belt also hosts a wide range of compositions, with the inner regions dominated by S‐type and other water‐poor asteroids and the outer regions dominated by C‐type and other primitive asteroids. We discuss a model of early inner solar system evolution whereby the gas‐driven migration of Jupiter and Saturn brings them inwards to 1.5 AU, truncating the disk of planetesimals in the terrestrial planet region, before migrating outwards toward their current locations. This model, informally titled “The Grand Tack,” examines the planetary dynamics of the solar system bodies during the final million years of the gaseous solar nebula lifetime—a few million years (Myr) after the formation of the first solids, but 20–80 Myr before the final accretion of Earth, and approximately 400–600 Myr before the Late Heavy Bombardment of the inner solar system. The Grand Tack attempts to solve some outstanding problems for terrestrial planet formation, by reproducing the size of Mars, but also has important implications for the asteroid population. The migration of Jupiter causes a very early depletion of the asteroid belt region, and this region is then repopulated from two distinct source regions, one inside the formation region of Jupiter and one between and beyond the giant planets. The scattered material reforms the asteroid belt, producing a population the appropriate mass, orbits, and with overlapping distributions of material from each parent source region.     

Origin of water in the terrestrial planets

Abstract— I examine the origin of water in the terrestrial planets. Late‐stage delivery of water from asteroidal and cometary sources appears to be ruled out by isotopic and molecular ratio considerations, unless either comets and asteroids currently sampled spectroscopically and by meteorites are unlike those falling to Earth 4.5 Ga ago, or our measurements are not representative of those bodies. However, the terrestrial planets were bathed in a gas of H, He, and O. The dominant gas phase species were H₂, He, H₂ O, and CO. Thus, grains in the accretion disk must have been exposed to and adsorbed H2 and water. Here I conduct a preliminary analysis of the efficacy of nebular gas adsorption as a mechanism by which the terrestrial planets accreted “wet.” A simple model suggests that grains accreted to Earth could have adsorbed 1‐3 Earth oceans of water. The fraction of this water retained during accretion is unknown, but these results suggest that examining the role of adsorption of water vapor onto grains in the accretion disk bears further study.     

Planetary growth with collisional fragmentation and gas drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 10⁷ years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.     

The Leonard Award Address Presented 1998 July 27, Dublin,Ireland: On the difficulties of making Earth-like planets

Abstract— Here I discuss the series of events that led to the formation and evolution of our planet to examine why the Earth is unique in the solar system. A multitude of factors are involved: These begin with the initial size and angular momentum of the fragment that separated from a molecular cloud; such random factors are crucial in determining whether a planetary system or a double star develops from the resulting nebula. Another requirement is that there must be an adequate concentration of heavy elements to provide the 2% “rock” and “ice” components of the original nebula. An essential step in forming rocky planets in the inner nebula is the loss of gas and depletion of volatile elements, due to early solar activity that is linked to the mass of the central star. The lifetime of the gaseous nebula controls the formation of gas giants. In our system, fine timing was needed to form the gas giant, Jupiter, before the gas in the nebula was depleted. Although Uranus and Neptune eventually formed cores large enough to capture gas, they missed out and ended as ice giants. The early formation of Jupiter is responsible for the existence of the asteroid belt (and our supply of meteorites) and the small size of Mars, whereas the gas giant now acts as a gravitational shield for the terrestrial planets. The Earth and the other inner planets accreted long after the giant planets, from volatile-depleted planetesimals that were probably already differentiated into metallic cores and silicate mantles in a gas-free, inner nebula. The accumulation of the Earth from such planetesimals was essentially a stochastic process, accounting for the differences among the four rocky inner planets—including the startling contrast between those two apparent twins, Earth and Venus. Impact history and accretion of a few more or less planetesimals were apparently crucial. The origin of the Moon by a single massive impact with a body larger than Mars accounts for the obliquity (and its stability) and spin of the Earth, in addition to explaining the angular momentum, orbital characteristics, and unique composition of the Moon. Plate tectonics (unique among the terrestrial planets) led to the development of the continental crust on the Earth, an essential platform for the evolution of Homo sapiens. Random major impacts have punctuated the geological record, accentuating the directionless course of evolution. Thus a massive asteroidal impact terminated the Cretaceous Period, resulted in the extinction of at least 70% of species living at that time, and led to the rise of mammals. This sequence of events that resulted in the formation and evolution of our planet were thus unique within our system. The individual nature of the eight planets is repeated among the 60-odd satellites—no two appear identical. This survey of our solar system raises the question whether the random sequence of events that led to the formation of the Earth are likely to be repeated in detail elsewhere. Preliminary evidence from the “new planets” is not reassuring. The discovery of other planetary systems has removed the previous belief that they would consist of a central star surrounded by an inner zone of rocky planets and an outer zone of giant planets beyond a few astronomical units (AU). Jupiter-sized bodies in close orbits around other stars probably formed in a similar manner to our giant planets at several astronomical units from their parent star and, subsequently, migrated inwards becoming stranded in close but stable orbits as “hot Jupiters”, when the nebula gas was depleted. Such events would prevent the formation of terrestrial-type planets in such systems.     

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.     

Mean Motion Resonances, Gas Drag, and Supersonic Planetesimals in the Solar Nebula

We examine the orbital evolution of planetesimals under the influence of Jupiter's perturbations and nebular gas drag, under the assumption that gas persisted in the asteroid region for some time after Jupiter attained its final mass. Two distinct mechanisms, associated with the 2 : 1 and 3 : 2 mean motion resonances, can excite eccentricities to high values, despite the damping effect of drag. If Jupiter's eccentricity was comparable to its present value, planetesimals can be temporarily trapped in the 2 : 1 resonance. Bodies crossing the 3 : 2 resonance can enter a region of phase space with overlapping high-order resonances. Both mechanisms can produce eccentricities greater than 0.5 for asteroid-sized planetesimals. The combination of resonant perturbations and drag causes secular decay of semimajor axes, resulting in migration of bodies from the outer to inner belt. Inclinations remain low, implying significant collisional evolution during this migration. Velocities of resonant bodies relative to the gas are highly supersonic; these would have been a source of shock waves in the solar nebula.This revised version was published online in October 2005 with corrections to the Cover Date.     

Planetesimal dissolution in the envelopes of the forming,giant planets

We have assessed the ability of planetesimals to penetrate through the envelopes of growing giant planets that form by a “core-instability” mechanism. According to this mechanism, a core grows by the accretion of solid bodies in the solar nebula and the growing core becomes progressively more effective in gravitationally concentrating gas from the surrounding solar nebula in an envelope until a “runaway” accretion of gas occurs. In performing this assessment, we have considered the ability of gas drag to slow down a planetesimal; the effectiveness of gas dynamical pressure in fracturing and ultimately finely fragmenting it; the ability of its strength and self-gravity to resist such fracturing; and the degree to which it is evaporated due to heating by the surrounding envelope, including shock heating that develops during the supersonic portion of its trajectory. We also consider what happens if the planetesimal is able to reach the core at free-fall velocity and the ability of the envelope to convectively mix dissolved materials to different radial distances. These calculations were performed for various epochs in the growth of a giant planet with the model envelopes derived by Bodenheimer and Pollack (1986,67, 391–408). As might have been anticipated, our results vary significantly with the size of the planetesimal, its composition, and the stage of growth of the giant planet and hence the mass of its envelope. Over much of the growth phase of the core, prior to its reaching its critical mass for runaway gas accretion, icy planetesimals less than about 1 m in size dissolve in the outer region of the envelope, ones larger than about 1 m and smaller than about 1 km dissolve in the middle region of the envelope, ones larger than 1 km either reach the core interface or dissolve in the deeper regions of the envelope. Similarly rocky planetesimals smaller than about a kilometer dissolve in the middle portion of the envelope, while larger ones can penetrate more deeply. Furthermore, the convection zones of the envelopes during this stage are confined to localized regions and hence dissolved materials experience little radial mixing then. Thus, if much of the accreted mass is contained in planetesimals larger than about a kilometer, the critical core mass for runaway accretion is not expected to change significantly when planetesimal dissolution is taken into account. After accretion is terminated and the planet contracts toward its present size, the convection zone grows until it encompasses the entire envelope. Therefore, dissolved material should eventually become well mixed through the envelope. We proposed that the envelopes of the giant planets should contain significant enhancements above solar proportions in the abundances of virtually all elements relative to that of hydrogen, with the magnitude of the enhancement increasing approximately linearly with the ratio of the high Z mass to the (H, He) mass for the bulk of the planet. This prediction is in accord both qualitatively and quantitatively with the systematic increase in the atmospheric C/H ratio from Jupiter to Saturn to Uranus and Neptune and semiquantitatively with the results of recent interior models of the giant planets. It is not clear whether it is consistent with the abundances of H₂O and NH₃ in the atmospheres of some of the outer planets. Finally, the complete reduction of some dissolved materials, especially C containing compounds, is expected to consume some of the H₂ in the envelopes. Consequently, the He/H₂ ratios in the atmospheres of Uranus and Neptune may be slightly enhanced over the solar ratio. We estimate that the He/H₂ ratios for Uranus' and Neptune's atmospheres should be about 6 and 15% larger, respectively, than the solar ratio.     

Oligarchic growth of giant planets

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as “oligarchic growth.” Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk (∼10 × minimum-mass) is required to produce giant planet core-sized bodies (∼10 M_⊕) within the lifetime of the nebular gas (?10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of noninteracting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System’s two outer giant planets at their current locations on any timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.     

Formation of gas giant planets: core accretion models with fragmentation and planetary envelope

We have calculated formation of gas giant planets based on the standard core accretion model including effects of fragmentation and planetary envelope. The accretion process is found to proceed as follows. As a result of runaway growth of planetesimals with initial radii of ∼10 km, planetary embryos with a mass of ∼10²⁷ g (∼ Mars mass) are found to form in ∼10⁵ years at Jupiter's position (5.2 AU), assuming a large enough value of the surface density of solid material (25 g/cm²) in the accretion disk at that distance. Strong gravitational perturbations between the runaway planetary embryos and the remaining planetesimals cause the random velocities of the planetesimals to become large enough for collisions between small planetesimals to lead to their catastrophic disruption. This produces a large number of fragments. At the same time, the planetary embryos have envelopes, that reduce energies of fragments by gas drag and capture them. The large radius of the envelope increases the collision rate between them, resulting in rapid growth of the planetary embryos. By the combined effects of fragmentation and planetary envelope, the largest planetary embryo with 21M_⊕ forms at 5.2 AU in 3.8×10⁶ years. The planetary embryo is massive enough to start a rapid gas accretion and forms a gas giant planet.     

Formation of terrestrial planets in disks with different surface density profiles

We present the results of an extensive study of the final stage of terrestrial planet formation in disks with different surface density profiles and for different orbital configurations of Jupiter and Saturn. We carried out simulations in the context of the classical model with disk surface densities proportional to \({r^{-0.5}}, {r^{-1}}\) and \({r^{-1.5}}\), and also using partially depleted, non-uniform disks as in the recent model of Mars formation by Izidoro et al. (Astrophys J 782:31, 2014). The purpose of our study is to determine how the final assembly of planets and their physical properties are affected by the total mass of the disk and its radial profile. Because as a result of the interactions of giant planets with the protoplanetary disk, secular resonances will also play important roles in the orbital assembly and properties of the final terrestrial planets, we will study the effect of these resonances as well. In that respect, we divide this study into two parts. When using a partially depleted disk (Part 1), we are particularly interested in examining the effect of secular resonances on the formation of Mars and orbital stability of terrestrial planets. When using the disk in the classical model (Part 2), our goal is to determine trends that may exist between the disk surface density profile and the final properties of terrestrial planets. In the context of the depleted disk model, results of our study show that in general, the \(\nu _5\) resonance does not have a significant effect on the dynamics of planetesimals and planetary embryos, and the final orbits of terrestrial planets. However, \(\nu _6\) and \(\nu _{16}\) resonances play important roles in clearing their affecting areas. While these resonances do not alter the orbits of Mars and other terrestrial planets, they strongly deplete the region of the asteroid belt ensuring that no additional mass will be scattered into the accretion zone of Mars so that it can maintain its mass and orbital stability. In the context of the classical model, the effects of these resonances are stronger in disks with less steep surface density profiles. Our results indicate that when considering the classical model (Part 2), the final planetary systems do not seem to show a trend between the disk surface density profile and the mean number of the final planets, their masses, time of formation, and distances to the central star. Some small correlations were observed where, for instance, in disks with steeper surface density profiles, the final planets were drier, or their water contents decreased when Saturn was added to the simulations. However, in general, the final orbital and physical properties of terrestrial planets seem to vary from one system to another and depend on the mass of the disk, the spatial distribution of protoplanetary bodies (i.e., disk surface density profile), and the initial orbital configuration of giant planets. We present results of our simulations and discuss their implications for the formation of Mars and other terrestrial planets, as well as the physical properties of these objects such as their masses and water contents.     

Terrestrial planet formation with strong dynamical friction

We have performed 8 numerical simulations of the final stages of accretion of the terrestrial planets, each starting with over 5× more gravitationally interacting bodies than in any previous simulations. We use a bimodal initial population spanning the region from 0.3 to 4 AU with 25 roughly Mars-mass embryos and an equal mass of material in a population of ∼1000 smaller planetesimals, consistent with models of the oligarchic growth of protoplanetary embryos. Given the large number of small planetesimals in our simulations, we are able to more accurately treat the effects of dynamical friction during the accretion process. We find that dynamical friction can significantly lower the timescales for accretion of the terrestrial planets and leads to systems of terrestrial planets that are much less dynamically excited than in previous simulations with fewer initial bodies. In addition, we study the effects of the orbits of Jupiter and Saturn on the final planetary systems by running 4 of our simulations with the present, eccentric orbits of Jupiter and Saturn (the EJS simulations) and the other 4 using a nearly circular and co-planar Jupiter and Saturn as predicted in the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465] (the CJS simulations). Our EJS simulations provide a better match to our Solar System in terms of the number and average mass of the final planets and the mass-weighted mean semi-major axis of the final planetary systems, although increased dynamical friction can potentially improve the fit of the CJS simulations as well. However, we find that in our EJS simulations, essentially no water-bearing material from the outer asteroid belt ends up in the final terrestrial planets, while a large amount is delivered in the CJS simulations. In addition, the terrestrial planets in the EJS simulations receive a late veneer of material after the last giant impact event that is likely too massive to reconcile with the siderophile abundances in the Earth's mantle, while the late veneer in the CJS simulations is much more consistent with geochemical evidence.     

Terrestrial planet and asteroid formation in the presence of giant planets. I. Relative velocities of planetesimals subject to Jupiter and Saturn perturbations.

We investigate the orbital evolution of 10(13)- to 10(25) -g planetesimals near 1 AU and in the asteroid belt (near 2.6 AU) prior to the stage of evolution when the mutual perturbations between the planetesimals become important. We include nebular gas drag and the effects of Jupiter and Saturn at their present masses and in their present orbits. Gas drag introduces a size-dependent phasing of the secular perturbations, which leads to a pronounced dip in encounter velocities (Venc) between bodies of similar mass. Plantesimals of identical mass have Venc approximately 1 and approximately 10 m s-1 (near 1 and 2.6 AU, respectively) while bodies differing by approximately 10 in mass have Venc approximately 10 and approximately 100 m s-1 (near 1 and 2.6 AU, respectively). Under these conditions, growth, rather than erosion, will occur only by collisions of bodies of nearly the same mass. There will be essentially no gravitational focusing between bodies less than 10(22) to 10(25) g, allowing growth of planetary embryos in the terrestrial planet region to proceed in a slower nonrunaway fashion. The environment in the asteroid belt will be even more forbidding and it is uncertain whether even the severely depleted present asteroid belt could form under these conditions. The perturbations of Jupiter and Saturn are quite sensitive to their semi-major axes and decrease when the planets' heliocentric distances are increased to allow for protoplanet migration. It is possible, though not clearly demonstrated, that this could produce a depleted asteroid belt but permit formation of a system of terrestrial planet embryos on a approximately 10(6)-year timescale, initially by nonrunaway growth and transitioning to runaway growth after approximately 10(5) years. The calculations reported here are valid under the condition that the relative velocities of the bodies are determined only by Jupiter and Saturn perturbations and by gas drag, with no mutual perturbations between planetesimals. If, while subject to these conditions, the bodies become large enough for their mutual perturbations to influence their velocity and size evolution significantly, the problem becomes much more complex. This problem is under investigation.     

Making other earths: dynamical simulations of terrestrial planet formation and water delivery

We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×10²⁴ g), and vary in mass between 0.23M_⊕ and 3.85M_⊕. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M_⊕ protoplanets which can serve as the accretion seeds for 2+M_⊕ planets with large water contents.     

Meteoritical and dynamical constraints on the growth mechanisms and formation times of asteroids and Jupiter

Thermal models and radiometric ages for meteorites show that the peak temperatures inside their parent bodies were closely linked to their accretion times. Most iron meteorites come from bodies that accreted <0.5 Myr after CAIs formed and were melted by ²⁶Al and ⁶⁰Fe, probably inside 2 AU. Rare carbon-rich differentiated meteorites like ureilites probably also come from bodies that formed <1 Myr after CAIs, but in the outer part of the asteroid belt. Chondrite groups accreted intermittently from diverse batches of chondrules and other materials over a 4 Myr period starting 1 Myr after CAI formation when planetary embryos may already have formed at ∼1 AU. Meteorite evidence precludes accretion of late-forming chondrites on the surface of early-formed bodies; instead chondritic and non-chondritic meteorites probably formed in separate planetesimals. Maximum metamorphic temperatures in chondrite groups are correlated with mean chondrule age, as expected if ²⁶Al and ⁶⁰Fe were the predominant heat sources. Because late-forming bodies could not accrete close to large, early-formed bodies, planetesimal formation may have spread across the nebula from regions where the differentiated bodies formed. Dynamical models suggest that the asteroids could not have accreted in the main belt if Jupiter formed before the asteroids. Therefore Jupiter probably reached its current mass >3-5 Myr after CAIs formed. This precludes formation of Jupiter via a gravitational instability <1 Myr after the solar nebula formed, and strongly favors core accretion. Jupiter probably formed too late to make chondrules by generating shocks directly, or indirectly by scattering Ceres-sized bodies across the belt. Nevertheless, shocks formed by gravitational instabilities or Ceres-sized bodies scattered by planetary embryos may have produced some chondrules. The minimum lifetime for the solar nebula of 3-5 Myr inferred from the total spread of CAI and chondrule ages may exceed the median lifetime of 3 Myr for protoplanetary disks, but is well within the 1-10 Myr observed range. Shorter formation times for extrasolar planets may help to explain their unusual orbits compared to those of solar giant planets.     

Accretion rates of planetesimals by protoplanets embedded in nebular gas

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals’ eccentricity, inclination, planet’s radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals’ kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with 1-10 Earth masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.     

Model consideration of the bombardment event of the asteroidal belt by the planetesimals scattered from the Jupiter zone

The temporal evolutions of the planetesimals scattered from the Jupiter zone for different masses of the proto-Jupiter [(a) 0.1 and (b) 1.0 of the present mass] are investigated. Due to the combined effects of the orbital evolution of the planetesimals and the elimination of these projectiles either via impact capture or injection into escape velocity by the outer planets, the whole scattering process lasts about 10⁸ yr for case (a) and about 10⁷ yr for case (b). The longer time scale may be a good estimate for the accretion time interval of Jupiter while the shorter one (10⁷) gives the upper time limit of the late heavy-bombardment epoch of the terrestrial planets due to planetesimals scattered from the Jupiter zone. The limiting value of the encounter velocity U at the end of the scattering process is ≈0.6. Consideration of the collisional interaction of these projectiles with the asteroids indicates that the corresponding bombardment effect could be rather appreciable. Also, the asteroids on the inner edge of the main asteroid belt would have been bombarded more severely than those on the outer edge. From this point of view, the structure of the asteroidal belt could be affected significantly not only by Jupiter's gravitational perturbation effect but also by its early scattering process.     

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.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region