1

Geodsicas en la mtrica de Schwarzschild y Kerr.

Tratamiento numrico
Zamora Ramos, Rafael
Licenciatura de Ciencias Fsicas. Facultad de Fsica. Avd. Reina Mercedes, s/n.
Universidad de Sevilla. E-41012 Sevilla, Spain.

En este trabajo, se pretende estudiar las rbitas de un planeta ligero inmerso en un es-
pacio deformado por una estrella masiva esfrica y esttica, haciendo uso de la Teora
General de la Relatividad mediante resolucin numrica de las ecuaciones correspon-
dientes. Se analizarn los posibles casos y sus diferencias con el caso no relativista. El
mismo procedimiento se aplicar a la trayectoria de un rayo de luz en el marco relativis-
ta. Finalmente veremos, sin profundizar mucho, qu efectos se siguen de una estrella en
rotacin. Se har nfasis en sealar que, para que estos efectos sean observables de for-
ma evidente a los sentidos del ser humano es necesario que la estrella sea en realidad un
agujero negro.

En todo el trabajo vamos a utilizar unidades geometrizadas, en las que la velocidad de la luz y la
constante de la gravitacin valen la unidad. En la Teora General de la Relatividad el espacio, el
tiempo y la masa estn totalmente ligados, por lo que puede usarse la misma unidad para las tres.
Este sistema tiene una gran ventaja en el tratamiento de las ecuaciones. De todas formas, siempre
podemos cambiar al sistema de unidades MLT, usando la tabla 2 de conversin que aparece al final
del trabajo.

Supondremos un espacio vaco en el que hay un cuerpo masivo esfrico esttico. Resolviendo las
ecuaciones de Einstein para el vacio se obtiene la solucin de Schwarzschild que cumple la siguien-
te mtrica en coordenadas polares:

= 1 + + 1 .

En esta mtrica se puede comprobar que existen varias simetras que van asociadas a la existencia
una magnitud invariante. Las dos magnitudes invariantes que en nuestra discusin nos interesan
son, el momento angular l y la energa e.

Nos interesa buscar la trayectoria de partculas con masa en ese espacio. Para eso, buscamos la geo-
dsica de esta mtrica, que sera la trayectoria de esa partcula, teniendo siempre presente que la
masa de la partcula es tan pequea respecto a la de la estrella que se puede anular su aportacin
a la mtrica del espacio.
Tenemos pues que:

= = 1 , (1)

que para nuestra mtrica corresponde con:

= 1 + + 1 = 1 . (2)
 2

Las ecuaciones de Lagrange asociadas son:

=0

Teniendo en cuenta que:

=
=
=
=

De las ecuaciones de Lagrange obtenemos 4 ecuaciones, de las que usaremos las tres ltimas.

2 2 =0 2 1 =0 =0 = (3)

=0 = (4)

1 =0 1 = . (5)

(Las derivadas son respecto al tiempo propio .)

De (3) se deduce que las geodsicas y, por tanto, las trayectorias son planas. Podemos escoger, sin
prdida de generalidad, = /2 .
Por tanto, (4) quedara de la forma:

= . (6)

Sustituyendo (5) y (6) en (2) :

2 2
1 + 1 = 1
2
1

1 + = 1.

Multiplicando por 1 :

2 2
+ 1 = 1

= + + . (7)

Podemos interpretar como la energa y + como un potencial efectivo ( ).

Por lo tanto,
3

= + ( ). (8)

Y podemos compararlo con el caso Newtoniano que vamos a recordar brevemente.

Caso Newtoniano

Partimos de la conservacin de la energa total de la partcula.

= = + , (9)

y teniendo en cuenta la conservacin del momento angular:

= . (10)

Sustituyendo (10) en (9):

= + . (11)

Definimos el potencial efectivo:

( )= , (12)

quedando (11) de la forma:

= + ( ). (13)

Si nos fijamos en la funcin ( ), vemos que los puntos de corte con la recta E nos dan puntos
de retorno, ya que en ellos = 0, como se deduce de (13).

En la figura 1 representamos la funcin ( ) en rojo, en azul, en negro y la ener-

ga total en lnea discontinua. Podemos distinguir cuatro situaciones dependiendo del valor de E:

El primer caso sera cuando la energa corte en el mnimo del potencial y la rbita sera de
radio constante y, por tanto una circunferencia, que es el lmite del caso siguiente.
En este segundo caso hay dos puntos de corte y la rbita estar ligada entre el valor de los
dos radios en que corta la energa con el potencial efectivo. Como veremos posteriormente,
la rbita corresponde a una elipse.
El tercer caso corresponde a energa cero, en la que hay un solo punto de retorno y la part-
cula se para en el infinito. Este caso coincide con una trayectoria parablica.
El cuarto caso es el de energas positivas, en el que hay un solo punto de retorno y la rbita
es no ligada, coincidiendo con una hiprbola.
 4

Fig. 1. Representamos en color rojo el potencial efectivo Newtoniano frente a r y

los dos trminos del potencial efectivo. En lneas discontinuas se representan
tres valores de la energa: negativo, nulo y positivo.
El caso B-B corresponde con trayectorias elpticas, el A con parbolas y el C
con hiprbolas.

Para obtener la ecuacin diferencial que gobierna las trayectorias partimos de (11) y (10). Adems,
teniendo en cuenta que = , obtenemos:

= + , (14)

que es ms fcil de resolver si usamos el siguiente cambio de variable: = 1/ .

(1/ ) 1
= =

De modo que (14) pasa a tener la siguiente forma:

= + . (15)

Si derivamos (15) respecto de :

+ = 0. (16)

Supondremos que 0, es decir, que las trayectorias no son circulares; y que 0, trayectorias
no rectas (no cada libre), por lo que podemos dividir , obteniendo:

+ = , (17)
5

Ecuacin inhomognea1 de fcil solucin analtica:

( )= + cos , (18)

siendo e la excentricidad de la rbita2 y cumpliendose las siguientes condiciones iniciales:

(0) = 0
( )
(0) = = . (19)

Deshaciendo el cambio de variables y definiendo = , tenemos:

( )= , (20)
con
(0) = 0
(0) = = (perihelio)
( )= = (afelio) (21)

La ecuacin (20) es una cnica en polares en la que distinguimos:

= 0 Circunferencia
< 1 Elipse
= 1 Parbola
> 1 Hiprbola

