Hi there,(adsbygoogle = window.adsbygoogle || []).push({});

First post, so be gentle. I'm currently doing a project on Kerr black holes, part of which is based on a project in Edwin F. Taylor's book "Exploring Black Holes" (the chapter is found here) as well as some of Hartle's book (chapter 15). Part of my project, and as a way of better understanding some principles, as well as getting up to speed with MATLAB again, is making a little program that can compute and plot the orbits of particles.

Basically I'm using the terms from Hartle (eqs. (15.18b) and (15.19) (I won't bother with my derivations of them here) derived from the Kerr metric in the equatorial plane:

[tex]

d\tau^2 = \left(1 - \frac{2M}{r}\right) dt^2

+ \frac{4Ma}{r} dt d\phi

- \frac{r^2}{\Delta} dr^2

- R_a^2 d\phi^2,

[/tex]

where [itex]a[/itex] is the spin parameter, [itex]R_a^2 = r^2 + a^2 + \frac{2Ma^2}{r}[/itex] is the reduced circumference, and [itex]\Delta \equiv r^2 - 2Mr + a^2[/itex]. Calculations give (15.18b): [tex]

\frac{d\phi}{d\tau} = \frac{1}{\Delta}\left[\left(1 - \frac{2M}{r}\right)\ell + \frac{2Ma}{r}e\right]

[/tex]

and (15.19): [tex]

\frac{e^2 - 1}{2} = \frac{1}{2}\left(\frac{dr}{d\tau}\right)^2 + V_\text{eff}(r,e,\ell)

[/tex]

where

[tex]

V_\text{eff}(r,e,\ell) = -\frac{2M}{r} + \frac{\ell^2 - a^2(e^2 - 1)}{r^2} - \frac{2M(\ell - a\cdot e)^2}{r^3}.

[/tex]

Now, following Hartle's example (I think it's in Taylor as well), one finds the orbit of a particle, moving radially inwards, by [itex](dr/d\tau)/(d\phi/d\tau)[/itex], which is somewhat simple enough. I haven't written an expression for it, as it isn't very pretty, but by solving that differential equation, I get a nice plot of what I wanted. That is

http://dicktaid.com/kerrorbit.jpg [Broken]

(here [itex]a=M[/itex], [itex]e = 1[/itex] and [itex]\ell = 0[/itex]). The same can be achieved by using the program Taylor supplies on his webpage (the GRorbits.jar file) which he uses in his drafts of the second edition of his book (chapter 16). I want to be able to plot bound orbits as depicted in fig. 16.9-10. Inserting those values in my MATLAB code gives an error (the values he uses for angular momentum, energy and spin parameter implies [itex](dr/d\tau)^2[/itex] is negative).

I'm trying to understand the coding done in the GRorbits.jar file, that is, I'm trying to understand how those orbits are achieved. I've looked around in the books, as well as the net, to find something that helps me make sense of it. Perhaps I'm searching for the wrong things, but I'm not able to locate what I'm looking for. Any pointers you guys can give me is appreciated.

EDIT: Fig. 16.9 from Taylor's drafts so you don't have to download it:

http://dicktaid.com/boundorbit.jpg [Broken]

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

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

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

# Computer plot of particle orbits around Kerr black holes

**Physics Forums | Science Articles, Homework Help, Discussion**