Dynamic fission instability of dissipative protoplanets

All theories of fission require a catastrophic, dynamic phase in order to produce two separate bodies. We have used nonlinear numerical and linear analytical calculations to show that the dynamic fission instability probably does not occur in dissipative protoplanets. The numerical calculations were performed with a three-spatial-dimension hydrodynamical code, with the proto-planet represented by a fluid with a Murnaghan equation of state. The kinetic energy in the protoplanet (other than rigid body rotation) is dissipated throughout the evolution in order to simulate the effects of viscous dissipation. Protoplanets rotating above the limit for dynamic instability were given initial asymmetric density perturbations; in each case the asymmetry did not grow during a time on the order of the rotational period. This dynamical stability has been verified by including the dissipative terms in the tensor-virial equation analysis for the stability of a Maclaurin spheroid: the dynamic instability vanishes when the dissipative terms are included, while the secular instability (with a growth time much larger than the rotational period) remains. The result applies to bodies of radius R with a kinematic viscosity ν? 4 × 10¹³ (R/6400 km)²cm²sec^?1, and hence may be applicable to any terrestrial protoplanet which is not totally molten. Current thermal histories for the Earth predict a partially molten mantle with a viscosity greater than this critical value. Depending on the detailed rheology of the early Earth, our results appear to rule out the possibility of forming the Earth-Moon system through a dynamic fission instability.     

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.     

A model for planetesimal meltdown by 26Al and its implications for meteorite parent bodies

Abstract— The melting of planetesimals heated by ²⁶Al has been modelled using a new finite difference method that incorporates convection. As an example, we consider a planetesimal with a radius of 64 km, which accretes instantaneously at t = 0.75 Myr (after the formation of calcium‐aluminum‐rich inclusions) from cold (250 K) nebular dust with 50% porosity. At t = 0.9 Myr (T = 700 K), the planetesimal shrinks to a radius of 50 km due to sintering. At t = 1.2 Myr (T= 1425 K), the fully insulated interior, deeper than a few kilometers, starts to melt, and at t = 1.5 Myr (T = 1725 K), with 50% melting, convection starts. By t = 2 Myr, the planetesimal is a globe of molten, convecting slurry inside a thin residual crust. From about t = 2.5 Myr, the crust thickens rapidly as the power of ²⁶Al fades. Planetesimals probably melt in this manner when they accrete before t = 1.3 Myr and are large enough to insulate themselves (R >20 km for accretion at t = 0, rising to >80 km at t = 1.3 Myr). Melting behavior will also be affected by the level of ⁶⁰Fe in nebular dust, by the extent of devolatilization reactions and basalt segregation during heating, and by gradual accretion. The model suggests that a) the parent bodies of differentiated meteorites had accreted before about t = 1.5 to 2 Myr and before most chondritic parent bodies had formed, and b) that molten planetesimals may be a source for chondrule melt droplets.     

Planetesimal capture in the disk instability model

We follow the contraction and evolution of a typical Jupiter-mass clump created by the disk instability mechanism, and compute the rate of planetesimal capture during this evolution. We show that such a clump has a slow contraction phase lasting ∼3×¹⁰⁵ yr. By following the trajectories of planetesimals as they pass through the envelope of the protoplanet, we compute the cross-section for planetesimal capture at all stages of the protoplanet's evolution. We show that the protoplanet can capture a large fraction of the solid material in its feeding zone, which will lead to an enrichment of the protoplanet in heavy elements. The exact amount of this enrichment depends upon, but is not very sensitive to the size and random speed of the planetesimals.     

Tidal disruptions: A continuum theory for solid bodies

