ASTRONOMY • THE SOLAR SYSTEM

Planetesimals & Planetary Formation — Describe planetesimals, accretion, and differentiation and what they imply about planetary interiors.

How dust grains became worlds with layered interiors through accretion and gravitational differentiation.

Historical Context & Motivation

The question of how planets form from diffuse clouds of gas and dust is one of the oldest in natural philosophy. Long before telescopes could resolve protoplanetary disks around distant stars, thinkers speculated about the origins of the solar system's architecture. The modern scientific narrative begins with the nebular hypothesis, first articulated in the eighteenth century, which proposed that a rotating cloud of material collapsed under its own gravity to form the Sun and its attendant planets. Over the following centuries, advances in physics, meteoritics, and space exploration refined this picture into the detailed framework of planetesimal accretion and planetary differentiation that we study today.

1755
Kant's Nebular Hypothesis
Immanuel Kant proposed that the solar system condensed from a rotating cloud of gas, laying philosophical groundwork for modern planet-formation theory.
1796
Laplace's Refinement
Pierre-Simon Laplace independently developed a mathematical version of the nebular hypothesis, suggesting concentric rings of material cooled to form planets.
1969
Apollo Samples & Safronov's Model
Viktor Safronov published his landmark work on planetesimal accretion, while Apollo lunar samples provided direct geochemical evidence of differentiation in planetary bodies.
1995
First Exoplanet Around a Sun-like Star
The detection of 51 Pegasi b by Mayor and Queloz demonstrated that planet formation is a universal process, spurring comparative planetology across stellar systems.
2014–present
ALMA Disk Imaging
The Atacama Large Millimeter Array revealed ring-and-gap structures in protoplanetary disks such as HL Tauri, providing direct observational evidence of ongoing planetesimal formation.

Despite centuries of progress, a central question persists: how do micrometer-scale dust grains grow by more than thirteen orders of magnitude in mass to become full-fledged planets, and how does the internal structure of those planets reflect the thermal and compositional history of their formation? Answering this question requires understanding three interconnected processes—planetesimal formation, accretion, and differentiation—each of which leaves indelible fingerprints on the worlds we observe today.

Core Principles & Definitions

The formation of planets from a protoplanetary disk proceeds through a sequence of stages, each governed by distinct physical mechanisms. Dust grains collide and stick via electrostatic and van der Waals forces to build centimeter-scale aggregates. These aggregates must then overcome the so-called meter-size barrier—a regime where collisional fragmentation and rapid radial drift toward the central star threaten to destroy growing bodies—before reaching kilometer scales. Once bodies attain sizes of roughly one to ten kilometers, they are termed planetesimals, and their subsequent growth is dominated by gravitational interactions rather than surface forces.

1

Planetesimals

Rocky or icy bodies approximately 1–100 km in diameter that serve as the fundamental building blocks of planets. They form via coagulation or gravitational instability within protoplanetary disks and grow through mutual collisions.
2

Accretion

The process by which planetesimals collide and merge under gravitational attraction, progressively building larger bodies. Runaway and oligarchic growth phases lead to planetary embryos and eventually full-sized planets.
3

Differentiation

The internal separation of a planetary body into compositionally distinct layers—typically an iron-rich core, a silicate mantle, and a lighter crust—driven by gravitational potential energy release and radiogenic heating.
4

Protoplanetary Disk

A circumstellar disk of gas and dust orbiting a young star, typically persisting for 1–10 Myr. Temperature and composition gradients within the disk set the initial conditions for planetary chemistry.
5

Gravitational Focusing

An enhancement of the effective collision cross-section of a growing body due to its gravitational field. Larger bodies sweep up material from a region much wider than their physical size, accelerating growth.
KEY TAKEAWAY
Think of planetary formation like building a snowman during a snowstorm. Initially, you pack a small snowball by hand (surface forces dominating at grain scales). Once the ball is large enough, rolling it across the ground causes snow to stick gravitationally and by compression (gravitational accretion). Finally, if you brought the snowman indoors and it partially melted, the denser dirt and gravel would sink to the bottom while water pooled at the top—an analog for differentiation, where dense materials settle toward the center of a molten body.

Visual Explanation — From Dust to Differentiated Planet

The diagram traces four major stages of planet formation from left to right: micrometer-scale dust grains coagulate into kilometer-scale planetesimals, which accrete into protoplanets, and ultimately become differentiated planets with layered interiors (core, mantle, crust). The lower panel summarizes the dominant physical mechanism at each transition.

The diagram above encapsulates the central narrative of this lesson. Notice how each stage involves a qualitative shift in the dominant physics: surface forces yield to gravity as bodies grow, and internal thermal processes (radioactive decay, gravitational potential energy release, and impact heating) eventually drive the differentiation that produces the layered interiors we infer from seismology and moment-of-inertia measurements. The characteristic timescales span five orders of magnitude, from roughly 105 years for initial planetesimal assembly (in optimistic models invoking the streaming instability) to 108 years for the final assembly of terrestrial planets through giant impacts.

