Blog entry by Bill_K(adsbygoogle = window.adsbygoogle || []).push({});

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 vectorWat 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: DW/Ds = 0. For example if the particle is nonaccelerating, then the velocity vector is parallel-transported: Dv/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 tov. As the particle accelerates, the rest frame changes. Consequently forWto remain orthogonal it must undergo a Lorentz transformation in thev-aplane. The definition of Fermi-Walker transport is

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

Note that using F-W transport, ifW·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 wayWchanges is a Lorentz boost applied in the direction of the particle's acceleration.

The spin vectorSis such a quantity. The fact thatSevolves 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= (cos ωt, sin ωt, 0, 0)_{r}

e= (-sin ωt, cos ωt, 0, 0)_{φ}

e= (0, 0, 0, 1)_{t}

Thusepoints radially outward, and_{r}epoints tangentially along the circumference._{φ}

Let τ be the particle's proper time. Then

De/Dτ = γω_{r}e_{φ}

De/Dτ = - γω_{φ}e_{r}

where γ = dt/dτ.

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

x= re+ t_{r}e_{t}

v= Dx/Dt = rγωe+ γ_{φ}e_{t}

a= Dv/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-vectorSorthogonal tov. In generalScan be written as some time-dependent linear combination ofe,_{r}eand_{φ}e:_{t}

S= a(τ)e+ b(τ)_{r}e+ c(τ)_{φ}e_{t}

S·v= 0 ⇒ c(τ) = rω b(τ)

S·a= - rγ^{2}ω^{2}a

The time derivative ofScomes from derivatives of a, b, c and also from the fact that the basis vectorse,_{r}eare rotating:_{φ}

DS/Dτ = (da/dτe+ db/dτ_{r}e+ dc/dτ_{φ}e) + (γωa_{t}e- γωb_{φ}e)_{r}

Finally, this must be inserted into the Fermi-Walker transport equation:

DS/Dτ = (S·a)v= - r^{2}γ^{3}ω^{3}ae- rγ_{φ}^{3}ω^{2}ae_{t}

and the coefficients matched:

coefficients ofe: da/dτ - γωb = 0_{r}

coefficients ofe: db/dτ + γωa = - r_{φ}^{2}γ^{3}ω^{3}a

coefficients ofe: 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,Sexhibits 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 beforeS·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}

**Physics Forums - The Fusion of Science and Community**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Precession in Special and General Relativity

Tags:

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

Loading...

Similar Threads - Precession Special General | Date |
---|---|

A Perihelion precession in GR using Robertson expansion | May 25, 2017 |

I Mercury precession | Dec 13, 2016 |

Insights Orbital Precession in the Schwarzschild and Kerr Metrics - Comments | Aug 11, 2016 |

I GR and Mercury precession | Aug 8, 2016 |

Insights Precession in Special and General Relativity - Comments | Aug 8, 2016 |

**Physics Forums - The Fusion of Science and Community**