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

Simplifying the geodesic equation

Measure the path by its own length and the geodesic equation sheds its square roots; solving for the acceleration brings in the inverse metric and the Christoffel symbols.

You need this first

Full reading path (23 pages)
  1. Vectors (8.01)
  2. Position and velocity in one dimension (8.01)
  3. The Doppler effect and special relativity (8.286)
  4. Acceleration in one dimension (8.01)
  5. Motion in two dimensions and projectiles (8.01)
  6. Newton's laws of motion (8.01)
  7. The universal law of gravitation (8.01)
  8. The homogeneous expanding universe (8.286)
  9. Why a static universe cannot work (8.286)
  10. The expanding sphere and the equation for a(t) (8.286)
  11. The Friedmann equation (8.286)
  12. Open, closed and flat universes (8.286)
  13. The closed universe as a cycloid (8.286)
  14. The age of the universe and the Big Bang singularity (8.286)
  15. The age of a closed universe from H and Ω (8.286)
  16. The open universe and the age of any matter universe (8.286)
  17. Non-Euclidean geometry (8.286)
  18. The metric of a sphere and of a closed universe (8.286)
  19. Curvature and expansion — the Robertson–Walker metric (8.286)
  20. The geometry of open universes (8.286)
  21. The spacetime interval (8.286)
  22. The spacetime metric of the universe (8.286)
  23. The geodesic equation (8.286)
On this page · 6 sections
  1. Recap: what ds² means, and why geodesics
  2. Choose the length itself as the parameter
  3. Solving for the second derivative
  4. The Christoffel symbols
  5. Spacetime notation
  6. Summary

Lecture 14 finishes the geodesic equation. The version derived in Lecture 13 is correct but clumsy: the square root A\sqrt A sits in two denominators, and on the left it is even differentiated. One sensible choice of parameter removes it, and a little matrix algebra turns the result into the form found in most relativity books.

Recap: what ds² means, and why geodesics#

Guth's review of Lecture 13 (slide 1) adds one word: the observers who give ds2ds^2 its meaning are local free-falling observers. The equivalence principle cancels gravity only in a small region, so the special-relativity measurements behind ds2ds^2 are made by an observer arbitrarily close to the events. A special case of "ds2<0ds^2 < 0" is someone reading their own wristwatch, which is always at the same place relative to them.

Asked whether time always enters as just −c2dt2-c^2dt^2, Guth says no: that relied on isotropy. To describe the real, lumpy universe rather than the ideal one, you would probably need dt drdt\,dr terms (transcript p. 3).

Why care about geodesics? Because according to general relativity, particles with no forces on them except gravity move along geodesics of spacetime. Applying this to the Robertson–Walker metric is an 8.286 homework problem; later in the lecture (p. 13) Guth says it shows how massive particles slow down as they travel through the expanding universe, just as photons lose energy.

Guth adds some history (p. 5). The method of varying a path, the calculus of variations, dates from 1696, when Johann Bernoulli posed the brachistochrone problem: the curve along which a frictionless bead slides between two points in the least time. The answer is a cycloid, the same curve as the closed matter-dominated universe. Newton is said to have solved it overnight after a day at the Royal Mint. [Guth has him coming home at 4 a.m.; the usual version of the story has him home at 4 p.m. and done by 4 a.m.]

Choose the length itself as the parameter#

The parameter λ\lambda was arbitrary, and the equation holds for every choice. So pick the most convenient one: let λ\lambda be the arc length ss along the path. Since ds2=A dλ2ds^2 = A\,d\lambda^2, this makes A=1A = 1. Guth stresses that what really matters is that AA is constant: then it is not differentiated, and it cancels between the two sides. The square roots disappear:

Geodesic equation, parametrised by arc length
dds[gij dxjds]=12 ∂gjk∂xi dxjdsdxkds\frac{d}{ds}\left[g_{ij}\,\frac{dx^j}{ds}\right] = \frac12\,\frac{\partial g_{jk}}{\partial x^i}\,\frac{dx^j}{ds}\frac{dx^k}{ds}
In words

d over ds [g sub ij dx to the power j over ds ] equals one half partial g sub jk over partial x to the power i dx to the power j over ds dx to the power k over ds

