Consider a space station moving in a circular orbit around some massive body (a planet, a star, or perhaps a black hole), “tidally locked” so that it always keeps the same face towards whatever it’s orbiting. At the centre of the station objects will be weightless, but away from the centre how will they move? Or if they are kept from moving by appropriately oriented “floors”, what weight will they experience? These seem like very simple questions, but it turns out that they involve a lot of interesting physics, and can even reveal whether or not the station is in an orbit where relativistic effects come into play, without the need to make observations of anything outside the station itself.

An overview of the situation, with lots of illustrations, is given in this companion page. Here, we give a detailed mathematical treatment.

Since the station is much smaller than the size of its orbit, we can treat all forces as *linear* in the distances
from the centre. This means it’s possible to completely describe the motion of all test particles in terms of
the tidal stretching or squeezing in each of three orthogonal directions, along with the rate of spin experienced by the station.
The tidal accelerations will grow linearly with
distance from the centre of the station, so we will describe them as accelerations per unit distance. The centrifugal force due to spin also grows linearly with distance
from the centre, within the plane orthogonal to the spin axis (which will be the same as the plane of the orbit). At first glance it might then seem that we can simply
absorb the effects of spin into two of the tidal accelerations, but that would only account for centrifugal force; there is also the Coriolis force that is experienced by particles
moving within the same plane, and which depends on their velocity in that plane rather than their position.

We will label the directions within the station as follows:

Label for direction | Meaning | Name in Incandescence | Coordinate for measuring distance from centre |
---|---|---|---|

r | Outwards, away from central body | sard | x |

–r | Inwards, towards central body | garm | –x |

φ | Along direction of orbital motion | rarb | y |

–φ | Opposite direction to orbital motion | sharq | –y |

z | “Above” plane of orbit (orbit looks counter-clockwise from above) | shomal | z |

–z | “Below” plane of orbit (orbit looks clockwise from below) | junub | –z |

We have used *r* and φ to label the radial and orbital directions, because these letters make the meaning of the directions in the context of the orbit very clear,
but because *r* is used for distance from the central massive body and φ is an angular coordinate, we need to introduce coordinates *x* and *y* which measure distances
in the *r* and φ directions from the centre of the station itself.

We will call the pure tidal accelerations *a*_{r}, *a*_{z} and *a*_{φ}, with the sign convention that tidal stretching is positive, and tidal
squeezing is negative. In other words, given that moving away from the centre of the station in the positive *r* direction gives rise to an acceleration in the same
direction, *a*_{r} will be positive. *In the absence of spin*, this definition would mean d^{2}*x*/d*t*^{2} = *a*_{r} *x*, etc., but
of course the equations for *x* and *y* will need to be modified to take account of spin.

We will call the angular velocity of the station’s spin (which is around the *z*-axis), ω_{s}. In Newtonian mechanics this can really only have one
meaning, but in general relativity we need to be clear and state that it is the rate (measured within the station) at which directions fixed with respect to the
structures of the station rotate, compared to **parallel transport**. This can be discovered experimentally by looking at the rate at which a gyroscope
precesses relative to the station walls: the gyroscope is taken as fixed, and the station is taken to be rotating with the same rate in the opposite direction.

We will call the net accelerations per unit distance, including the effects of centrifugal force,
*A*_{r}, *A*_{z} and *A*_{φ}, where:

A_{r}=a_{r}+ ω_{s}^{2}

A_{z}=a_{z}

A_{φ}=a_{φ}+ ω_{s}^{2}

Now, it’s clear that particles that share exactly the same circular orbit as the station’s centre of mass, but are simply displaced a short distance along the orbit, will remain fixed relative to the station. So we will always have:

A_{φ}=a_{φ}+ ω_{s}^{2}= 0

a_{φ}= – ω_{s}^{2}

The remaining, non-zero net accelerations (*A*_{r} > 0, *A*_{z} < 0), along with the Coriolis force, will determine the motion of anything free-falling within the station,
and the weight felt by anything that is not in free fall. For a particle in free fall we have the differential equations:

d^{2}x/dt^{2}=A_{r}x+ 2 ω_{s}dy/dt

d^{2}y/dt^{2}= –2 ω_{s}dx/dt

