Schwarzschild orbit

Note

I’m not sure “Schwarzschild orbit” is a real thing, that’s just what I’m calling a body orbiting r=0r=0 in a Schwarzschild spacetime.

The Schwarzschild metric has some nice conserved quantities. In particular:

E^=(12GMr) ⁣dt ⁣dτ,and L^z=r2 ⁣dϕ ⁣dτ \begin{align*} \hat E &= \left(1-\frac{2GM}{r}\right) \frac{\dt}{\d\tau}, \\ \text{and } \hat L_z &= r^2 \frac{\d\phi}{\d\tau} \end{align*}

We will use the invariant uu=1\fv u \cdot \fv u = -1 and the above conserved quantities to explore orbits in the θ=π/2\theta=\pi/2 plane.

1=(12GM/r)t˙2+r˙212GM/r+r2ϕ˙2. \begin{align*} -1 = - (1-2GM/r) \dot t^2 + \frac{\dot r^2}{1-2GM/r} + r^2 \dot \phi^2. \end{align*}

We rewrite the above as

r˙2=E^2Veff(r)Veff(r)=(12GMr)(1+L^z2r2). \begin{align*} \dot r^2 &= \hat E^2 - V_\eff(r) \\ V_\eff(r) &= \left(1-\frac{2GM}{r}\right)\left(1 + \frac{\hat L_z^2}{r^2}\right). \end{align*}

Since r˙\dot r must be real, we require E^>Veff\hat E > V_\eff.

Plotting VeffV_\eff, we see the physically valid regions for some E^\hat E.

When E^\hat E is greater than the trough of VeffV_\eff, an eccentric orbit is possible. When it is less, the only solution is a body crashing into r=0r=0. Various energy levels lead to various stable or unstable orbits.

The effective potential depends on the angular momentum. As the angular momentum becomes small, stable orbits are no longer possible.

Light

Light can orbit in a Schwarzschild spacetime, too. The equations are a little different.

Consider a photon moving in from r=r=\infty towards the center of the Schwarzschild spacetime, separated by a transverse distance bb.

We write the Lagrangian in terms of λ\lambda as we do for light in strong gravity (assuming θ=π/2\theta = \pi/2)

L=12[(12GM/r)t˙2+(12GM/r)1r˙2+r2ϕ˙2], \begin{align*} L &= \frac12 \left[-(1-2GM/r) \dot t^2 + (1-2GM/r)^{-1} \dot r^2 + r^2 \dot\phi^2 \right], \\ \end{align*}

We see that there are two conserved quantities:

E^=(12GMr)t˙=gttpt=ptL^z=r2ϕ˙=bE^ \begin{align*} \hat E &= -\left(1 - \frac{2GM}{r}\right) \dot t = g_{tt} p^t = p_t \\ \hat L_z &= r^2 \dot \phi = b \hat E \end{align*}

From pp=0\fv p \cdot \fv p = 0 we find

r˙2=E^2L^z2r2(12GMr)=L^z2b2L^z2r2(12GMr)L^z2r˙2=1b21r2(12GMr)=:1b2Vlight(r). \begin{align*} \dot r^2 &= \hat E^2 - \frac{\hat L_z^2}{r^2} \left(1-\frac{2GM}{r}\right) \\ &= \frac{\hat L_z^2}{b^2} - \frac{\hat L_z^2}{r^2} \left(1-\frac{2GM}{r}\right) \\ \hat L_z^2 \dot r^2 &= \frac{1}{b^2} - \frac{1}{r^2} \left(1-\frac{2GM}{r}\right) \\ &=: \frac{1}{b^2} - V_\text{light}(r). \end{align*}

Again, we get certain valid radii given some bb and the light potential.

There are three regimes light can be in:

  1. b>33GMb>3\sqrt3GM: light comes in from r=r=\infty, bends as it passes

around the center, and then goes out to \infty again.

  1. v<33GMv<3\sqrt3GM: light comes in from r=r=\infty and crashes into

r=0r=0. Light disappears forever.

  1. v=33GMv=3\sqrt3GM: light comes in from r=r=\infty, then reaches r=3GMr=3GM

and orbits there forever! This is an unstable orbit, though, so the light is likely to crash in or fly away given any perturbation.