Free coursewareRewritten from MIT OpenCourseWareCC BY-NC-SA 4.0
Open CosmologyPhysics & cosmology: four years, plus a graduate year

The expanding sphere and the equation for a(t)

Newton's gravity applied to a uniformly expanding ball of matter gives the equation that governs the scale factor — and a conserved quantity that decides the fate of the universe.

You need this first

Full reading path (8 pages)
  1. Vectors (8.01)
  2. Position and velocity in one dimension (8.01)
  3. Acceleration in one dimension (8.01)
  4. Motion in two dimensions and projectiles (8.01)
  5. Newton's laws of motion (8.01)
  6. The universal law of gravitation (8.01)
  7. The homogeneous expanding universe (8.286)
  8. Why a static universe cannot work (8.286)
On this page · 6 sections
  1. The model
  2. Following the shells
  3. The cancellation
  4. Density and deceleration
  5. A conserved quantity
  6. Summary

Gravity cannot leave a uniform universe standing still, as the previous page showed. So how does it change? In Lecture 6 Guth answers with a model that uses only Newton's law of gravity — and promises that, for ordinary matter, it gives exactly the same answer as general relativity.

The model#

To keep infinity under control, start with something finite and take a limit at the end:

  • At some initial time tit_i, take a sphere of radius Rmax⁡,iR_{\max,i}, filled uniformly with matter of density ρi\rho_i. A sphere is the simplest shape that keeps the model the same in every direction.
  • Give it a Hubble flow: every particle moves away from the centre with velocity proportional to its distance.
Initial Hubble flow
v⃗i=Hi r⃗\vec v_i = H_i\,\vec r
In words

vector v sub i equals H sub i vector r

Source: 8.286 Lecture 6 slides, slide 8
  • Treat the matter as "dust": particles with no pressure, moving slowly compared with light. This describes the real universe from about 50,000 years after the Big Bang (when matter took over from radiation) until a few billion years ago (when dark energy took over).

Following the shells#

Divide the sphere into thin spherical shells. Label each by its initial radius rir_i, and let r(ri,t)r(r_i, t) be its radius at time tt. All velocities and forces point radially, so each shell stays a shell.

By the shell theorem, a shell feels only the mass inside it, and as long as shells don't cross, that mass never changes:

Mass inside a shell
M(ri)=4π3 ri3 ρiM(r_i) = \frac{4\pi}{3}\,r_i^3\,\rho_i
In words

M(r sub i ) equals 4 pi over 3 r sub i cubed rho sub i

Initially every shell moves outward faster than the ones inside it, so none are about to cross. Guth's argument: the equations, valid up to any first crossing, would have to predict that crossing — and they never do. So there is none.

Newton's second law for the shell is then

Equation of motion of a shell
r¨=−G M(ri)r2=−4π3 G ri3 ρir2,r(ri,ti)=ri,r˙(ri,ti)=Hi ri\ddot r = -\frac{G\,M(r_i)}{r^2} = -\frac{4\pi}{3}\,\frac{G\,r_i^3\,\rho_i}{r^2}, \qquad r(r_i, t_i) = r_i,\quad \dot r(r_i, t_i) = H_i\,r_i
In words

r double dot equals minus G M(r sub i ) over r squared equals minus 4 pi over 3 G r sub i cubed rho sub i over r squared , r(r sub i , t sub i ) equals r sub i , r dot (r sub i , t sub i ) equals H sub i r sub i

The cancellation#

Now the key step. Measure each shell's radius relative to where it started, u(ri,t)≡r(ri,t)/riu(r_i, t) \equiv r(r_i, t)/r_i. Dividing the equation by rir_i and writing r=u rir = u\,r_i:

Rescaled equation
u¨=−4π3 G ri3 ρiri u2 ri2=−4π3 Gρiu2,u(ti)=1,u˙(ti)=Hi\ddot u = -\frac{4\pi}{3}\,\frac{G\,r_i^3\,\rho_i}{r_i\,u^2\,r_i^2} = -\frac{4\pi}{3}\,\frac{G\rho_i}{u^2}, \qquad u(t_i) = 1,\quad \dot u(t_i) = H_i
In words