d^{2}z/dt^{2}=A_{z}z

In general, solutions to this set of equations will be unbounded, but there are two families of bounded, periodic solutions. One is fairly obvious:

x= 0

y=B

z=Ccos(ω_{z}t)

whereB,Care any constants and ω_{z}^{2}= –A_{z}

This describes a particle oscillating back and forth across the plane of the orbit. We can make a particle execute this motion by placing
it at rest at *x*=0, *y*=*B*, *z*=*C*, and releasing it. The other solution is a little less obvious:

x=Dsin(ω_{r}t)

y= 2D(ω_{s}/ ω_{r}) cos(ω_{r}t) +B

z= 0

whereB,Dare any constants and ω_{r}^{2}= 4 ω_{s}^{2}–A_{r}

This describes a particle moving in an ellipse that is centred on the line *x*=0, the line of the station’s orbit (the “Null Line” in *Incandescence*).
How do we get a particle to move like this? When *t*=0, we have *x*=0, d*x*/d*t* = *D* ω_{r}, while *y* can be anything by choice of *B*, and d*y*/d*t*=0. So by firing a
particle in the *r*-direction from the Null Line, we can guarantee that it executes this kind of cyclic motion.

From the perspective of orbital motion, the particle is travelling along a non-circular orbit which periodically returns it to the same relationship with the station; however, in
general (in the relativistic cases) that period won’t be the same as the period of the station’s orbit.
The particle’s orbit represents a slight perturbation to the station’s orbit
that nonetheless remains close to it, an indication that the orbit is stable. But note that if *A*_{r} ≥ 4 ω_{s}^{2}, the motion
ceases to be periodic, implying that the station’s orbit ceases to be stable!

The qualitative behaviour of various test particle trajectories and the relationships between their periods will depend on *ratios* between the three parameters, *A*_{r}, *A*_{z}
and ω_{s}, so if we are concerned with such relationships rather than absolute time scales the number of parameters is reduced to two.
Furthermore, there is an equation, *a*_{r} + *a*_{z} + *a*_{φ} = 0, that arises from the gravitational physics of tidal accelerations,
and reduces the number of independent ratios to just one. So every *dimensionless* fact that these test particle motions can reveal by themselves (without reference to clocks or rulers calibrated by
some non-gravitational process) can be encapsulated in a single number. For example, we have seen that if *A*_{r} / ω_{s}^{2} ≥ 4,
the station’s orbit is not stable. But we can answer the same question about the orbit’s stability by looking at another ratio, such as *A*_{r} / *A*_{z}, the ratio
between the net accelerations. It’s easy to show that:

A_{r}/ ω_{s}^{2}= 2 (A_{r}/A_{z}) / [1 +A_{r}/A_{z}]

A_{r}/A_{z}= (A_{r}/ ω_{s}^{2}) / [2 –A_{r}/ ω_{s}^{2}]

and so the orbit will become unstable if *A*_{r} / *A*_{z} ≥ –2.

We will now find explicit formulas for *A*_{r}, *A*_{z}, and ω_{s} in three cases: Newtonian gravity, the Schwarzschild geometry around
a non-rotating black hole, and the Kerr geometry around a rotating black hole. In everything that follows, we will adopt units where the speed of light, *c*,
and the gravitational constant, *G*, are equal to 1.

The station orbits the central mass *M* with an orbital angular velocity, ω_{orbit}, such that orbital centrifugal acceleration balances
gravitational acceleration:

ω_{orbit}^{2}r=M/r^{2}

ω_{orbit}^{2}=M/r^{3}

But in the Newtonian case, the fact that the rotation is tidally locked implies that the station’s rate of spin is exactly the same as its orbital angular velocity:

ω_{s}= ω_{orbit}= √(M/r^{3})

We can find *A*_{z} = *a*_{z} easily by either of two approaches. The first is to note that by the spherical symmetry of the gravitational field, *a*_{z} = *a*_{φ},
and since *a*_{φ} = – ω_{s}^{2} we have:

A_{z}=a_{z}=a_{φ}= – ω_{s}^{2}= –M/r^{3}