Una vez recordado el caso Newtoniano, volvamos a nuestra ecuacin para el caso relativista. En
este caso porque nos interesa compararlo con el caso Newtoniano, vamos a usar unidades MLT, el
cambio es fcil teniendo en cuanta las tablas 1 y 2 del final del trabajo:

= + + . (22)

Podemos encontrar mayor similitud si definimos

=
y
= .

Sustituyendo en (22):

= + + (23)

que salvo por el ltimo trmino 3, que sera una correccin de la Teora General de la Relatividad, es
comparable al caso Newtoniano. Pero debemos siempre tener presente una diferencia entre (23) y

1
La solucin general de la ecuacin inhomognea es la suma de la solucin general de la homognea ms una solu-
cin particular de la inhomognea.
2
No confundir con la energa definida en (5).
 6

(11): en la ecuacin relativista estamos derivando con respecto al tiempo propio y r no es la distan-
cia a ningn centro, ya que la mtrica no es la mtrica plana en coordenadas esfricas y, por tanto,
no coincide con la distancia del centro al planeta.

A pesar de las diferencias, an podemos clasificar el tipo de rbitas de la misma forma que haca-
mos en el caso clsico.
Volvemos a la ecuacin (7).

= + + . (24)

Dependiendo de los valores de l y M tendremos un tipo de potencial efectivo, y el valor de nos

dar el tipo de rbita posible, de forma anloga al caso Newtoniano.

Si definimos las variables adimensionales, = = , el potencial efectivo puede escri-
birse como:

( ) = + . (25)

Por lo que el potencial efectivo depende slo del parmetro . Podemos encontrar varios tipos:
V

0.15

0.10

0.05

r M
20 40 60 80

0.05

0.10

0.15

0.20
Fig. 2. En este grfico representamos en el eje de ordenadas el potencial efectivo y en
el de abscisas r / M. En rojo representamos un potencial con = 2; en azul un po-
tencial con un punto de inflexin con = 12, cuyo punto de inflexin est en
= 6; en negro, con = 4, el potencial presenta un mximo en cero; y en verde
= 5, en el que el mximo del potencial es ya positivo.

Las posibles rbitas dependen, por tanto, del tipo de potencial y del valor de .
Para estudiar el potencial efectivo, es interesante encontrar sus mximos y mnimos. Para ello deri-
vamos el potencial efectivo:
3
Este trmino solo es importante para M muy grandes o para r muy pequeas. En el caso del Sol el trmino llegara a
ser importante para distancia de pocos kilmetros, y teniendo en cuenta que el radio del Sol es 696000 Km, esas dis-
tancias no son alcanzables y por esa razn la teora Newtoniana funciona tan bien en el sistema solar.
7

( )
= + ,
y lo igualamos a cero:

1 3


+
=0

+3 = 0,

y obtenemos dos soluciones posibles:

= 1 1 , (26)

que son los mximos y mnimos de potencial.

Cuando = 12 hay una solucin doble y tenemos un punto de inflexin, en el que = 6,
que est representado en color azul en la figura 2.
Un caso interesante del potencial efectivo se da cuando el valor mximo es nulo. Lo podemos ver
en la figura 2, en negro, y sucede cuando = 4.
Todo lo anterior podramos realizarlo con la variable r en vez de , obteniendo:

= 1 1 12
2

=6 con = 12 =
=4 con =4 y ( ) = 0.

El caso ms simple lo encontramos para l =0, en el que tenemos el problema de cada libre, ya que
al ser = 0 la geodsica seria radial. Este caso lo trataremos de forma diferenciada ms adelante.

Casos posibles para energas negativas

Claramente la geometra queda definida por las races de la ecuacin ( ) = 0, que son tres.

Esta situacin que es ms inmediata de comprobar si realizamos el cambio = , quedndonos:


+ + = 0 , (27)

que es un polinomio de tercer orden y, por tanto, tiene tres races, pudiendo ser las tres reales o una
real y dos complejas conjugadas. Pero como es mas intuitivo usando la variable r, seguiremos la
busqueda de las races de

+ + =0. (28)

Podemos ver en la figura 3 los posibles casos en los que nos podemos encontrar:
 8

Races reales Caso

, , a
= , b
, = c
= = d
e
f

Suponemos que el orden de las races reales, que son positivas, es:

. (29)

Caso a
Hay dos posibles rbitas:

Una rbita ligada que oscilara entre r2 y r3. Tendra su analoga Newtoniana en las rbitas
elpticas, aunque en este caso no son realmente elpticas, como podremos observar en las fi-
guras 8, 9, 10 y 11.
Una rbita que empieza en r1 y termina cayendo en el origen.
Caso b
Tenemos en este caso tambin dos tipos de rbitas:
Una rbita circular inestable en r1. Sabemos que es inestable porque al ser un mximo del
potencial cualquier perturbacin en la energa lo alejara de una rbita ligada, esto podemos
comprobarlo en el clculo numrico dejando que pasen suficientes iteraciones para que el
error de clculo nos aleje del mximo y vemos que es inestable como podemos observar en
la figura 13.
Una rbita cuya forma recuerda a la espiral logartmica, desde r3 asinttica a r1.
Caso c
De nuevo podemos considerar dos tipos de rbitas:
Una rbita circular estable en r3 que tiene un radio mayor que la anterior rbita circular ines-
table y que coincide con el valor mnimo del potencial.
La segunda es una rbita que empieza en r1 y cae al origen.

Caso d
En este caso las tres races coinciden y tenemos una rbita estable circular en el valor de la raz, que
es la ltima rbita circular estable posible. Tambin hay rbitas de cada al origen. El radio de esta
rbita es 6M y el punto triple aparece cuando = 12 en el potencial efectivo. Esta raz coincide
con el punto de inflexin r = 6M, que normalmente en la literatura se denomina r , haciendo
referencia a Innermost Stable Circular Orbit.

Caso e
Slo hay una raz real positiva y la nica rbita posible es cada en espiral a la singularidad.

Caso f
En este caso solo son posibles rbitas de cada en el origen. El potencial efectivo en este caso no
posee mximos ni mnimos relativos, esta situacin se da cuando < 12.
 9

V
rM
V 10 20 30 40 50
rM
5 10 15 20 25 30 35

0.05
0.05

0.10
0.10

0.15
0.15

0.20 0.20

(a) (b)
V V
rM rM
5 10 15 20 25 30 5 10 15 20 25 30