Source: 8.286 Lecture 14 transcript, p. 9 (board)

Guth boxes this: it is the final geodesic equation, and for solving problems it is almost always the best form.

Solving for the second derivative#

Most books instead give a formula for d2xi/ds2d^2x^i/ds^2. To get it, expand the left side with the product rule. The metric depends on ss only through the position xk(s)x^k(s), so the chain rule gives ddsgij=∂kgij dxkds\frac{d}{ds}g_{ij} = \partial_k g_{ij}\,\frac{dx^k}{ds}, where ∂k\partial_k is shorthand for ∂/∂xk\partial/\partial x^k. Moving that term to the right:

Expanded form
gij d2xjds2=(12 ∂igjk−∂kgij)dxjdsdxkdsg_{ij}\,\frac{d^2x^j}{ds^2} = \left(\frac12\,\partial_i g_{jk} - \partial_k g_{ij}\right)\frac{dx^j}{ds}\frac{dx^k}{ds}
In words

g sub ij d squared x to the power j over ds squared equals ( one half partial sub i g sub jk minus partial sub k g sub ij ) dx to the power j over ds dx to the power k over ds

Source: 8.286 Lecture 14 transcript, pp. 10–11 (board)

(On the board this took several attempts — students caught a misplaced bracket and two wrong indices before it came out right.)

The left side is a matrix gijg_{ij} times a vector d2xj/ds2d^2x^j/ds^2. To isolate the vector, multiply by the inverse matrix. In general relativity the inverse metric is written with upper indices, gijg^{ij}:

The inverse metric
gil glj=δjig^{il}\,g_{lj} = \delta^i_j
In words

g to the power il g sub lj equals delta to the power i sub j

Source: 8.286 Lecture 14 transcript, p. 11 (board)

Summing over the neighbouring index ll is exactly matrix multiplication, and δji\delta^i_j (1 if i=ji = j, 0 otherwise) is the identity matrix.

The Christoffel symbols#

Multiplying through by gilg^{il} and renaming dummy indices gives the textbook form. The coefficients of the two first derivatives get their own symbol:

Geodesic equation, textbook form
d2xids2=− Γjki dxjdsdxkds,Γjki=12 gil(∂jglk+∂kglj−∂lgjk)\frac{d^2x^i}{ds^2} = -\,\Gamma^i_{jk}\,\frac{dx^j}{ds}\frac{dx^k}{ds}, \qquad \Gamma^i_{jk} = \frac12\,g^{il}\left(\partial_j g_{lk} + \partial_k g_{lj} - \partial_l g_{jk}\right)
In words

d squared x to the power i over ds squared equals minus gamma to the power i sub jk dx to the power j over ds dx to the power k over ds , gamma to the power i sub jk equals one half g to the power il ( partial sub j g sub lk plus partial sub k g sub lj minus partial sub l g sub jk )

Source: 8.286 Lecture 14 transcript, p. 12 (board)

