The Shell Theorem and Spherical Mass Distributions

Every gravitation calculation you have done so far contains a quiet act of faith. Computing your weight, you plugged in $r = 6371$ km — the distance to Earth's centre — as if all $5.97 \times 10^{24}$ kg of planet sat compressed at that one point. But the rock under your feet is only metres away and pulls ferociously; the rock under New Zealand is a planet-width away — some twelve million times farther, and so weaker per kilogram by the square of that. Why should this lopsided tug-of-war average out to exactly the point-mass answer, rather than roughly? Newton asked the same question and could not answer it for years — by some accounts this single gap delayed the Principia. The answer, when he found it, is the shell theorem, really two theorems: a uniform spherical shell attracts any external mass exactly as if the shell were a point at its centre; and the same shell exerts zero net force on a mass anywhere inside it — not just at the middle, anywhere. Both are miracles of the inverse square: change the exponent and both fail. This article proves the pair honestly — one integral, done in full — and then spends the winnings. A solid sphere is nested shells, so theorem one licenses every "centre-to-centre" step in the track. Theorem two answers a stranger question: what does gravity do inside a planet? The answer (it fades linearly to zero at the centre, in the uniform model) turns a tunnel through the Earth into a giant harmonic oscillator with a 42-minute one-way ride, and hands stellar and planetary science its basic working tool.

Key takeaways

Shell theorem 1 (outside): a uniform spherical shell of mass $M$ attracts an external mass exactly as a point mass $M$ at its centre — $F = GMm/r^2$ , $U = -GMm/r$ , no approximation.
Shell theorem 2 (inside): everywhere inside the shell, the net gravitational force is exactly zero — near-wall and far-wall pulls cancel patch by patch.
A solid sphere is nested shells, so for an external point only $r$ to the centre matters — the licence behind every $g = GM/R^2$ in the track.
Inside a body only the mass interior to your radius counts: for uniform density, $g(r) = GMr/R^3$ — gravity fades linearly to zero at the centre.
A straight tunnel through a uniform Earth gives simple harmonic motion: about 42 minutes to fall through, same period as a surface-skimming orbit.
Both theorems are exclusive perks of the $1/r^2$ law — the same geometry that gives Gauss's law in electromagnetism.
Key idea: for inverse-square forces, spheres act as points from outside and shells vanish from inside — geometry and the force law cancel exactly.

The core insight

For an inverse-square force, a uniform spherical shell gravitates like a point mass at its centre for exterior points and exerts zero net force on interior points; consequently a spherically symmetric body of any layered density pulls as a point from outside, while inside it only the mass beneath your feet counts.

Where this lives today

Equation (5) is the working equation of planetary interiors. Seismologists invert earthquake travel times for the density profile $\rho(r)$ , and the shell decomposition turns that profile into the interior gravity, pressure, and moment of inertia that models must match — this is how we know Earth's outer core is liquid iron and how the Cassini and Juno missions, tracking tiny Doppler shifts in spacecraft radio, weigh the layers inside Saturn and Jupiter without touching them. The same theorem runs stellar structure: a star is spherical gravity feeding on its own enclosed mass, and $g(r) = GM(r)/r^2$ appears on the first page of every stellar-interior text. The failures of spherical symmetry are an industry of their own. Once the shell theorem tells you what a layered sphere would do, every measured deviation is signal: satellite gravimetry (GRACE, GOCE) maps mountain roots, ocean currents, and vanishing aquifers as parts-per-million wrinkles on the shell-theorem baseline. And the interior cancellation has a laboratory career — inside a spherical (or well-designed) shell, external-style gravity gradients cancel, a principle exploited in gravitational-wave detector environmental modelling and in the "drag-free" test masses of missions like LISA Pathfinder, which float in a shielded pocket where the only surviving field is the one being measured.

The physical picture