The second is to note that, again by spherical symmetry, a particle displaced in the *z* direction will move in an inclined orbit with the same period as that of the centre of the station,
and for small displacements this will look just like simple harmonic motion with a “spring constant” equal to minus the square of the angular frequency of the motion.
In other words, if:

z=Ccos(ω_{orbit}t)

then

d^{2}z/dt^{2}= – ω_{orbit}^{2}z

But by definition, d^{2}*z*/d*t*^{2} = *A*_{z} *z*, giving us the same result for *A*_{z}.

We find *a*_{r} by taking the derivative of the gravitational acceleration:

a_{r}= d(–M/r^{2}) / dr= 2M/r^{3}

A_{r}=a_{r}+ ω_{s}^{2}= 3M/r^{3}

Note that *a*_{r} + *a*_{z} + *a*_{φ} = 0, a result that turns out to hold in general relativity as well as it does in Newtonian physics.
In the novel *Incandescence*, this is known as “Zak’s principle”.

We have a ratio for the net accelerations of:

A_{r}/A_{z}= –3

If we compute ω_{r}, the angular frequency for the periodic motion we get if we fire a particle radially from the Null Line, we find:

ω_{r}^{2}= 4 ω_{s}^{2}–A_{r}= 4M/r^{3}– 3M/r^{3}=M/r^{3}

This is the same as the other angular frequencies: that of the orbit, and the station’s spin. Another way we can think of
ω_{r} is in
relation to the “effective potential” for a particle orbiting the central mass. A particle in a circular orbit of radius *r*_{0} has angular momentum per unit mass of:

J= ω_{orbit}r_{0}^{2}= √(Mr_{0})

The angular velocity of a particle with the same angular momentum, but not necessarily in a circular orbit, is:

ω(r) = (J/r^{2}) = √(Mr_{0}) /r^{2}

The energy per unit mass of this particle is:

E= –M/r+ (1/2) (dr/dt)^{2}+ (1/2) (rω(r))^{2}

= –M/r+ (1/2) (dr/dt)^{2}+ (1/2) (Mr_{0}) /r^{2}

The second term is the kinetic energy (per unit mass) due to radial motion, *K*_{r}, so we have *E* = *K*_{r} + *V*(*r*), where:

V(r) = –M/r+ (1/2) (Mr_{0}) /r^{2}

So the radial motion of the particle will correspond to motion in this potential *V*(*r*). The first two derivatives of *V*(*r*), evaluated at *r*=*r*_{0}, are:

dV(r)/dr|_{r=r0}= 0

d^{2}V(r)/dr^{2}|_{r=r0}=M/r_{0}^{3}

So locally the potential well will look like a harmonic oscillator potential with a “spring constant” of *M* / *r*_{0}^{3}, and a particle perturbed slightly
from a circular orbit at *r*=*r*_{0} will undergo oscillations with a frequency whose square is equal to *M* / *r*_{0}^{3}.

In the Schwarzschild geometry around a non-rotating black hole, the Schwarzschild *r* coordinate is defined so that the circumference of a circle of constant *r*, as measured
by someone at rest with respect to the black hole, is 2πr. Far from the hole, changes in the *r* coordinate
correspond to proper distances,
but in the highly curved geometry close to the hole, the actual distance between, say, circular orbits with two different
*r* values *r*_{1} and *r*_{2}, will be greater
than |*r*_{1} – *r*_{2}|.

The hole’s event horizon is at *r*=2M. For objects with rest mass, circular orbits only exist for *r* > 3*M* (though photons can orbit at *r* = 3*M*), and those orbits are only
stable for *r* > 6M.

The proper angular velocity of a circular orbit, as measured by a clock moving with the orbiting body, is given by:

ω_{orbit}^{2}=M/ [r^{2}(r– 3M)]

which is faster than the Newtonian angular velocity for the same r. Surprisingly, though, the rate at which a gyroscope carried on the orbiting body precesses takes exactly the same form as in the Newtonian case, so we have:

ω_{s}^{2}=M/r^{3}

anda_{φ}= –M/r^{3}

