Dust to planetesimals: Settling and coagulation in the solar nebula

The behavior of solid particles in a low-mass solar nebula during settling to the central plane and the formation of planetesimals is examined. Gravitational instability in a dust layer and collisional accretion are considered as possible mechanisms of planetesimal formation. Non-Keplerian rotation of the nebula results in shear between the gas and a dust layer. This shear produces turbulence within the layer which inhibits gravitational instability, unless the mean particle size exceeds a critical value, ～1 cm at 1 AU. The size requirement is less stringent at larger heliocentric distances, suggesting a possible difference in planetesimal formation mechanisms between the inner and outer nebula. Coagulation of grains during settling is expected in the solar nebula environment. Van der Waals forces appear adequate to produce centimeter-sized aggregates. Growth is primarily due to sweepup of small particles by larger ones due to size-dependent settling velocities. A numerical model for computing simultaneous coagulation and settling is described. Relative velocities are determined by gas drag and the non-Keplerian rotation of the nebula. The settling is very nonhomologous. Most of the solid matter reaches the central plane as centimeter-sized aggregates in a few times 10³ revolutions, but some remains suspended in the form of fine dust. Drag-induced relative velocities result in collisions. The growth of bodies in the central plane is initially rapid. After sizes reach ～10³ cm, relative velocities decrease and the growth rate declines. Gas drag rapidly damps the out-of-plane motions of these intermediate-sized bodies. They settle into a thin layer which is subject to gravitational instability. Kilometer-sized planetesimals are formed by this composite process.     

Nebular shock waves generated by planetesimals passing through Jovian resonances: Possible sites for chondrule formation

Abstract— The primordial asteroid belt contained at least several hundred and possibly as many as 10,000 bodies with diameters of 1000 km or larger. Following the formation of Jupiter, nebular gas drag combined with passage of such bodies through Jovian resonances produced high eccentricities (e = 0.3‐0.5), low inclinations (i < 0.5°), and, therefore, high velocities (3–10 km/s) for “resonant” bodies relative to both nebular gas and non‐resonant planetesimals. These high velocities would have produced shock waves in the nebular gas through two mechanisms. First, bow shocks would be produced by supersonic motion of resonant bodies relative to the nebula. Second, high‐velocity collisions of resonant bodies with non‐resonant bodies would have generated impact vapor plume shocks near the collision sites. Both types of shocks would be sufficient to melt chondrule precursors in the nebula, and both are consistent with isotopic evidence for a time delay of ?1‐1.5 Myr between the formation of CAIs and most chondrules. Here, initial simulations are first reported of impact shock wave generation in the nebula and of the local nebular volumes that would be processed by these shocks as a function of impactor size and relative velocity. Second, the approximate maximum chondrule mass production is estimated for both bow shocks and impact‐generated shocks assuming a simplified planetesimal population and a rate of inward migration into resonances consistent with previous simulations. Based on these initial first‐order calculations, impact‐generated shocks can explain only a small fraction of the minimum likely mass of chondrules in the primordial asteroid belt (?10²⁴‐10²⁵g). However, bow shocks are potentially a more efficient source of chondrule production and can explain up to 10–100 times the estimated minimum chondrule mass.     

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.     

Eccentricity growth of planetesimals in a self-gravitating protoplanetary disc

We investigate the orbital evolution of planetesimals in a self-gravitating circumstellar disc in the size regime (∼1–5000 km) where the planetesimals behave approximately as test particles in the disc's non-axisymmetric potential. We find that the particles respond to the stochastic, regenerative spiral features in the disc by executing large random excursions (up to a factor of 2 in radius in ∼1000 yr), although typical random orbital velocities are of the order of one tenth of the Keplerian speed. The limited time frame and small number of planetesimals modelled do not permit us to discern any net direction of planetesimal migration. Our main conclusion is that the high eccentricities (∼0.1) induced by interaction with spiral features in the disc is likely to be highly unfavourable to the collisional growth of planetesimals in this size range while the disc is in the self-gravitating regime. Thus if , as recently argued by Rice et al., the production of planetesimals gets under way when the disc is in the self-gravitating regime (either at smaller planetesimal size scales, where gas drag is important, or via gravitational fragmentation of the solid component), the planetesimals thus produced would not be able to grow collisionally until the disc ceases to be self-gravitating. It is unclear, however, given the large amplitude excursions undergone by planetesimals in the self-gravitating disc, whether they would be retained in the disc throughout this period, or whether they would instead be lost to the central star.     

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.     

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.     