Start with the theorem that sounds false. Stand inside a spherical shell, well off-centre, close to one wall. The near wall looms — surely it wins? Draw a narrow double cone through your position: it cuts a small patch of shell nearby and a large patch far away. The near patch is close (strong pull per kilogram) but small; the far patch is weak but big. For a cone of small opening angle, patch area grows as the square of the distance $s$ to the wall, while the pull of each kilogram falls as $1/s^2$ — the two squares cancel exactly, and the two patches tug with equal and opposite force. Sweep the double cone over every direction: the whole shell cancels, patch by patch, wherever you stand inside. The balance is perfect only for the exponent 2; a $1/r^{2.01}$ force would pull you to the near wall.

Figure

Inside a shell — every double cone cancels exactly

Now the outside theorem. Slice a shell into rings centred on the line to the external mass. Each ring is a set of points all at the same distance $s$ , pulling along the axis by symmetry (the sideways components of opposite ring-elements cancel). Adding rings means adding simple axial pulls — an honest one-variable integral, done below. The miraculous part is the answer: the sum over rings, near-side rocks and far-side rocks and everything between, collapses to exactly $GMm/r^2$ , as if the shell had packed itself into its centre. And once one shell behaves this way, any sphere built of concentric shells — uniform or not, so long as density depends only on radius — inherits the point-mass property, layer by layer. Earth's iron core, rocky mantle, thin crust: from outside, one point. The two theorems together produce the strangest picture in Newtonian gravity: the journey to the centre of the Earth. Descend a shaft, and every shell above your radius switches off (theorem 2) while every shell below still pulls like a central point (theorem 1). Gravity does not grow toward the centre — in the uniform model it weakens, linearly, reaching exactly zero at the core, where all shells surround you and all cancel.

Figure

Gravity of a uniform sphere — linear inside, inverse-square outside

The math: one integral, honestly

Everything rests on a single computation: the potential energy of a point mass $m$ at distance $r$ from the centre of a uniform shell (mass $M$ , radius $R$ ). We work with the potential rather than the force for a practical reason worth naming: potential contributions are scalars — they add without direction — whereas force contributions are vectors whose axial projections drag a cosine factor through the integral. Same physics, half the trigonometry; the force at the end is one derivative away.

Derivation AThe shell integral — both theorems in one computation

We are deriving the potential energy of a mass $m$ at distance $r$ from the centre of a uniform spherical shell of mass $M$ and radius $R$ , and reading both shell theorems from the result. The strategy: slice the shell into rings centred on the axis through $m$ , add each ring's scalar contribution $-Gm\,dm/s$ , and convert the awkward angular integral into a clean integral over the ring distance $s$ itself. The one step where all the magic lives is the change of variables in Step 3 — after it, the inverse-square law's signature has quietly arranged for the answer to depend on the limits alone, and the limits are where inside differs from outside.

Step 1. Slice into rings. The ring at polar angle $\theta$ (measured from the axis to $m$ ) with angular width $d\theta$ has circumference $2\pi R\sin\theta$ and width $R\,d\theta$ ; since the shell is uniform, its mass is the area fraction of the total: $dm = \dfrac{M}{4\pi R^2}\,(2\pi R^2 \sin\theta\,d\theta) = \dfrac{M}{2}\sin\theta\,d\theta$ . Why rings: every point of a ring sits at the same distance from $m$ , so a ring contributes one clean scalar term.

Step 2. Write the distance to a ring by the law of cosines: $s^2 = r^2 + R^2 - 2rR\cos\theta$ , and the total potential energy as

U = -Gm\int \frac{dm}{s} = -\frac{GmM}{2}\int_0^{\pi} \frac{\sin\theta\,d\theta}{s} \qquad (1)

Why the potential: contributions add as numbers; the force route would need each ring's axial projection $\cos\alpha$ and a second application of the cosine law.

Step 3. Change variables from $\theta$ to $s$ . Differentiating the cosine law at fixed $r, R$ : $2s\,ds = 2rR\sin\theta\,d\theta$ , so $\sin\theta\,d\theta = \dfrac{s\,ds}{rR}$ . Substituting, the integrand collapses:

U = -\frac{GmM}{2rR}\int s\,\frac{ds}{s}\cdot 1 = -\frac{GmM}{2rR}\int_{s_{\min}}^{s_{\max}} ds \qquad (2)