The Schwarzschild solution has the same perfect spherical symmetry as the gravitational field around a spherical body in Newtonian gravity. However, relativistic effects
mean that we *can’t* assume, as we did in the Newtonian case, that *a*_{z} = *a*_{φ}; the fact that the station is moving
in the φ direction breaks the symmetry. What we can do is use the spherical symmetry to note that an orbit inclined with respect to the station’s
orbit will have an identical period, and so we can use the same argument as in the Newtonian case to connect *a*_{z} to
ω_{orbit}^{2}, giving us:

A_{z}=a_{z}= – ω_{orbit}^{2}= –M/ [r^{2}(r– 3M)]

If we trust in Zak’s principle, *a*_{r} + *a*_{z} + *a*_{φ} = 0 (which is really the vacuum Einstein equation in disguise), we
can now get everything we need:

a_{r}= –a_{z}–a_{φ}=M(2r– 3M) / [r^{3}(r– 3M)]

A_{r}=a_{r}+ ω_{s}^{2}= 3M(r– 2M) / [r^{3}(r– 3M)]

This value for *a*_{r} agrees with an explicit calculation of the Riemann curvature tensor in a suitable basis. Note that *a*_{r} and *A*_{r}
are greater than the Newtonian values.

The ratio of net accelerations is:

A_{r}/A_{z}= – 3 (1 – 2M/r)

which will lie between –3, for very large orbits, and –1 for the smallest possible orbits at *r*=3*M*, with a value of –2 for the smallest *stable* orbits at *r*=6M.
So if we measure a ratio of tidal accelerations
other than the Newtonian value of –3, that tells us immediately that we are in a relativistic orbit of some kind, without any need to make external observations.

To confirm that *r*=6*M* is the transition between stable and unstable orbits, as we claimed earlier, we can
compute:

ω_{r}^{2}= 4 ω_{s}^{2}–A_{r}=M(r– 6M) / [r^{3}(r– 3M)]