Although the theory of Roche 1847 for the tidal disruption limits of orbiting satellites assumes a fluid body, a length to diameter of exactly 2.07:1, and a particular body orientation, the theory is commonly applied to the satellites of the Solar System and to small asteroids and comets passing nearby a planet. Clearly these bodies are neither fluid nor generally are that elongated, so a more appropriate theory is needed. Here we present exact analytical results for the distortion and disruption limits of solid spinning ellipsoidal bodies subjected to tidal forces, using the Drucker-Prager strength model with zero cohesion. It is the appropriate model for dry granular materials such as sands and rocks, for rubble-pile asteroids and comets, and for larger satellites, asteroids and comets where the cohesion can be ignored. This study uses the same approach as the studies of spin limits for solid ellipsoidal bodies given in [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303]. It is a static theory that predicts conditions for breakup and predicts the nature of the deformations at the limit state, but does not track the dynamics of the body as it comes apart. The strength is characterized by a single material parameter associated with an angle of friction, which can range from zero to 90°. The case with zero friction angle has no shear strength whatsoever, so it is then the model of a fluid or gas. The case of 90° represents a material that cannot fail in shear, but still has zero tensile strength. Typical dry soils have angles of friction of 30°-40°. Since the static fluid case is included in the theory as a special case, the classical results of Roche [Roche, E.A., 1847. Acad. Sci. Lett. Montpelier. Mem. Section Sci. 1, 243-262] and Jeans [Jeans, J.H., 1917. Mem. R. Astron. Soc. London 62, 1-48] are included and re-derived in their entirety; but the general solid case has much more variety and applicability. We consider both the spin-locked case, appropriate for most satellites of the Solar System; and the zero spin case, a possible case for a passing stray body. Detailed plots of many special cases are presented, in terms of shape, orientation and mass densities. A very typical result gives a closest approach d=1.5(ρ/ρP)^1/3R in terms of the planet radius R, and the satellite and planet mass densities ρ and ρP. We also use the theory to distinguish between conditions allowing global shape changes leading to new equilibrium states, or those leading to complete disruption. We apply the theory to the potentially hazardous Asteroid 99942 Apophis due to pass very near the Earth in 2029, and conclude it is extremely unlikely to experience any tidal readjustments during its passage. The states of many of the satellites of the Solar System are compared to the theory, and we find that all are well within their tidal disruption limits for expected values of the internal friction.     

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.     

Studies of microparticle impact phenomena leading to the development of a highly sensitive micrometeoroid detector

A highly sensitive device has been developed for the detection of micrometeoroids with masses extending down to and below the expected classical solar radiation pressure limit. This limit occurs at a mass of about 10⁻¹⁶ kg for a particle density of 3 × 10³ kg m⁻³. The detector operates on the principle that charge (Q) is released when a hypervelocity microparticle impacts on a solid surface. This process has previously been investigated at particle velocities (ν) greater than 1 km sec⁻¹ and the empirical relationship Q = Km^αν^β obtained where m is the particle mass and K, α and β are constants. In the present study the validity of this relationship has been demonstrated for iron particles with masses from 10⁻¹⁶ to 10⁻¹³ kg in the velocity range 0.05–1.4 km sec⁻¹ impacting on a molybdenum target. A value of α ≈ 1 is indicated and a value of β of 3.2 ± 0.4 has been obtained.     

Against all odds? Forming the planet of the HD 196885 binary

HD 196885 Ab is the most ??extreme?? planet-in-a-binary discovered to date, whose orbit places it at the limit for orbital stability. The presence of a planet in such a highly perturbed region poses a clear challenge to planet-formation scenarios. We investigate this issue by focusing on the planet-formation stage that is arguably the most sensitive to binary perturbations: the mutual accretion of kilometre-sized planetesimals. To this effect we numerically estimate the impact velocities dv amongst a population of circumprimary planetesimals. We find that most of the circumprimary disc is strongly hostile to planetesimal accretion, especially the region around 2.6 AU (the planet??s location) where binary perturbations induce planetesimal-shattering dv of more than 1 kms^?1. Possible solutions to the paradox of having a planet in such accretion-hostile regions are (1) that initial planetesimals were very big, at least 250 km (2) that the binary had an initial orbit at least twice the present one, and was later compacted due to early stellar encounters (3) that planetesimals did not grow by mutual impacts but by sweeping of dust (the ??snowball?? growth mode identified by Xie et al., in Astrophys J 724:1153, 2010b), or (4) that HD 196885 Ab was formed not by core-accretion but by the concurrent disc instability mechanism. All of these 4 scenarios remain however highly conjectural.     

Numerical simulations of the differentiation of accreting planetesimals with 26Al and 60Fe as the heat sources