Mathematical Framework

Several key equations underpin our quantitative understanding of accretion and differentiation. These expressions connect observable properties—such as a planet's mass, radius, and moment of inertia—to the physics of growth and internal structure.

GRAVITATIONAL FOCUSING CROSS-SECTION
σ = π R² (1 + v_esc² / v_∞²)
Here σ is the effective collisional cross-section, R is the physical radius of the accreting body, vesc = √(2GM/R) is the escape velocity, and v is the velocity of incoming planetesimals at infinity. When vesc ≫ v, the factor (1 + vesc²/v²) becomes very large, enabling runaway accretion.
ACCRETION RATE
dM/dt = ρ_s σ v_∞ = π R² ρ_s v_∞ (1 + v_esc² / v_∞²)
The mass growth rate dM/dt depends on the surface mass density of planetesimals in the swarm (ρs), the effective cross-section σ, and the encounter velocity v. This expression shows that more massive bodies grow faster—the hallmark of runaway growth.
MOMENT OF INERTIA FACTOR
C / (M R²) = (2/5) for uniform sphere; < 2/5 for differentiated body
The dimensionless moment-of-inertia factor C/(MR²) equals 0.4 for a homogeneous sphere. Earth's measured value of ≈ 0.3307 confirms that mass is concentrated toward the center, providing direct evidence for a dense metallic core.
GRAVITATIONAL ENERGY OF DIFFERENTIATION
ΔE_grav = −(3 GM²) / (5 R) × [1 − (ρ̄/ρ_core)^(2/3) × f(geometry)]
During differentiation, the release of gravitational potential energy as dense iron sinks to the center provides substantial heating—on the order of 1030 J for Earth-sized bodies. This energy, combined with radiogenic heating from 26Al and 60Fe, is sufficient to melt the entire body and drive core–mantle separation.

These equations collectively tell a powerful story: gravitational focusing accelerates the growth of the largest bodies in a swarm, the accretion process deposits enormous kinetic energy as heat, and the resulting temperature rise—augmented by short-lived radioactive isotopes—drives the global melting that allows density-driven stratification. The measurable consequence—a moment-of-inertia factor less than 0.4—provides the observational test that confirms differentiation has occurred.

Differentiation & Planetary Interiors

Once a protoplanet has accreted enough mass—and generated enough internal heat—it undergoes differentiation, the wholesale reorganization of its interior into concentric layers of decreasing density from center to surface. Iron and siderophile (iron-loving) elements such as nickel, cobalt, and platinum-group metals sink to form a metallic core. Silicate minerals rich in magnesium and iron (olivine, pyroxene) constitute the mantle, and the lightest aluminosilicate minerals and incompatible elements float to form the crust. This layering is not merely academic—it governs a planet's magnetic field generation, tectonic activity, volcanic outgassing, and long-term thermal evolution.

Cross-sections of four differentiated bodies drawn to qualitatively illustrate the relative sizes of their cores, mantles, and crusts. The moment-of-inertia factor C/(MR²) listed beneath each body quantifies the degree of mass concentration toward the center: lower values indicate stronger differentiation and larger metallic cores.
Comparison of differentiated bodies in the solar system
BodyCore CompositionCore Fraction (by mass)C/(MR²)Evidence for Differentiation
EarthFe–Ni (solid inner + liquid outer)~32.5%0.3307Seismic P & S waves, magnetic dynamo, siderophile-depleted mantle
MarsFe–Ni–S (partially liquid)~20–25%0.3644InSight seismology, SNC meteorites, crustal magnetism
MoonSmall Fe core~2–4%0.3931Apollo seismometers, laser ranging, anorthosite crust
VestaFe–Ni~18%~0.33HED meteorites, Dawn gravity data, basaltic surface

Several lines of evidence converge to establish that differentiation occurred early in solar system history—within the first few tens of millions of years. Hafnium-tungsten isotope systematics (182Hf → 182W, half-life ≈ 8.9 Myr) demonstrate that Earth's core segregated within approximately 30–60 Myr of solar system formation. The asteroid 4 Vesta, with a diameter of only ~525 km, also differentiated—showing that even relatively small bodies can reach the requisite temperatures if short-lived radionuclides such as 26Al (half-life ≈ 0.72 Myr) were present in sufficient abundance at the time of accretion.

Worked Example — Interpreting a Planet's Interior

