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

The geodesic equation

How to find the paths of stationary length in any metric — index notation, a varied path, one integration by parts, and the equation that comes out.

You need this first

Full reading path (22 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)
On this page · 8 sections
  1. Geodesics: stationary paths
  2. A better notation
  3. The length of a path
  4. Varying the path
  5. Differentiating the length
  6. The key step: integration by parts
  7. The geodesic equation
  8. Summary

With the metric of the universe in hand, the second half of Lecture 13 asks how things move in a given metric. The tool is the geodesic: the path of stationary length. (That free particles move on geodesics is taken up in Lecture 14.) The method is general: it works for any metric, including the Schwarzschild metric of stars and black holes, which Guth mentions as a later example. In the lecture he showed the long equations from his lecture notes on screen rather than writing them on the board; the same equations appear on slides 2–5 of Lecture 14, which are cited here.

Geodesics: stationary paths#

Guth starts with the simplest case: a two-dimensional space with coordinates xx and yy and some metric, the kind of space Gauss, Bolyai and Lobachevsky studied. Everything carries over to spacetime unchanged.

A geodesic between two points is a path whose length is stationary: change the path slightly, keeping the endpoints fixed, and the length does not change to first order. In an ordinary space a geodesic is usually the shortest path. In spacetime, where the metric is not positive definite, geodesics can also be maxima or saddle points. The derivation only uses "the first-order change vanishes", so it covers all of these at once.

[Guth suggests that in two-dimensional spaces the stationary paths are always minima. Not quite: on a sphere the long way round a great circle between two points is also stationary, but it is not the shortest path: it is a saddle point.]

A better notation#

Two improvements make the metric easier to handle.

Number the coordinates. Write x1,x2x^1, x^2 instead of x,yx, y. The superscripts are labels, not powers — an unfortunate clash, tolerated because general relativity needs both upper and lower indices.

Sum with indices. The metric becomes ∑i=12∑j=12gij dxidxj\sum_{i=1}^{2}\sum_{j=1}^{2} g_{ij}\,dx^i dx^j. The (1,1)(1,1) and (2,2)(2,2) terms are the old dx2dx^2 and dy2dy^2 terms; g12g_{12} and g21g_{21} are equal and both multiply dx dydx\,dy. Then comes Einstein's summation convention: an index that appears twice, once up and once down, is summed automatically, so the sum signs can be dropped.

Metric in index notation
ds2=gij(xk) dxi dxjds^2 = g_{ij}(x^k)\,dx^i\,dx^j
In words

ds squared equals g sub ij (x to the power k ) dx to the power i dx to the power j

In four dimensions this one line stands for 16 terms (10 different ones, since gij=gjig_{ij} = g_{ji}).

The length of a path#

Describe a path from AA to BB by two functions xi(λ)x^i(\lambda) of a parameter λ\lambda, which runs from 00 to λf\lambda_f:

A path from A to B
xi(λ),xi(0)=xAi,xi(λf)=xBix^i(\lambda), \qquad x^i(0) = x^i_A, \qquad x^i(\lambda_f) = x^i_B
In words

x to the power i ( lambda ), x to the power i (0) equals x to the power i sub A , x to the power i ( lambda sub f ) equals x to the power i sub B

Cut the path into small pieces from λ\lambda to λ+dλ\lambda + d\lambda. Over one piece the coordinates change by dxi=(dxi/dλ) dλdx^i = (dx^i/d\lambda)\,d\lambda; putting this into the metric gives

One small piece
ds2=gij(xk(λ)) dxidλdxjdλ dλ2ds^2 = g_{ij}\big(x^k(\lambda)\big)\,\frac{dx^i}{d\lambda}\frac{dx^j}{d\lambda}\,d\lambda^2
In words

ds squared equals g sub ij (x to the power k ( lambda ) ) dx to the power i over d lambda dx to the power j over d lambda d lambda squared

For a space, ds2ds^2 is positive, so the square root is the length of the piece. Adding up all pieces:

Length of the path
S[xi(λ)]=∫0λfgij(xk(λ)) dxidλdxjdλ  dλS\big[x^i(\lambda)\big] = \int_0^{\lambda_f} \sqrt{g_{ij}\big(x^k(\lambda)\big)\,\frac{dx^i}{d\lambda}\frac{dx^j}{d\lambda}}\;d\lambda
In words

