# Precession in Special and General Relativity

**The Absolute Derivative**

In relativity we typically deal with two types of quantities – fields, which are defined everywhere – and particle properties, which are defined only along a curve or world line. The familiar covariant derivative is appropriate when we need to differentiate a field. A field is a function of all four coordinates, and the covariant derivative of φ(x) consists of the four partial derivatives ∂φ/∂x^{µ} plus correction terms involving the Christoffel symbols, one for each tensor index on φ.

A particle property φ(s), on the other hand, is a function only of a single parameter s running along the curve. In this situation, the partial derivatives of φ with respect to the four coordinates do not exist. (Unfortunately many references miss this point!) Writing partial derivatives would require that φ be defined everywhere in a neighborhood of the curve, which is not the case. The correct derivative to use for a particle property is called the absolute derivative, and is written either δφ/δs or Dφ/Ds. Similar in form to the covariant derivative, it consists of an ordinary derivative dφ/ds plus a correction term for each tensor index. For example, if the world line is given parametrically as x^{µ}(s), the absolute derivative of a contravariant vector φ^{µ} is Dφ^{µ}/Ds = dφ^{µ}/ds + φ^{ν} Γ^{µ}_{νσ} dx^{σ}/ds.

For a given timelike curve x^{µ}(s) we define the vector v^{µ} = Dx^{µ}/Ds. This 4-vector is clearly tangent to the curve. If the parameter s is chosen to measure proper time, then v^{µ} will be a unit vector, **v·v** = -1. Next we define an acceleration 4-vector a^{µ} = Dv^{µ}/Ds. Since D(v_{µ}v^{µ})/Ds = 0 = v_{µ}a^{µ}, a^{µ} is orthogonal to v^{µ} and therefore lies in the instantaneous rest frame of the particle.

**Fermi-Walker Transport**

Given a vector **W** at a single point on the world line, we can define it at other points along the curve by specifying how it is to be propagated or transported along the curve. The simplest method is parallel transport: D**W**/Ds = 0. For example if the particle is nonaccelerating, then the velocity vector is parallel-transported: D**v**/Ds = 0. But a generally more useful concept is Fermi-Walker transport. We often need to deal with a previously 3-dimensional quantity. In 4 dimensions it will be confined to the instantaneous rest frame of the particle and represented by a vector orthogonal to **v**. As the particle accelerates, the rest frame changes. Consequently for **W** to remain orthogonal it must undergo a Lorentz transformation in the **v-a** plane. The definition of Fermi-Walker transport is

DW^{µ}/Ds = W^{ν}(v^{µ}a_{ν} – a^{µ}v_{ν})

Note that using F-W transport, if **W·v** = 0 initially, it will remain so. Physically, Fermi-Walker transport defines what we mean by a nonnrotating instantaneous rest frame. To say it again, the only way **W** changes is a Lorentz boost applied in the direction of the particle’s acceleration.

The spin vector **S** is such a quantity. The fact that **S** evolves according to Fermi-Walker transport is the source of precession effects in both special and general relativity.

**Precession in Flat Space**

Consider a particle moving in a circular orbit of radius r in the equatorial plane. Define a comoving orthonormal basis of 4-vectors:

**e _{r}** = (cos ωt, sin ωt, 0, 0)

**e**= (-sin ωt, cos ωt, 0, 0)

_{φ}**e**= (0, 0, 0, 1)

_{t}Thus **e _{r}** points radially outward, and

**e**points tangentially along the circumference.

_{φ}Let τ be the particle’s proper time. Then

D**e _{r}** /Dτ = γω

**e**

_{φ}D

**e**/Dτ = – γω

_{φ}**e**

_{r}where γ = dt/dτ.

The particle’s position, velocity and acceleration 4-vectors are:

**x** = r **e _{r}** + t

**e**

_{t}**v**= D

**x**/Dt = rγω

**e**+ γ

_{φ}**e**

_{t}**a**= D

**v**/Dt = – rγ

^{2}ω

^{2}

**e**

_{r}From this we can derive that the magnitude of the 3-velocity is v = rω, and that γ^{2} = 1/(1 – r^{2}ω^{2}).

Now let’s suppose that the particle carries a spin: a unit 4-vector **S** orthogonal to **v**. In general **S** can be written as some time-dependent linear combination of **e _{r}**,

**e**and

_{φ}**e**:

_{t}**S** = a(τ) **e _{r}** + b(τ)

**e**+ c(τ)

_{φ}**e**

_{t}**S·v**= 0 ⇒ c(τ) = rω b(τ)

**S·a**= – rγ

^{2}ω

^{2}a

The time derivative of **S** comes from derivatives of a, b, c and also from the fact that the basis vectors **e _{r}**,

**e**are rotating:

_{φ}D**S**/Dτ = (da/dτ **e _{r}** + db/dτ

**e**+ dc/dτ

_{φ}**e**) + (γωa

_{t}**e**– γωb

_{φ}**e**)

_{r}Finally, this must be inserted into the Fermi-Walker transport equation:

D**S**/Dτ = (**S·a**)**v** = – r^{2}γ^{3}ω^{3}a **e _{φ}** – rγ

^{3}ω

^{2}a

**e**

_{t}and the coefficients matched:

coefficients of **e _{r}**: da/dτ – γωb = 0

coefficients of

**e**: db/dτ + γωa = – r

_{φ}^{2}γ

^{3}ω

^{3}a

coefficients of

**e**: dc/dτ = – rγ

_{t}^{3}ω

^{2}a

Which simplify to:

da/dτ = γω b

db/dτ = – γ^{3}ω a

(The dc/dτ equation is identically satisfied by our earlier result, c(τ) = rω b(t).)

The first two equations may be combined into a harmonic oscillator equation:

d^{2}a/dτ^{2} + γ^{4}ω^{2} a = 0,

with solution a ~ cos(Ωτ) where Ω = γ^{2}ω.

The important thing to note is the extra factor of γ. Although the orbital position varies like (sin, cos)(ωt) = (sin, cos)(γωτ), the spin vector S varies like (sin, cos)(γ^{2}ωτ), that is, more slowly. Thus, relative to a stationary frame, **S** exhibits a net precession in a retrograde sense, an effect known as Thomas precession.

**Precession in the Schwarzschild Metric**

Now consider the same particle moving in a circular orbit in a Schwarzschild field, again with coordinate angular velocity ω = dφ/dt. The same procedure is to be followed as we did in flat space, and we’ll mostly just compare results. The velocity 4-vector is

v^{μ} = (v^{r}, v^{φ}, v^{t}) = (0, γω, γ)

with γ determined by the normalization condition

v·v = 1 = γ^{2}(g_{tt} + ω^{2} g_{φφ}) ⇒ γ^{2} = (1 – 2M/r – r^{2}ω^{2})^{-1}. Before it was just (1 – r^{2}ω^{2})^{-1}.

The acceleration 4-vector is

a^{μ} = Dv^{μ}/Dτ = dv^{μ}/dt + Γ^{μ}_{νσ} v^{ν}v^{σ}

The only nonzero component of a^{μ} is the radial component:

a^{r} = -½ g^{rr}(g_{tt,r} v^{t}v^{t} + g_{φφ,r} v^{φ}v^{φ}) = -γ^{2}(r – 2M)(ω^{2} – M/r^{3}), where before we had a^{r} = – γ^{2}rω^{2}.

(Note that for a free particle following a geodesic, the acceleration is zero, a^{r} = 0, and the orbital velocity is given by ω^{2} = M/r^{3}. We recognize this as Kepler’s Law, “period squared goes as distance cubed.” It’s remarkable that in terms of the coordinate angular velocity, the circular orbits in the Schwarzschild field obey Kepler’s Law exactly!)

Now consider the spacelike vector S^{μ} orthogonal to vμ:

S^{μ} = (S^{r}, S^{φ}, S^{t})

**S·v** = 0 = g_{tt}S^{t}v^{t} + g_{φφ}S^{φ}v^{φ} ⇒ S^{t} = r^{3}ω/(r – 2M) S^{φ}, whereas before c(t) = rω b(t)

**S·a** = g_{rr}S^{r}a^{r} = rγ^{2}(ω^{2} – M/r^{3})S^{r}, whereas before **S·a** = rγ^{2}ω^{2} a

Supposing that S^{μ} is Fermi-Walker transported, we get evolution equations for each component:

dS^{r}/dτ = γω(r – 3M)S^{φ} (before it was da/dτ = γω b)

dS^{φ}/dτ = -γ^{3}ω(r – 3M)/r^{2} (before it was db/dτ = -γ^{3}ω a)

Combining these:

d^{2}S^{r}/dτ^{2} + γ^{4}ω^{2}(r – 3M)/r^{2} S^{r} = 0,

results in a harmonic oscillator equation with solutions S^{r} ~ A cos(Ωτ) where Ω = γ^{2}ω(r – 3M)/r. The precession rate Ω for Schwarzschild may be thought of as the combination of two effects:

– a retrograde Thomas precession at a rate γ^{2}ω

– a prograde de Sitter precession at a rate 3γ^{2}ωM/r

**Precession in the Kerr Metric**

For Kerr, results only.

γ^{2} = (1 – (r^{2} + a^{2})ω^{2} – 2M(1 – aω)^{2}/r)^{-1}

a^{r} = – Δγ^{2}[ω^{2} – M(1 – aω)^{2}/r^{3}], where Δ = r^{2} – 2Mr + a^{2}

S^{t} = S^{φ}[rω(r^{2} + a^{2}) – 2Ma(1 – aω)]/[r – 2M(1 – aω)]

**S·a** = rγ^{2}[ω^{2} – M(1 – aω)^{2}/r^{3}]

d^{2}S^{r}/dτ^{2} + Ω S^{r} = 0 where Ω = γ^{2}ω[r – 3M(1 – aω)]/r + γ^{2}Ma(1 – aω)^{2}/r^{3}

— This article was originally part of Physics Forums member Bill_K‘s PF blog. He may not respond to comments.

Great article!

In the section "Precession in Flat Space", the velocity and acceleration 4-vectors should be defined as v = Dx / Dτ and a = Dv / Dτ. (The article as it stands has Dt instead of Dτ.) The actual expressions for v and a are correct.

thanks for your postdoes the Fermi-Walker transport could be determined by fiber bundle?