Consider a hypothetical terrestrial exoplanet for which radial velocity and transit observations have yielded a mass M = 5.0 × 1024 kg and a radius R = 6.0 × 106 m. Rotational oblateness measurements give a moment-of-inertia factor C/(MR²) = 0.315. We wish to determine whether this planet is differentiated and estimate the gravitational focusing factor for a 50-km planetesimal accreting in a swarm with v = 100 m/s.

Determining Differentiation & Gravitational Focusing
1
Step 1 — Assess Differentiation from the MOI FactorA uniform-density sphere has C/(MR²) = 2/5 = 0.400. Our planet's measured value is 0.315, which is significantly less than 0.400. This tells us that mass is concentrated toward the center, confirming the planet is strongly differentiated—even more so than Earth (0.331). We conclude a substantial iron-rich core exists.
C/(MR²) = 0.315 < 0.400 → Planet is differentiated with a large dense core.
2
Step 2 — Compute Mean DensityThe mean density is ρ̄ = M / (4/3 π R³). Substituting: ρ̄ = 5.0 × 1024 / (4/3 × π × (6.0 × 106)³) = 5.0 × 1024 / (9.047 × 1020) ≈ 5527 kg/m³. This is comparable to Earth's mean density of 5514 kg/m³, consistent with a similar bulk composition.
ρ̄ ≈ 5527 kg/m³
3
Step 3 — Calculate Escape Velocity of a 50-km PlanetesimalFor a 50-km radius body with density ≈ 3000 kg/m³ (silicate-rich), mass m = (4/3)π(5 × 10⁴)³ × 3000 ≈ 1.57 × 1018 kg. The escape velocity is vesc = √(2Gm/R) = √(2 × 6.674 × 10⁻¹¹ × 1.57 × 1018 / 5 × 10⁴) = √(4.19 × 10³) ≈ 64.7 m/s.
vesc64.7 m/s
4
Step 4 — Compute the Gravitational Focusing FactorThe focusing factor F = 1 + vesc²/v² = 1 + (64.7)²/(100)² = 1 + 4188/10000 = 1 + 0.419 ≈ 1.42. The effective cross-section is 42% larger than the geometric cross-section, providing a modest but real enhancement to the accretion rate. For larger bodies where vesc ≫ v, this factor becomes enormous—for a 500-km body with vesc ≈ 450 m/s, F ≈ 21.3.
F ≈ 1.42 (for 50-km body); demonstrates onset of gravitational focusing
5
Step 5 — Interpret Physical ImplicationsThe low moment-of-inertia factor (0.315) combined with the Earth-like mean density tells us this planet experienced thorough melting and iron–silicate separation early in its history. The gravitational focusing calculation shows that at 50-km scales, accretion is already beginning to accelerate beyond simple geometric cross-section estimates—this is the incipient stage of runaway growth that eventually builds planetary embryos.
Conclusion: The planet is strongly differentiated with a proportionally larger core than Earth, and the 50-km planetesimal is at the threshold of runaway accretion.

Strengths & Limitations of the Accretion–Differentiation Model

The planetesimal accretion and differentiation framework is remarkably successful at explaining a wide range of solar system observations, from the chemical zonation of meteorite parent bodies to the seismically determined structure of Earth's interior. However, several outstanding problems remain, particularly regarding the earliest stages of solid growth and the diversity of planetary architectures revealed by exoplanet surveys.

Strengths and limitations of the planetesimal accretion–differentiation model
StrengthLimitation / Open Question
Explains layered interiors (core–mantle–crust) via density-driven segregationThe 'meter-size barrier' remains poorly understood—how exactly do cm-sized pebbles grow to km-sized planetesimals?
Predicts siderophile element depletion in silicate mantles, confirmed by geochemistryTimescale tension: Hf–W data suggest rapid core formation (~30 Myr), but classical accretion models predict slower assembly
Gravitational focusing naturally produces runaway growth, explaining the dominance of a few large bodiesHot Jupiters and super-Earths challenge simple in-situ formation; migration must be invoked
Isotopic chronometry (²⁶Al, Hf–W) provides absolute timing of differentiation eventsPartial differentiation of some asteroids (e.g., Psyche) is difficult to model in detail
Consistent with ALMA disk observations showing ring/gap structures at expected planetesimal-forming locationsPebble accretion (an alternative/complementary mechanism) may dominate in some regimes, complicating the classical picture
KEY TAKEAWAY
The accretion–differentiation model is analogous to a successful engineering blueprint that explains the finished product (a layered planet) extremely well, but whose manufacturing instructions (the earliest assembly steps) still have some missing pages. Recent advances—particularly the streaming instability mechanism, which concentrates solids via aerodynamic coupling with disk gas—are helping to fill those gaps, but a complete, self-consistent model from dust to planet remains an active frontier of research.

Connections to Advanced Planetary Science

The principles of planetesimal accretion and differentiation extend naturally into several advanced topics that are at the forefront of current research. Understanding these connections deepens appreciation for how the basic framework scales to diverse contexts—from the formation of giant planet cores to the interpretation of exoplanetary mass–radius relationships.

