The formation of the atmospheres of the terrestrial planets by impact

It is generally supposed that the atmospheres of the terrestrial planets were formed by secondary degassing processes. We propose, instead, that they are of primary origin, forming as an immediate and necessary consequence of the final stages of planetary accretion. Once the planetary embryo reached a critical size, the impacting material began to vaporize. The atmosphere, so created, then decelerated other impacting material, thus limiting the rate of atmospheric growth. We show that, given reasonable assumptions concerning the chemical composition of the impacting material, an acceptable model for the early atmosphere of the Earth, and the present atmospheres of Venus and Mars results.A discussion of the noble gas data for the terrestrial atmosphere indicates that these can be readily reconciled with an impact origin.     

Capture of planetesimals by the giant planets

We investigate the possibility of gravitational capture of planetesimals as temporary or permanent satellites of Uranus and Neptune during the process of planetary growth. The capture mechanism is based in the enhancement of the Hill's sphere of action not only due to the mass acquired by the planet, but also by the variation of the planet-Sun distance as a consequence of the scattering of planetesimals by the planets of the outer solar system. Our calculations indicate that satellite capture was very important, specially during the first stages of the accretion process, contributing in a significant way to the planetary growth.     

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 formation surrounding close binary stars

Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M_⊙, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities.     

系外类地行星的内部结构及大气成分研究进展

系外类地行星是目前搜寻地外生命的主要目标.随着观测仪器的发展,现在已经能探测到低于10个地球质量的系外行星.该文简要回顾了系外类地行星的形成与演化,介绍了当前研究它们内部结构的模型和方法,以及由此得出的类地行星质量-半径关系.同时,对应不同的行星初始物质成分,讨论了各种可能的大气结构.最后介绍了未来的空间任务在相关方面的工作.     

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.     

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.     

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.     

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.     

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.     

Solar wind origin of 36Ar on Venus

A hypothesis is considered in which the ³⁶Ar found on Venus is of solar origin. This possibility is quantitatively discussed within the framework of present theories of planetary accumulation by sweep up of planetesimals under gas-free conditions. Solar wind implantation of ³⁶Ar would take place by irradiation of accumulating material during the first ≈10⁵ years of planetary growth, provided that the flux of solar wind was enhanced by a factor of ≈100 at that time. Enrichment of Venus in implanted gas would be a consequence of the irradiated material being initially confined to the innermost edge of the radially opaque circusolar planetesimal disk predicted by these theories. The observed atmospheric data require a Ne/Ar fractionation by a factor of ≈100 during the planetesimal stage. It is also necessary that there be very little mixing of irradiated planetesimals from the inner edge of disk to the distance (≈1 AU) at which the Earth formed. The hypothesis can be tested by measurement of the abundance of Kr and Xe in the Venus atmosphere. Venera data indicate a terrestrial ³⁶Ar/Kr ratio, in disagreement with the solar wind hypothesis. In contrast, the Pioneer experiments find a lower limit to this ratio, well above the terrestrial value, that is compatible with the hypothesis. These experiments also show that Venus' ³⁶Ar/Xe ratio does not correspond to the so-called “planetary” trapped inert gas composition. The inert of Venus could be related to result of admixture of gas with solar composition. The inert gas on Venus could be related to that found in enstatite chondrites.     

N-Body Simulations of Late Stage Planetary Formation with a Simple Fragmentation Model

We present results of two-dimensional gravitationalN-body simulations of the late stage of planetary formation. This stage is characterized by the direct accretion of hundreds of lunar-sized planetesimals into planetary bodies. Our simulation code is based on the Hermite Individual Timestep integration algorithm, and gravitational interactions among all bodies are included throughout the simulations. We compare our simulation with earlier works that do not include all interactions, and we find very good agreement. A previously published collisional fragmentation model is included in our simulation to study the effects of the production of fragments on the subsequent evolution of the larger planetary bodies. It is found that for realistic two-body collisions that, according to this model, both bodies will suffer fragmentation, and that the outcome of the collision will be a relatively large core containing most of the mass and a few small fragments. We present the results of simulations that include this simple fragmentation model. They indicate that the presence of small fragments have only a small effect on the growth or orbital evolution of the large planet-sized bodies.     

Computer simulation of planet formation in a binary star system: Terrestrial planets