0.05
0.05

0.10
0.10

0.15

0.15

0.20

(c) 0.20
(d)
V V
r M rM
10 20 30 40 50 60 10 20 30 40 50 60

0.05 0.05

0.10 0.10

0.15
0.15

0.20
0.20
(e) (f)

Fig. 3. Observamos diferentes soluciones para la ecuacin = 0 con < 0 :

(a) Hay tres races reales negativas. Coincide con el caso equivalente a rbitas Keplerianas.
(b) Una de las races es doble y coincide con el mximo del potencial efectivo. Nos lleva a una rbita cir-
cular inestable.
(c)Tenemos una raz doble, que en este caso coincide con el mnimo del potencial efectivo y da una rbita
circular estable.
(d) Es precisamente cuando el potencial tiene un punto de inflexin y la raz triple coincide con este punto,
dndose una rbita circular estable, la de menor r posible.
(e)Slo hay una raz real.
( f) Es para valores de < 12 , en donde no hay ni mximos ni mnimos relativos.
 10

Casos posibles para energas positivas

Consideramos potenciales efectivos con > 4 y, por tanto, que tiene valores en el primer cua-
drante, como podemos ver en la figura 4.
Vemos que en este caso puede haber races negativas y no hay posibilidad de rbitas circulares es-
tables, ni de rbitas ligadas distintas de la circular inestable.
Tenemos, pues, tres casos posibles, que son:

V
0.5

0.4

0.3

0.2

0.1

rM
40 20 20 40

0.1
Fig. 4. El potencial tiene un mximo positivo. Como vemos en la grfica tambin
tiene valores negativos, lo que implica que existir una raz negativa para
= 0.

V
0.5

0.4

0.3

0.2

0.1

r M
40 20 20 40

0.1
Fig. 5. (A) Hay tres races reales: 2 positivas y una negativa, que no tiene sentido fsico
11

V
0.5

0.4

0.3

0.2

0.1

rM
40 20 20 40

0.1
Fig. 6. (B) Hay tres races reales: 2 positivas coincidentes y una negativa que
no tiene sentido fsico
.
V
0.5

0.4

0.3

0.2

0.1

rM
40 20 20 40

0.1
Fig. 7. (C) Hay una raz real negativa que no tiene sentido fsico y dos com-
plejas conjugadas.

Caso A

La ecuacin ( ) = 0 tiene dos races reales positivas y una negativa que carece de sentido
fsico.
Las rbitas posibles son: una rbita de dispersin, que equivaldra a las trayectorias hiperblicas
Newtonianas, y una rbita de cada en la singularidad.

Caso B

Las dos races reales coinciden en un mximo del potencial efectivo y llevan a una rbita circular
inestable. Adems aparecen rbitas de cada en el origen.

Caso C

Hay dos races complejas conjugadas adems de la consabida solucin negativa. Las rbitas posi-
bles son de cada en la singularidad.
 12

Estudio numrico

Para poder representar las rbitas anteriormente comentadas, podemos resolver las ecuaciones dife-
renciales analticamente usando las aproximaciones oportunas o intentar resolverlas numricamente.
ste ltimo es el mtodo que me planteo en este proyecto, para lo cual parto de las siguientes ecua-
ciones:

Mtodo 1

Partimos de = + +

y derivamos respecto al tiempo propio

= + . (30)

Despejando de (6):

= . (31)

De (2):

1 + 1 = 1.

Deducimos:
1
= + +
2 2 2
1 1 1

= + +( )
. (32)

Las ecuaciones (30), (31) y (32) las he usado con las condiciones iniciales:

(0) =
(0) = 04
(0) = 0.

Este sistema de ecuaciones presenta el problema de divergencia para = 2 , es decir, para = ,

siendo = 2 el radio de Schwarzschild, que es el lmite del horizonte de sucesos y el tamao
mximo que debe tener una estrella para considerarse agujero negro.

4
El tomar esta condicin inicial nos restringe simplemente a un punto concreto de la rbita, facilitando el clculo
numrico, sin menoscabar la generalizacin de las soluciones.
13

Mtodo 2
Otra forma de integracin numrica es partir de = + + y usar = y
(31). Obtenemos:

= + + . (33)

Pero es ms conveniente para la resolucin de la ecuacin pasar a la variable = , de donde

( / )
= = .

(33) queda de la forma:

= + + .

Reordenando trminos:

+ 2 =0, (34)

con la condicin inicial

(0) = .

Mtodo 3
Partimos de (34) y derivamos respecto de :

2 +2 6 = 0.

Suponiendo que 0, es decir, que la rbita no es circular, dividimos por 2 :

+ 3 = , (35)

con las condiciones iniciales

(0) =
(0) = 0

Los tres mtodos los he usado para representar la trayectoria de la partcula de prueba en una mtri-
ca de Schwarzschild, dando resultados idnticos a excepcin del primer mtodo en = .
Para la integracin numrica me he ayudado de la herramienta informtica Wolfram Mathemtica
8.0.

Particularidades de la mtrica de Schwarzschild

Hemos de tener en cuenta las posibles singularidades del problema, en nuestro caso de la propia
mtrica:
 14

= 1 + + 1 ,

que se puede reescribir:

= 1 + + 1 (36)

Con = 2 , que es el llamado radio de Schwarzschild, vemos que la mtrica tiene dos singulari-
dades: una en r = 0 y otra en = . Para comprobar si son singularidades fsicas o aparecen por la
eleccin del sistema de coordenadas, podemos buscar un escalar, que no depende de la eleccin de
las coordenadas, y ver si es singular respecto a este punto. sta es una condicin necesaria para que
un punto sea singular, pero no suficiente [21]. En nuestro caso, para la mtrica de Schwarzschild po-
demos calcular el escalar:

= , (37)
(Vase bibliografa [21])

que nos muestra que r = 0 es una singularidad fsica. Pero, sin embargo, se puede comprobar que
ninguno de los invariantes de curvatura diverge en = , por lo que esta singularidad en ese punto
no es realmente tal sino que es debida a la eleccin del sistema de coordenadas y con otra eleccin
se puede salvar.