S [x to the power i ( lambda ) ] equals the integral from 0 to lambda sub f of the square root of g sub ij (x to the power k ( lambda ) ) dx to the power i over d lambda dx to the power j over d lambda d lambda

Varying the path#

To test whether a path is a geodesic, compare it with every nearby path between the same endpoints. This is the calculus of variations, which slide 4 dates to 1696. Write a nearby path as

A varied path
x~i(λ)=xi(λ)+α wi(λ),wi(0)=0,wi(λf)=0\tilde x^i(\lambda) = x^i(\lambda) + \alpha\,w^i(\lambda), \qquad w^i(0) = 0, \quad w^i(\lambda_f) = 0
In words

x tilde to the power i ( lambda ) equals x to the power i ( lambda ) plus alpha w to the power i ( lambda ), w to the power i (0) equals 0, w to the power i ( lambda sub f ) equals 0

Here wi(λ)w^i(\lambda) is any ordinary function (with values like 2 or 5) and α\alpha is a small number that controls how far the new path strays. The conditions on wiw^i keep the endpoints fixed — otherwise you could always find a shorter path by moving them. The path xi(λ)x^i(\lambda) is a geodesic if

Condition for a geodesic
dS[x~i(λ)]dα∣α=0=0for all wi(λ)\left.\frac{d S\big[\tilde x^i(\lambda)\big]}{d\alpha}\right|_{\alpha = 0} = 0 \quad \text{for all } w^i(\lambda)
In words

d S [ x tilde to the power i ( lambda ) ] over d alpha evaluated at alpha equals 0 equals 0 for all w to the power i ( lambda )

Paths from (0, 0) to (1, 1) in the flat plane, varied with w = x(1 − x), which vanishes at both ends. The length, computed exactly from the integral above, is S(α) = [F(1 + α) − F(1 − α)]/(2α) with F(u) = ½(u√(1 + u²) + asinh u). It is flat at α = 0, the straight line (√2 ≈ 1.414), and grows on both sides: 1.429 at α = ½, 1.479 at α = ±1, 1.700 at α = ±2.

Differentiating the length#

Name the integrand of the varied path's length:

The integrand A
S[x~i(λ)]=∫0λfA(λ,α)  dλ,A(λ,α)=gij(x~k(λ)) dx~idλdx~jdλS\big[\tilde x^i(\lambda)\big] = \int_0^{\lambda_f} \sqrt{A(\lambda, \alpha)}\;d\lambda, \qquad A(\lambda, \alpha) = g_{ij}\big(\tilde x^k(\lambda)\big)\,\frac{d\tilde x^i}{d\lambda}\frac{d\tilde x^j}{d\lambda}
In words

S [ x tilde to the power i ( lambda ) ] equals the integral from 0 to lambda sub f of the square root of A( lambda , alpha ) d lambda , A( lambda , alpha ) equals g sub ij ( x tilde to the power k ( lambda ) ) d x tilde to the power i over d lambda d x tilde to the power j over d lambda

α\alpha appears in three places. The metric does not depend on α\alpha directly, but it is evaluated at x~\tilde x, which does; by the chain rule

Chain rule for the metric
∂gij∂α=∂gij∂xk ∂x~k∂α=∂gij∂xk wk\frac{\partial g_{ij}}{\partial\alpha} = \frac{\partial g_{ij}}{\partial x^k}\,\frac{\partial\tilde x^k}{\partial\alpha} = \frac{\partial g_{ij}}{\partial x^k}\,w^k
In words

partial g sub ij over partial alpha equals partial g sub ij over partial x to the power k partial x tilde to the power k over partial alpha equals partial g sub ij over partial x to the power k w to the power k

Source: 8.286 Lecture 13 transcript, pp. 23–24 (on screen from the lecture notes)

and each dx~i/dλd\tilde x^i/d\lambda gives dwi/dλdw^i/d\lambda. The derivative of a square root brings 1/(2A)1/(2\sqrt A). After differentiating, everything else is set to α=0\alpha = 0, so the tildes disappear. The two terms with dw/dλdw/d\lambda are equal once the summed ("dummy") indices ii and jj are renamed, using gij=gjig_{ij} = g_{ji}, which gives a factor 2:

Derivative of the length
dSdα∣α=0=12∫0λf1A(λ,0)[∂gij∂xk wk dxidλdxjdλ+2 gij dwidλdxjdλ]dλ\left.\frac{dS}{d\alpha}\right|_{\alpha = 0} = \frac12\int_0^{\lambda_f}\frac{1}{\sqrt{A(\lambda, 0)}}\left[\frac{\partial g_{ij}}{\partial x^k}\,w^k\,\frac{dx^i}{d\lambda}\frac{dx^j}{d\lambda} + 2\,g_{ij}\,\frac{dw^i}{d\lambda}\frac{dx^j}{d\lambda}\right]d\lambda
In words

dS over d alpha evaluated at alpha equals 0 equals one half the integral from 0 to lambda sub f of 1 over the square root of A( lambda , 0) [ partial g sub ij over partial x to the power k w to the power k dx to the power i over d lambda dx to the power j over d lambda plus 2 g sub ij dw to the power i over d lambda dx to the power j over d lambda ]d lambda

The key step: integration by parts#

The first term is proportional to ww, the second to its derivative, so it is not yet clear when the sum vanishes for every ww. Guth calls this the only real difficulty in the derivation; the fix is to integrate the second term by parts, ∫u dv=−∫v du+[uv]\int u\,dv = -\int v\,du + [uv], with

  • u=1A gij dxjdλu = \dfrac{1}{\sqrt A}\,g_{ij}\,\dfrac{dx^j}{d\lambda} and dv=dwidλ dλdv = \dfrac{dw^i}{d\lambda}\,d\lambda, so v=wiv = w^i.

The boundary term [uv][uv] contains wiw^i at the two endpoints, where it vanishes: no boundary term survives. Now ww appears undifferentiated in both terms. Renaming dummy indices so it carries the same index ii in each, it can be taken outside:

After integrating by parts
dSdα∣α=0=∫0λf{12A ∂gjk∂xidxjdλdxkdλ−ddλ[1A gij dxjdλ]}wi(λ) dλ\left.\frac{dS}{d\alpha}\right|_{\alpha = 0} = \int_0^{\lambda_f}\left\{\frac{1}{2\sqrt A}\,\frac{\partial g_{jk}}{\partial x^i}\frac{dx^j}{d\lambda}\frac{dx^k}{d\lambda} - \frac{d}{d\lambda}\left[\frac{1}{\sqrt A}\,g_{ij}\,\frac{dx^j}{d\lambda}\right]\right\} w^i(\lambda)\,d\lambda
In words

dS over d alpha evaluated at alpha equals 0 equals the integral from 0 to lambda sub f of 1 over 2 the square root of A partial g sub jk over partial x to the power i dx to the power j over d lambda dx to the power k over d lambda minus d over d lambda [ 1 over the square root of A g sub ij dx to the power j over d lambda ] w to the power i ( lambda ) d lambda

The geodesic equation#

This must vanish for every wiw^i. One integral can be zero by luck — positive in some places, negative in others. But if the expression in curly brackets were not zero, we could choose wiw^i equal to it; the integrand would then be a perfect square and the integral could not vanish. So the brackets must vanish identically:

The geodesic equation
ddλ[1A gij dxjdλ]=12A ∂gjk∂xi dxjdλdxkdλ\frac{d}{d\lambda}\left[\frac{1}{\sqrt A}\,g_{ij}\,\frac{dx^j}{d\lambda}\right] = \frac{1}{2\sqrt A}\,\frac{\partial g_{jk}}{\partial x^i}\,\frac{dx^j}{d\lambda}\frac{dx^k}{d\lambda}
In words

d over d lambda [ 1 over the square root of A g sub ij dx to the power j over d lambda ] equals 1 over 2 the square root of A partial g sub jk over partial x to the power i dx to the power j over d lambda dx to the power k over d lambda

There is one equation for each value of ii. Both sides depend on the metric and on the path, through dx/dλdx/d\lambda. A continuous path is a geodesic if and only if it satisfies this equation. Lecture 14 simplifies it and tries it on examples.

Worked example

Straight lines in the flat plane

For the flat plane in Cartesian coordinates, g11=g22=1g_{11} = g_{22} = 1 and g12=g21=0g_{12} = g_{21} = 0. What does the geodesic equation say? Then compare the straight path from (0,0)(0,0) to (1,1)(1,1) with the parabola x=λx = \lambda, y=λ2y = \lambda^2, 0≤λ≤10 \le \lambda \le 1.

Right-hand side

The metric is constant, so ∂gjk/∂xi=0\partial g_{jk}/\partial x^i = 0 and the right-hand side vanishes.

