A semi-analytic model for oligarchic growth

A new semi-analytic model for the oligarchic growth phase of planetary accretion is developed. The model explicitly calculates damping and excitation of planetesimal eccentricities e and inclinations i due to gas drag and perturbations from embryos. The effects of planetesimal fragmentation, enhanced embryo capture cross sections due to atmospheres, inward planetesimal drift, and embryo-embryo collisions are also incorporated. In the early stages of oligarchic growth, embryos grow rapidly as e and i fall below their equilibrium values. The formation of planetesimal collision fragments also speeds up embryo growth as fragments have low-e, low-i orbits, thereby optimizing gravitational focussing. At later times, the presence of thick atmospheres captured from the nebula aids embryo growth by increasing their capture cross sections. Planetesimal drift due to gas drag can lead to substantial inward transport of solid material. However, inward drift is greatly reduced when embryo atmospheres are present, as the drift timescale is no longer short compared to the accretion timescale. Embryo-embryo collisions increase embryo growth rates by 50% compared to the case where growth is solely due to accretion of planetesimals. Formation of 0.1-Earth-mass protoplanets at 1 AU and 10-Earth-mass cores at 5 AU requires roughly 0.1 and 1 million years respectively, in a nebula where the local solid surface density is 7 g cm⁻² at each of these locations.     

Radial diffusion rate of planetesimals due to gravitational encounters

We study the rate of radial diffusion of planetesimals due to mutual gravitational encounters under Hill’s approximations in the three-body problem. Planetesimals orbiting a central star radially migrate inward and outward as a result of mutual gravitational encounters and transfer angular momentum. We calculate the viscosity in a disk of equal-sized planetesimals due to their mutual gravitational encounters using three-body orbital integrations, and obtain a semianalytic expression that reproduces the numerical results. We find that the viscosity is independent of the velocity dispersion of planetesimals when the velocity dispersion is so small that Kepler shear dominates planetesimals’ relative velocities. On the other hand, in high-velocity cases where random velocities dominate the relative velocities, the viscosity is a decreasing function of the velocity dispersion, and is found to agree with previous estimates under the two-body approximation neglecting the solar gravity. We also calculate the rate of radial diffusion of planetesimals due to gravitational scattering by a massive protoplanet. Using these results, we discuss a condition for formation of nonuniform radial surface density distribution of planetesimals by gravitational perturbation of an embedded protoplanet.     

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.     

Envelope instability in giant planet formation

We compute the growth of isolated gaseous giant planets for several values of the density of the protoplanetary disk, several distances from the central star and two values for the (fixed) radii of accreted planetesimals. Calculations were performed in the frame of the core instability mechanism and the solids accretion rate adopted is that corresponding to the oligarchic growth regime. We find that for massive disks and/or for protoplanets far from the star and/or for large planetesimals, the planetary growth occurs smoothly. However, notably, there are some cases for which we find an envelope instability in which the planet exchanges gas with the surrounding protoplanetary nebula. The timescale of this instability shows that it is associated with the process of planetesimals accretion. The presence of this instability makes it more difficult the formation of gaseous giant planets.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

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.     

Protoatmospheres and Surface Environment of Protoplanets

Protoatmospheres and surface environment of terrestrial protoplanets during the oligarchic accretion phase and the giant impacts phase are discussed from theoretical points of view. Mars-sized protoplanets form during the stage of the oligarchic growth. Since protoplanets are formed from more or less ‘local’ planetesimals, the surface environment of the accreting protoplanets depends on availability of volatile material in planetesimals. Even if no volatile-bearing planetesimals are available, a gravitationary captured solar composition atmosphere is formed during accretion. In such cases the surface temperature is always kept under the melting temperature of mantle silicate and only a subsurface magma ocean is formed. Core formation proceeds under dry conditions, and volatile elements are not partitioned into metallic iron. Accretion of water-bearing planetesimals results in impact degassing. A surface hydrous magma ocean forms in response to the thermal blanketing effect of the proto-atmosphere. Then, some volatile materials dissolve into the magma ocean. If we consider reaction with metallic iron, the proto-atmosphere is likely to be rich in hydrogen. In addition, a large amount of hydrogen may be partitioned into metallic iron under high pressure, and delivered to the core. In the stage of giant impacts, both dry and water-bearing protoplanets collide on the proto-Earth. Substantial amount of proto-atmosphere (including water vapor) survives giant impacts. Moreover, giant impacts on protoplanets with oceans result in relative concentration of water against other gases.     