El lmite = es un horizonte de sucesos que indica que lo que suceda tras ese punto no es obser-
vable por puntos en el exterior. Esto ha querido ampliarse en general con la llamada conjetura de
censura csmica, formulada por Roger Penrose en 1969, en la que se supone que un horizonte de
sucesos ocultara las singularidades fsicas a los observadores exteriores al mismo para evitar el
problema de causalidad que pueden crear tales singularidades. Pero hasta el momento es una conje-
tura y no se ha demostrado en caso general.

Vamos a mostrar un conjunto de rbitas representativas tanto para energas negativas como positi-
vas. Todas ellas han sido obtenidas mediante integracin numrica con alguno de los tres mtodos
mostrados anteriormente. En algunas grficas, donde la trayectoria no atraviesa el radio de Sch-
warzschild, el interior del crculo de radio ha sido pintado en negro.
15

Empezamos mostrando las rbitas con energa negativa.

rbitas para < 0
1000

40

500
20

1000 500 500 1000 40 20 20 40

20
500

40

1000

Fig. 8. l=10, r0=1000, M = 0.5, = 0.

Representa un campo dbil5 con M= 0.5 y con Fig. 9. l=4, r0=35, M = 0.5, = 0.
r0=1000. En ella podemos observar que las rbitas se Observamos el mismo efecto que en fig. 8 para una rbita
acercan a las Keplerianas pero con una precesin6 menos excntrica.
del perihelio de la rbita.

40
30

20 20

10

40 20 20 40

30 20 10 10 20 30
20

10

40
20

Fig. 10. l=4, r0=35, M = 0.8, = 0.

El campo no es dbil por lo que observamos que la 30

rbita no se asemeja a una elipse con una precesin Fig. 11. l=19/5, r0=6, M = 1, = 0
de su perihelio. En el caso del Sistema Solar para
poder sentir un campo tan fuerte deberamos estar Vemos de forma ms llamativa la rbita ligada de tipo
en una rbita menor al radio solar, cosa en todo (a) de la figura 3, en la que existe un radio mnimo y
imposible. La rbita representada se observara en otro mximo entre los cuales discurre la trayectoria.
estrellas cuya relacin entre su masa y radio fue-
7
ra suficientemente pequea .

5
An con estos datos, el campo es mucho ms fuerte que en el caso del Sol y Mercurio. M = 2.96 Km y r0 = 58000000 Km.
6
En el caso de Mercurio la precesin es de 42.9 segundos de arco por siglo.
7
En este caso = 2 = 1.6 y = 35 = 21.875 . Es decir, en estrellas cuyo radio sea inferior a 22 pueden existir rbitas
de ese tipo, en las que se sienta este campo gravitatorio tan fuerte o, en trminos ms precisos, que el espacio tiempo se haya de-
formado fuertemente (recordemos que en Teora General de la Relatividad no existe campo gravitatorio).
 16

6
4

6 4 2 2 4 6 4 2 2 4

4
6

Fig. 12. l=19/5, r0=rmax, M = 1, = 0.

Podemos observar una trayectoria circular inesta- Fig. 13. l=19/5, r0=rmax, M = 1, = 0 .
Trayectoria circular inestable en la que una ligera per-
ble, que sera representativa del caso (b). Represen- turbacin, en este caso de clculo numrico, hace que
tamos por un disco negro el crculo interior del ra- la rbita se precipite al origen.
dio de Schwarzschild.
15

10 5

15 10 5 5 10 15 5 5

10 5

15

Fig. 15. l= 12, r0=rinf=6, M = 1, = 0.

Fig. 14. l=19/5, r0=rmin, M = 1, = 0.
Observamos otra vez una rbita circular estable, pero
Representa una rbita circular estable cuyo radio
en este caso es la de radio8 ms pequeo posible y
es el mnimo del potencial efectivo. Por tanto, nos
coincide concretamente con el punto de inflexin del
encontramos en el caso (c).
potencial efectivo para = 12, que es el caso (d).

8
Ya mencionamos que a este radio se le denomina rISCO. El poner entre comillas la palabra radio se debe a que real-
mente en la mtrica de Schwarzschild r no coincide con la distancia al centro, aunque s en el lmite Newtoniano.
17

30 30

20 20

10 10

30 20 10 10 20 30 30 20 10 10 20 30

10 10

20 20

30 30

Fig. 16. l=19/5, r0=40, M = 1, = 0. Fig. 17. l=19/5, r0=40, M = 1, = 0.

Las figuras 16 y 17 corresponden con el caso (e), en el que la energa es mayor que el valor mximo del potencial.
En este caso no hay rbitas ligadas y aparece una rbita que cae a la singularidad.
40 20

20 10

40 20 20 40 20 10 10 20

20 10

40 20

Fig. 18. l=2, r0=40, M = 1, = 0. Fig. 19 l=2, r0=40, M = 1, = 0.

Estas dos figuras representan el caso (f), en el que el potencial efectivo no tiene ni mximos ni mnimos y no exis-
ten rbitas ligadas.
6 6

4 4

2 2

6 4 2 2 4 6 6 4 2 2 4 6

2 2

4 4

6 6

Fig. 20. l=19/5, r0= 3.50467, M = 1, = 0. Fig. 21. l=19/5, r0=3.5046, M = 1, = 0.

Estas dos figuras representan el caso (a) con rbitas no ligadas de cada en el centro.
 18

15

30

10
20

5
10

30 20 10 10 20 30 15 10 5 5 10 15

10
5

20
10

30

15

Fig. 22. l=19/5, r0=r3=33.7328, M = 1, = 0. Fig. 23. l=19/5, r0=r3=33.7328, M = 1, = 0.

15

10

15 10 5 5 10 15

10

15

Fig. 24. l=19/5, r0=r3=33.7328, M = 1, = 0.

En las figuras 22, 23 y 24 observamos trayectorias no ligadas, que tambin ocurren en el caso (b).
19

rbitas para > 0

Para energas positivas, como ya vimos, existen tres tipos de rbitas que quedan representadas por
las figuras de abajo:

30 4

20

10

30 20 10 10 20 30 4 2 2 4

10

20

30 4

Fig. 25. l=6, r0=5.45203 ,E=0.2 M=1. Fig. 26. l=6, r0=2.541,E=0.2 M=1.
Representa una rbita del tipo (A), que consideramos de Representa otra rbita del tipo (A), en la que la part-
dispersin. Es la equivalente a las rbitas hiperblicas cula cae en la singularidad. En este grfico aparece en
Newtonianas. negro el interior del horizonte.
30
4

20

2
10

30 20 10 10 20 30
4 2 2 4