A model of planetary formation in a binary system with a small relative mass of primary is computed on the assumption of a mass transfer from the less massive component to the more massive one with no mass and angular momentum carried away from the system under consideration. At the last stage of mass transfer the condensed Moon-like objects (planetoids) are ejected through the inner Lagrange point of the primary Roche lobe with the outflow of gaseous matter.The whole system is considered in the plane of binary star rotation. Newtonian equations of motion are integrated with the initial conditions for the planetoids referred to as the coordinates and velocity of the inner Lagrangian point at the moments of planetoid ejections, all the pairwise gravitational interactions being included in computations but without a gas-drag. The mass transfer ceases at the primary relative mass 10^–3 which corresponds to the present Sun-Jupiter system. The total mass of planetoids approximates that of the terrestrial planets. Those are formed through coagulation of the planetoids with the effective radius of capture cross-section as an input parameter in the computer simulation. When the minimum separation between the pair of bodies becomes less than this radius they coalesce into a single body with their masses and momenta summed. If the effective radius value is under a certain limit the computer simulation yields the planetary system like that of terrestrial planets of the present Sun system.Numerical computations reveal the division of the planetoids into 4 groups along their distances from the Sun. Further, each group forms a single planet or a planet and a less massive body at the nearest orbits. The parameters of simulated planet orbits are close to the present ones and the interplanetary spacings are in accord with the Titius-Bode law.     

Core instability models of giant planet accretion – II. Forming planetary systems

We develop a simple model for computing planetary formation based on the core instability model for the gas accretion and the oligarchic growth regime for the accretion of the solid core. In this model several planets can form simultaneously in the disc, a fact that has important implications especially for the changes in the dynamic of the planetesimals and the growth of the cores since we consider the collision between them as a source of potential growth. The type I and type II migration of the embryos and the migration of the planetesimals due to the interaction with the disc of gas are also taken into account. With this model we consider different initial conditions to generate a variety of planetary systems and analyse them statistically. We explore the effects of using different type I migration rates on the final number of planets formed per planetary system such as on the distribution of masses and semimajor axis of extrasolar planets, where we also analyse the implications of considering different gas accretion rates. A particularly interesting result is the generation of a larger population of habitable planets when the gas accretion rate and type I migration are slower.     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.     

Making the Earth: Combining dynamics and chemistry in the Solar System

No terrestrial planet formation simulation completed to date has considered the detailed chemical composition of the planets produced. While many have considered possible water contents and late veneer compositions, none have examined the bulk elemental abundances of the planets produced as an important check of formation models. Here we report on the first study of this type. Bulk elemental abundances based on disk equilibrium studies have been determined for the simulated terrestrial planets of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58]. These abundances are in excellent agreement with observed planetary values, indicating that the models of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58] are successfully producing planets comparable to those of the Solar System in terms of both their dynamical and chemical properties. Significant amounts of water are accreted in the present simulations, implying that the terrestrial planets form “wet” and do not need significant water delivery from other sources. Under the assumption of equilibrium controlled chemistry, the biogenic species N and C still need to be delivered to the Earth as they are not accreted in significant proportions during the formation process. Negligible solar photospheric pollution is produced by the planetary formation process. Assuming similar levels of pollution in other planetary systems, this in turn implies that the high metallicity trend observed in extrasolar planetary systems is in fact primordial.     

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.     

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.     

Galactic civilizations: Population dynamics and interstellar diffusion

The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential and difference equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the planetary environments. Dispersal is fundamentally a diffusion process; a directed density-dependent diffusivity describes interstellar emigration. We concentrate on two models, the first describing zero population growth (ZPG) and the second which also includes local growth and saturation of a planetary population, and for which we find an asymptotic travelling wave solution. For both models the colonization wavefront expands slowly and uniformly, but only the frontier worlds are sources of further expansion. For nonlinear diffusion with growth and saturation, the colonization wavefront from the nearest independently arisen galactic civilization can have reached the Earth only if its lifetime exceeds 2.6 × 10⁶ years. If discretization can be neglected, the critical lifetime is 2.0 × 10⁷ years. For ZPG the corresponding number is 1.3 × 10¹⁰ years. These numerical results depend on our choices for the specific emigration rate, the distribution of colonizable worlds, and, in the second model, the population growth rate; but the dependence on these parameters is entrancingly weak. We conclude that the Earth is uncolonized not because interstellar spacefaring societies are rare, but because there are too many worlds to be colonized in the plausible lifetime of the colonization phase of nearby galactic civilizations. This phase is, we contend, eventually outgrown. We also conclude that, except possibly early in the history of the Galaxy, there are no very old galactic civilizations with a consistent policy of conquest of inhabited worlds; there is no Galactic Empire. There may, however, be abundant groups of ～10⁵ to 10⁶ worlds linked by a common colonial heritage. The radar and television announcement of an emerging technical society on Earth may induce a rapid response by nearby civilizations, thus newly motivated to reach our system directly rather than by diffusion.