Abstract— Numerical simulations have been performed for the differentiation of planetesimals undergoing linear accretion growth with ²⁶Al and ⁶⁰Fe as the heat sources. Planetesimal accretion was started at chosen times up to 3 Ma after Ca‐Al‐rich inclusions (CAIs) were formed, and was continued for periods of 0.001–1 Ma. The planetesimals were initially porous, unconsolidated bodies at 250 K, but became sintered at around 700 K, ending up as compact bodies whose final radii were 20, 50, 100, or 270 km. With further heating, the planetesimals underwent melting and igneous differentiation. Two approaches to core segregation were tried. In the first, labelled A, the core grew gradually before silicate began to melt, and in the second, labelled B, the core segregated once the silicate had become 40% molten. In A, when the silicate had become 20% molten, the basaltic melt fraction began migrating upward to the surface, carrying ²⁶Al with it. The ⁶⁰Fe partitioned between core and mantle. The results show that the rate and timing of core and crust formation depend mainly on the time after CAIs when planetesimal accretion started. They imply significant melting where accretion was complete before 2 Ma, and a little melting in the deep interiors of planetesimals that accreted as late as 3 Ma. The latest melting would have occurred at <10 Ma. The effect on core and crust formation of the planetesimal's final size, the duration of accretion, and the choice of (⁶⁰Fe/⁵⁶Fe)_initial were also found to be important, particularly where accretion was late. The results are consistent with the isotopic ages of differentiated meteorites, and they suggest that the accretion of chondritic parent bodies began more than 2 or 3 Ma after CAIs.     

The penetration limit of thin films

One sensor of the Helios micrometeoroid experiment is covered by a thin film consisting of 3000 Å parylene and 750 Å aluminium. Micrometeoroids must penetrate this film before they are detected. In order to study the effects of the film on the detection of micrometeoroids simulation experiments were performed with iron, aluminium, glass and polyphenylene projectiles in the mass range of 5 × 10^?13g < m < 2 × 10^?10g and in the speed range of 1.5 km/sec <ν < 13 km/sec. The bulk densities of the projectiles ranged from 1.25 g/cm³ (polyphenylene) to 7.9 g/cm³ (iron). By measuring the speed of the projectiles before and after the film penetration the speed loss Δν caused by the film was determined. The angle of incidence was varied in three steps (0°, 30° and 60°). This deceleration strongly depends on the projectiles' densities: Vertically impacting iron projectiles of m = 10^?11g and ν₁ = 3 km/sec were subject to a relative speed loss of Δν/ν₁ = 4%, aluminium projectiles of the same mass and speed showed Δν/ν₁ = 8%, glass projectiles Δν/ν₁ = 9% and polyphenylene projectiles Δν/ν₁ = 14%. The total charge of the plasma produced upon impact on a gold target of a projectile which had penetrated the film before that was compared with the plasma produced by a projectile without a penetration. For iron projectiles these two signals did not differ significantly even at an angle of incidence of 60°. Whereas polyphenylene projectiles showed an attenuation of the charge signal by a factor of 10 after the penetration already at an angle of incidence of 0°. Polyphenylene projectiles impacting the film at an angle of incidence of 60° could no longer be detected behind the film. This experiment defined the penetration limit of the Helios film. Comparison with other penetration data yielded a penetration formula which is applicable to projectiles with diameters in the submicron to centimeter range. This penetration formula gives the penetration limit of a film as a function of the projectile's mass, speed and density.     

Cumulative subject index volumes 27–29

It is premature to establish a chronology for Mars and Mercury, relative to the known lunar chronology, to better than an order of magnitude. Lunar evidence neither requires nor excludes a “cataclysmic” episode of bombardment about 4.0 b.y. ago. Such a cataclysm might have resulted naturally from tidal disruption by a planet or collisional fragmentation in the asteroid belt of either a Uranus/Neptune-scattered planetesimal or a large asteroid, in which case any lunar cataclysm would have occurred as well on other planets. There is no independent evidence in Mariner 10 imagery for (or against) an early episodic bombardment on Mercury. Crater densities on plains units of the Moon, Mars, and Mercury have not been shown to be “strikingly similar” and do not imply, in the absence of definitive dynamical calculations of planetary impact rates of plausible populations of planetesimals, any similarity in the geological chronologies for those planets. Photogeological studies alone cannot determine absolute chronologies for planets. In combination with dynamical analyses, they can help us date to no better than a factor of 3 to 10 the formation of the Caloris Basin or the epoch when the Martian rivers ran.     

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.     

Formation of the galilean satellites in a gaseous nebula