Because ω_{r} is *slower* than the orbital angular frequency
ω_{orbit} by a factor of √(1 – 6*M*/*r*), the periastrons of a non-circular orbit (the points closest to the central body) will occur after the orbiting body has
completed more than one orbit, and so these points will advance around the central body in the same direction as the orbit.

We can cross-check this value by noting that the effective potential for a particle moving in the Schwarzschild geometry, with the same angular momentum as it would have in a circular
orbit of radius *r*_{0}, is ^{[1]}:

V(r) = (1/2) (1 – 2M/r) (1 +Mr_{0}^{2}/ [r^{2}(r_{0}– 3M)])

The first two derivatives of *V*(*r*), evaluated at *r*=*r*_{0}, are:

dV(r)/dr|_{r=r0}= 0

d^{2}V(r)/dr^{2}|_{r=r0}=M(r_{0}– 6M) / [r_{0}^{3}(r_{0}– 3M)]

This tells us that for *r* > 6*M* a particle perturbed slightly from a circular orbit at *r*=*r*_{0} will experience oscillations with an angular frequency whose square is equal
to this second derivative.

The Kerr geometry ^{[2]} describes a black hole with angular momentum, which is usually parameterised by a quantity, *a*, that lies between –*M* and *M*. We’ll define the positive
direction to be whichever way our station is orbiting, so *a* > 0 means the station is in a “co-rotating” orbit, and *a* < 0 means it is in a “counter-rotating” orbit,
with respect to the black hole’s angular momentum.
We will also take it as given that our station is orbiting in the black hole’s equatorial plane, the plane orthogonal to its angular momentum vector.

We’ll use a generalisation of the Schwarzschild coordinates known as Boyer-Lindquist coordinates. The Boyer-Lindquist *r* coordinate no longer corresponds to the circumference of circles divided
by 2π, but we’ll see shortly that it does have one very nice property.

The event horizon is located at *r* = *M* + √(*M*^{2} – *a*^{2}). For general values of *a*, the locations of the smallest orbit and the
smallest stable orbit are complicated expressions, but for the extremal values *a*=±*M* they become simple: for *a*=*M* there are stable orbits all the way down to the event horizon at *r*=*M*,
while for *a*=–*M* the smallest stable orbit is at 9*M*, and the smallest orbit at 4*M*.

The nice property of *r* that we alluded to is that the rate at which a gyroscope carried on the station precesses *still* takes exactly the same form as in the Newtonian case,
so we again have:

ω_{s}^{2}=M/r^{3}

anda_{φ}= –M/r^{3}

The proper angular velocity for the orbital motion is given by:

ω_{orbit}^{2}=M/ (rd)

whered=r(r– 3M) + 2a√(Mr)

Because the space-time is no longer spherically symmetrical, there is no guarantee that orbits inclined to the equatorial plane will be circular,
or even form closed curves,
let alone possess the same period as orbits of the same
size lying in the plane. So the link we used previously between *a*_{z} and ω_{orbit} no longer applies, and we need to derive
*a*_{z} from first principles: by computing the Riemann curvature tensor in a suitable basis, to find all the tidal accelerations. The results are:

A_{z}=a_{z}= –M[r^{2}– 4a√(Mr) + 3a^{2}] / (r^{3}d)

a_{r}=M[r(2r– 3M) – 2a√(Mr) + 3a^{2}] / (r^{3}d)

along with *a*_{φ} which we’ve already noted. If *a* > 0, ω_{z}
(defined by ω_{z}^{2} = –*a*_{z}) will be less than ω_{orbit},
and the nodes where an ocillating test particle crosses the plane of the unperturbed equatorial orbit will be separated by more than 180 degrees; if *a* < 0 then
ω_{z} > ω_{orbit}, and the nodes will be separated by less than 180 degrees.

We also have:

A_{r}= 3M[r(r– 2M) +a^{2}] / (r^{3}d)

A_{r}/A_{z}= – 3 (1 – 2 [a– √(Mr)]^{2}/ [r^{2}– 4a√(Mr) + 3a^{2}])

ω_{r}^{2}= 4 ω_{s}^{2}–A_{r}=M[r(r– 6M) + 8a√(Mr) – 3a^{2}] / (r^{3}d)

We can re-express the quantity *d* that appears in the denominators of *A*_{r}, *A*_{z} and
ω_{orbit}
in terms that include the ratio *A*_{r} / *A*_{z}:

d= 3 (1 +A_{r}/A_{z}) (a^{2}+r(r– 2M)) / (2A_{r}/A_{z})

From this, we can see that if *A*_{r} / *A*_{z} = –1, *d* = 0, and ω_{orbit}
becomes singular. So
the value of –1 for this ratio marks the smallest possible circular orbit for a given black hole; this is a limiting case that is only
possible for orbiting photons. As we’ve previously noted, orbits become unstable when *A*_{r} / *A*_{z} = –2.

In the first section of this page, we discussed the fact that test particle motions on their own only tell us a single dimensionless number, and now that we’re considering
orbits around a rotating black hole, which are parameterised by *two* dimensionless numbers (*r*/*M*, the size of the orbit, and *a*/*M*, the angular momentum of the hole), the
consequences of that become clear: without an external observation, we can never hope to determine both *r*/*M* and *a*/*M*.

In the diagram below, each point represents a possible circular orbit around a rotating black hole, with *r*/*M* and *a*/*M* shown on the coordinate axes. If we only know about the motions of
test particles inside the station, we can only determine a single dimensionless ratio, such as the ratio *A*_{r} / *A*_{z}. The grey curves
are curves of constant *A*_{r} / *A*_{z}, and knowing which one the station lies on will tell us
how close the orbit is to the transition from stability to instability, but by itself it can’t fix the values of *r*/*M* and *a*/*M*.
However, by making an external observation that reveals the ratio of the station’s orbital frequency,
ω_{orbit}, to the frequency of a test particle oscillating back and forth across the orbital plane,
ω_{z}, we can identify which green curve the station lies on, and deduce *r*/*M* and *a*/*M* from its intersection
with the grey curve determined by *A*_{r} / *A*_{z}.

[1] *Gravitation* by Charles Misner, Kip Thorne and John Wheeler,
W.H. Freeman, San Francisco, 1973. Adapted from equation 25.16; what we use as *V*(*r*) is half the square of Misner, Thorne & Wheeler’s *V*(*r*), in order to get the oscillation
rate from the second derivative by making a precise analogy with a non-relativistic particle in a potential well.

[2] Misner, Thorne & Wheeler, *op. cit.*, Chapter 33.