Install this package!

The KerrGeodesics package for Mathematica provides functions for computing bound timelike geodesics and their properties in Kerr spacetime.

Example usage

As a quick example, the figure at the top of this page is made using the simple commands:

orbit = KerrGeoOrbit[0.998, 3, 0.6, Cos[π/4]];
{t, r, θ, φ} = orbit["Trajectory"];

Followed by the plot command:

 ParametricPlot3D[{r[λ] Sin[θ[λ]] Cos[φ[λ]], r[λ] Sin[θ[λ]] Sin[φ[λ]], r[λ] Cos[θ[λ]]}, {λ, 0, 20}, 
  ImageSize -> 700, Boxed -> False, Axes -> False, PlotStyle -> Red, PlotRange -> All],
 Graphics3D[{Black, Sphere[{0, 0, 0}, 1 + Sqrt[1 - 0.998^2]]}]

Orbital parametrization

The orbits are parameterized by the following

$a$ - the black hole spin
$p$ - the semi-latus rectum
$e$ - the eccentricity
$x_\text{inc} = \cos\theta_\text{inc}$ - the orbital inclination.

The parametrization $\{a,p,e,\theta_\text{inc}\}$ is described in, e.g., Sec. II of arXiv:gr-qc/0509101

Orbital constants and frequencies

The constants of the motion can be computed using


The above three can be computed together using KerrGeoConstantsOfMotion[a,p,e,x].

The orbital frequencies (w.r.t Boyer-Lindquist time $t$) are computed using KerrGeoFrequencies[a,p,e,x]. For this function you can pass the option Time->"Mino" to compute the frequencies w.r.t. Mino time.

Special orbits

The package allows you compute a variety of special orbits including the innermost stable circular/spherical orbit (ISCO/ISSO), innermost bound spherical orbit (IBSO), the photon orbit and the location of the separatrix between stable and plunging orbits. The relevant functions are:


Further examples

See the Documentation Centre for a tutorial and documentation on individual commands. Example notebooks can also be found in the Mathematica Toolkit Examples repository.

Authors and contributors

Niels Warburton, Maarten van de Meent, Zach Nasipak, Thomas Osburn, Charles Evans, Leo Stein, Philip Lynch, Oliver Long