The origin of the Kuiper Belt high-inclination population

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

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.     

The planetesimal bow shock model for chondrule formation: A more quantitative assessment of the standard (fixed Jupiter) case

Abstract– One transient heating mechanism that can potentially explain the formation of most meteoritic chondrules 1–3 Myr after CAIs is shock waves produced by planetary embryos perturbed into eccentric orbits via resonances with Jupiter following its formation. The mechanism includes both bow shocks upstream of resonant bodies and impact vapor plume shocks produced by high‐velocity collisions of the embryos with small nonresonant planetesimals. Here, we investigate the efficiency of both shock processes using an improved planetesimal accretion and orbital evolution code together with previous simulations of vapor plume expansion in the nebula. Only the standard version of the model (with Jupiter assumed to have its present semimajor axis and eccentricity) is considered. After several hundred thousand years of integration time, about 4–5% of remaining embryos have eccentricities greater than about 0.33 and shock velocities at 3 AU exceeding 6 km s^?1, currently considered to be a minimum for melting submillimeter‐sized silicate precursors in bow shocks. Most embryos perturbed into highly eccentric orbits are relatively large—half as large as the Moon or larger. Bodies of this size could yield chondrule‐cooling rates during bow shock passage compatible with furnace experiment results. The cumulative area of the midplane that would be traversed by highly eccentric embryos and their associated bow shocks over a period of 1–2 Myr is <1% of the total area. However, future simulations that consider a radially migrating Jupiter and alternate initial distributions of the planetesimal swarm may yield higher efficiencies.     

Orbital resonances in the solar nebula: Implications for planetary accretion

The influence of gas drag and gravitational perturbations by a planetary embryo on the orbit of a planetesimal in the solar nebula was examined. Non-Keplerian rotation of the gas causes secular decay of the orbit. If the planetesimal's orbit is exterior to the perturber's, resonant perturbations oppose this drag and can cause it to be trapped in a stable orbit at a commensurability of order j/(j + 1), where j is an integer. Numerical and analytical demonstrations show that resonant trapping occurs for wide ranges of perturbing mass, planetesimal size, and j. Induced eccentricities are large, causing overlap of orbits for bodies in different resonances with j > 2. Collisions between planetesimals in different resonances, or between resonant and nonresonant bodies, result in their disruption. Fragments smaller than a critical size can pass through resonances under the influence of drag and be accreted by the embryo. This effect speeds accretion and tends to prevent dynamical isolation of planetary embryos, making gas-rich scenarios for planetary formation more plausible.     

Towards initial mass functions for asteroids and Kuiper Belt Objects

Our goal is to understand primary accretion of the first planetesimals. Some examples are seen today in the asteroid belt, providing the parent bodies for the primitive meteorites. The primitive meteorite record suggests that sizeable planetesimals formed over a period longer than a million years, each of which being composed entirely of an unusual, but homogeneous, mixture of millimeter-size particles. We sketch a scenario that might help explain how this occurred, in which primary accretion of 10-100 km size planetesimals proceeds directly, if sporadically, from aerodynamically-sorted millimeter-size particles (generically “chondrules”). These planetesimal sizes are in general agreement with the currently observed asteroid mass peak near 100 km diameter, which has been identified as a “fossil” property of the pre-erosion, pre-depletion population. We extend our primary accretion theory to make predictions for outer Solar System planetesimals, which may also have a preferred size in the 100 km diameter range. We estimate formation rates of planetesimals and explore parameter space to assess the conditions needed to match estimates of both asteroid and Kuiper Belt Object (KBO) formation rates. For parameters that satisfy observed mass accretion rates of Myr-old protoplanetary nebulae, the scenario is roughly consistent with not only the “fossil” sizes of the asteroids, and their estimated production rates, but also with the observed spread in formation ages of chondrules in a given chondrite, and with a tolerably small radial diffusive mixing during this time between formation and accretion. As previously noted, the model naturally helps explain the peculiar size distribution of chondrules within such objects. The optimum range of parameters, however, represents a higher gas density and fractional abundance of solids, and a smaller difference between Keplerian and pressure-supported orbital velocities, than “canonical” models of the solar nebula. We discuss several potential explanations for these differences. The scenario also produces 10-100 km diameter primary KBOs, and also requires an enhanced abundance of solids to match the mass production rate estimates for KBOs (and presumably the planetesimal precursors of the ice giants themselves). We discuss the advantages and plausibility of the scenario, outstanding issues, and future directions of research.     

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.     

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