Why this is the miracle step: the $1/s$ of the potential and the $s\,ds$ of the geometry cancel, leaving the integrand equal to one. The entire answer now sits in the limits.

Step 4. Read the limits for an external point ( $r > R$ ): the nearest ring is at $s_{\min} = r - R$ ( $\theta = 0$ ), the farthest at $s_{\max} = r + R$ ( $\theta = \pi$ ). The integral is $(r+R) - (r-R) = 2R$ :

U_{\rm out} = -\frac{GmM}{2rR}\,(2R) = -\frac{GMm}{r} \qquad (3)

Why this is theorem 1: this is precisely the potential of a point mass $M$ at the centre — and differentiating, $F = -\dfrac{dU}{dr} = -\dfrac{GMm}{r^2}$ , attractive, exactly the point-mass force. No approximation anywhere.

Step 5. Read the limits for an internal point ( $r < R$ ): now the nearest ring is at $s_{\min} = R - r$ , the farthest still at $R + r$ . The integral is $(R+r) - (R-r) = 2r$ :

U_{\rm in} = -\frac{GmM}{2rR}\,(2r) = -\frac{GMm}{R} \qquad (4)

Why this is theorem 2: the answer does not contain $r$ — the potential is constant throughout the interior, so the force $F = -dU/dr = 0$ everywhere inside, off-centre included. Constant potential is the calculus form of the double-cone cancellation.

Step 6. Stack shells. A spherically symmetric body is concentric shells; each obeys (3)–(4) independently. From outside, all shells act as central points, so the whole body does. From radius $r$ inside, shells with radius greater than $r$ contribute nothing, shells beneath contribute as a central point of mass $M(r)$ , the mass interior to $r$ :

g(r) = \frac{G\,M(r)}{r^2} \qquad (5)

Why this is the working tool: (5) holds for any radial density profile — uniform, layered, realistic Earth — which is what makes it the backbone of planetary and stellar structure. *As a consistency check: equations (3) and (4) agree at $r = R$ , both giving $U(R) = -GMm/R$ — the potential is continuous at the shell surface, as it must be. Dimensions balance throughout, and the results correspond to Newton's Principia Book I, Propositions LXX–LXXI.*

Go deeper · Derivation B · 2 minGravity inside a uniform planet, and the tunnel through the Earth

Derivation BGravity inside a uniform planet, and the tunnel through the Earth

The goal is the interior field of a uniform-density sphere, $g(r) = GMr/R^3$ , and its famous corollary: a mass dropped into a straight tunnel through the centre executes simple harmonic motion with the same period as a surface-skimming orbit, about 84 minutes round trip. The strategy is equation (5) with the uniform-density mass function, then a comparison with the SHM template $a = -\omega^2 x$ from the oscillations chapters. Watch the cancellation in Step 1 — the $r^3$ of enclosed mass against the $r^2$ of the force law — which is where "linear in $r$ " comes from.

Step 1. Enclosed mass at radius $r$ : with uniform density $\rho = \dfrac{M}{\tfrac43\pi R^3}$ , $\ M(r) = \tfrac43\pi\rho\,r^3 = M\dfrac{r^3}{R^3}$ . Insert into (5):

g(r) = \frac{G}{r^2}\cdot M\frac{r^3}{R^3} = \frac{GM}{R^3}\,r \qquad (6)

Why linear: the pull grows like the enclosed mass ( $r^3$ ) but is diluted by distance squared — net one power of $r$ . At $r = R$ this meets the exterior branch $GM/r^2$ continuously.

Step 2. Drop a mass into a diametric tunnel (idealized: no drag, no rotation). At displacement $x$ from the centre, the force is $F = -m\,g(x) = -\dfrac{GMm}{R^3}\,x$ . Why the minus sign: the pull always aims at the centre — a restoring force proportional to displacement, the harmonic-oscillator signature.

Step 3. Match to SHM: $a = -\omega^2 x$ with $\omega^2 = \dfrac{GM}{R^3}$ , so the period is

