Schwarzschild Innermost Orbit Stability

In general relativity, gravity is modeled as the curvature of spacetime. For a non-rotating, uncharged black hole, the spacetime geometry is described by the Schwarzschild metric. Test particles follow geodesics in this curved spacetime. One intriguing consequence is that circular orbits can exist but might be unstable if they are too close to the event horizon, even outside it.

Effective Potential and Circular Orbits

By utilizing the symmetries of the Schwarzschild metric (in natural units G=M=c=1), we find conserved quantities for a test particle: specific energy E and specific angular momentum L. For motion restricted to the equatorial plane, the radial equation reduces to a one-dimensional problem governed by an effective potential V(r):

V(r) = (1 - 2/r)(1 + L^2/r^2)

Circular orbits occur at the extrema of the effective potential, where V'(r) = 0 and E² = V(r). Solving these conditions yields the requisite energy and angular momentum for a circular orbit at radius r:

L² = r² / (r - 3)\nE² = (r - 2)² / (r(r - 3))

These equations imply that circular orbits are only possible for r \gt 3. The stability of these orbits is determined by the second derivative of the effective potential, related to the epicyclic frequency κ²:

κ² = (r - 6) / r^4

For r \gt 6, the potential has a local minimum (κ² \gt 0), so the orbit is stable against small radial perturbations. For 3 \lt r \lt 6, it's a local maximum (κ² \lt 0), making the orbit unstable. The boundary r = 6 marks the Innermost Stable Circular Orbit (ISCO).

Calculating Stability for r = 8

Let's consider a test particle in a circular orbit at a radius r = 8. We can calculate its specific energy, specific angular momentum, and the coordinate radial frequency squared to determine its stability.

First, find the specific angular momentum squared (L²):
L² = 8² / (8 - 3) = 64 / 5 = 12.8
L = √(12.8) ≈ 3.5777
Next, find the specific energy squared (E²):
E² = (8 - 2)² / (8(8 - 3)) = 36 / 40 = 0.9
E = √(0.9) ≈ 0.9487
Finally, calculate the radial epicycle frequency squared (κ²) to check stability:
κ² = (8 - 6) / 8^4 = 2 / 4096 ≈ 0.00048828

Since κ² \gt 0, the orbit is stable. The effective potential at this minimum is V(8) = (1 - 2/8)(1 + 12.8/64) = 0.75 × 1.2 = 0.9, which matches E².

Interactive Experiment

Use the interactive simulation to explore circular geodesics and the effective potential in the Schwarzschild metric.

Select the Circular radius to set the orbit of the test particle. Notice how the shape of the effective potential curve and the position of the stationary point change.
Adjust the Potential probe radius to evaluate the effective potential V(x) at different distances for the current orbit's angular momentum.
Observe the calculated values for Specific energy and Coordinate radial frequency squared. Verify that when the radius is below 6, the frequency squared becomes negative, indicating an unstable orbit.

Try it out in the module: Schwarzschild Innermost Orbit Stability.

Model Limitations

This model assumes a non-rotating (Schwarzschild) black hole. For spinning black holes (Kerr metric), the ISCO radius depends on the black hole's angular momentum and whether the orbit is prograde or retrograde. Additionally, this treats the orbiting body as a "test particle" whose mass does not back-react on the spacetime geometry. The analysis also strictly applies to timelike geodesics, not null geodesics (photon orbits).

References

The following reference provides a comprehensive mathematical derivation of geodesics in the Schwarzschild metric, including the effective potential and the Innermost Stable Circular Orbit.

David Tong: Lectures on General Relativity (Section 6: The Schwarzschild Solution)

Schwarzschild Innermost Orbit Stability