u double dot equals minus 4 pi over 3 G r sub i cubed rho sub i over r sub i u squared r sub i squared equals minus 4 pi over 3 G rho sub i over u squared , u(t sub i ) equals 1, u dot (t sub i ) equals H sub i

The ri3r_i^3 from the mass has cancelled against the ri3r_i^3 from rescaling the force. Neither the equation nor its starting values depend on rir_i any more, so every shell follows the same function of time. That function is the scale factor, u=a(t)u = a(t), and

Uniform expansion stays uniform
r(ri,t)=a(t) rir(r_i, t) = a(t)\,r_i
In words

r(r sub i , t) equals a(t) r sub i

Guth emphasises that this works because gravity is exactly inverse-square: with a force law of most other forms, the cancellation would fail and the expansion would not stay uniform. (A force growing in proportion to distance would also preserve uniformity — that is how a cosmological constant behaves, a later topic.) The radius Rmax⁡,iR_{\max,i} never appeared, so letting the sphere become infinite changes nothing.

Density and deceleration#

The mass inside a shell is fixed while its volume grows as a3a^3, so the density falls as

Matter density
ρ(t)=ρia3(t)\rho(t) = \frac{\rho_i}{a^3(t)}
In words

rho (t) equals rho sub i over a cubed (t)

Putting this back into the rescaled equation gives the equation of motion of the universe (for matter only):

Deceleration of the scale factor
a¨=−4π3 G ρ(t) a\ddot a = -\frac{4\pi}{3}\,G\,\rho(t)\,a
In words

a double dot equals minus 4 pi over 3 G rho (t) a