10

2
20

30
4

Fig. 27. l=6, r0=2.541,E=0.2 M=1. Fig. 28. l=6, r0= 18 1 = ,E= 0.348, M=1.
Representa una rbita del tipo (A), en la que la partcula Est representado el caso (B) de rbita circular ines-
cae en la singularidad. table, en la que acaba escapando de la rbita cerrada
debido a perturbaciones numricas.
 20

4
4

2 2

4 2 2 4 4 2 2 4

2 2

4
4

Fig. 29. l=6, r0= 18 1 = ,E= 0.348, M=1. Fig. 30. l = 6, M =1, r0 = 5,E = 0.4.
Se representa exclusivamente la rbita circular inestable Esta figura corresponde al caso (C), de cada a la
y, en negro, el interior del horizonte de sucesos. singularidad.

Geodsica radial

Cuando el momento angular es cero, lo que aparece es una geodsica radial, que solo depende de r.
Las ecuaciones que gobiernan estas geodsicas son (24) y (5) que en nuestro caso quedan:

= = + +

1 = .

Al eliminar l:

= + 1 (38)
= . (39)

Consideraremos las trayectorias de las partculas que empiezan desde el reposo a una distancia ri y
caen al origen. La distancia inicial est relacionada con e por:

= ( = cuando = 0).

Para obtener las ecuaciones de movimiento es ms conveniente realizar el cambio de variables:

= (1 + cos ) = cos = cos . (40)

Claramente = 0 cuando = . Los valores de al cruzar el horizonte =2 y en la singulari-

dad = 0 son respectivamente:

= =2 cuando =2
21

y
= cuando = 0.

Las ecuaciones (38) y (39) con ese cambio de variable son:

1
= (1 )
2

= , (41)

que juntas dan:

= . (42)

Para partculas en cada libre nos quedamos con el signo negativo

= (1 ) = (43)

Usando las ecuaciones (42) y (43) obtenemos:

= = (1 + cos ) , (44)

de donde:

= ( + sin ) , (45)

donde hemos supuesto que = 0 en = 0. De la ecuacin se sigue que la partcula cruza el hori-
zonte y alcanza la singularidad en un tiempo propio finito:

= ( + sin ) y = . (46)

La situacin en muy diferente cuando consideramos la coordenada temporal t. La ecuacin que

puede ser integrada para obtener t se deduce de (41) y (44), y es:

= . (47)

La integracin de esta ecuacin da:

= ( + sin ) + (1 ) +2 (48)

Conforme a esta ecuacin,

lim =, (49)
 22

que tiene un fuerte contraste con la conducta del tiempo propio . As para un observador estaciona-
rio en el infinito una partcula tarda un tiempo infinito en llegar al horizonte, mientras que para el
tiempo propio la partcula llegara al horizonte en un tiempo finito e incluso llegara a la singulari-
dad tambin en un tiempo finito. Estos hechos se ilustran en la figura 31.

Fig. 31 Representamos la variacin de la coordenada tiempo (t) y

tiempo propio a lo largo de una geodsica radial.

Geodsicas nulas

Vamos a estudiar las rbitas de los rayos luminosos en la mtrica de Schwarzschild, clculo que es
parecido al anterior pero con diferencias importantes.
En el caso de la luz, sta sigue una geodsica nula respecto a algn parmetro :

1 + + 1 =0. (49)

La nica diferencia con (2) es que a la derecha aparece 0 en vez de -1.

Tenemos:
= 1 (50)
= (51)

Anlogamente la rbita es plana y podemos escoger = 2, por lo que (49) puede escribirse:

1 + 1 =0 . (52)

Usando (45) y (46) para eliminar y , tenemos:

1 + 1 =0. (53)
23

1
Multiplicando por , podemos ponerlo en la forma:

= + ( ) , (54)

donde

= (55)

( )= 1 (56)

que es el potencial efectivo para rbitas de fotones.

La ecuacin (54) tiene la forma de una integral de energa para movimiento radial con ( ) ju-
-2
gando el papel de un potencial efectivo y b jugando el papel de la energa. Esta ecuacin puede
usarse para hallar las rbitas de las partculas luminosas de la misma forma que usamos (8) para
analizar las rbitas de las partculas con masa, pero a diferencia del caso anterior las rbitas de los
rayos luminosos dependen slo de la relacin .
El signo de l indica el camino que sigue el rayo luminoso alrededor del centro de atraccin.
El parmetro b es el parmetro de impacto.

Wefe

0.04

0.02

rM
30 20 10 10 20 30

0.02

Fig. 32. Representacin de respecto a r/M. He representado tam-

bin la parte negativa para observar que = tiene una raz nega-
tiva y dos positivas o complejas conjugadas.

Las soluciones de (54) dependen de las races de ( ) = . Como podemos ver en la figura 32
las races posibles son tres: una siempre negativa y las otras dos pueden ser reales positivas o com-
plejas conjugadas. Por lo tanto los casos posibles los podemos resumir en la figura 33.
 24

Wefe

0.04

0.02

rM
30 20 10 10 20 30

0.02

Wefe

0.04

0.02

rM
30 20 10 10 20 30

0.02

Wefe

0.04

0.02

rM
30 20 10 10 20 30

0.02

Fig. 33. En las tres grficas de arriba se representan y

y podemos observar los puntos de corte. De arriba a abajo te-
nemos:
Tres races reales: de ellas una negativa y dos positivas.
Una raz real negativa y una raz doble real positiva. ste
es un caso lmite que nos llevar a una rbita inestable
circular.
Una raz real negativa y dos complejas conjugadas.

Para encontrar las rbitas posibles podemos partir de (54), teniendo en cuenta que = y
(51):

1 1 1 2
= + 1

1 1 1 2
= + 1

= + . (57)

Si usamos el cambio = 1/ :
25

(1/ ) 1
= =
nos queda:

=2 + . (58)

Si derivamos con respecto a tenemos:

=3 . (59)

Tanto la ecuacin (58) como la (59) las podemos utilizar para calcular numricamente las rbitas de
los rayos luminosos. Los casos posibles sern tres y estn caracterizados por las races de ( )=
, que vemos representadas en la figura 33. Debemos pues encontrar el mximo de la funcin
( ), ya que es un punto crtico de las soluciones.

= + =0, (60)

que da como solucin:

=3 . (61)