Left-hand side

gij dxj/dλ=dxi/dλg_{ij}\,dx^j/d\lambda = dx^i/d\lambda, so the equation says ddλ[1Adxidλ]=0\dfrac{d}{d\lambda}\left[\dfrac{1}{\sqrt A}\dfrac{dx^i}{d\lambda}\right] = 0. Since A dλ=ds\sqrt A\,d\lambda = ds, this is d2xi/ds2=0d^2x^i/ds^2 = 0: the coordinates change at a constant rate with distance. The geodesics are straight lines, as they should be.

The parabola is not one

Along the parabola A=1+4λ2A = 1 + 4\lambda^2. For i=1i = 1: ddλ[(1+4λ2)−1/2]=−4λ (1+4λ2)−3/2\dfrac{d}{d\lambda}\left[(1+4\lambda^2)^{-1/2}\right] = -4\lambda\,(1+4\lambda^2)^{-3/2}, which is not zero for λ>0\lambda > 0. So the parabola fails the equation.

Lengths

Straight line: 2≈1.414\sqrt 2 \approx 1.414. Parabola: ∫011+4λ2 dλ≈1.479\int_0^1\sqrt{1+4\lambda^2}\,d\lambda \approx 1.479, longer. It is exactly the α=−1\alpha = -1 path of the family in the graph; its mirror image, y=2x−x2y = 2x - x^2 at α=+1\alpha = +1, has the same length.

Summary#

  • A geodesic is a path of stationary length between fixed endpoints: usually a minimum in space, possibly a maximum or saddle in spacetime.
  • Index notation: ds2=gij dxidxjds^2 = g_{ij}\,dx^i dx^j, summed over repeated upper–lower indices (Einstein convention).
  • A path xi(λ)x^i(\lambda) has length S=∫gij dxidλdxjdλ dλS = \int\sqrt{g_{ij}\,\frac{dx^i}{d\lambda}\frac{dx^j}{d\lambda}}\,d\lambda. Compare it with xi+αwix^i + \alpha w^i, wi=0w^i = 0 at the ends, and demand dS/dα=0dS/d\alpha = 0 at α=0\alpha = 0 for all wiw^i.
  • Integration by parts, with no boundary term, leaves ∫{⋯ } wi dλ=0\int\{\cdots\}\,w^i\,d\lambda = 0 for all wiw^i, so {⋯ }=0\{\cdots\} = 0: the geodesic equation ddλ[1Agijdxjdλ]=12A ∂igjk dxjdλdxkdλ\frac{d}{d\lambda}\left[\frac{1}{\sqrt A}g_{ij}\frac{dx^j}{d\lambda}\right] = \frac{1}{2\sqrt A}\,\partial_i g_{jk}\,\frac{dx^j}{d\lambda}\frac{dx^k}{d\lambda}.

Symbols on this page

Show the 10 symbols used here
α\alpha alpha
The small number that sizes a variation of the path; the geodesic condition is taken at α = 0. Say it: “alpha”.
λ\lambda lambda
The parameter along a path; in the spacetime formulas λ is also an index running over 0–3. Say it: “lambda”.
≈\approx approximately equal
Is approximately equal to. Say it: “is approximately”.
≤\le less than or equal
Is less than or equal to. Say it: “is less than or 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”.
∫\int integral
Integral: the sum of infinitely many small pieces, such as the area under a curve. Say it: “the integral of”.
∑\sum sum
Sum: add up the terms over the range written below and above. Say it: “the sum of”.
a~\tilde a tilde
A modified version of a quantity (ã, k̃); its meaning is given where it is introduced. Say it: “tilde”.
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. What is a geodesic?
  2. In ds2=gij dxidxjds^2 = g_{ij}\,dx^i dx^j with the Einstein convention, what does the expression mean?
  3. Why do the varied paths satisfy wi(0)=wi(λf)=0w^i(0) = w^i(\lambda_f) = 0?
  4. In the flat plane, the paths y=x+α x(1−x)y = x + \alpha\,x(1-x) run from (0,0)(0,0) to (1,1)(1,1). What is the length at α=0\alpha = 0 and what is dS/dαdS/d\alpha there?
  5. Why can the quantity in curly brackets be set to zero once its integral against wiw^i vanishes for every wiw^i?
0 of 5 answered

Flashcards

1 / 5