0
What Is a Planet?
From five wandering lights to a Galaxy of a hundred billion worlds
For most of history "planet" meant five points of light that refused to stay put among the fixed stars. Today the IAU definition requires a body that orbits a star, is massive enough for gravity to pull it round, and has cleared its orbital neighborhood — the third clause famously demoting Pluto in 2006. But the deeper shift is quantitative: since the first exoplanet around a Sun-like star was confirmed in 1995, the census has grown to roughly six thousand confirmed planets, and Kepler's statistics imply that planets outnumber stars in the Milky Way. Planet formation is not a fluke; it is the default byproduct of star formation.
The discoveries broke every solar-system prejudice. Jupiters orbiting in four days, planets denser than iron and puffier than cotton candy, worlds orbiting two suns, chains of resonant super-Earths — a class of planet (1–4 R⊕) that is the most common outcome of formation yet has no representative in the Solar System at all. The physics on this sheet is what turned those surprises into a discipline:
Celestial mechanics — the oldest exact science, now stability theory for planetary systems. Disk & formation physics — how micron dust becomes 10⁴-km planets in a few Myr. Detection methods — five independent ways to find what cannot be seen. Interiors & atmospheres — hydrostatic bodies read through mass, radius, and spectra. Habitability — where liquid water, and perhaps biology, is physically permitted.
As on the companion gravity, stellar, and galaxies sheets, every equation is paired with a plain-language reading of what it physically asserts, a Use in Research column with key references, and each section ends with the open unknowns. Newton's law, Kepler's third law, and the Roche limit live on the gravity sheet; stellar masses and radii — the reference frame every exoplanet measurement is ratioed against — live on the stellar-astrophysics sheet. Toggle the Dark theme at top-right for a dark background.
I
Orbital Mechanics & the Two-Body Problem
6 equationsThe two-body problem is the only exactly solved problem in gravitational dynamics, and planetary science is built on the solution. These six lines carry spacecraft across the Solar System, weigh exoplanets, and define the sphere of influence within which a planet owns its moons. Kepler's third law itself — the mass-measurement engine — lives on the companion gravity sheet.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| The Conic-Section Orbit | \[ r(\theta) = \frac{a(1-e^2)}{1+e\cos\theta} \]
Every bound two-body orbit is an ellipse; every unbound one a parabola or hyperbola. One geometric family, indexed by a single number e, exhausts all possible trajectories — the miracle Newton extracted from the inverse-square law. |
a = semi-major axis; e = eccentricity; θ = true anomaly (angle from perihelion) |
The coordinate skeleton of every ephemeris, transit model, and RV fit; perihelion q = a(1−e) and aphelion Q = a(1+e) follow immediately.
Key referencesKepler (1609); Newton (1687); Murray & Dermott (1999).
|
| Vis-Viva Equation | \[ v^2 = GM\left(\frac{2}{r}-\frac{1}{a}\right) \]
Speed anywhere on an orbit from just the current distance and the orbit size — energy conservation in its most useful disguise. Fall inward and speed up; the semi-major axis alone fixes the total energy. |
r = current distance; a = semi-major axis; M = central mass |
The working equation of mission design: every ΔV budget, gravity assist, and orbit-insertion burn is a vis-viva calculation before it is anything else.
Key referencesNewton (1687); Bate, Mueller & White (1971).
|
| Kepler's Equation | \[ M = E - e\sin E \]
The link between clock time (mean anomaly M) and position on the ellipse (eccentric anomaly E). It has no closed-form solution — the first famous transcendental equation in physics, solved by iteration for four centuries. |
M = 2πt/P, mean anomaly; E = eccentric anomaly; e = eccentricity |
Solved (typically by Newton–Raphson) millions of times inside every RV-fitting code, transit ephemeris, and N-body integrator initializer.
Key referencesKepler (1609); Colwell (1993); Fulton et al. (2018, RadVel).
|
| Hill Sphere | \[ r_{\rm H} = a\left(\frac{m}{3M}\right)^{1/3} \]
The radius within which a planet's gravity beats the star's tide — the planet's sphere of ownership. Moons live inside it; anything outside is on loan to the star. |
m = planet mass; M = stellar mass; a = orbital distance |
Sets where satellites are stable, the feeding zone of a growing planet, and the truncation radius of circumplanetary disks; stable prograde orbits fill roughly the inner half of r_H.
Key referencesHill (1878); Hamilton & Burns (1992).
|
| Restricted Three-Body Problem & Lagrange Points | \[ C_J = x^2+y^2 + \frac{2(1-\mu)}{r_1} + \frac{2\mu}{r_2} - v^2 \]
Add a third, massless body and orbits stop closing — but one quantity, the Jacobi constant, survives. Its level surfaces carve space into allowed and forbidden zones, with five equilibrium points where a test particle co-rotates with the two masses forever. |
μ = m₂/(m₁+m₂); r₁, r₂ = distances to the two masses; rotating-frame coordinates |
The framework of space-mission staging (halo orbits at L1/L2), Trojan asteroids (L4/L5), and the zero-velocity curves that decide whether a moon can escape or an impactor be captured.
Key referencesLagrange (1772); Jacobi (1836); Szebehely (1967).
|
| Hohmann Transfer | \[ \Delta v_1 = \sqrt{\frac{GM}{r_1}}\left(\sqrt{\frac{2r_2}{r_1+r_2}}-1\right) \]
The minimum-energy path between two circular orbits is half an ellipse touching both — burn to leave, coast, burn to stay. Patience is the cheapest propellant. |
r₁, r₂ = inner/outer orbit radii; second burn Δv₂ at apoapsis completes the transfer |
The baseline against which every interplanetary trajectory is costed; combined with the synodic period 1/P_syn = 1/P₁ − 1/P₂ it sets launch windows.
Key referencesHohmann (1925); Bate, Mueller & White (1971).
|
Open unknowns · Orbital Mechanics
Is the Solar System Stable?
Will the planets stay on their orbits for the Sun's remaining lifetime?
The N-body problem has no general solution, and numerical integrations show the inner Solar System is chaotic with a ~5-Myr Lyapunov time. There is a ~1% chance Mercury's orbit destabilizes before the Sun dies — stability is probabilistic, not guaranteed.
Planet Nine
Does an unseen planet shepherd the distant Kuiper Belt?
The clustering of extreme trans-Neptunian orbits suggests a ~5–10 M⊕ planet at hundreds of AU — pure celestial-mechanics inference, exactly how Neptune was found. A decade of searching (and Vera Rubin Observatory's survey) has yet to settle it.
Free-Floating Planets
How many planets orbit no star at all?
Microlensing surveys suggest unbound or very-wide planets may rival or outnumber stars — ejected by the same dynamical chaos these equations describe. Roman Space Telescope should measure their mass function this decade.
II
Tides, Resonances & Long-Term Dynamics
6 equationsOver millions of orbits, effects invisible on any single orbit take over: tides drain energy, resonances trade angular momentum, and secular terms slowly reshape whole systems. This is the physics that decides which architectures survive. The tidal force and Roche limit themselves are on the gravity sheet; here is what tides and resonances do over Gyr.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Mean-Motion Resonance | \[ \phi = (p+q)\,\lambda' - p\,\lambda - q\,\varpi \]
When orbital periods sit near a ratio of small integers, the resonant angle φ librates instead of circulating — the two bodies exchange angular momentum coherently, orbit after orbit. Resonance is a phase lock, and it can protect or destabilize. |
λ, λ′ = mean longitudes; ϖ = periapse longitude; p, q = integers (q = order) |
Explains Pluto's survival (3:2 with Neptune, protected from close approaches), the Kirkwood gaps, and the resonant chains that fossilize disk migration in exoplanet systems.
Key referencesPeale (1976); Murray & Dermott (1999); Luger et al. (2017).
|
| Tidal Circularization Timescale | \[ \tau_{\rm circ} \simeq \frac{4}{63}\,\frac{Q_p'}{n}\,\frac{M_p}{M_*}\left(\frac{a}{R_p}\right)^{5} \]
Tides raised on an eccentric planet flex it every orbit; internal friction (parameterized by Q′) converts orbital energy to heat until the orbit is a circle. The (a/R)⁵ makes this a switch: violent up close, negligible farther out. |
Q′_p = 3Q/2k₂, reduced tidal quality factor; n = mean motion; R_p = planet radius |
Why hot-Jupiter orbits are circular and a sharp eccentricity–period boundary appears near P ≈ 3–5 d; measured e(P) distributions are the main empirical constraint on giant-planet Q′ ~ 10⁵–10⁶.
Key referencesGoldreich & Soter (1966); Hut (1981); Jackson et al. (2008).
|
| Tidal Heating | \[ \dot E = \frac{21}{2}\,\frac{k_2}{Q}\,\frac{G M_*^2 R^5\, n\, e^2}{a^6} \]
An eccentric orbit in a strong tidal field kneads a body like bread dough, at power set by how deformable it is (k₂) and how lossy (1/Q). Orbital eccentricity becomes geology. |
k₂ = tidal Love number; Q = quality factor; e = eccentricity; n = mean motion |
Powers Io's volcanism and Europa's and Enceladus's oceans; for exoplanets it can dominate internal heat budgets and may render eccentric M-dwarf habitable-zone planets volcanic waterworlds — or magma balls.
Key referencesPeale, Cassen & Reynolds (1979); Yoder (1979); Spencer et al. (2000).
|
| Tidal Despinning (Locking) Timescale | \[ \tau_{\rm lock} \sim \frac{\omega\, a^6\, I\, Q}{3\,G M_*^2\, k_2\, R^3} \]
The same tidal friction that circularizes orbits also brakes rotation, until one face points at the star forever. The a⁶ dependence means close-in worlds lock in a cosmic instant; distant ones never do. |
ω = initial spin; I ≈ 0.33 MR², moment of inertia; Q, k₂ = tidal parameters |
Decides which exoplanets are tidally locked — including essentially all habitable-zone planets of M dwarfs, forcing permanent day/night sides and reshaping their climates (and observability).
Key referencesGoldreich & Peale (1966); Gladman et al. (1996); Barnes (2017).
|
| Resonance Overlap & Chaos | \[ \Delta a_{\rm chaos} \simeq 1.3\,\mu^{2/7} a \]
Where neighboring resonances overlap, no orbit can remember which one it belongs to — motion becomes chaotic (Chirikov's criterion). Every planet is wrapped in a chaotic zone whose width grows as a weak power of its mass. |
μ = planet/star mass ratio; a = planet's semi-major axis |
Predicts which test particles survive near a planet — the edge of the Kuiper Belt, clearing of the Kirkwood gaps, and the packing limit ("dynamical spacing") of tightly packed exoplanet systems.
Key referencesChirikov (1979); Wisdom (1980); Deck et al. (2013).
|
| Secular Dynamics & Kozai–Lidov Cycles | \[ \sqrt{1-e^2}\,\cos i \;=\; {\rm const} \]
Averaged over orbits, planets trade eccentricity and inclination while keeping this combination fixed. A distant companion inclined by more than 39.2° drives huge coherent e–i oscillations — a slingshot hiding in plain algebra. |
e, i = eccentricity and inclination relative to the outer perturber's orbit |
The leading mechanism (with scattering) for making hot Jupiters: Kozai cycles pump e until tides grab the planet at periastron and shrink the orbit. Also drives Milankovitch-type climate cycles and satellite orbital decay.
Key referencesKozai (1962); Lidov (1962); Wu & Murray (2003); Naoz (2016).
|
Open unknowns · Tides & Resonances
What Is Q, Physically?
Why can't we predict any body's tidal dissipation from first principles?
Q is measured (Earth ~12 at tidal frequencies, Jupiter ~10⁵, hot Jupiters inferred ~10⁵–10⁶) but not derived — it hides unresolved physics of internal friction, oceans, and resonant fluid modes. Jupiter's measured dissipation (from Io's orbit) is ten times some theoretical expectations.
Why So Few Resonant Systems?
Migration predicts resonant chains — why are most Kepler pairs just wide of resonance?
Kepler multis pile up slightly outside 3:2 and 2:1 rather than in them. Tidal drift, stochastic disk turbulence, and obliquity tides are all proposed to break chains after formation; which mechanism dominates is unresolved.
Hot-Jupiter Origins
Disk migration, Kozai, scattering — what actually delivers giant planets to 0.05 AU?
All three channels work on paper; misaligned and retrograde hot Jupiters favor high-eccentricity migration, resonant warm Jupiters favor disks. The population is probably a blend, and its proportions remain one of the field's founding unsolved problems (since 51 Peg b, 1995).
III
Protoplanetary Disks
6 equationsPlanets are made from the leftovers of star formation: a rotating disk of gas and dust, ~99% gas by mass, alive for only a few Myr. The disk sets every initial condition of planet formation — how much solid material exists, where ices survive, and how long the clock runs. ALMA now images these disks directly, rings, gaps and all.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Minimum-Mass Solar Nebula | \[ \Sigma(r) \simeq 1700\left(\frac{r}{\rm AU}\right)^{-3/2}\ {\rm g\,cm^{-2}} \]
Smear the present planets back into a smooth disk of solar-composition gas, and this is the least surface density that could have built them. Not a theory — an accounting exercise that became the reference disk of the field. |
Σ = gas surface density; solids ~1% of this, boosted ~3–4× beyond the snow line |
The default initial condition of formation models for forty years; total implied disk mass ~0.01–0.02 M☉, consistent with the median dust masses ALMA measures in 1–3 Myr disks.
Key referencesWeidenschilling (1977); Hayashi (1981); Chiang & Laughlin (2013).
|
| Disk Temperature Profile | \[ T(r) \simeq 280\left(\frac{r}{\rm AU}\right)^{-1/2}\left(\frac{L_*}{L_\odot}\right)^{1/4}\,{\rm K} \]
A dust grain in the disk equilibrates where stellar heating balances its own radiation, so temperature falls as the inverse square root of distance. The disk is a thermometer laid across the future planetary system. |
L_* = stellar luminosity; flared, optically thick disks run somewhat shallower (T ∝ r^{−3/7}) |
Sets which solids can exist where — the condensation sequence from refractory calcium–aluminum minerals inward to ices outward — and hence the composition of everything that later accretes.
Key referencesHayashi (1981); Chiang & Goldreich (1997); Lecar et al. (2006).
|
| Snow Line | \[ r_{\rm snow} \simeq 2.7\left(\frac{L_*}{L_\odot}\right)^{1/2}\,{\rm AU} \]
The orbit beyond which water freezes onto grains. Crossing it multiplies the solid mass available for planet-building severalfold — the single most important boundary in a forming system. |
L_* = stellar luminosity; moves inward as the disk cools and accretion heating fades |
The classic explanation for why giant planets live in the outer Solar System: only beyond the snow line could cores grow big enough, fast enough, to capture gas before the disk died.
Key referencesHayashi (1981); Kennedy & Kenyon (2008); Agol et al. (2021).
|
| α-Viscosity & Disk Accretion | \[ \nu = \alpha c_s H, \qquad \dot M = 3\pi\nu\Sigma \]
Whatever turbulence transports angular momentum in the disk, hide it in one dimensionless number α times the largest eddy speed (sound speed) and size (scale height). Ignorance, efficiently parameterized — and it works disturbingly well. |
c_s = sound speed; H = c_s/Ω, scale height; α ~ 10⁻⁴–10⁻² observed |
The engine of disk evolution models: sets accretion rates onto the star, disk lifetimes, gap-opening thresholds, and migration speeds. Measuring α is a central goal of ALMA turbulence studies (so far: smaller than theorists hoped).
Key referencesShakura & Sunyaev (1973); Lynden-Bell & Pringle (1974); Hartmann et al. (1998).
|
| Toomre Stability Criterion | \[ Q = \frac{c_s\,\kappa}{\pi G \Sigma} > 1 \]
A rotating disk fights self-gravity with pressure (c_s) and shear (κ). Where Q drops below unity, the disk fragments under its own weight — the express lane from gas directly to giant planets or brown dwarfs. |
c_s = sound speed; κ ≈ Ω, epicyclic frequency; Σ = surface density |
Decides whether gravitational instability can form planets: only cold, massive, young disks at wide separations qualify. The same Q governs spiral structure in galaxy disks — one criterion, twenty orders of magnitude apart.
Key referencesToomre (1964); Boss (1997); Kratter & Lodato (2016).
|
| Disk Dispersal by Photoevaporation | \[ \dot M_{\rm wind} \sim 10^{-10}\!-\!10^{-9}\ M_\odot\,{\rm yr^{-1}} \quad (r > r_g \simeq \frac{GM_*}{c_s^2}) \]
High-energy photons heat the disk surface until gas beyond the gravitational radius r_g simply evaporates. When the wind outruns viscous resupply, the disk dies from the inside out in ~10⁵ yr — a two-timescale death: long life, sudden end. |
r_g = radius where sound speed beats escape; driven by stellar EUV/X-ray flux |
Explains the observed disk half-life of ~2–3 Myr and the rarity of "transition" disks (caught mid-clearing) — and sets the deadline every gas-giant core must beat.
Key referencesClarke et al. (2001); Haisch et al. (2001); Alexander et al. (2014).
|
Open unknowns · Protoplanetary Disks
What Drives Disk Accretion?
Turbulent viscosity or magnetized winds?
The classic driver (magnetorotational instability) is quenched in the poorly ionized planet-forming zone ("dead zone"), and ALMA measures turbulence weaker than α-models assumed. Magnetothermal disk winds may carry the angular momentum instead — a live paradigm shift.
The Rings Everywhere
Are ALMA's ubiquitous rings and gaps all carved by planets?
If so, Neptune-mass planets at tens of AU are common at age ~1 Myr — demanding planet formation far faster and farther out than any standard model. Alternatives (snow lines, MHD zonal flows) remain viable; the answer rewrites formation timelines.
The Mass Budget Problem
Measured disk dust masses seem too small to build the observed exoplanets.
Median Class II disk solids fall short of the median planetary system's solid mass. Either dust masses are underestimated (optically thick emission), or planet formation is already largely done by ~1 Myr — pushing the action into the embedded protostellar phase.
IV
From Dust to Planets
7 equationsThirteen orders of magnitude in size separate an interstellar grain from a giant-planet core, and no single process spans them. Sticking builds pebbles; a still-debated instability builds km-scale planetesimals; gravity does the rest — if it can beat the disk's ~3-Myr clock. These are the growth laws that must run to completion.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Radial Drift (the Meter-Size Barrier) | \[ v_{\rm drift,max} \simeq \eta\, v_{\rm K} \sim 50\ {\rm m\,s^{-1}} \]
Pressure support makes the gas orbit slightly slower than Keplerian (by ηv_K); solids feel no pressure, so they orbit faster and plough into a headwind. Meter-scale boulders drift toward the star in ~100 years — growth must sprint through this size or die. |
η ≈ ½|d ln P/d ln r|(c_s/v_K)² ~ 2×10⁻³, pressure-support parameter; v_K = Kepler speed |
The central obstacle of planetesimal formation; conversely, drift concentrates pebbles at pressure bumps — which is why rings in disks may be planet factories, and why pebble fluxes power outer-planet growth.
Key referencesWhipple (1972); Weidenschilling (1977); Birnstiel et al. (2012).
|
| Streaming Instability | \[ Z = \frac{\Sigma_{\rm solid}}{\Sigma_{\rm gas}} \gtrsim Z_{\rm crit}(\tau_s) \sim 0.02 \]
The headwind problem contains its own cure: drifting pebbles drag the gas, which reduces the headwind, which lets pebbles pile up — a traffic-jam instability. Where the local solid fraction crosses a threshold, clumps collapse gravitationally straight into ~100-km planetesimals. |
Z = local solids-to-gas ratio; τ_s = dimensionless stopping time (pebble aerodynamic size) |
The leading solution to the meter-size barrier; predicts planetesimals born big (~100 km) and in binaries — both confirmed in the Kuiper Belt (size distribution turnover; Arrokoth's contact-binary shape).
Key referencesYoudin & Goodman (2005); Johansen et al. (2007); McKinnon et al. (2020).
|
| Gravitationally Focused Growth | \[ \frac{dM}{dt} = \pi R^2 \rho\, v_{\rm rel}\left(1+\frac{v_{\rm esc}^2}{v_{\rm rel}^2}\right) \]
A big body's gravity bends incoming trajectories onto it, multiplying its geometric cross-section by the focusing factor in parentheses. Once v_esc exceeds the swarm's random speeds, the biggest body eats fastest — runaway growth, the rich getting richer. |
ρ = planetesimal space density; v_rel = velocity dispersion; v_esc = escaper's escape speed |
The classical (Safronov) growth mode: runaway, then oligarchic growth of embryos. Sets the ~10⁵–10⁶ yr emergence of Moon-to-Mars-mass oligarchs in the inner disk, followed by ~10⁷–10⁸ yr of giant impacts to finish Earth.
Key referencesSafronov (1969); Wetherill & Stewart (1989); Kokubo & Ida (1998); Kleine et al. (2009).
|
| Isolation Mass | \[ M_{\rm iso} \simeq \frac{(4\pi b\, a^2 \Sigma_{\rm s})^{3/2}}{(3M_*)^{1/2}} \]
An embryo can only eat what its gravitational reach (a few Hill radii, b ≈ 10 r_H/a) contains. When the feeding zone is empty, growth stalls at the isolation mass — local solids decide the largest body that neighborhood can grow. |
Σ_s = solid surface density; b ≈ 10, feeding-zone width in Hill units |
Explains why the inner Solar System made Mars-mass embryos needing later mergers, while beyond the snow line isolation masses reach the ~10 M⊕ needed to trigger gas capture directly.
Key referencesLissauer (1987); Kokubo & Ida (2002); Morbidelli et al. (2012).
|
| Pebble Accretion | \[ \dot M_{\rm peb} \simeq 2\, r_{\rm H}\, v_{\rm H}\, \Sigma_{\rm peb}\,\tau_s^{2/3} \]
For cm-sized pebbles, gas drag kills the relative velocity during an encounter, so a core captures everything drifting through a large fraction of its Hill sphere — a cross-section vastly beyond gravitational focusing. Drag turns the drift problem into a delivery service. |
r_H, v_H = Hill radius and Hill speed Ωr_H; Σ_peb = pebble surface density; τ_s = stopping time |
Grows 10-M⊕ cores at 5–30 AU within disk lifetimes — the fix for core accretion's timescale crisis, and the current default theory for giant-planet cores, ice giants, and possibly super-Earths.
Key referencesOrmel & Klahr (2010); Lambrechts & Johansen (2012); Johansen & Lambrechts (2017).
|
| Critical Core Mass & Runaway Gas Accretion | \[ M_{\rm core} \gtrsim M_{\rm crit} \approx 10\, M_\oplus \;\Rightarrow\; \dot M_{\rm gas} \to {\rm runaway} \]
A core's gas envelope can sit in hydrostatic balance only while the core dominates the mass. Past ~10 M⊕ no equilibrium exists: the envelope's own weight pulls in gas ever faster, and a rocky core becomes a gas giant in ~10⁵ yr. |
M_crit depends on envelope opacity and accretion luminosity (~5–20 M⊕) |
The defining threshold of core accretion: it splits the planet population into those that crossed it (gas giants) and those that missed the deadline (super-Earths, ice giants) — a bimodality visible across 6,000 exoplanets.
Key referencesMizuno (1980); Pollack et al. (1996); Wahl et al. (2017).
|
| Type I Migration & Gap Opening | \[ \tau_{\rm I} \simeq \frac{1}{2.7+1.1\beta}\,\frac{M_*}{M_p}\,\frac{M_*}{\Sigma a^2}\left(\frac{H}{r}\right)^{2}\Omega^{-1} \]
A planet's spiral wakes in the gas torque it — almost always inward, faster for heavier planets. An Earth at 1 AU migrates in ~10⁵ yr, far shorter than the disk's life: planets do not, in general, stay where they form. |
β = −d ln Σ/d ln r; H/r = disk aspect ratio; gap opening (→ slower Type II) once M_p ≳ 10⁻⁴ M_* for typical α, H/r |
The reason system architecture is a dynamical outcome, not a birth record: resonant chains, hot super-Earths, and warm Neptunes are all migration fossils. Stopping migration (at disk inner edges, pressure bumps) is half of modern formation theory.
Key referencesGoldreich & Tremaine (1980); Tanaka et al. (2002); Crida et al. (2006); Kley & Nelson (2012).
|
Open unknowns · From Dust to Planets
How Do Planetesimals Actually Form?
Is streaming instability the whole story — and can it start from realistic dust?
The instability needs pre-enriched solids and pebbles of the right size; whether real disks reach those thresholds everywhere planetesimals are needed (e.g. at 0.4 AU for Mercury) is unproven. Forty years after the meter-size barrier was named, the km-scale step is still the least constrained in all of planet formation.
Why Isn't Every System Solar?
The most common planets (close-in super-Earths) are the ones we lack.
Roughly half of Sun-like stars host 1–4 R⊕ planets inside Mercury's orbit; we have none. Whether the Solar System's inner emptiness reflects Jupiter's interference ("Grand Tack"), early drift-starvation, or luck is a live question with anthropic stakes.
Where Did Earth's Water Come From?
Local accretion, snow-line pebbles, or late cometary delivery?
Earth formed inside the snow line, yet is 0.02% water by mass with D/H matching carbonaceous chondrites, not most comets. The mix of dry local material and hydrated outer-belt material — and when it arrived — is still being read out of isotopes.
Did the Giant Planets Wander?
How much did Jupiter and Saturn move — and did they reshape everything?
The Grand Tack (Jupiter in to 1.5 AU, then back out) explains Mars's small mass; the Nice model's later instability explains the Kuiper Belt and Trojans. Both are inference from dynamical scars; neither is directly testable, and their timing is actively contested.
V
Detection I: Radial Velocities & Transits
7 equationsNearly every known exoplanet was found without being seen — inferred from what it does to its star. Radial velocities weigh planets; transits size them; together they give density, and hence composition, for worlds hundreds of parsecs away. These seven equations are the workhorses of the two methods that built the census.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Radial-Velocity Semi-Amplitude | \[ K = \left(\frac{2\pi G}{P}\right)^{1/3}\frac{M_p \sin i}{(M_*+M_p)^{2/3}}\frac{1}{\sqrt{1-e^2}} \]
The star orbits the system barycenter too, and its line-of-sight wobble is Doppler-shifted starlight. The planet's mass — times the unknowable sin i — is written in the star's spectrum, at meters per second. |
K = stellar velocity semi-amplitude; i = orbital inclination; P, e = period, eccentricity |
The method of the first discoveries and still the mass-measurement engine; modern spectrographs (ESPRESSO, EXPRES, NEID) reach ~0.2–0.5 m/s, where stellar surface jitter — not photons — is the wall.
Key referencesMayor & Queloz (1995); Anglada-Escudé et al. (2016); Wright (2018).
|
| Transit Depth | \[ \delta = \left(\frac{R_p}{R_*}\right)^{2} \]
A planet crossing its star blocks light in proportion to the area ratio. The simplest equation in exoplanet science — and the one that made planet-finding an industry, because photometry scales to 100,000 stars at once. |
R_p, R_* = planet and stellar radii; limb darkening shapes the ingress/egress |
Gives the planet's radius (given the star's — every transit radius is really a stellar-astrophysics measurement); with an RV mass, bulk density and composition follow.
Key referencesCharbonneau et al. (2000); Henry et al. (2000); Winn (2010).
|
| Transit Probability | \[ p_{\rm tr} \simeq \frac{R_*}{a} \]
Alignment is luck: the chance a randomly oriented orbit carries the planet across the stellar disk is just the star's angular size as seen from the planet. Close-in planets are found because geometry loves them. |
a = semi-major axis; (1+e cos…) corrections for eccentric orbits |
The completeness correction at the heart of every occurrence-rate calculation — each detected transiting planet stands in for ~a/R_* undetected twins with unlucky inclinations.
Key referencesBorucki & Summers (1984); Winn (2010); Bryson et al. (2021).
|
| Transit Duration | \[ T \simeq \frac{P}{\pi}\arcsin\!\left(\frac{R_*}{a}\sqrt{1-b^2}\right) \]
How long the crossing lasts encodes the orbit's size relative to the star and where the chord cuts (impact parameter b). Duration is geometry's second, independent message in the light curve. |
b = (a/R_*)cos i, impact parameter; P = period |
Combined with depth and period, duration yields the stellar density directly from photometry (Seager–Mallén-Ornelas) — a transit is a stellar-density meter, used to vet false positives and to spot eccentric orbits ("photoeccentric" effect).
Key referencesSeager & Mallén-Ornelas (2003); Winn (2010); Dawson & Johnson (2012).
|
| Transit-Timing Variations | \[ \delta t \;\propto\; \frac{M_{\rm pert}}{M_*}\,\frac{P}{\Delta} \]
In a multi-planet system transits don't tick like a clock — mutual pulls advance and delay them, hugely amplified near resonance (small Δ). The deviations weigh planets with no spectrograph at all. |
M_pert = perturber mass; Δ = fractional distance from exact resonance |
Mass measurements for systems too faint for RVs — the method that weighed the TRAPPIST-1 planets to a few percent and revealed a population of surprisingly low-density sub-Neptunes.
Key referencesAgol et al. (2005); Holman & Murray (2005); Agol et al. (2021).
|
| Rossiter–McLaughlin Effect | \[ \Delta v_{\rm RM} \simeq \delta\; v\sin i_* \]
A transiting planet alternately hides the approaching and receding limbs of its rotating star, imprinting a red-then-blue (or blue-then-red) squiggle in the RVs. The squiggle's shape reveals whether the planet orbits with, against, or over the star's spin. |
δ = transit depth; v sin i_* = star's projected rotation speed; shape → sky-projected obliquity λ |
The obliquity meter: measures spin–orbit alignment, the key dynamical fossil separating quiet disk migration (aligned) from high-eccentricity migration (misaligned, retrograde, polar).
Key referencesRossiter (1924); McLaughlin (1924); Winn et al. (2005); Albrecht et al. (2022).
|
| Stellar Jitter Floor | \[ \sigma_{\rm RV} \sim 1\ {\rm m\,s^{-1}}\ \; ({\rm granulation + activity}) \]
The star is not a rigid reference: convective granulation, oscillations, and spots move its spectral lines by ~m/s on every timescale from minutes to magnetic cycles. Below this floor, the noise is astrophysics, not instrumentation. |
Granulation ~0.4–1 m/s; p-modes ~0.5 m/s (Sun); spots up to many m/s on active stars |
The obstacle between current instruments (~0.3 m/s) and the 0.09 m/s Earth–Sun signal; beating it — via line-by-line diagnostics, ML activity models, and long baselines — is the defining challenge of extreme-precision RV.
Key referencesSaar & Donahue (1997); Dumusque et al. (2015); Crass et al. (2021, EPRV report).
|
Open unknowns · RV & Transits
Can Jitter Be Beaten?
Is 10 cm/s precision — an Earth twin — physically reachable?
No fundamental law forbids modeling granulation and activity to the cm/s level, but no demonstrated technique gets there yet. The answer decides whether RV mass measurements exist for the planets future flagship imagers will target.
η⊕
How common are Earth-size planets in habitable zones of Sun-like stars?
Kepler ran out of photons exactly at Earth-analog parameter space; extrapolations span η⊕ ~ 0.1–0.5. The number sets the required aperture of any mission designed to image and characterize an exo-Earth — a multi-billion-dollar uncertainty.
The Obliquity Record
Are misaligned orbits primordial, dynamical, or stellar in origin?
Stars themselves may be born tilted relative to their disks, and hot stars retain misalignment while cool stars tidally realign — so obliquity may say as much about stars as about migration. Disentangling the channels is ongoing.
VI
Detection II: Imaging, Microlensing & Astrometry
5 equationsTransits and RVs are blind to wide orbits and unlucky geometries. The other three methods fill in the census: direct imaging sees young giants far out, microlensing weighs planets (even starless ones) at kiloparsecs, and astrometry — with Gaia's data releases — is now delivering the wide-orbit giants around nearby stars. Plus the method that started it all: pulsar timing.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Planet–Star Contrast | \[ \frac{F_p}{F_*}\bigg|_{\rm refl} = A_g\left(\frac{R_p}{a}\right)^{2}\Phi(\alpha) \]
In reflected light a planet is its star diluted by albedo times the tiny solid angle it intercepts — parts per billion for an Earth twin. Direct imaging is the art of suppressing a star by ten orders of magnitude one wavefront correction at a time. |
A_g = geometric albedo; Φ(α) = phase function; young self-luminous giants instead have F_p/F_* ~ 10⁻⁴–10⁻⁶ in the IR |
Sets the coronagraph + adaptive-optics requirements of every imaging instrument; current facilities reach ~10⁻⁶ (young Jupiters in the IR), while an exo-Earth at 10⁻¹⁰ defines the Habitable Worlds Observatory's design contrast.
Key referencesTraub & Oppenheimer (2010); Marois et al. (2008, 2010); The LUVOIR/HabEx reports (2019).
|
| Einstein Radius (Microlensing) | \[ \theta_{\rm E} = \sqrt{\frac{4GM}{c^2}\,\frac{D_S-D_L}{D_S D_L}} \]
A foreground mass bends background starlight into (unresolved) images separated by ~θ_E, brightening the source as they align. A planet near the lens star adds an hours-to-days spike on the weeks-long stellar event: gravity betrays mass, with no light from the lens required. |
M = lens mass; D_L, D_S = lens/source distances; planetary anomalies scale as √(M_p/M_*)·t_E |
The only method sensitive to cold, low-mass planets at 1–10 AU around distant (and unbound) hosts — the demographics of the Galaxy's typical planets, not just the bright nearby ones. Roman's bulge survey will find thousands.
Key referencesEinstein (1936); Mao & Paczyński (1991); Gould & Loeb (1992); Beaulieu et al. (2006).
|
| Astrometric Signature | \[ \alpha = \frac{M_p}{M_*}\,\frac{a}{d} \]
The star's barycentric wobble, seen as a tiny ellipse on the sky. Unlike RV it grows with orbital distance and yields the full 3-D orbit — inclination included, so masses come out true, not "×sin i". |
a = planet's semi-major axis; d = distance; α in radians (→ µas) |
Gaia's epoch astrometry (DR4, 2026) is delivering the long-promised harvest: thousands of true-mass giant planets at 1–5 AU around nearby stars — the Jupiter-analog census, and the inclinations that break RV degeneracies.
Key referencesPerryman et al. (2014); Gaia Collaboration; Sozzetti et al. (2023).
|
| Pulsar Timing | \[ \Delta t = \frac{M_p}{M_*}\,\frac{a}{c}\,\sin i \]
A millisecond pulsar is a clock stable to parts in 10¹⁵; an orbiting planet makes the clock's ticks arrive early and late as the pulsar circles the barycenter. Light-travel time across a wobble of kilometers is measurable. |
Δt = timing residual amplitude; a = pulsar's barycentric orbit; c = light speed |
Historic more than demographic — but the precision (sub-lunar-mass sensitivity) remains unmatched, and the same residual analysis now underpins pulsar-timing-array gravitational-wave astronomy.
Key referencesWolszczan & Frail (1992); Wolszczan (1994).
|
| Secondary Eclipse & Phase Curves | \[ \frac{\Delta F}{F} \simeq \left(\frac{R_p}{R_*}\right)^{2}\frac{B_\lambda(T_p)}{B_\lambda(T_*)} \]
When the planet passes behind the star, its own light vanishes from the sum — measure the dip and you have measured the planet's thermal emission. Around the orbit, the waxing and waning dayside traces the planet's weather map in longitude. |
B_λ = Planck function; T_p = dayside brightness temperature |
The gateway to exoplanet meteorology: dayside temperatures, hotspot offsets (winds), day–night contrasts (heat redistribution), and — via eclipse mapping — the first crude 2-D images of any exoplanet.
Key referencesDeming et al. (2005); Charbonneau et al. (2005); Knutson et al. (2007); Bell et al. (2024).
|
Open unknowns · Imaging, Lensing & Astrometry
10⁻¹⁰ Contrast
Can a coronagraph + wavefront control actually hold ten orders of suppression?
Roman's coronagraph (flying now as a technology demo) targets ~10⁻⁸–10⁻⁹; an exo-Earth needs another order of magnitude, held stable for days, plus exozodiacal dust cooperation. This is the make-or-break technology of the Habitable Worlds Observatory.
The Free-Floater Mass Function
What is the mass spectrum of unbound planets — and what ejected them?
Short-timescale microlensing events suggest free-floating Earth-to-Neptune-mass objects are abundant, but single events give degenerate masses. Roman + Euclid parallax observations should break the degeneracy and read out the Galaxy's ejection history.
Joining the Censuses
Do transit, RV, lensing, imaging, and astrometry populations knit into one distribution?
Each method samples a different (mass, separation, age, host) slice with different biases. A single occurrence surface across 0.01–100 AU is the field's demographic grail — and current stitched estimates still disagree at factor-of-few level in the overlap regions.
VII
Planetary Interiors
6 equationsA planet's inside is forever invisible; everything we know arrives through mass, radius, gravity harmonics, tides, seismology, and magnetic fields. The same hydrostatic scaffolding that holds up stars (see the stellar sheet) holds up planets — with pressure supplied by rock, ice, and degenerate hydrogen instead of fusion-fed gas.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Hydrostatic Equilibrium & Central Pressure | \[ \frac{dP}{dr} = -\frac{G m(r)\,\rho}{r^2} \;\Rightarrow\; P_c = \frac{3GM^2}{8\pi R^4}\,({\rm uniform}\ \rho) \]
Every layer's weight is carried by the pressure gradient beneath it. Integrated crudely, a planet's central pressure follows from mass and radius alone — no seismometer required for the first estimate. |
m(r) = enclosed mass; P_c = central pressure; real profiles need an equation of state P(ρ, T) |
The first of the four structure equations every interior model integrates; paired with a rock/ice/H–He equation of state it turns one (M, R) point into a constrained set of internal layerings.
Key referencesDziewonski & Anderson (1981, PREM); Guillot (2005); Militzer et al. (2016).
|
| Mass–Radius Relations | \[ R \propto M^{0.27}\ ({\rm rocky}); \qquad R \approx {\rm const} \approx R_{\rm J}\ ({\rm giants}) \]
Rock compresses gently, so terrestrial radii creep up with mass; in giants, electron degeneracy fights gravity to a draw and radius goes flat from Saturn's mass to the brown-dwarf line. A planet's location in the M–R plane is a composition X-ray. |
Rocky slope ~0.27 for Earth-like mix; water/H–He layers lift R at fixed M; degeneracy turnover near ~4 M_J |
The decoder ring for the density census: given M and R, the space of allowed iron/rock/water/gas mixes follows — degenerate, but powerfully exclusionary (a 1.8 R⊕, 2 M⊕ planet cannot be bare rock).
Key referencesZapolsky & Salpeter (1969); Fortney et al. (2007); Zeng et al. (2019); Thorngren & Fortney (2018).
|
| Moment of Inertia & Darwin–Radau | \[ \bar C = \frac{C}{MR^2}: \; 0.4\ {\rm uniform};\; {\rm less\ if\ centrally\ condensed} \]
How hard a planet is to spin up depends on where its mass sits. Measured from precession or (via Darwin–Radau) from rotational flattening, one dimensionless number reveals whether a world has a core — no drilling required. |
C = polar moment of inertia; Darwin–Radau links C̄ to J₂ and the flattening for a body in hydrostatic balance |
The classical constraint on differentiation: combined with gravity harmonics (J₂, J₄…) from orbiting spacecraft it fixes core sizes throughout the Solar System — and is the template for what exoplanet Love numbers (below) may someday do.
Key referencesRadau (1885); Darwin (1899); Hubbard (1984); Iess et al. (2010).
|
| Tidal Love Number k₂ | \[ \Phi_{\rm induced} = k_2\,\Phi_{\rm tide}\quad (k_2 \le 3/2) \]
How much a body's own gravity field responds to being tidally squeezed. A fluid, uniform body answers with k₂ = 1.5; a rigid or centrally condensed one barely whispers. The tide is a free seismology experiment run by the companion. |
Φ_tide = perturbing tidal potential; measured via orbiting spacecraft or transit-shape/apsidal precession for exoplanets |
Distinguishes interior models with identical (M, R): Juno's k₂ for Jupiter, Cassini's for Saturn (whose rings double as a seismograph), and first exoplanet k₂ estimates from apsidal precession of ultra-hot Jupiters.
Key referencesLove (1911); Durante et al. (2020); Yee et al. (2020).
|
| Mantle Convection (Rayleigh Number) | \[ Ra = \frac{\rho g \alpha \Delta T\, d^3}{\kappa \mu} \gtrsim 10^3 \Rightarrow {\rm convection} \]
Solid rock flows — on million-year timescales — if the buoyancy stored in a temperature contrast beats diffusion and viscosity. Ra measures the contest; planetary mantles win it by four orders of magnitude, so geology is convection's surface expression. |
d = layer depth; μ = dynamic viscosity (~10²¹ Pa·s, Earth's mantle); κ = thermal diffusivity |
Governs how planets lose heat, whether they sustain plate tectonics or stagnant lids, and — through volcanic outgassing and the carbon cycle — couples interior physics to atmospheric habitability.
Key referencesRayleigh (1916); Schubert, Turcotte & Olson (2001); Valencia et al. (2007).
|
| Dynamo Criterion (Magnetic Reynolds & Elsasser) | \[ Re_m = \frac{u L}{\eta} \gtrsim 10\!-\!100; \qquad \Lambda = \frac{\sigma B^2}{\rho \Omega} \sim 1 \]
A conducting fluid stirred fast enough regenerates its own magnetic field; the field then grows until magnetic forces rival Coriolis (Λ ~ 1). Planetary magnetism is convection plus rotation plus conductivity — lose any one and the field dies. |
u, L = flow speed/scale; η = magnetic diffusivity; σ = conductivity; Ω = rotation rate |
Explains which worlds hold magnetospheres (Earth, Jupiter's 4-gauss giant, Ganymede) and which don't (Venus — too slow a heat loss? Mars — dynamo died 4 Gyr ago); shields atmospheres from stellar-wind stripping, with contested habitability stakes.
Key referencesElsasser (1946); Stevenson (2003); Christensen (2010); Jakosky et al. (2018).
|
Open unknowns · Interiors
Jupiter's Dilute Core
Why is the core fuzzy — smeared across half the planet's radius?
Juno's gravity field rules out a compact core; the heavy elements grade smoothly into the envelope. A giant impact, eroding convection, or primordial layered accretion could do it — and Cassini found Saturn's core is dilute too. Formation models didn't predict this.
Super-Earth Degeneracy
Is a 2.5 R⊕ planet a water world or a rocky core under hydrogen?
Radically different planets occupy identical (M, R) points. JWST atmospheres break some ties (K2-18b's spectrum launched the contested "hycean" debate), but for most sub-Neptunes the interior remains formally unknowable from bulk density alone.
Why Is Venus Not Earth?
Same size, same star — no plates, no dynamo, no water. What diverged?
Whether Venus lost its water early (runaway greenhouse) or late, and whether plate tectonics requires water in the mantle, are open. DAVINCI and EnVision fly this decade; the answer calibrates the fate of billions of Venus-zone exoplanets.
VIII
Planetary Atmospheres
7 equationsThe atmosphere is the only part of an exoplanet we can actually do spectroscopy on — a skin of gas, often a millionth of the planet's mass, forced to reveal composition, temperature, winds, and chemistry. The same seven equations run Earth's weather models and JWST's transmission-spectrum fits.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Pressure Scale Height | \[ H = \frac{k_{\rm B} T}{\mu m_{\rm H}\, g} \]
Hydrostatic air thins by e every scale height: hotter and lighter means puffier. One number sets the thickness of the atmospheric skin — and, crucially, the size of every spectral signature we can hope to see. |
T = temperature; μ = mean molecular weight; g = surface gravity |
The first number computed for any planet, Solar System or exo-; the amplitude unit of transmission spectroscopy (below) and the vertical grid of every atmosphere model.
Key referencesChamberlain & Hunten (1987); Seager (2010).
|
| Equilibrium Temperature | \[ T_{\rm eq} = T_*\sqrt{\frac{R_*}{2a}}\,\big[f(1-A_B)\big]^{1/4} \]
Absorbed starlight in, thermal radiation out; balance gives the temperature a planet would have with no greenhouse and no internal heat. Not the real surface temperature — the baseline every real atmosphere is measured against. |
A_B = Bond albedo; f = redistribution factor (1 uniform, 2 dayside-only) |
The universal first-order label of the census (every catalog lists T_eq) and the input to chemistry regimes: above ~1500 K rock vaporizes; near 350–2000 K clouds of everything from water to corundum condense.
Key referencesSagan & Mullen (1972); Seager (2010); Fortney et al. (2021).
|
| Gray Greenhouse (Two-Stream) | \[ T^4(\tau) = T_{\rm eq}^4\left(1+\tfrac{3}{4}\tau\right) \]
An atmosphere transparent to sunlight but opaque (optical depth τ) to infrared forces the surface warmer than equilibrium — radiation must random-walk out. The simplest honest model of why greenhouse gases heat planets. |
τ = thermal-IR optical depth at the surface; gray = wavelength-averaged |
The pedagogical spine of climate physics and the zeroth-order temperature–pressure profile in retrieval codes; full models replace gray τ with line-by-line radiative transfer but inherit this structure.
Key referencesSchwarzschild (1906); Chandrasekhar (1960); Pierrehumbert (2010).
|
| Adiabatic Lapse Rate | \[ \Gamma_{\rm dry} = -\frac{dT}{dz} = \frac{g}{c_p} \]
Convecting air cools as it expands on the way up, at a rate set only by gravity and heat capacity. Wherever radiation would demand a steeper gradient, convection takes over and pins the profile to this slope — tropospheres are adiabats. |
c_p = specific heat; latent heat of condensing species (water, ammonia, silicates) flattens it to a moist adiabat |
Fixes the deep temperature structure of every convective atmosphere from Earth's troposphere to Jupiter's envelope to brown dwarfs; the moist version powers hurricanes and sets the water content of rising exoplanet air.
Key referencesChamberlain & Hunten (1987); Seiff et al. (1998); Pierrehumbert (2010).
|
| Jeans Escape Parameter | \[ \lambda_{\rm esc} = \frac{G M m}{k_{\rm B} T\, r_{\rm exo}} = \left(\frac{v_{\rm esc}}{v_{\rm th}}\right)^{\!2}\!\Big/ 2 \]
At the exobase, molecules on the fast tail of the Maxwell distribution simply leave. λ compares gravitational binding to thermal energy per molecule — heavy gases on big cold planets stay (λ ≫ 10); light gases on warm small ones leak away over Gyr. |
m = molecular mass; r_exo = exobase radius; escape flux ∝ e^{−λ}(1+λ) |
Explains the Solar System's atmospheric roll call — why Earth kept N₂ but lost H, why Mars is thin, why Titan (cold) keeps N₂ at Moon-like gravity, why Jupiter keeps everything — and which exoplanets can hold what.
Key referencesJeans (1925); Hunten (1973); Catling & Kasting (2017).
|
| Energy-Limited Hydrodynamic Escape | \[ \dot M \simeq \frac{\epsilon\, \pi F_{\rm XUV} R_p^3}{G M_p} \]
When stellar X-ray/UV heating is strong, the upper atmosphere stops leaking molecule-by-molecule and flows off as a bulk transonic wind — a planet with a comet's tail. Mass loss is then set by the energy budget, not the Maxwell tail. |
F_XUV = X-ray+EUV flux; ε ≈ 0.1–0.3, heating efficiency; enhancement by the low-density factor (R_XUV/R_p)³ folded into R_p³ |
The sculptor of the small-planet population: drives the radius valley (Section XI), strips hot Neptunes (the "desert"), and turns young sub-Neptunes into bare super-Earths during their star's active first Gyr.
Key referencesWatson et al. (1981); Vidal-Madjar et al. (2003); Ehrenreich et al. (2015); Owen (2019).
|
| Transmission-Spectroscopy Signal | \[ \Delta\delta \simeq \frac{2\, n_H\, H\, R_p}{R_*^2} \]
In transit, starlight filters through the planet's limb; at wavelengths where the atmosphere absorbs, the planet looks a few scale heights bigger. The wavelength-dependence of a shadow — that is the entire data stream of exoplanet atmospheric science. |
n_H ≈ 1–5, scale heights of modulation; H = scale height; signal in fractional transit depth |
The master equation of JWST exoplanet observing: it dictates target choice (big H: hot, light, low-g atmospheres; small stars) and explains every "flat spectrum" heartbreak (clouds or high μ crush H).
Key referencesSeager & Sasselov (2000); Charbonneau et al. (2002); Kreidberg et al. (2014); JWST Transiting Exoplanet ERS Team (2023).
|
Open unknowns · Atmospheres
Clouds
Can we ever see through — or credibly model — exoplanet aerosols?
Clouds and hazes mute or erase spectral features on hot Jupiters and sub-Neptunes alike, and their microphysics (nucleation, sedimentation, composition) is essentially uncalibrated at exotic conditions. They are the systematic error on every abundance JWST reports.
Do M-Dwarf Rockies Have Air?
Did the TRAPPIST-1 worlds keep atmospheres through their star's violent youth?
JWST thermal measurements show TRAPPIST-1 b and c lack thick atmospheres; the habitable-zone members (e, f) remain open. Since M dwarfs host most of the Galaxy's temperate rocky planets, "atmosphere or airless rock?" is the current billion-dollar binary.
Chemical Disequilibrium
Why do measured abundances defy equilibrium chemistry — and what does that hide?
Vertical mixing, photochemistry (JWST's SO₂ detection was a triumph here), and interior fluxes all push atmospheres off equilibrium. Retrieval models disagree on basic quantities (WASP-39b's C/O; K2-18b's CH₄/DMS controversy) — the inference layer is now the bottleneck.
IX
Climate & Habitability
5 equations"Habitable" is a physical claim: surface liquid water, maintained over geological time. That requires the right stellar flux, a working greenhouse, a thermostat against runaway in both directions, and geological luck. These equations define the habitable zone — and remind us how conditional it is.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Habitable-Zone Boundaries | \[ d_{\rm HZ} = \left(\frac{L_*/L_\odot}{S_{\rm eff}}\right)^{1/2}\,{\rm AU} \]
Scale the Sun–Earth flux by climate-model limits: inside the inner edge water is lost to a runaway or moist greenhouse; outside the outer edge, even a maximum CO₂ greenhouse cannot keep water liquid. A belt, not a promise. |
S_eff ≈ 1.02–1.78 (inner, model-dependent) to ~0.32–0.36 (outer, maximum greenhouse) |
The target-selection function for every "potentially habitable planet" claim and for mission design (η⊕ is defined over this belt); conservative limits for the Sun: ~0.95–1.67 AU.
Key referencesKasting, Whitmire & Reynolds (1993); Kopparapu et al. (2013); Gillon et al. (2017).
|
| Zero-Dimensional Energy Balance & Ice–Albedo Feedback | \[ \frac{S}{4}\big(1-A(T)\big) = \sigma T_{\rm eff}^4 \]
Make albedo a function of temperature — ice is bright — and the simplest climate model becomes nonlinear: cooling grows ice, which reflects more light, which cools. Multiple equilibria appear, including a fully frozen one, and transitions between them are abrupt. |
S = stellar constant; A(T) = temperature-dependent albedo (ice ~0.6, ocean ~0.1) |
The origin of climate bistability and hysteresis: snowball states, deglaciation thresholds, and the warning that "in the habitable zone" does not mean "habitable now" — initial conditions matter.
Key referencesBudyko (1969); Sellers (1969); Hoffman et al. (1998); Pierrehumbert et al. (2011).
|
| Runaway Greenhouse Limit | \[ F_{\rm out} \le F_{\rm SN} \simeq 282\ {\rm W\,m^{-2}} \]
A moist atmosphere cannot radiate more than a ceiling flux: pile in more sunlight and the surface cannot cool by warming — oceans boil into a steam greenhouse, and water is then photolyzed and lost forever. The one-way door of planetary climate. |
F_SN = Simpson–Nakajima limit (for an Earth-like water world); Earth currently radiates ~240 W/m² |
Defines the habitable zone's hard inner edge, explains Venus's fate, and bounds Earth's own future: the brightening Sun crosses the limit in ~1–2 Gyr, and no thermostat can veto it.
Key referencesSimpson (1927); Ingersoll (1969); Nakajima et al. (1992); Goldblatt & Watson (2012).
|
| The Carbonate–Silicate Thermostat | \[ {\rm CO_2 + CaSiO_3 \rightleftharpoons CaCO_3 + SiO_2} \quad (\tau \sim 10^5\!-\!10^6\ {\rm yr}) \]
Weathering of silicate rock consumes CO₂ faster when it's warmer and wetter; volcanoes resupply it steadily. The loop is a planetary thermostat with a ~Myr response time — the reason Earth's climate survived a 30% dimmer young Sun. |
Weathering rate rises with T and rainfall; requires liquid water, exposed rock, and volcanic recycling |
The mechanism that widens the habitable zone (more CO₂ where colder) and the resolution of the faint young Sun paradox; its dependence on tectonics is why "habitability" quietly assumes geology, not just orbit.
Key referencesWalker, Hays & Kasting (1981); Sagan & Mullen (1972); Catling & Kasting (2017).
|
| The Drake Equation | \[ N = R_*\, f_p\, n_e\, f_l\, f_i\, f_c\, L \]
The number of communicating civilizations, factored into astrophysical, biological, and sociological probabilities. Not a predictive equation — an agenda: each factor is a research program, and the first three have been measured since 1961. |
R_* = star-formation rate; f_p, n_e = planet occurrence terms; f_l, f_i, f_c, L = life/intelligence/communication/lifetime |
The astrophysical factors are now data: f_p ≈ 1, n_e ~ 0.1–0.5 (this sheet's Sections V and XI). Everything past f_l is a prior — which is exactly what biosignature missions and technosignature searches exist to change.
Key referencesDrake (1961); Burchell (2006); Frank & Sullivan (2016).
|
Open unknowns · Climate & Habitability
f_l
Does life start easily — or was Earth a fluke?
One data point, with a survivorship bias. Life appeared fast on Earth (<0.5 Gyr after cooling), weakly suggesting ease — but abiogenesis has never been reproduced, and no second genesis (Mars, Enceladus, exoplanet biosignature) has been found. This is the largest unknown in science expressible in one symbol.
What Counts as a Biosignature?
Is any single gas — O₂, CH₄, DMS — ever conclusive?
Every proposed biosignature now has proposed abiotic sources (O₂ from water photolysis on M-dwarf planets is the cautionary tale). The field has retreated to "disequilibrium pairs + context"; the K2-18b DMS episode showed how fraught single-molecule claims are at JWST signal-to-noise.
Tidally Locked Climates
Can an eyeball planet — eternal day, eternal night — stay habitable?
GCMs say yes: ocean and atmospheric heat transport can prevent nightside atmospheric collapse, and the substellar "eyeball" ocean stays open. But cloud feedbacks at the substellar point differ by model, and the answer governs most of the Galaxy's HZ real estate.
Does Habitability Need a Moon, Plates, or a Jupiter?
Which of Earth's contingencies are requirements?
Obliquity stabilization by the Moon, tectonic thermostats, and Jupiter as impact shield are each argued to be essential and each argued to be irrelevant (Jupiter also sends impactors in). Rare-Earth versus mediocrity remains unsettled for want of a second inhabited example.
X
Small Bodies, Rings & Debris
6 equationsAsteroids, comets, Kuiper Belt objects, rings, and dust are the crumbs of planet formation — dynamically alive, collisionally evolving, and drifting under forces a planet never feels. Small bodies are where radiation pressure and thermal recoil compete with gravity, and where the Solar System's construction records survive unmelted.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Radiation Pressure (β) | \[ \beta = \frac{F_{\rm rad}}{F_{\rm grav}} \simeq \frac{0.57\, Q_{\rm pr}}{\rho\, s_{\rm \mu m}} \]
Sunlight pushes; gravity pulls; both fall as 1/r², so their ratio is a property of the grain alone. For sub-micron dust β exceeds ½ and the Sun effectively repels it — the smallest solids are blown out of the system entirely. |
s_μm = grain radius in microns; ρ = grain density (g/cm³); Q_pr ≈ 1, radiation-pressure efficiency |
Sets the blowout size (~0.5 µm for the Sun) — the lower cutoff of every debris-disk and zodiacal-dust grain distribution — and shapes comet dust tails into their curved syndyne fans.
Key referencesBurns, Lamy & Soter (1979); Wyatt (2008).
|
| Poynting–Robertson Drag | \[ t_{\rm PR} \simeq 400\,\frac{(r/{\rm AU})^2}{\beta}\ {\rm yr} \]
An orbiting grain sees sunlight arriving slightly head-on (aberration), so absorbing and re-emitting it costs orbital angular momentum. Dust quietly spirals into the Sun — centuries for micron grains at 1 AU. |
r = orbital distance; β = radiation-pressure ratio above |
Proves the zodiacal cloud is replenished (its dust dies in ~10³–10⁴ yr; comets and asteroid collisions resupply ~10⁴ kg/s) and cleans debris disks inside their birth rings — why inner holes are expected even without planets.
Key referencesPoynting (1904); Robertson (1937); Burns et al. (1979); Nesvorný et al. (2010).
|
| Yarkovsky Effect | \[ \frac{da}{dt} \;\propto\; \frac{\cos\gamma}{\rho\, s\, }\quad (\sim 10^{-4}\ {\rm AU/Myr\ at\ km\ size}) \]
A rotating asteroid radiates its afternoon heat from its evening side — a photon exhaust pipe pointed slightly along the orbit. Prograde rotators drift outward, retrograde inward: sunlight steers rocks. |
γ = spin obliquity; s = radius; strongest for ~0.1–10 km bodies; the torque cousin (YORP) spins them up/down |
The delivery mechanism of meteorites and near-Earth asteroids (drift into resonances → chaos → inner Solar System), the age-dater of asteroid families, and a term planetary-defense trajectory predictions cannot skip.
Key referencesYarkovsky (1901); Bottke et al. (2006); Chesley et al. (2014); Lauretta et al. (2019).
|
| Comet Sublimation Balance | \[ \frac{S_\odot(1-A)}{4\,r^2} = \epsilon\sigma T^4 + L_s\, Z(T) \]
Inside a critical distance, absorbed sunlight can no longer be balanced by radiation alone — ice sublimates, and the excess energy goes into throwing gas and dust off the nucleus. A comet is this equation switching branches. |
Z(T) = sublimation rate; L_s = latent heat; r = heliocentric distance (AU) |
Explains why water-ice comets switch on near ~3 AU (CO and CO₂ comets much farther out), sets production rates Q(H₂O) measured by spectroscopy, and underlies activity models for interstellar visitors.
Key referencesWhipple (1950); Cowan & A'Hearn (1979); Hansen et al. (2016); Bodewits et al. (2020).
|
| Collisional Cascade | \[ \frac{dN}{ds} \propto s^{-3.5} \]
Grind a population of bodies against itself long enough and the size distribution forgets its initial condition, settling to a universal power law (Dohnanyi's −3.5 for self-similar strength). Most of the mass sits in the biggest bodies; most of the surface area — and all the observable dust — in the smallest. |
N(s) = number per size bin; slope steepens/flattens where material strength or gravity regime changes |
The engine of debris disks: unseen Pluto-sized stirrers keep km-belts grinding, and the infrared excess we detect around ~20% of AFG stars is the cascade's dust floor. Deviations from −3.5 in the asteroid belt date its collisional history.
Key referencesDohnanyi (1969); Wyatt (2008); Hughes, Duchêne & Matthews (2018).
|
| Impact-Crater Scaling | \[ D_{\rm crater} \;\propto\; E^{\,1/3.4} \]
Crater size grows as a weak power of impact energy (gravity-regime π-scaling), so a modest energy range spans the Moon's face. Run backwards, every cratered surface is a particle detector with a 4-Gyr exposure time. |
E = impact kinetic energy; full π-group scaling adds gravity, velocity, angle, target strength |
Crater counting is the only chronometer for surfaces we can't sample: calibrated on Apollo samples, exported (with debated corrections) to Mars, Mercury, and the outer-planet moons. Also the sizing law for planetary defense.
Key referencesHolsapple (1993); Melosh (1989); Neukum et al. (2001); Daly et al. (2023).
|
Open unknowns · Small Bodies & Debris
The Interstellar Visitors
Why were 'Oumuamua, Borisov, and 3I/ATLAS all so different?
Three interstellar objects, three phenomenologies: one dry and non-gravitationally accelerating, one a hyperactive CO comet, one massive and ancient-seeming. Each is a core sample of another system's small-body population; Rubin Observatory should turn anecdotes into a distribution.
Rubble All the Way Down?
Bennu and Ryugu barely hold together — is anything under 1 km monolithic?
Both sampled asteroids are shockingly weak rubble piles (Bennu's surface behaved like a ball pit). Strength assumptions propagate directly into deflection strategies and cascade models; the answer may vary by taxonomic class, and only sample-return tells.
How Old Are Saturn's Rings?
Primordial, or younger than the dinosaurs?
Cassini's finale measured a low ring mass and high pollution-based youth (~10⁸ yr) — but "young rings" demand destroying a Mimas-sized icy body recently, which is dynamically awkward. Recent work argues appearances can deceive; the debate is genuinely open.
XI
Demographics: The Exoplanet Census
5 equationsWith thousands of planets, the science graduates from objects to populations: occurrence rates, distributions, and the features carved into them. The census is where formation theory faces its exam — every mechanism on this sheet must reproduce these statistics or die.
| Name | Equation | Variables | Use in Research |
|---|---|---|---|
| Occurrence Rate | \[ \eta = \frac{1}{N_*}\sum_{\rm det} \frac{1}{p_{\rm tr}\; C} \]
Count detections, then repay the debts: divide by the geometric transit probability and the pipeline completeness C for each. What survey selection taketh away, statistics giveth back — the true abundance of planets per star. |
N_* = searched stars; C = detection completeness (injection–recovery); vetting reliability corrects false alarms |
The framework of every Kepler demographic result; its headline outputs — planets outnumber stars; ~half of Sun-like stars host a 1–4 R⊕ planet within 1 AU — are the boundary conditions of formation theory.
Key referencesHoward et al. (2012); Petigura et al. (2013); Bryson et al. (2021).
|
| The Radius Valley | \[ R_{\rm valley} \simeq 1.8\left(\frac{P}{10\,{\rm d}}\right)^{-0.1} R_\oplus \]
The small-planet population is bimodal: super-Earths (~1.3 R⊕) and sub-Neptunes (~2.4 R⊕), separated by a genuine desert whose location slides down with period. A gap this clean is a physical process caught in the act. |
P = orbital period; slope ~ −0.1 favors photoevaporation/core-powered loss over gas-poor formation (+ slope) |
The census's sharpest feature and the direct fingerprint of atmospheric escape (Section VIII): sub-Neptunes are cores that kept ~1% H/He envelopes; super-Earths are the same cores stripped bare. Its slope, width, and stellar-mass dependence discriminate the stripping mechanism.
Key referencesOwen & Wu (2013, 2017); Fulton et al. (2017); Ginzburg et al. (2018); Van Eylen et al. (2018).
|
| Giant Planet–Metallicity Relation | \[ P({\rm giant}) \propto 10^{\,2.0\,[{\rm Fe/H}]} \]
Double a star's metal content and its odds of hosting a gas giant roughly quadruple. Planets are made of metals; giants are made of planets that grew fast — a straight line from stellar spectroscopy to core accretion. |
[Fe/H] = log iron abundance relative to solar; relation strong for giants, nearly flat for small planets |
The first population-level confirmation of core accretion (more solids → cores beat the disk clock more often) and a live target-selection tool; its absence for small planets says their formation is not solid-starved anywhere.
Key referencesGonzalez (1997); Santos et al. (2004); Fischer & Valenti (2005); Johnson et al. (2010).
|
| Eccentricity Distributions | \[ p(e)\, \sim\, {\rm Beta}(a,b): \;\langle e\rangle_{\rm giants} \approx 0.2\!-\!0.3 \gg \langle e\rangle_{\rm SS} \]
Solar System orbits are nearly circular; exoplanet giants are not. The eccentricity distribution is a dynamical thermometer — reading the violence (scattering, Kozai) each population survived after the disk's damping died. |
Beta/Rayleigh forms fit populations; Kepler compact multis are cold (σ_e ~ 0.03), single-transit systems warmer |
Constrains the frequency of planet–planet scattering (it reproduces the giant e-distribution almost too well) and flags system classes: dynamically cold flat multis versus hot, mutually inclined singles — the "Kepler dichotomy."
Key referencesJurić & Tremaine (2008); Chatterjee et al. (2008); Van Eylen & Albrecht (2015); Xie et al. (2016).
|
| Planet Mass Function | \[ \frac{dN}{d\log M} \;\uparrow\; {\rm toward\ low\ } M \quad (\propto M^{-0.5\ \rm to\ -1}) \]
Nature makes small planets prodigiously and giants rarely, with breaks in the power law marking physical thresholds — runaway gas accretion (~0.1–0.3 M_J deserts) and the stripping valley. The mass function is the formation process's final answer sheet. |
Slope from RV + microlensing + transit surveys jointly; Neptune desert and brown-dwarf desert are carved features |
The synthesis product of the whole census — the function population-synthesis models must hit; microlensing extends it to cold orbits and (via short events) below Earth mass, where bound and free-floating branches meet.
Key referencesCumming et al. (2008); Mayor et al. (2011); Suzuki et al. (2016); Fulton et al. (2021).
|
Open unknowns · Demographics
Is the Solar System Rare?
Cold Jupiter, empty inner zone, circular orbits — how unusual is our architecture?
Jupiter analogs orbit ~10% of Sun-like stars; inner super-Earths are absent here but present around ~half; our eccentricities are unusually low. Multiply the factors and the Solar System looks ~1-in-100 — but the factors' independence is untested, and Gaia DR4's cold-giant census will sharpen it.
Sub-Neptunes: Gas, Water, or Both?
What is the most common type of planet in the Galaxy actually made of?
The 2–3 R⊕ population fits either rocky cores with H/He skins or ice-rich "water worlds" from beyond the snow line. JWST atmospheres (K2-18b, TOI-270d, GJ 1214b's high-metallicity air) hint at volatile-rich interiors for some — resolving the mix rewrites the water budget of the Galaxy.
Planets of the Majority
M dwarfs are 75% of stars — do their compact systems follow different rules?
M dwarfs host more small planets per star (~2.5 with P < 200 d) but far fewer giants — core accretion starved of time and mass, as theory expects. Whether their HZ planets kept water and air (Section VIII) decides if the Galaxy's typical planetary system is habitable or a cinder field.
Planetary reference values: 1 AU = 1.496×10¹¹ m; M⊕ = 5.97×10²⁴ kg, R⊕ = 6,371 km; M_J = 318 M⊕ = 9.5×10⁻⁴ M☉, R_J = 11.2 R⊕; Earth: v_esc 11.2 km/s, H = 8.5 km, T_eq 255 K (A_B = 0.29), P_c 364 GPa, C̄ = 0.331; Jupiter: r_H 0.35 AU, k₂ ≈ 0.56, B_eq ≈ 4.2 G; Sun's RV wobble: 12.5 m/s (Jupiter), 0.09 m/s (Earth); transit depths: 1.05% (Jupiter–Sun), 84 ppm (Earth–Sun); snow line ≈ 2.7 AU (young Sun); conservative HZ 0.95–1.67 AU; runaway-greenhouse ceiling ≈ 282 W/m² (Earth radiates ~240); disk half-life ~2–3 Myr; MMSN Σ(1 AU) ≈ 1,700 g/cm² gas; blowout grain size ~0.5 µm; Io tidal power ~10¹⁴ W; radius valley ≈ 1.8 R⊕; hot-Jupiter occurrence ~1%; ~6,000 confirmed exoplanets (2026).