Si sustituimos este valor en ( )= obtenemos:

1 1 2
= 1
(3 ) 3

= 27 , (62)

que es el valor crtico del parmetro de impacto, pudindose resumir las rbitas en tres tipos:

TIPO
Para > 27 las rbitas son de dispersin.

TIPO
Para = 27 las rbitas son circunferencias inestables.

TIPO
Para < 27 las rbitas son de cada en el centro.
 26

Fig. 34. rbita para el tipo , en el que el rayo luminoso se

ve dispersado por el agujero negro.

4 2 2 4

Fig. 35. rbita inestable que corresponde al caso crtico .

6

6 4 2 2 4 6

Fig. 36. Caso que corresponde a cada en la singularidad, es de-

cir, tipo .
27

Caso de agujeros negros rotantes, mtrica de Kerr

Debemos tener presente que la mtrica de Schwarzschild es vlida para estrellas esfricas y est-
ticas. En realidad esto es una idealizacin, ya que las estrellas no son estticas, sino que giran sobre
s mismas y adems suelen estar achatadas por los polos debido a este giro. El achatamiento no es lo
suficientemente importante para ser necesario considerarlo. Pero no es as en el caso de la rotacin,
que en algunos casos puede ser muy importante, por lo que debemos considerar estrellas esfricas
que giran sobre s mismas. La rotacin de la estrella implica una energa y como en la Teora Gene-
ral de la Relatividad la energa y la masa son en realidad equivalentes, la energa tambin afecta a la
mtrica, por lo que la mtrica en este caso resuelta por Kerr en 1963 cumple el diferencial de lnea:

= 1 + + + + + sin (63)

= (64)
= + cos (65)
= 2 + (66)

Podemos ver que la mtrica es invariante respecto de t e independiente de . Por lo tanto: es esta-
cionario y tiene simetra axial por lo que tendremos dos magnitudes que se conservan.

Adems es invariable ante una reflexin en el plano = 2, lo que significa que va desde a .

Podemos observar que cuando a = 0 obtenemos la mtrica de Schwarzchild, como cabra esperar.

En cuanto a las singularidades de esta mtrica tenemos que se observan cuando 0 y 0.

En el caso de = 0 tenemos que + cos = 0. Por lo tanto, = 0 = 2, que es una

singularidad equivalente a la singularidad r = 0 de la mtrica de Schwarzchild (cuando a = 0).
Por otra parte tenemos que si = 0 entonces:

2 + =0 = (67)
Supondremos que .

Estos radios no son singularidades fsicas, al igual que suceda con r = 2M en la mtrica de Sch-
warzchild, es decir, que podemos encontrar transformaciones de coordenadas para que estos radios
no sean singulares.
Concretamente nos da el horizonte de sucesos del agujero negro de Kerr, y dentro de este hori-
zonte es donde se da la singularidad de = 0.
En la mtrica de Kerr el horizonte r = r+ gira con una velocidad angular = . Adems, como
no es una mtrica esfrica, r = r+ no corresponde con una esfera. Sustituyendo r = r+ en el elemen-
to queda:

= ( ) + sin , (68)
( )

que no es una esfera. Concretamente para el caso extremo a = M, la distancia en el ecuador ( = )

es 4 y la distancia en los polos es < 7.6 M.
 28

Fig. 35. El horizonte de un agujero negro rotante. La figura

muestra una superficie en el espacio tridimensional que
tiene la misma geometra para t constante que el horizonte
de un agujero negro de Kerr. En este caso tomamos a/M =
86. La superficie est aplastada a lo largo del eje de rota-
cin.

Vemos, pues, que a diferencia del caso de Schwarzchild el horizonte de sucesos no es una esfera.
Otra diferencia importante es que las rbitas no son planas en general, aunque s lo son las rbitas
en el plano ecuatorial. En la figura 36 podemos observar dos rbitas ligadas en un agujero de Kerr,
que claramente no son planas; han sido obtenidas directamente por integracin numrica de las
ecuaciones de movimiento.

rbitas en el plano ecuatorial

En general, como ya mencionamos, las rbitas no estn confinadas a un plano porque la magnitud
que se conserva es la componente del momento angular paralela al eje de simetra.

Fig. 36. rbitas ligadas en el caso de un agujero negro de Kerr. Podemos observar que no son planas y que
son ms complicadas que el caso de Schwarzchild. Estn obtenidas con clculo numrico usando Wolfram
Mathemtica 8.0
29

Sin embargo, en el plano ecuatorial = 2 s lo son. Veamos este caso concreto con un poco ms
de detalle. La mtrica restringida a este plano es:

= 1 + + + + . (69)

Las dos magnitudes que se conservan en esta mtrica son e que se puede interpretar como la ener-
ga y l como la componente del momento angular paralela al eje de simetra por unidad de masa en
reposo. En este caso, al ser rbitas ecuatoriales, l es el momento angular total y no una componente.

Concretamente son (vase para su deduccin [1]):

= + (70)
= + . (71)

Y obtenemos:

= + + (72)
= 1 + (73)

= + ( , , ) (74)
( )
( , , )= + . (75)

Para la luz tenemos:

= ( , , ) (76)
( , , )= 1 1 (77)

= = ( ).

Los potenciales efectivos (75) y (77) tienen la misma dependencia con r que en el caso de la geo-
metra de Schwarzschild. Una diferencia importante es que los potenciales son dependientes de la
energa y el momento angular. Por ejemplo, las partculas de prueba o rayos luminosos que vienen
desde el infinito rotando en la misma direccin que el agujero negro (valores positivos de l o ) se
mueven en un diferente potencial efectivo que partculas que giran al contrario (valores negativos
de l o ). Estas diferencias reflejan, en parte, el arrastre espaciotemporal del agujero negro rotante.
Las partculas son arrastradas por esta rotacin.

Vemos, pues, que otra diferencia con la mtrica de Schwarzchild es que la rbita es diferente si gira
en el mismo sentido que el agujero negro, a la que llamaremos rbita corrotante, que si gira en sen-
tido opuesto, llamada contrarrotante.

Otra diferencia importante con la mtrica no rotante, es la aparicin de la llamada ergoesfera, la

cual se define como la zona donde ningn observador, independientemente de su capacidad propul-
sora, puede permanecer inmvil. La deduccin del lmite de la Ergoesfera podemos encontrarlo en
la bibliografa y es:

( )= + . (78)
 30