Γjki\Gamma^i_{jk} is the affine connection (Weinberg's name), also called the Christoffel symbol, and — according to Sean Carroll's book, which Guth quotes — sometimes the Riemann or Levi-Civita connection. Lots of names, Guth jokes, presumably because lots of people invented it. Everyone agrees on the definition up to a sign, and books differ on that sign. [The transcript does not make the sign on the board audible; with Γ\Gamma defined as above, the derivation gives the minus sign shown, which is the most common convention.]

For Guth this form is "window dressing" to connect with other books; in practice the arc-length form above is usually easier to solve.

Worked example

Which paths on a sphere are geodesics?

On a sphere of unit radius, ds2=dθ2+sin⁡2θ dϕ2ds^2 = d\theta^2 + \sin^2\theta\,d\phi^2. Write down the geodesic equations, and check two paths: the equator, and the circle θ=60∘\theta = 60^\circ, each travelled at unit speed. (θ\theta is the polar angle, measured from the north pole, so θ=60∘\theta = 60^\circ is the circle of latitude 30∘30^\circ.)

The θ equation

Take i=θi = \theta. Only gθθ=1g_{\theta\theta} = 1 contributes on the left; on the right only gϕϕ=sin⁡2θg_{\phi\phi} = \sin^2\theta depends on θ\theta, with ∂θgϕϕ=2sin⁡θcos⁡θ\partial_\theta g_{\phi\phi} = 2\sin\theta\cos\theta. So d2θds2=sin⁡θcos⁡θ(dϕds)2\dfrac{d^2\theta}{ds^2} = \sin\theta\cos\theta\left(\dfrac{d\phi}{ds}\right)^2.

The φ equation

Take i=ϕi = \phi. No metric component depends on ϕ\phi, so the right side vanishes: dds[sin⁡2θ dϕds]=0\dfrac{d}{ds}\left[\sin^2\theta\,\dfrac{d\phi}{ds}\right] = 0. The quantity in brackets is conserved along a geodesic.

Christoffel symbols

From the definition, the non-zero ones are Γϕϕθ=−sin⁡θcos⁡θ\Gamma^\theta_{\phi\phi} = -\sin\theta\cos\theta and Γθϕϕ=Γϕθϕ=cot⁡θ\Gamma^\phi_{\theta\phi} = \Gamma^\phi_{\phi\theta} = \cot\theta. The textbook form then gives the same θ\theta equation: d2θ/ds2=−Γϕϕθ(dϕ/ds)2d^2\theta/ds^2 = -\Gamma^\theta_{\phi\phi}(d\phi/ds)^2.

Testing the two circles

On a circle θ=θ0\theta = \theta_0 at unit speed, dϕ/ds=1/sin⁡θ0d\phi/ds = 1/\sin\theta_0 and d2θ/ds2=0d^2\theta/ds^2 = 0. The θ\theta equation demands sin⁡θ0cos⁡θ0/sin⁡2θ0=cot⁡θ0\sin\theta_0\cos\theta_0/\sin^2\theta_0 = \cot\theta_0. At the equator, cot⁡90∘=0\cot 90^\circ = 0: a geodesic. At θ0=60∘\theta_0 = 60^\circ, cot⁡60∘≈0.577≠0\cot 60^\circ \approx 0.577 \neq 0: not a geodesic. The φ equation holds for both, since sin⁡2θ0⋅(1/sin⁡θ0)\sin^2\theta_0 \cdot (1/\sin\theta_0) is constant.

The graph shows what the latitude circle loses. Between two points on the same latitude, the geodesic — a great circle — bulges towards the pole and is shorter.

Both paths join two points at the same latitude β₀, at longitudes ±D (both sliders in degrees). The great circle is computed from tan(latitude) = tan(latitude of ends) · cos φ / cos D. For ends at 45° latitude and D = 60°, the great circle peaks at 63.4° latitude and has length 1.318 radii, against 1.481 radii along the circle of latitude.

Spacetime notation#

A particle moves along a timelike path. So instead of ss it is natural to use proper time τ\tau, defined by ds2=−c2dτ2ds^2 = -c^2d\tau^2, and — a near-universal convention — Greek indices μ,ν,λ,σ\mu, \nu, \lambda, \sigma for spacetime (with 00 for time), keeping Latin i,j,ki, j, k for space alone. The geodesic equation becomes

Geodesic equation in spacetime
ddτ[gμν dxνdτ]=12 ∂gλσ∂xμ dxλdτdxσdτ\frac{d}{d\tau}\left[g_{\mu\nu}\,\frac{dx^\nu}{d\tau}\right] = \frac12\,\frac{\partial g_{\lambda\sigma}}{\partial x^\mu}\,\frac{dx^\lambda}{d\tau}\frac{dx^\sigma}{d\tau}
In words

d over d tau [g sub mu nu dx superscript nu over d tau ] equals one half partial g sub lambda sigma over partial x superscript mu dx superscript lambda over d tau dx superscript sigma over d tau

Source: 8.286 Lecture 14 transcript, p. 17 (board)

[Guth says "the partial of gλσg_{\lambda\sigma} with respect to xνx^\nu"; the free index must be μ\mu, as written here.]

Nothing in the derivation used a positive-definite metric, so it goes through unchanged. The difference is that spacetime geodesics can be maxima or other stationary points as well as minima.

Summary#

  • Choosing the arc length as parameter makes A=1A = 1 (constant is what matters), giving dds[gijdxjds]=12∂igjkdxjdsdxkds\frac{d}{ds}\left[g_{ij}\frac{dx^j}{ds}\right] = \frac12\partial_i g_{jk}\frac{dx^j}{ds}\frac{dx^k}{ds}, usually the best form for solving problems.
  • The inverse metric gijg^{ij} satisfies gilglj=δjig^{il}g_{lj} = \delta^i_j; with it, d2xids2=−Γjkidxjdsdxkds\frac{d^2x^i}{ds^2} = -\Gamma^i_{jk}\frac{dx^j}{ds}\frac{dx^k}{ds}, where Γjki=12gil(∂jglk+∂kglj−∂lgjk)\Gamma^i_{jk} = \frac12 g^{il}(\partial_j g_{lk} + \partial_k g_{lj} - \partial_l g_{jk}) is the Christoffel symbol (sign conventions vary).
  • For particles, use proper time τ\tau and Greek indices: ddτ[gμνdxνdτ]=12∂μgλσdxλdτdxσdτ\frac{d}{d\tau}\left[g_{\mu\nu}\frac{dx^\nu}{d\tau}\right] = \frac12\partial_\mu g_{\lambda\sigma}\frac{dx^\lambda}{d\tau}\frac{dx^\sigma}{d\tau}.

Symbols on this page

Show the 21 symbols used here
β\beta beta
Latitude of the two end points (in the plot). Say it: “beta” (sounds like “bay-ta or bee-ta”).
δ\delta delta
The Kronecker delta: δ^i_j is 1 when i = j and 0 otherwise. Say it: “delta”.
θ\theta theta
Polar angle: measured from the north pole (0) to the south pole (π). Say it: “theta” (sounds like “thay-ta or thee-ta”).
λ\lambda lambda
The parameter along a path; in the spacetime formulas λ is also an index running over 0–3. Say it: “lambda”.
μ\mu mu
A spacetime index, running over 0 (time) and 1, 2, 3 (space). Say it: “mu” (sounds like “myoo”).
ν\nu nu
A spacetime index, running over 0–3. Say it: “nu” (sounds like “new”).
σ\sigma sigma
A spacetime index, running over 0–3. Say it: “sigma”.
τ\tau tau
Proper time: the time measured by a clock moving with the object. Say it: “tau” (sounds like “rhymes with cow”).
ϕ\phi phi
Azimuthal angle (longitude), running from 0 to 2π around the axis. Say it: “phi” (sounds like “fye”).
Γ\Gamma capital gamma
A Christoffel symbol, Γ^i_jk, built from derivatives of the metric. Say it: “gamma”.
∘\circ degree
Degrees of angle (as in 45°). Say it: “degrees”.
≈\approx approximately equal
Is approximately equal to. Say it: “is approximately”.
≠\neq not equal
Is not equal to. Say it: “is not equal to”.
⋅\cdot dot
Times (multiplication). Between two vectors: the dot product. Say it: “times”.
∂\partial partial
Partial derivative: the rate of change with one variable, holding the others fixed. Say it: “partial”.
sin⁡\sin sine
Sine of an angle. Say it: “sine”.
cos⁡\cos cosine
Cosine of an angle. Say it: “cos”.
cot⁡\cot cotangent
Cotangent of an angle, cos/sin. Say it: “cot”.
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”.
cc c
The speed of light, c ≈ 2.998 × 10⁸ m/s. Say it: “c”.
kk k
An index running over the coordinates. Say it: “k”.

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. Why does choosing λ=s\lambda = s, the arc length, simplify the geodesic equation?
  2. What is gijg^{ij}?
  3. On a unit sphere (ds2=dθ2+sin⁡2θ dϕ2ds^2 = d\theta^2 + \sin^2\theta\,d\phi^2), a path runs along the circle θ=60∘\theta = 60^\circ (latitude 30∘30^\circ, since θ\theta is measured from the pole) at unit speed. Is it a geodesic?
  4. Why do relativists switch to τ\tau and Greek indices for particle paths?
0 of 4 answered

Flashcards

1 / 4