The physically meaningful combination is a¨/a\ddot a/a: the deceleration rate depends only on the density at that moment. (In Guth's convention, set up formally in Lecture 7, aa carries units of metres per notch, and these units cancel only in a¨/a\ddot a/a.)

Worked example

How hard is gravity braking the expansion?

Ignore dark energy and take today's matter density as ρ=2.7×10−27 kg/m3\rho = 2.7\times10^{-27}\ \text{kg/m}^3 and H=70H = 70 km/s/Mpc =2.27×10−18 s−1= 2.27\times10^{-18}\ \text{s}^{-1}. Compute a¨/a\ddot a/a and compare it with H2=(a˙/a)2H^2 = (\dot a/a)^2.

Deceleration rate

a¨/a=−4π3Gρ=−4.19×6.67×10−11×2.7×10−27=−7.5×10−37 s−2\ddot a/a = -\tfrac{4\pi}{3}G\rho = -4.19 \times 6.67\times10^{-11} \times 2.7\times10^{-27} = -7.5\times10^{-37}\ \text{s}^{-2}.

Expansion rate squared

H2=(2.27×10−18)2=5.1×10−36 s−2H^2 = (2.27\times10^{-18})^2 = 5.1\times10^{-36}\ \text{s}^{-2}.

Compare

The ratio is about 0.150.15: in one Hubble time, gravity from matter alone would slow the recession speed a˙\dot a of any given galaxy by roughly 15%. (The Hubble rate H=a˙/aH = \dot a/a falls faster, because aa also grows.) In reality dark energy more than cancels this braking today — see the acceleration page and the later lectures on the cosmological constant.

A conserved quantity#

The equation a¨=−4π3Gρi/a2\ddot a = -\tfrac{4\pi}{3}G\rho_i/a^2 can be integrated once. Multiply by a˙\dot a and notice that the result is a total time derivative:

Conserved quantity E
E=12 a˙2−4π3 Gρia,dEdt=a˙[a¨+4π3Gρia2]=0E = \tfrac12\,\dot a^2 - \frac{4\pi}{3}\,\frac{G\rho_i}{a}, \qquad \frac{dE}{dt} = \dot a\left[\ddot a + \frac{4\pi}{3}\frac{G\rho_i}{a^2}\right] = 0
In words

E equals one half a dot squared minus 4 pi over 3 G rho sub i over a , dE over dt equals a dot [ a double dot plus 4 pi over 3 G rho sub i over a squared ] equals 0

Multiplied by m ri2m\,r_i^2, EE is the kinetic plus potential energy of a test particle of mass mm on the outer edge of the sphere. (For a particle inside the sphere this is not its true potential energy, because the outer shells also contribute; Guth notes the interpretation is cleanest at the edge.) So, as Lecture 7 develops, the universe behaves like a ball thrown upward:

  • E>0E > 0: the expansion slows but never stops.
  • E<0E < 0: the expansion stops at a maximum size and the universe recollapses.
  • E=0E = 0: the borderline case, just barely expanding forever.
The speed of expansion as the universe grows, in units where 8πGρᵢ/3 = 1, computed from the conserved E. Negative E: the curve hits zero at a = −1/(2E), the largest the universe gets.

The next lecture works out these three behaviours in time and names them closed, open and flat; what those names mean for the geometry of space comes later in the course.

Summary#

  • A uniformly expanding sphere of dust, followed shell by shell with Newton's gravity, stays uniform: r=a(t) rir = a(t)\,r_i.
  • This works only because gravity is exactly inverse-square; the result does not depend on the sphere's size.
  • Matter density falls as ρ∝1/a3\rho \propto 1/a^3, and the scale factor obeys a¨=−4π3Gρ a\ddot a = -\tfrac{4\pi}{3}G\rho\,a.
  • E=12a˙2−4π3Gρi/aE = \tfrac12\dot a^2 - \tfrac{4\pi}{3}G\rho_i/a is conserved; its sign decides between eternal expansion and recollapse.

Symbols on this page

Show the 13 symbols used here
π\pi pi
The ratio of a circle's circumference to its diameter, π ≈ 3.14159. Say it: “pi”.
ρ\rho rho
Density: mass (or energy divided by c²) per unit volume. ρ_c is the critical density, ρ̄ the mean density. Say it: “rho”.
∝\propto proportional to
Is proportional to: doubling one doubles the other. Say it: “is proportional to”.
≡\equiv identical to
Is defined as (or is identically equal to). Say it: “is defined as”.
×\times times
Times (multiplication). Between two vectors: the cross product. Say it: “times”.
max⁡\max maximum
Maximum: the largest value. Say it: “max”.
v⃗\vec v vector arrow
A vector: a quantity with a size and a direction. Say it: “vector”.
a˙\dot a dot
Rate of change with time: ȧ = da/dt. Say it: “dot”.
a¨\ddot a double dot
Second rate of change with time: ä = d²a/dt². Say it: “double dot”.
aa a
The scale factor a(t): how much the universe has expanded. Physical distance = a × comoving distance (Guth's units: metres per notch). Say it: “a”.
HH H
The Hubble parameter H = ȧ/a: the expansion rate of the universe at a given time. Say it: “H”.
GG G
Newton's gravitational constant, G ≈ 6.674 × 10⁻¹¹ m³ kg⁻¹ s⁻². Say it: “G”.
ρi\rho_{i} initial density
The density when a = 1, so that ρ = ρᵢ/a³ for matter. Say it: “rho i”.

Tap or click a symbol in a formula to see the same explanation. All symbols: symbols and notation.

Where this leads

Pages that build on this

Check your understanding

  1. If the scale factor doubles, what happens to the density of matter?
  2. Why does every shell of the sphere follow the same function a(t)a(t)?
  3. Which feature of gravity makes the rir_i cancel?
  4. E=12a˙2−4π3Gρi/aE = \tfrac12\dot a^2 - \tfrac{4\pi}{3}G\rho_i/a is constant. If E<0E < 0, the expansion…
0 of 4 answered

Flashcards

1 / 4