Claramente ( ) > , por lo que la superficie ( ) que limita la ergoesfera est fuera del hori-
zonte, como muestra la figura 37.

Fig. 37. En esta grfica representamos el lmite de la

ergoesfera ( ) y el del horizonte de sucesos . La
representacin no es fidedigna ya que como podemos
observar el horizonte est representado por un crcu-
lo cuando sabemos que est achatado. Grafico sacado
de [1].

winding number : 47.03045 initial radius : 40.90

winding number: 52.00690 initial radius : 44.02

40
40

20
20

0
0

20
20

40 40

40 20 0 20 40 40 20 0 20 40

Fig. 38. En estas grficas representamos dos rbitas ligadas en el plano ecuatorial de un agujero negro de
Kerr. Ambas han sido obtenidas computacionalmente.
31

Tabla 1 de Unidades Masa-Longitud y Masa-Longitud-Tiempo

Unidades Unidades Conversin MLT

Cantidad Smbolo tpico
ML (c=1) MLT ML
Masa m M M M
Longitud L L L l
Tiempo t L T ct
Distancia espacio-
s L L s
tiempo
Tiempo propio L T c
Energa E M M(L/T)2 Ec2
Momento p M M(L/T) pc
Velocidad V adimensional L/T V/c

Tabla 2 de Unidades geometrizadas y Masa-Longitud-Tiempo

Unidad geometri-
Unidades Conversin
Cantidad Smbolo tpico zada
MLT MLT geometrizada
(c=G=1)
Masa M L M GM/c2
Longitud L L L L
Tiempo t L T ct
Distancia espacio-
s L L s
tiempo
Tiempo propio L T c
Energa E L M(L/T)2 GE/c4
Cantidad de mo-
p L M(L/T) Gp/c3
vimiento
Momento angular J L2 M(L2/T) GJ/c3
Potencia (lumni-
L Adimensional M(L2/T 3) GL/c5
nosa)
Densidad de
L-2 M/(LT2) G /c4
energa
Densidad de mo-
mento (flujo de L-2 M/(L2T) G /c3
energa)
Presin (tensin) p L-2 M/(LT2) Gp/c4
Energa de una
rbita por unidad e Adimensional (L/T)2 e/c2
de masa
Momento angular
de una rbita por L L2/T /
unidad de masa
Constante de
L2 M(L2/T) /
Planck

Para convertir masa en kilogramos a masa en metros usamos la primera fila de la tabla 2 y encontramos que M (en m) =
(G/c2)M (en kilogramos).
Para la conversin inversa a unidades MTL desde cualquiera de los otros dos sistemas, reemplazar las cantidades por la
expresin en la ltima columna con c y G restauradas. Por ejemplo, la ecuacin dada para la velocidad de escape de una
partcula desde una coordenada de radio r de Schwarzschild en el exterior de un agujero negro esfrico es =
(2 / ) / en unidades geometrizadas. Para encontrar la misma relacin en unidades MLT, buscamos en la tabla 1:
debe ser reemplazada por / , y en la tabla 2: M debe ser reemplazado por GM/c2. Por lo que queda:
/ /
= = .
 32

Bibliografa

[1] James B. Hartle, Gravity An Introduction to Einstein's General Relativity, 2003, Pearson Education.
[2] E. Castillo, A. Iglesias, J.M. Gutirrez, E. lvarez y A. Cobo, Mathematica, Paraninfo, 1993.
[3] S. Chandrasekhar, The Mathematical theory of black holes, Oxford University Press, 1998.
[4] M.L. Abell y J.P. Braselton, Differential Equations with Mathematica, Academic Press, 1993.
[5] R.M. Wald, General Relativity, The University of Chicago, 1984.
[6] L. D. Landau y E.M. Lifshitz, Teora Clsica de los Campos, volumen 2 del curso de Fsica Terica, Editorial Revert, S.A., 1992.
[7] A. Einstein, Zur Elektrodynamik bewegter Korper, Annalen der Physik.
[8] A. Einstein, Erklrung der Perihelbewegung des Merkur aus der Allgemeinen Relativittstheorie, Sitzungsberichte der Preussischen
Akademie der Wissenschaften, 1915.
[9] A. Einstein, Kosmologische Betrachtungen zur Allgemeinen Relativitatstheorie , Sitzungsberichte der Preussischen Akademie der
Wissenschaften.
[10] A. Einstein and N. Rosen, The Particle Problem in the General Theory of Relativity, July 1, 1935, Physical Review Volume 48.
[11] K. Schwarzschild, ber das Gravitationfeld eines Massenpunktes nach der Einsteinschen Theorie, Sitzungsberichte der Kniglichen
Preussischen Akademie der Wissenschaften, 1916.
[12] J. Foster, J. D . Nightingale, A short course in General Relativity, Longman Group Limited, 1979.
[13] Wolfgang Rindler, Relativity Special, General, and Cosmological, Oxford University Press, 2006.
[14] C. Misner, K. Thorn, A. Wheeler, Gravitation, Freeman, 1973.
[15] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity, Wiley. 1972.
[16] Frank Y.-H. Wang, Relativistic orbits with computer algebra, American Journal of Physics, August 2004.
[17] D.C. Wilkins, Bound Geodesics in the Kerr Metric, Physical Review D Volume 5, Number 4, 15 February 1972.
[18] H. Goldstein, Mecnica Clsica, Aguilar SA, 1977.
[19] J. Garca-Bellido, Astrophysics and Cosmology, Theoretical Physics Group, Blackett Laboratory, Imperial College of Science.
[20] A.V. Pogorelov, Geometra Diferencial, Ed. Mir, 1984.
[21] S. Carroll, Lecture Notes on General Relativity, http://fr.arxiv.org/PS_cache/gr-qc/pdf/9712/9712019v1.pdf.
[22] R.P. Kerr, Gravitational field of a spinning mass as an example of algebraically special metrics Phys. Rev. Lett. 11 237.
[23] F. de Felice, Bound states in the Kerr metric, Physical Review D Volume 19, Number 2, 15 January 1979.

Otras referencias Bibliogrficas relacionadas con el tema son:

[24] F. De Felice, C. J. S. Clarke, Relativity on curved manifolds, Cambridge University Press, 1990.
[25] E. M. Vargas a y J. M. Tejeiro, Simulacin Numrica de un Agujero Negro de Schwarzschild, Revista Colombiana de Fsica, Vol. 40,
No.2, Julio 2008.
[26] P.A.M. Dirac, General Theory of Relativity, John Wiley & Sons,1975.
[27] B.F. Schutz, A first course in general relativity, Cambridge University Press, 1990.
[28] JeanAlain Marck, Shortcut method of solution of geodesic equations for Schwarzschild black hole, Class.Quant.Grav. 13, 1996.
[29] Janna Levin and Gabe Perez-Giz, A periodic table for black hole orbits, Physical Review D 77, 103005, 2008.
[30] J.V. Narlikar, An introduction to cosmology, Third edition, Cambridge Univ. Press, Cambridge, UK, 2002.
[31] William C. Daywitt, The Newton and Schwarzschild Geodesics Compared.
[32] Hans-Jrgen Schmidt, Perihelion advance for orbits with large eccentricities in the Schwarzschild black hole, Physical Review D 83,
124010, 2011.
[33] G V Kraniotis, Precise relativistic orbits in Kerr and Kerr(anti) de Sitter spacetimes, Class. Quantum Grav. 21, 2004.
[34] W Robin, Schwarzschild relativistic orbits, Eur. J. Phys. 12 204, 1991.
[35] J.V. Narlikar, Introduction to cosmology, Cambrige University Press, 1993.
[36] J, L. SYNGK, Angular Momentum, Mass-Center and the Inverse Square Law in Special Relativity, Physical Review, May 15, 1935.
[37] C. Fronsdal, Completion and Embedding of the Schwarzschild Solution, Physical Review volume 116, number 3, 1959.
[38] Arthur Komar, Covariant Conservation Laws in General Relativity, Physical Review Volume 113, Number 3, 1959.
[39] Karel V. Kuchar, Geometrodynamics of Schwarzschild black holes, Physical Review D Volume 50, Number 6, 1994.
[40] S.W.Hawking & G. F.R. Ellis, The large scale structure of space-time, Cambridge University Press, 1973.
[41] M. D. Kruskal, Maximal Extension of Schwarzschild Metric, Physical Review Volume 119, Number 5, 1960.
[42] L. McFadden, P.R. Weissman y T.V. Johnson, Encyclopedia of the Solar System, 2007.
[43] G. Szekeres, New Formulation Of The General Theory of Relativity, Physical Review Volume 97, Number 1, 1955.
[44] David Finkelstein, Past-Future Asymmetry of the Gravitational Field of a Point Particle, Physical Review Volume 110, Number 4,
1958.
[45] B. P. Jensen and P. Candeias, Schwarzschild Radial Functions, Physical Review D Volume 33, Number 6, 1986.
[46] Z. Stuchlk and S. Hledk, Some properties of the Schwarzschildde Sitter and Schwarzschildanti-de Sitter spacetimes, Physical Re-
view D, Volume 60, 044006.
[47] Akira Tomimatsu, New Exact Solution for the Gravitational Field of a Spinning Mass, Physical Review Letters volume 29, number
19, 1972.
[48] Bert Janssen, Teora de la Relatividad General, Dpto de Fsica Terica y del Cosmos, Universidad de Granada, 2010.
[49] Fernando Chamizo Lorente, Seminario relatividad: La odisea continua.
33

[50] Frank Y.-H. Wanga, Relativistic orbits with computer algebra, American Association of Physics Teachers, 2004.
[51] Matthew R. Francis and Arthur Kosowsky, Geodesics in the Generalized Schwarzschild Solution, Am.J.Phys. 72, 2004.
[52] Eduardo Rodrguez S., La Solucin de Schwarzschild: Agujeros Negros no rotantes, Universidad de Concepcin, 6 de septiembre de
2002.
[53] Javier Bussons Gordo, http://webs.um.es/bussons/indexGyA2009_10.html.
[54] H. Falcke y F.W. Hehl,The.galactic.black.hole. Lectures on general relativity and astrophysics, IOP, 2003.
[55] N. Straumann, General relativity and relativistic astrophysics, Spriger-Verlag, 1984.
[56] E.F. Taylor y J. A. Wheeler, Exploring Black Holes Introduction to General Relativity, Addison Wesley Longman, 2000.
[57] Kraniotis G V and Whitehouse, Compact calculation of the perihelion precession of Mercury in general relativity, the cosmological
constant and Jacobis inversion problem Class. Quantum Grav. 20 481735.
[58] Boyer R H and Lindquist R W, Maximal analytic extension of the Kerr metric, J. Math. Phys. 8 26581.
[59] Boyer R H and Price T G, An interpretation of the Kerr metric in general relativity, Proc. Camb. Phil. Soc. 61 531.
[60] Stuchlik Z and Calvani M, Null geodesics in black hole metrics with nonzero cosmological constant, Gen. Rel. Grav. 23 50719.
[61] Carter B, HamiltonJacobi and Schrodinger separable solutions of Einsteins equations Commun, Math. Phys. 10 280310.
[62] Carter B, Global structure of the Kerr family of gravitational fields, Phys. Rev. 174 155971.
[63] Wilkins D C, Bound geodesics in the Kerr metric, Phys. Rev. D 5 81422.
[64] Wiles A, Modular elliptic curves and Fermats last theorem, Ann. Math. 141 443551.
[65] Byrd Gene G, Chernin A.D., Valtonen M.J., Cosmology, URSS, 2007.
[66] Andrew Liddle, An Introduction to Modern Cosmology, Wiley, 2003.
[67] M.A. Hernndez Cifre y J. A. Pastor Gonzlez ,Un curso de Geometra diferencial , Consejo Superior de Investigaciones Cientficas,
Madrid, 2010.
[68] L.P. Hughston, K.P. Tod, An introduction to General Relativity. London Mathematical Society Student Texts.
[69] Stefan Waner, Introduction to Differential Geometry and General Relativity,
http://people.hofstra.edu/Stefan_Waner/RealWorld/pdfs/DiffGeom.pdf.
[70] Duncombe, R. L., Relativity effects for the three inner planets, Astronom. J., 61, 174-5.
[71] Frank Wang, rbita relativista con lgebra computacional, American Journal of Physics, 72, 2004.
[72] R. W. Hamming. Mtodos Numricos para cientficos e ingenieros, Segunda Edicin, Dover Publications, Nueva York, 1973.
[73] N. A. Sharp, General Relativity and Gravitation, 10, 659-70, 1979.