A model for Galilean satellite formation was analyzed in which the satellites accrete in the presence of a dense, gaseous disk-shaped nebula and rapidly form optically thick, gravitationally bound primordial atmospheres. Upper-bound temperatures expected during accretion lead to partially differentiated structures for both Ganymede and Callisto, although with Ganymede much more differentiated than Callisto. When allowance is made for the aerodynamic breaking of infalling planetesimal fragments, lower surface temperatures result, and the amount of partial differentiation of Callisto is small, possibly approaching zero for a narrow size distribution of infalling planetesimals. The model is chosen to be consistent with the observed densities of the Galilean satellites and our current understanding of Jupiter formation. The retention of ices more volatile than H₂O is considered but not modeled in detail. A nominal nebula of ～0.1 Jupiter masses is constructed by consideration of likely surface density profiles and existing Jupiter collapse calculations. This nebula is optically thick (even if grain opacity is ignored) in both radial and vertical directions and has a temperature profile T ～ 3600 (R_J/R), where R_J is Jupiter's radius and R is the radial distance in the disk midplane. Satellites accrete very rapidly (dynamical time scales being 10²–10⁴ years) and their optically thick gaseous envelopes are unable to eliminate the heat of accretion by radiation. Water-saturated, convective, adiabatic envelopes form, through which planetesimals fall, break up, and partially disseminate their mass. The resulting satellite surface temperatures during accretion are calculated. Possible implications of these models for the subsequent evolution of Ganymede and Callisto are explored and it is suggested that the extensive differentiation undergone by Ganymede may provide the right environment for subsequent resurfacing, whereas the relative lack of extensive differentiation for Callisto may explain the inferred absence of endogenic tectonism.     

Tidal interactions and disruptions of giant planets on highly eccentric orbits

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.     

Tidal disruptions: II. A continuum theory for solid bodies with strength, with applications to the Solar System

We present a comprehensive theory for the breakup conditions for ellipsoidal homogeneous secondary bodies subjected to the tidal forces from a nearby larger primary: for materials ranging from purely fluid ones, to granular rubble-pile gravel-like ones, and to those with either cohesive or granular strength including cohesive rocks and metals. The theory includes but greatly extends the classical analyses given by Roche in 1847, which dealt only with fluids, and also our previous analysis [Holsapple, K.A., Michel, P., 2006. Icarus 183, 331-348], which dealt only with solid but non-cohesive bodies. The results here give the distance inside of which breakup must occur, for both a steadily orbiting satellite and for a passing or impacting object. For the fluid bodies there is a single specific shape (a “Roche Ellipsoid”) that can be in equilibrium at any given distance from a primary, and especially only one shape that can exist at the overall minimum distance (d/R)(ρ/ρp)^1/3=2.455, the classical well-known “Roche limit.” In contrast, solid bodies can exist at a given distance from a primary with a range of shapes. Here we give multiple plots of the minimum distances for various important combinations of body shape, spin, mass density, and the strength parameters characterized by an angle of friction and cohesive strength. Such results can be used in different ways. They can be used to estimate limits on strengths and mass densities for orbiting bodies at a known distance and shape. They can be used to determine breakup distances for passing bodies with an assumed strength and shape. They can be used to constraint physical properties such as bulk density of bodies with a known shape that were known to breakup at a given distance. A collection of approximately 40 satellites of the Solar System is used for comparison to the theory. About half of those bodies are closer than the Roche fluid limit and must have some cohesion and/or friction angle to exist at their present orbital distance. The required solid strength for those states is determined. Finally, we apply the theory to the break up of the SL9 comet at close approach with Jupiter. Our results make clear that the literature estimates of its bulk density depend markedly on unknown parameters such as shape, orientation and spin, and most importantly, material strength characterization.     

The origin of the Moon and the single-impact hypothesis,II

This is the second paper devoted to the numerical study of planetary collisions as a possible scenario for forming the Moon. We present a series of nine simulations of a collision between the protoearth and an impactor of various sizes. The mass ratio between the protoearth and the impactor ranged from 0.1 to 0.25. We were able to model both planets with iron cores, having modified our smoothed particle hydrodynamics code to allow the inclusion of up to 10 different material types. Two different relative velocities at infinity for the impactor were considered: ν_∞ = 0 km/sec and ν_∞ = 10 km/sec. We show that for a low-velocity collision and an impactor in the mass range 6.5 × 10²⁶ ≤ M_impactor ≤ 8.2 × 10²⁶ g, more than a lunar mass of iron-poor material is thrown into orbit. For an impactor with a mass within this range, the ejected mass that goes into orbit is for the most part divided comparably into material orbiting inside the Roche limit and into material orbiting outside the Roche limit. This material is either spread out in the form of a disk, or, for a relatively narrow range of masses toward the lower end of the range, clumped into an object of about lunar mass beyond the Roche limit. For impactors more massive than about 8.2 × 10²⁶ g we found that there is too little mass thrown into orbit. For very small mass impactors well over a lunar mass is placed in orbit, but a large amount of it is iron. In the high-velocity range we did not find a possible mass range for the impactor that would lead to the formation of an iron-poor disk massive enough to form the Moon.     

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.     

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.