From classical planetesimal theory to modern planetary science
Classical ConceptAdvanced ExtensionKey Implication
Planetesimal accretionPebble accretion — aerodynamically coupled cm-sized particles are swept up efficiently by protoplanetsExplains rapid core growth of gas giants within disk lifetimes (~3–5 Myr)
Gravitational focusingN-body simulations — full dynamical modeling of planetesimal swarms with thousands of interacting bodiesReveals stochastic outcomes: final planet counts and orbital architectures are probabilistic
Differentiation (core formation)Magma ocean dynamics — modeling iron droplet rain-out through a fully molten silicate mantleSets initial conditions for dynamo onset and early atmospheric outgassing
Moment-of-inertia factorExoplanet interior modeling — using mass–radius data with equations of state to infer core mass fractionsEnables classification of rocky vs. water-world vs. iron-enriched exoplanets
Disk temperature gradientGrand Tack / Nice Model — giant planet migration reshuffles planetesimal reservoirsExplains the small mass of Mars and the compositional structure of the asteroid belt

A particularly exciting frontier is the use of mass–radius diagrams for transiting exoplanets. When both mass (from radial velocity) and radius (from transit depth) are known, one can compute the mean density and compare it against theoretical models assuming different core mass fractions. A planet with the same mass as Earth but a higher density likely has a proportionally larger iron core—perhaps the remnant of a giant impact that stripped much of the silicate mantle, analogous to Mercury's formation scenario. Conversely, a low-density planet of the same mass might possess a thick water or ice layer, implying formation beyond the snow line followed by inward migration. In every case, the logic connects directly back to the accretion and differentiation principles developed in this lesson.

Practice Problems

PROBLEM 1CONCEPTUAL
Explain why the moment-of-inertia factor C/(MR²) for a differentiated planet is always less than 0.400. What physical process produces this departure from the uniform-sphere value, and what observable consequence does it have for a planet's rotational dynamics?
PROBLEM 2BASIC CALCULATION
A planetesimal has a radius of 25 km and a bulk density of 3500 kg/m³. Calculate its escape velocity and the gravitational focusing factor if the velocity dispersion of the surrounding swarm is v = 50 m/s. Use G = 6.674 × 10⁻¹¹ N m² kg⁻².
PROBLEM 3INTERMEDIATE
Two protoplanets in the same orbital zone have radii of 100 km and 500 km, respectively, with similar bulk densities of 3000 kg/m³. Both encounter planetesimals at v = 80 m/s. Calculate the ratio of their mass accretion rates dM/dt, accounting for gravitational focusing. What does this imply about the nature of growth in a planetesimal swarm?
PROBLEM 4APPLIED
An exoplanet has a measured mass of 8.0 × 10²⁴ kg and a radius of 7.2 × 10⁶ m. Its moment-of-inertia factor is determined to be 0.285. (a) Calculate the mean density. (b) Compare C/(MR²) with Earth's value and explain what this implies about the exoplanet's core mass fraction. (c) Propose a formation scenario that could produce such a planet.
PROBLEM 5CRITICAL THINKING
The asteroid 4 Vesta (diameter ≈ 525 km) is differentiated, while the asteroid 1 Ceres (diameter ≈ 940 km) appears only partially differentiated and retains significant water ice in its interior. Both formed in the asteroid belt, yet they have fundamentally different internal structures. Using the concepts of accretion timing, radiogenic heating, and volatile retention, construct a coherent argument explaining why a smaller body can be more thoroughly differentiated than a larger one.

Summary — Planetesimals, Accretion, and Differentiation

Planets form through a multi-stage process within protoplanetary disks. Micrometer-scale dust grains coagulate and—through mechanisms such as the streaming instability—overcome the meter-size barrier to form kilometer-scale planetesimals. These bodies grow via gravitational focusing, which amplifies their effective collision cross-section (σ = πR²[1 + vesc²/v²]), driving runaway accretion in which the largest body in each feeding zone dominates mass acquisition.

As protoplanets grow, kinetic energy from impacts combined with radiogenic heating (primarily from ²⁶Al and ⁶⁰Fe) raises internal temperatures above the iron–silicate melting point, initiating differentiation: dense iron sinks to form a metallic core while lighter silicates float to form the mantle and crust. The resulting internal layering is directly observable through the moment-of-inertia factor C/(MR²), which falls below 0.400 for any body with mass concentrated toward its center. Earth's value of 0.331, Mars's 0.364, and the Moon's 0.394 encode the degree to which each body melted and separated internally—providing a direct window into the thermal and accretionary histories written into planetary interiors billions of years ago.

Varsity Tutors • Astronomy • Planetesimals & Planetary Formation