Origin of scattered disk resonant TNOs: Evidence for an ancient excited Kuiper belt of 50 AU radius

Resonance occupation of trans-neptunian objects (TNOs) in the scattered disk (>48 AU) was investigated by integrating the orbits of 85 observed members for 4 Gyr. Twenty seven TNOs were locked in the 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4 resonances. We then explored mechanisms for the origin of the resonant structure in the scattered disk, in particular the long-term 9:4, 5:2, and 8:3 resonant TNOs (median 4 Gyr), by performing large scale simulations involving Neptune scattering and planetary migration over an initially excited planetesimals disk (wide range of eccentricities and inclinations). To explain the formation of Gyr-resident populations in such distant resonances, our results suggest the existence of a primordial planetesimal disk of at least 45-50 AU radius that suffered a dynamical perturbation leading to 0.1-0.3 or greater eccentricities and a range of inclinations up to ∼20° during early stages of the Solar System history, before planetary migration.     

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.     

Particles in the nebular midplane: Collective effects and relative velocities

Abstract– In the absence of global turbulence, solid particles in the solar nebula tend to settle toward the midplane, forming a layer with enhanced solids/gas ratio. Shear relative to the surrounding pressure‐supported gas generates turbulence within the layer, inhibiting further settling and preventing gravitational instability. Turbulence and size‐dependent drift velocities cause collisions between particles. Relative velocities between small grains and meter‐sized bodies are typically about 50 m s^?1 for isolated particles; however, in a dense particle layer, collective effects alter the motion of the gas near the midplane. Here, we develop a numerical model for the coupled motions of gas and particles of arbitrary size, based on the assumption that turbulent viscosity transfers momentum on the scale of the Ekman length. The vertical distribution of particles is determined by a balance between settling and turbulent diffusion. Self‐consistent distributions of density, turbulent velocities, and radial fluxes of gas and particles of different sizes are determined. Collective effects generate turbulence that increases relative velocities between small particles, but reduce velocities between small grains and bodies of decimeter size or larger by bringing the layer’s motion closer to Keplerian. This effect may alleviate the “meter‐size barrier” to collisional growth of planetesimals.     

The rotation rates of accreting planetesimals

There are obtained upper limits for the relative velocity at infinity of accreting planetesimals for a nearly constant mass of the largest accreting planetesimal and also in the case of variable mass. We conclude, that while the larger planets cannot be brought to the stage of rotational instability by stochastic collisions, the asteroids could be brought. provided that the relative velocities in the asteroid belt were larger than about 2 km s^–1.     

Sweeping of the Jovian resonances and the evolution of the asteroids

During the formation of the solar system the variation of the gravitational field produced by removal of a nebula with its moderately massive accretion disk led to sweeping of the Jovian commensurability resonances through the asteroid zone. This process produced increased eccentricities and random velocities of the early planetesimals which resulted in collisional comminution rather than accretion. The existence of the asteroids, their low mass density, and their high relative velocities are interpreted as due to disruption of the accretion processes of the terrestrial planets by the influence of Jupiter.     

Chondrule formation by repeated evaporative melting and condensation in collisional debris clouds around planetesimals

Abstract– A synthesis of previous work leads to a model of chondrule formation that involves periodic melting of dispersed dust in debris clouds that were generated by collisions between chondritic planetesimals. I suggest that chondrules formed by the passage of nebular shock waves through these dust clumps, which temporarily surrounded disrupted planetesimals. Type I chondrules formed by more intense evaporative heating of fewer particles in tenuous clumps, or at the edges of dense clumps, and type II chondrules formed by less intense evaporative heating of more particles deeper within dense clumps. Chondrules reaccreted by self‐gravity into the planetesimals, mixing with less heated dust and rock. This process of disruption, melting, and reaccretion could have repeated many times. In this way, chondrite components of various origins and thermal histories could remain preserved in planetesimals as a distinctive mix of materials for extended periods of time, while still allowing for a repetitive melting process that converted some of the planetesimal debris into chondrules. I also suggest that during chondrule formation, the inner solar nebula gas was evolving by the gradual incorporation and heating of icy bodies depleted in ¹⁶O, causing a general increase in gaseous Δ¹⁷O with time in most places, especially close to the “snow line.” In this model, early formed type I chondrules in C chondrites with lower Δ¹⁷O values were produced inside the snow line, and later formed type I and type II chondrules in C and O chondrites with higher Δ¹⁷O values were created nearer the snow line after it had moved closer to the young Sun.