Gaseous drag and planetary formation by accretion

We have made numerical experiments of the collisional and gravitational interaction of a planetesimal swarm in the early Solar System. In particular we study the dynamical evolution of an initial population of kilometer-size planetesimals subject to collisions (accretion, rebound, cratering, and catastrophic fragmentation). This study is based on a Monte-Carlo statistical method and provides the mass and velocity distributions of the planetesimal swarm as a function of time as well as their distribution in heliocentric distance. Several experiments have been performed and three of them are presented here. They simulate the accretional growth of numerous planetesimals in the absence (or presence) of gaseous drag, with (or without) one larger embryo among them, and with (or without) a size gradient. The results show that (i) for a population of planetesimals submitted to a negative gradient in size as the heliocentric distance increases, the outer planetesimals spiral toward the Sun faster than inner ones, leading after some time to an accumulation of bodies inside the cloud which allows the formation of an embryo; (ii) the growth of one embryo among a population of planetesimals is accelerated by the presence of gas and is warranted as long as its feeding zone is fed by the inward flow of planetesimals due to gas drag. These results offer some complementary new insights in the understanding of the accretional formation of 4–5 terrestrial planets instead of the numerous Moon-size planets generally found in numerical experiments.     

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.     

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.     

Cosmogony of the solar system

In Sections 1–6, we determine an approximate analytical model for the density and temperature distribution in the protoplanetary could. The rotation of the planets is discussed in Section 7 and we conclude that it cannot be determined from simple energy conservation laws.The velocity of the gas of the protoplanetary cloud is found to be smaller by about 5×10³ cm s^–1 in comparison to the Keplerian circular velocity. If the radius of the planetesimals is smaller than a certain limitr ₁, they move together with the gas. Their vertical and horizontal motion for this case is studied in Sections 8 and 9.As the planetesimals grow by accretion their radius becomes larger thanr ₁ and they move in Keplerian orbits. As long as their radius is betweenr ₁ and a certain limitr ₂ their gravitational interaction is negligible. In Section 10, we study the accretion for this case.In Section 11, we determine the change of the relative velocities due to close gravitational encounters. The principal equations governing the late stages of accretion are deduced in Section 12, In Section 13 there are obtained approximate analytical solutions.The effect of gas drag and of collisions is studied in Sections 14 and 15, respectively. Numerical results and conclusions concerning the last and principal stage of accretion are drawn in Section 16.     

Core formation in giant gaseous protoplanets

Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 MSaturn to 10 MJupiter. Giant protoplanets with masses of 5 MJupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ∼0.25 M_⊕ and ∼0.5 M_⊕, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M_⊕. Saturn's core mass is found to be ∼8 M_⊕.     

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.     

An efficient, low-velocity, resonant mechanism for capture of satellites by a protoplanet

Numerical simulations of the gravitational scattering of planetesimals by a protoplanet reveal that a significant fraction of scattered planetesimals can become trapped as so-called quasi-satellites in heliocentric 1:1 co-orbital resonance with the protoplanet. While trapped, these resonant planetesimals can have deep low-velocity encounters with the protoplanet that result in temporary or permanent capture onto highly eccentric prograde or retrograde circumplanetary orbits. The simulations include solar nebula gas drag and use planetesimals with diameters ranging from ∼1 to ∼1000 km. Initial protoplanet eccentricities range from ep=0 to 0.15 and protoplanet masses range from 300 Earth-masses (M_⊕) down to 0.1M_⊕. This mass range effectively covers the final masses of all planets currently thought to be in possession of captured satellites—Jupiter, Saturn, Neptune, Uranus, and Mars. For protoplanets on moderately eccentric orbits (ep?0.1) most simulations show from 5-20% of all scattered planetesimals becoming temporarily trapped in the quasi-satellite co-orbital resonance. Typically, 20-30% of the temporarily trapped quasi-satellites of all sizes came within half the Hill radius of the protoplanet while trapped in the resonance. The efficiency of the resonance trapping combined with the subsequent low-velocity circumplanetary capture suggests that this trapped-to-captured transition may be important not only for the origin of captured satellites but also for continued growth of protoplanets.     

Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals

The final stage of the accretion of Uranus and Neptune is numerically investigated. The four Jovian planets are considered with Jupiter and Saturn assumed to have reached their present sizes, whereas Uranus and Neptune are taken with initial masses 0.2 of their present ones. Allowance is made for the orbital variation of the Jovian planets due to the exchange of angular momentum with interacting bodies (“planetesimals”). Two possible effects that may have contributed to the accretion of Uranus and Neptune are incorporated in our model: (1) an enlarged cross section for accretion of incoming planetesimals due to the presence of extended gaseous envelopes and/or circumplanetary swarms of bodies; and (2) intermediate protoplanets in mid-range orbits between the Jovian planets. Significant radial displacements are found for Uranus and Neptune during their accretion and scattering of planetesimals. The orbital angular momentum budgets of Neptune, Uranus, and Saturn turn out to be positive; i.e., they on average gain orbital angular momentum in their interactions with planetesimals and hence they are displaced outwardly. Instead, Jupiter as the main ejector of bodies loses orbital angular momentum so it moves sunward. The gravitational stirring of planetesimals caused by the introduction of intermediate protoplanets has the effect that additional solid matter is injected into the accretion zones of Uranus and Neptune. For moderate enlargements of the radius of the accretion cross section (2–4 times), the accretion time scale of Uranus and Neptune are found to be a few 10⁸ years and the initial amount of solid material required to form them of a few times their present masses. Given the crucial role played by the size of the accretion cross section, questions as to when Uranus and Neptune acquired their gaseous envelopes, when the envelopes collapsed onto the solid cores, and how massive they were are essential in defining the efficiency and time scale of accretion of the two outer Jovian planets.     

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.     

Evolution of Planetesimal Velocities Based on Three-Body Orbital Integrations and Growth of Protoplanets

We obtain the viscous stirring and dynamical friction rates of planetesimals with a Rayleigh distribution of eccentricities and inclinations, using three-body orbital integration and the procedure described by Ohtsuki (1999, Icarus137, 152), who evaluated these rates for ring particles. We find that these rates based on orbital integrations agree quite well with the analytic results of Stewart and Ida (2000, Icarus 143, 28) in high-velocity cases. In low-velocity cases where Kepler shear dominates the relative velocity, however, the three-body calculations show significant deviation from the formulas of Stewart and Ida, who did not investigate the rates for low velocities in detail but just presented a simple interpolation formula between their high-velocity formula and the numerical results for circular orbits. We calculate evolution of root mean square eccentricities and inclinations using the above stirring rates based on orbital integrations, and find excellent agreement with N-body simulations for both one- and two-component systems, even in the low-velocity cases. We derive semi-analytic formulas for the stirring and dynamical friction rates based on our numerical results, and confirm that they reproduce the results of N-body simulations with sufficient accuracy. Using these formulas, we calculate equilibrium velocities of planetesimals with given size distributions. At a stage before the onset of runaway growth of large bodies, the velocity distribution calculated by our new formulas are found to agree quite well with those obtained by using the formulas of Stewart and Ida or Wetherill and Stewart (1993, Icarus106, 190). However, at later stages, we find that the inclinations of small collisional fragments calculated by our new formulas can be much smaller than those calculated by the previously obtained formulas, so that they are more easily accreted by larger bodies in our case. The results essentially support the previous results such as runaway growth of protoplanets, but they could enhance their growth rate by 10-30% after early runaway growth, where those fragments with low random velocities can significantly contribute to rapid growth of runaway bodies.     

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.     

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.     

Relative velocities among accreting planetesimals in binary systems: The circumprimary case

We investigate classical planetesimal accretion in a binary star system of separation ab?50 AU by numerical simulations, with particular focus on the region at a distance of 1 AU from the primary. The planetesimals orbit the primary, are perturbed by the companion and are in addition subjected to a gas drag force. We concentrate on the problem of relative velocities Δv among planetesimals of different sizes. For various stellar mass ratios and binary orbital parameters we determine regions where Δv exceed planetesimal escape velocities vesc (thus preventing runaway accretion) or even the threshold velocity vero for which erosion dominates accretion. Gaseous friction has two crucial effects on the velocity distribution: it damps secular perturbations by forcing periastron alignment of orbits, but at the same time the size-dependence of this orbital alignment induces a significant Δv increase between bodies of different sizes. This differential phasing effect proves very efficient and almost always increases Δv to values preventing runaway accretion, except in a narrow eb?0 domain. The erosion threshold Δv>vero is reached in a wide (ab,eb) space for small <10-km planetesimals, but in a much more limited region for bigger ?50-km objects. In the intermediate vesc<Δv<vero domain, a possible growth mode would be the type II runaway growth identified by Kortenkamp et al. [Kortenkamp, S., Wetherill, G., Inaba, S., 2001. Science 293, 1127-1129].