T = 2\pi\sqrt{\frac{R^3}{GM}} = 2\pi\sqrt{\frac{(6.371 \times 10^6)^3}{3.986 \times 10^{14}}} = 5060\ \mathrm{s} \approx 84\ \mathrm{min} \qquad (7)

Why the number matters twice: one-way, Pole to Pole, takes $T/2 \approx 42$ minutes — and (7) is exactly the period of a circular orbit at radius $R$ , since that orbit's $\omega^2 = GM/R^3$ too. The tunnel traveller and a surface-skimming satellite, released together, would meet again at the far side.

Step 4. Add the honest caveat: the real Earth is not uniform — density climbs from about $2700\ \mathrm{kg/m^3}$ in crustal rock to about $13{,}000$ in the inner core. Feeding the seismological density profile into (5), $g$ actually rises slightly with depth through the mantle (enclosed mass thins more slowly than $r^2$ shrinks), peaking near $10.7\ \mathrm{m/s^2}$ at the core–mantle boundary before falling to zero at the centre. Why this matters: (6) is a model result; (5) is the theorem. Only the uniform-density step used the model. *The arithmetic step by step: $(6.371 \times 10^6)^3 = 2.586 \times 10^{20}$ ; dividing by $3.986 \times 10^{14}$ gives $6.49 \times 10^5$ ; taking the square root gives $805\ \mathrm{s}$ ; multiplying by $2\pi$ gives $5060\ \mathrm{s} = 84.3\ \mathrm{min}$ , confirming the result in equation (7). The 42-minute half-trip and the coincidence with the orbital period are standard results.*

A remark that will pay dividends later in your physics: theorem-grade statements this clean usually signal a structural reason, and here it is flux. An inverse-square field spreads over spheres whose area grows as $r^2$ , so the "flow" of field lines through any closed surface counts only the mass inside — that is Gauss's law, which the electromagnetism track develops as its central tool. Equations (3)–(5) are Gauss's law for gravity, evaluated on spheres; the hard integral above is what the flux argument compresses into one line, and the E&M article on Gauss's law is the natural next place to see the machinery.

Limiting cases

Case	Result	Meaning
$r \gg R$	$U \to -GMm/r$	any compact blob looks like a point from far away — here it is exact at every $r > R$
$r \to R^{\pm}$	$U \to -GMm/R$ from both sides	potential continuous at the shell; force jumps (surface layer)
anywhere inside a shell	$F = 0$ , $U$ constant	perfect cancellation, off-centre included
centre of any spherical body	$g = 0$	all shells surround you; all cancel
uniform sphere, interior	$g = GMr/R^3$ , linear	tunnel SHM, $T \approx 84$ min
real Earth (PREM density)	$g$ peaks $\approx 10.7\ \mathrm{m/s^2}$ at core–mantle boundary	theorem (5) survives; the uniform model does not
force law $\ne 1/r^2$	both theorems fail	inside a shell you would drift wall-ward; spheres would not act as points

Common trap"Inside a shell, the nearer wall pulls you toward it"

The intuition is reasonable — closer means stronger — and it is exactly wrong, because it forgets that how much shell sits in each direction is also position-dependent. Off-centre, less shell is near you and more is far, in the precise proportion ( $\text{area} \propto s^2$ ) that cancels the strength advantage ( $\propto 1/s^2$ ). The net force inside is zero everywhere, not just at the centre. Two siblings of the error to refuse: first, "zero force means zero potential" — the interior potential is a nonzero constant, $-GMm/R$ , and it still costs energy to remove a mass from the shell's interior to infinity; second, "so gravity is zero inside the Earth" — the Earth is not a hollow shell but a stack of them, and only the shells above your radius switch off. At the bottom of the deepest mine, the overlying rock's pull is gone, but almost the entire planet still lies beneath you.

Worked examples

1 · Weight at the bottom of a deep shaft (confidence builder)

In the uniform-density model of the Earth ( $M = 5.97 \times 10^{24}\ \mathrm{kg}$ , $R = 6371\ \mathrm{km}$ ), what does a 60 kg student weigh at the bottom of a 3000 km shaft — nearly halfway to the centre?

Step 1. Locate the radius: $r = 6371 - 3000 = 3371\ \mathrm{km}$ , so $r/R = 0.529$ .

Step 2. Shells above $r$ contribute nothing (theorem 2); shells below act as a central point of mass $M(r) = M(r/R)^3$ (theorem 1). By equation (6), $g(r) = g_{\rm surface} \times (r/R)$ .

Step 3. Compute: $g = 9.81 \times 0.529 = 5.19\ \mathrm{m/s^2}$ , so the weight is $W = 60 \times 5.19 = \boxed{311\ \mathrm{N}}$ — down from 589 N at the surface, despite being 3000 km closer to the centre.

Step 4. Contrast with the exterior rule to see the theorems working: at 3000 km altitude ( $r = 9371$ km), the inverse square gives $g = 9.81 \times (6371/9371)^2 = 4.53\ \mathrm{m/s^2}$ . Above the surface, distance dilutes the pull; below it, disappearing shells do. (In the real, dense-cored Earth, $g$ at 3371 km radius is actually close to $10\ \mathrm{m/s^2}$ — the uniform model's one honest casualty, per Step 4 of the tunnel derivation.)

2 · The off-centre cavity — superposition meets the theorem (exam level)

A uniform sphere of density $\rho = 5500\ \mathrm{kg/m^3}$ contains a spherical cavity whose centre sits $d = 1000\ \mathrm{km}$ from the sphere's centre. Show that the gravitational field everywhere inside the cavity is uniform, and compute it.

Step 1. Superpose: the actual body = (full uniform sphere) + (a sphere of negative density $-\rho$ filling the cavity). Gravity adds linearly, so the field is the sum of the two spheres' fields.

Step 2. Write the interior field of a uniform sphere in vector form. From equation (6), at position $\vec r$ from the centre: $\vec g = -\dfrac{4}{3}\pi G\rho\,\vec r$ — inward, linear in the position vector. This vector form is the key trick and it is exact inside.

Step 3. Add the two contributions at a point $P$ inside the cavity. Let $\vec r$ run from the big sphere's centre to $P$ , and $\vec r\,'$ from the cavity's centre to $P$ : $\vec g = -\tfrac43\pi G\rho\,\vec r + \tfrac43\pi G\rho\,\vec r\,' = -\tfrac43\pi G\rho\,(\vec r - \vec r\,')$ .

Step 4. Recognize the geometry: $\vec r - \vec r\,' = \vec d$ , the fixed vector from the big centre to the cavity centre — independent of $P$ . So everywhere in the cavity,

\vec g = -\frac{4}{3}\pi G\rho\,\vec d

uniform in magnitude and direction, pointing from the cavity centre toward the sphere's centre.

Step 5. Numbers: $g = \tfrac43\pi \times 6.674 \times 10^{-11} \times 5500 \times 10^6 = 2.80 \times 10^{-10} \times 5500 \times 10^6 = \boxed{1.54\ \mathrm{m/s^2}}$ , the same at every point of the cavity. A dust grain released anywhere inside falls along the same parallel field lines — a pocket of perfectly uniform gravity manufactured out of two applications of the shell theorem.

Origins · 1 minIsaac Newton and how this was discovered

Isaac Newton

1643–1727 · Woolsthorpe & Cambridge

The shell theorem is the proof Newton owed himself. His 1666 Moon test treated the Earth as a point at its centre — a step he could justify at lunar distance but not for the apple, a few metres from the nearest rock and 12,700 km from the farthest. The gap gnawed at him; historians read his long delay in publishing gravitation partly as the wait for this proof, finally achieved with the Principia's geometric machinery as Propositions LXX and LXXI of Book I (1687). Proposition LXX is the interior cancellation, argued by exactly the double-cone construction of the hero figure; LXXI is the exterior point-equivalence. With them, "centre to centre" stopped being a hopeful approximation and became a theorem, and the apple and Moon calculation closed honestly. It is a fine irony that the modern one-page proof — the potential integral above — uses calculus Newton possessed but chose to hide behind Euclidean geometry, to keep the Principia in the mathematics his readers would trust.

History

1687

— Principia, Book I, Props. LXX–LXXI: interior force vanishes, exterior sphere acts as a point — the licence for the Moon test, proved at last.

1798

— Cavendish weighs the Earth: mean density 5.45 times water against surface rock's 2.7 — the interior must be dense, and equation (5) becomes a probe of what lies beneath.

1835

— Gauss formulates the flux law for inverse-square fields — the shell theorem's one-line generalization, later the backbone of electrostatics.

1909

— Mohorovičić spots the crust–mantle seismic discontinuity: layered-Earth models, all resting on the shell decomposition, begin in earnest.

1981

— Dziewonski and Anderson publish PREM, the reference density profile

\rho(r)

; run through equation (5), it shows

g

nearly constant to 2890 km depth, peaking at the core–mantle boundary.

Going further

Texts. Halliday, Resnick & Walker, Fundamentals of Physics, Ch. 13 · Young & Freedman, University Physics, Ch. 13 · Kleppner & Kolenkow, An Introduction to Mechanics, Ch. 9 · Newton, Principia, Book I, Props. LXX–LXXI for the geometric original · any E&M text's Gauss's-law chapter for the flux-side view. Key results to keep in mind: the ring mass element $dm = (M/2)\sin\theta\,d\theta$ ; the substitution $s\,ds = rR\sin\theta\,d\theta$ with limits $(r \pm R)$ for an exterior point and $(R \pm r)$ for an interior point; the tunnel period $T = 2\pi\sqrt{R^3/GM} = 5060\ \mathrm{s} \approx 84.3$ min; worked example 1 gives $g = 9.81 \times 0.529 = 5.19\ \mathrm{m/s^2}$ and $W = 311$ N; worked example 2 gives $\vec g = -\tfrac43\pi G\rho\,\vec d$ with magnitude $1.54\ \mathrm{m/s^2}$ for $\rho = 5500\ \mathrm{kg/m^3}$ and $d = 10^6\ \mathrm{m}$ .

Builds on: Newton's Law of Universal Gravitation · Gravitational Potential Energy and Escape Velocity · Work and the Line Integral of Force. Leads to: Central Force Fundamentals · The Two-Body Problem and Reduced Mass · Simple Harmonic Motion (the tunnel oscillator) · Gauss's Law in the E&M track.

Does the shell theorem require the shell to be thin, or uniform, or is it fully general?

The two theorems are proved for a thin, uniform shell — that is the honest statement. But the useful generality comes free by stacking: any body whose density depends only on radius, $\rho(r)$ , is a nest of thin uniform shells, each obeying both theorems independently. So layered planets qualify perfectly (crust, mantle, core — each layer is spherically symmetric), and the exterior point-mass rule is exact for them. What breaks the theorem is angular structure: mountains, tidal bulges, a lopsided core. Then the exterior field picks up correction terms (physicists organize them as multipole moments), which is exactly why satellite orbits wobble in ways that map Earth's lumpiness. The theorem is the clean baseline; the lumps are the science on top.

The interior force is zero but the interior potential is not. Doesn't zero force mean no gravity at all in there?

Zero force means no push — a mass inside a hollow shell floats indifferently, anywhere. But potential measures the energy cost of getting out, and that is real: lifting a mass from inside the shell to infinity costs $GMm/R$ , paid entirely in the climb from the shell's surface outward (the trip from the interior to the wall is free — constant potential means level ground). The distinction becomes vivid if you make the shell massive and compact: clocks deep inside a hollow massive shell tick slower than clocks at infinity (a general-relativity effect governed by the potential, not the force), even though a resident feels no pull whatsoever. Force is the slope of the potential; inside the shell the landscape is a flat plateau — but a plateau at the bottom of a well.

Why does the potential trick make the integral so much easier than integrating the force directly?

Because potential contributions are scalars and force contributions are vectors. Each mass element pulls toward itself, so summing forces means projecting every pull onto the symmetry axis — an extra factor of $\cos\alpha$ , where $\alpha$ is the angle between the element's pull and the axis, and that cosine must itself be written via the law of cosines. The integral remains doable (Newton effectively did its geometric equivalent) but doubles in trigonometric weight. Potentials just add as numbers; then one derivative at the very end, $F = -dU/dr$ , recovers the force with the direction supplied by symmetry. This "energy first, differentiate last" economy is a habit worth internalizing — it is the same reason the Lagrangian half of the track works with energies rather than forces.

Is the 42-minute tunnel ride real? Could a gravity train actually work?

The number is real physics, the train is not practical. The idealization ignores three fatal engineering facts: air drag (evacuating a planet-length tunnel), friction, and above all temperature and pressure — the route passes through rock at thousands of kelvin and millions of atmospheres. But the idealized result contains two lovely truths. First, the period is independent of the endpoints: a straight chord tunnel between any two cities (not through the centre) also gives SHM with the same 84-minute period — gravity's component along the chord is proportional to displacement along the chord, with the same $\omega^2 = GM/R^3$ — so every chord ride one-way is about 42 minutes, London to Paris or Pole to Pole. Second, the coincidence with the orbital period is exact and not accidental: a surface-skimming circular orbit is two perpendicular tunnel oscillations run 90° out of phase, exactly as circular motion decomposes into two SHMs.

Real Earth's gravity increases as you descend into the mantle. Doesn't that contradict the shell theorem?

It contradicts the uniform-density model, not the theorem. The theorem's working form, $g(r) = GM(r)/r^2$ , is agnostic about the density profile; feed it uniform $\rho$ and you get the linear law $g \propto r$ . But the real Earth stashes a third of its mass in an iron core occupying just 16% of its volume. Descending through the light mantle, you discard shells that were contributing little, while $r^2$ in the denominator shrinks fast — so $g$ rises, gently, from $9.81$ at the surface to about $10.7\ \mathrm{m/s^2}$ at the core–mantle boundary (2890 km down), and only then plunges toward zero through the dense core. The lesson generalizes: interior gravity profiles are density diagnostics, which is precisely how planetary scientists use them.

Is this the same as Gauss's law from electromagnetism?

Same theorem, different force. Gauss's law says the flux of an inverse-square field through any closed surface counts only the source enclosed — for gravity, $\oint \vec g \cdot d\vec A = -4\pi G M_{\rm enc}$ . Evaluate it on a sphere concentric with a spherically symmetric mass and both shell theorems fall out in one line: outside, all of $M$ is enclosed ( $g = GM/r^2$ ); inside a shell, nothing is ( $g = 0$ ). The hard ring integral in this article is what that one line is secretly summarizing — worth doing once honestly, precisely so that Gauss's law never feels like magic. The correspondence runs deep because both forces are inverse-square: swap $m \to q$ and $-G \to 1/4\pi\epsilon_0$ and every result here becomes electrostatics of charged shells, including the zero field inside a uniformly charged sphere's cavity. The E&M track's Gauss's law article picks the story up from the flux side.

Why exactly does the exponent 2 make everything work? What breaks for other force laws?

Two separate cancellations both demand the 2. Inside: the double-cone argument balances patch area (growing as $s^2$ ) against patch pull per kilogram (falling as $s^{-2}$ ) — any other exponent leaves a net force, dragging you wall-ward for steeper laws and centre-ward for shallower ones. Outside: in the shell integral, the substitution $\sin\theta\,d\theta = s\,ds/rR$ cancels the potential's $1/s$ perfectly, collapsing the integrand to a constant; for a potential $1/s^p$ with $p \ne 1$ the integral still runs but no longer resembles a point mass at the centre. There is a deeper way to say it: in three dimensions, the inverse square is the unique central force whose field lines spread without being created or destroyed in empty space (zero divergence), which is also why it alone has a Gauss's law. The shell theorem is 3D geometry recognizing its matching force law.