Geodésicas de Schwarzschild

Caminos de partículas en la solución de Schwarzschild a las ecuaciones de campo de Einstein

En relatividad general , las geodésicas de Schwarzschild describen el movimiento de partículas de prueba en el campo gravitacional de una masa fija central , es decir, el movimiento en la métrica de Schwarzschild. Las geodésicas de Schwarzschild han sido fundamentales en la validación de la teoría de la relatividad general de Einstein . Por ejemplo, proporcionan predicciones precisas de la precesión anómala de los planetas del Sistema Solar y de la desviación de la luz por la gravedad. ${\textstyle M,}$

Las geodésicas de Schwarzschild se refieren únicamente al movimiento de partículas de masas tan pequeñas que contribuyen poco al campo gravitacional. Sin embargo, son muy precisos en muchos escenarios astrofísicos, siempre que sean muchas veces más pequeñas que la masa central , por ejemplo, en el caso de los planetas que orbitan alrededor de su estrella. Las geodésicas de Schwarzschild también son una buena aproximación al movimiento relativo de dos cuerpos de masa arbitraria, siempre que la masa de Schwarzschild se establezca igual a la suma de las dos masas individuales y . Esto es importante para predecir el movimiento de las estrellas binarias en la relatividad general. ${\textstyle m}$ ${\textstyle M}$ ${\textstyle M}$ ${\textstyle m_{1}}$ ${\textstyle m_{2}}$

Contexto histórico

La métrica de Schwarzschild lleva el nombre de su descubridor Karl Schwarzschild , quien encontró la solución en 1915, sólo aproximadamente un mes después de la publicación de la teoría de la relatividad general de Einstein. Fue la primera solución exacta de las ecuaciones de campo de Einstein además de la solución trivial del espacio plano .

En 1931, Yusuke Hagihara publicó un artículo que mostraba que la trayectoria de una partícula de prueba en la métrica de Schwarzschild se puede expresar en términos de funciones elípticas . ^[1]

Samuil Kaplan en 1949 ha demostrado que existe un radio mínimo para que la órbita circular sea estable en la métrica de Schwarzschild. ^[2]

Métrica de Schwarzschild

Una solución exacta a las ecuaciones de campo de Einstein es la métrica de Schwarzschild , que corresponde al campo gravitacional externo de un cuerpo de masa esféricamente simétrico, no giratorio y sin carga . La solución de Schwarzschild se puede escribir como ^[3] ${\textstyle M}$

${\displaystyle ds^{2}=c^{2}{d\tau }^{2}=\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}dt^{2}-{\frac {dr^{2}}{1-{\frac {r_{\text{s}}}{r}}}}-r^{2}(d\theta ^{2}+\sin ^{2}\theta \,d\varphi ^{2})}$

dónde

${\displaystyle \tau }$, en el caso de una partícula de prueba de pequeña masa positiva, es el tiempo propio (tiempo medido por un reloj que se mueve con la partícula) en segundos,

${\displaystyle c}$es la velocidad de la luz en metros por segundo,

${\displaystyle t}$es, para , la coordenada de tiempo (tiempo medido por un reloj estacionario en el infinito) en segundos, ${\displaystyle r>r_{\text{s}}}$

${\displaystyle r}$es, para , la coordenada radial (circunferencia de un círculo centrado en la estrella dividida por ) en metros, ${\displaystyle r>r_{\text{s}}}$ ${\displaystyle 2\pi }$

${\displaystyle \theta }$es la colatitud (ángulo desde el norte) en radianes,

${\displaystyle \varphi }$es la longitud en radianes, y

${\displaystyle r_{\text{s}}}$es el radio de Schwarzschild del cuerpo masivo (en metros), que está relacionado con su masa por ${\textstyle M}$

${\displaystyle r_{\text{s}}={\frac {2GM}{c^{2}}},}$

¿Dónde está la constante gravitacional ? La teoría clásica newtoniana de la gravedad se recupera en el límite cuando la relación llega a cero. En ese límite, la métrica vuelve a la definida por la relatividad especial . ${\textstyle G}$ ${\textstyle {\frac {r_{\text{s}}}{r}}}$

En la práctica, esta proporción es casi siempre extremadamente pequeña. Por ejemplo, el radio de Schwarzschild de la Tierra es aproximadamente 9 mm ( 3 ⁄ 8  pulgadas); en la superficie de la Tierra, las correcciones a la gravedad newtoniana son sólo una parte entre mil millones. El radio de Schwarzschild del Sol es mucho mayor, aproximadamente 2953 metros, pero en su superficie la proporción es aproximadamente de 4 partes en un millón. Una estrella enana blanca es mucho más densa, pero incluso aquí la proporción en su superficie es de aproximadamente 250 partes en un millón. La proporción sólo aumenta cerca de objetos ultradensos como las estrellas de neutrones (donde la proporción es aproximadamente del 50%) y los agujeros negros . ${\textstyle r_{\text{s}}}$ ${\textstyle {\frac {r_{\text{s}}}{r}}}$

Órbitas de partículas de prueba.

Comparación entre la órbita de una partícula de prueba en el espacio-tiempo newtoniano (izquierda) y Schwarzschild (derecha); nótese la precesión absidal a la derecha.

Podemos simplificar el problema usando la simetría para eliminar una variable de la consideración. Dado que la métrica de Schwarzschild es simétrica respecto a , cualquier geodésica que comience a moverse en ese plano permanecerá en ese plano indefinidamente (el plano es totalmente geodésico ). Por lo tanto, orientamos el sistema de coordenadas de modo que la órbita de la partícula se encuentre en ese plano, y fijamos la coordenada para que la métrica (de este plano) se simplifique a ${\textstyle \theta ={\frac {\pi }{2}}}$ ${\textstyle \theta }$ ${\textstyle {\frac {\pi }{2}}}$

${\displaystyle c^{2}d\tau ^{2}=\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}dt^{2}-{\frac {dr^{2}}{1-{\frac {r_{\text{s}}}{r}}}}-r^{2}d\varphi ^{2}.}$

Se pueden identificar dos constantes de movimiento (valores que no cambian en el tiempo adecuado ) (consulte la derivación que se proporciona a continuación). Uno es la energía total : ${\displaystyle \tau }$ ${\textstyle E}$

${\displaystyle \left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}={\frac {E}{mc^{2}}}.}$

y el otro es el momento angular específico :

${\displaystyle h={\frac {L}{\mu }}=r^{2}{\frac {d\varphi }{d\tau }},}$

donde es el momento angular total de los dos cuerpos y es la masa reducida . Cuando , la masa reducida es aproximadamente igual a . A veces se supone que . En el caso del planeta Mercurio, esta simplificación introduce un error más del doble que el efecto relativista. Cuando se habla de geodésicas, se las puede considerar ficticias, y lo que importa son las constantes y . Para cubrir todas las geodésicas posibles, debemos considerar casos en los que sea infinita (dando trayectorias de fotones ) o imaginaria (para geodésicas taquiónicas ). Para el caso fotónico, también necesitamos especificar un número correspondiente a la relación de las dos constantes, es decir , que puede ser cero o un número real distinto de cero. ${\textstyle L}$ ${\textstyle \mu }$ ${\textstyle M\gg m}$ ${\textstyle m}$ ${\textstyle m=\mu }$ ${\textstyle m}$ ${\textstyle {\frac {E}{m}}}$ ${\textstyle h}$ ${\textstyle {\frac {E}{m}}}$ ${\textstyle {\frac {mh}{E}}}$

Sustituyendo estas constantes en la definición de la métrica de Schwarzschild

${\displaystyle c^{2}=\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}\left({\frac {dt}{d\tau }}\right)^{2}-{\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\left({\frac {dr}{d\tau }}\right)^{2}-r^{2}\left({\frac {d\varphi }{d\tau }}\right)^{2},}$

produce una ecuación de movimiento para el radio en función del tiempo propio : ${\textstyle \tau }$

${\displaystyle \left({\frac {dr}{d\tau }}\right)^{2}={\frac {E^{2}}{m^{2}c^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left(c^{2}+{\frac {h^{2}}{r^{2}}}\right).}$

La solución formal a esto es

${\displaystyle \tau =\int {\frac {dr}{\pm {\sqrt {{\frac {E^{2}}{m^{2}c^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left(c^{2}+{\frac {h^{2}}{r^{2}}}\right)}}}}.}$

Tenga en cuenta que la raíz cuadrada será imaginaria para las geodésicas taquiónicas.

Usando la relación más arriba entre y , también podemos escribir ${\textstyle {\frac {dt}{d\tau }}}$ ${\textstyle E}$

${\displaystyle t=\int {\frac {dr}{\pm c\left(1-{\frac {r_{\text{s}}}{r}}\right){\sqrt {1-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left(c^{2}+{\frac {h^{2}}{r^{2}}}\right){\frac {m^{2}c^{2}}{E^{2}}}}}}}.}$

Dado que asintóticamente el integrando es inversamente proporcional a , esto muestra que en el marco de referencia si se acerca lo hace exponencialmente sin llegar nunca a alcanzarlo. Sin embargo, en función de , sí alcanza . ${\textstyle r-r_{\text{s}}}$ ${\textstyle r,\theta ,\varphi ,t}$ ${\textstyle r}$ ${\textstyle r_{\text{s}}}$ ${\textstyle \tau }$ ${\textstyle r}$ ${\textstyle r_{\text{s}}}$

Las soluciones anteriores son válidas mientras el integrando sea finito, pero una solución total puede involucrar dos o una infinidad de piezas, cada una descrita por la integral pero con signos alternos para la raíz cuadrada.

Cuando y , podemos resolver para y explícitamente: ${\textstyle E=mc^{2}}$ ${\textstyle h=0}$ ${\textstyle t}$ ${\textstyle \tau }$

${\displaystyle {\begin{aligned}t&={\text{constant}}\pm {\frac {r_{\text{s}}}{c}}\left({\frac {2}{3}}\left({\frac {r}{r_{\text{s}}}}\right)^{\frac {3}{2}}+2{\sqrt {\frac {r}{r_{\text{s}}}}}+\ln {\frac {\left|{\sqrt {\frac {r}{r_{\text{s}}}}}-1\right|}{{\sqrt {\frac {r}{r_{\text{s}}}}}+1}}\right)\\\tau &={\text{constant}}\pm {\frac {2}{3}}{\frac {r_{\text{s}}}{c}}\left({\frac {r}{r_{\text{s}}}}\right)^{\frac {3}{2}}\end{aligned}}}$

y para geodésicas fotónicas ( ) con momento angular cero ${\textstyle m=0}$

${\displaystyle {\begin{aligned}t&={\text{constant}}\pm {\frac {1}{c}}\left(r+r_{\text{s}}\ln \left|{\frac {r}{r_{\text{s}}}}-1\right|\right)\\\tau &={\text{constant}}.\end{aligned}}}$

(Aunque el tiempo adecuado es trivial en el caso fotónico, se puede definir un parámetro afín y luego la solución a la ecuación geodésica es ). ${\textstyle \lambda }$ ${\textstyle r=c_{1}\lambda +c_{2}}$

Otro caso solucionable es aquel en el que y y son constantes. En el volumen donde esto da para el momento adecuado. ${\textstyle E=0}$ ${\textstyle t}$ ${\textstyle \varphi }$ ${\textstyle r<r_{\text{s}}}$

${\displaystyle \tau ={\text{constant}}\pm {\frac {r_{\text{s}}}{c}}\left(\arcsin {\sqrt {\frac {r}{r_{\text{s}}}}}-{\sqrt {{\frac {r}{r_{\text{s}}}}\left(1-{\frac {r}{r_{\text{s}}}}\right)}}\right).}$

Esto está cerca de soluciones pequeñas y positivas. El exterior de la solución es taquiónico y el "tiempo adecuado" es similar al espacio: ${\textstyle {\frac {E^{2}}{m^{2}}}}$ ${\textstyle r_{\text{s}}}$ ${\textstyle E=0}$

${\displaystyle \tau ={\text{constant}}\pm i{\frac {r_{\text{s}}}{c}}\left(\ln \left({\sqrt {\frac {r}{r_{\text{s}}}}}+{\sqrt {{\frac {r}{r_{\text{s}}}}-1}}\right)+{\sqrt {{\frac {r}{r_{\text{s}}}}\left({\frac {r}{r_{\text{s}}}}-1\right)}}\right).}$

Esto está cerca de otras soluciones taquiónicas con pequeñas y negativas. La geodésica taquiónica constante en el exterior no continúa con una geodésica constante en el interior , sino que continúa en una "región exterior paralela" (ver coordenadas Kruskal-Szekeres ). Otras soluciones taquiónicas pueden entrar en un agujero negro y volver a salir a la región exterior paralela. La solución constante dentro del horizonte de sucesos ( ) continúa con una solución constante en un agujero blanco . ${\textstyle {\frac {E^{2}}{m^{2}}}}$ ${\textstyle t}$ ${\textstyle r_{\text{s}}}$ ${\textstyle t}$ ${\textstyle r_{\text{s}}}$ ${\textstyle t}$ ${\textstyle r_{\text{s}}}$ ${\textstyle t}$

Cuando el momento angular no es cero, podemos reemplazar la dependencia del tiempo propio por una dependencia del ángulo usando la definición de ${\textstyle \varphi }$ ${\textstyle h}$

${\displaystyle \left({\frac {dr}{d\varphi }}\right)^{2}=\left({\frac {dr}{d\tau }}\right)^{2}\left({\frac {d\tau }{d\varphi }}\right)^{2}=\left({\frac {dr}{d\tau }}\right)^{2}\left({\frac {r^{2}}{h}}\right)^{2},}$

lo que produce la ecuación de la órbita

${\displaystyle \left({\frac {dr}{d\varphi }}\right)^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right)}$

donde, por brevedad, dos escalas de longitud, y , han sido definidas por ${\textstyle a}$ ${\textstyle b}$

${\displaystyle {\begin{aligned}a&={\frac {h}{c}},\\b&={\frac {cL}{E}}={\frac {hmc}{E}}.\end{aligned}}}$

Tenga en cuenta que en el caso taquiónico, será imaginario y real o infinito. ${\textstyle a}$ ${\textstyle b}$

La misma ecuación también se puede derivar utilizando un enfoque lagrangiano ^[4] o la ecuación de Hamilton-Jacobi ^[5] (ver más abajo). La solución de la ecuación de la órbita es

${\displaystyle \varphi =\int {\frac {dr}{\pm r^{2}{\sqrt {{\frac {1}{b^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {1}{a^{2}}}+{\frac {1}{r^{2}}}\right)}}}}.}$

Esto se puede expresar en términos de la función elíptica de Weierstrass . ^[6] ${\textstyle \wp }$

Velocidades locales y retardadas

A diferencia de la mecánica clásica, en Schwarzschild las coordenadas y no son las componentes radial y transversal de la velocidad local (en relación con un observador estacionario), sino que dan las componentes de la celeridad que están relacionadas con ${\textstyle {\frac {{\rm {d}}r}{{\rm {d}}\tau }}}$ ${\textstyle r\ {\frac {{\rm {d}}\varphi }{{\rm {d}}\tau }}}$ ${\textstyle v_{\parallel }}$ ${\textstyle v_{\perp }}$ ${\textstyle v}$ ${\textstyle v}$

${\displaystyle {\frac {{\rm {d}}r}{{\rm {d}}\tau }}=v_{\parallel }{\sqrt {1-{\frac {r_{\text{s}}}{r}}}}\ \gamma }$

para el radial y

${\displaystyle {\frac {{\rm {d}}\varphi }{{\rm {d}}\tau }}={\frac {v_{\perp }}{r}}\ \gamma }$

para la componente transversal del movimiento, con . El contable de coordenadas lejos de la escena observa la velocidad retardada de Shapiro , que está dada por la relación ${\textstyle v^{2}=v_{\parallel }^{2}+v_{\perp }^{2}}$ ${\textstyle {\hat {v}}}$

${\displaystyle {\hat {v}}_{\perp }=v_{\perp }{\sqrt {1-{\frac {r_{\text{s}}}{r}}}}}$y . ${\displaystyle {\hat {v}}_{\parallel }=v_{\parallel }\left(1-{\frac {r_{\text{s}}}{r}}\right)}$

El factor de dilatación del tiempo entre el contable y la partícula de prueba en movimiento también se puede expresar en la forma

${\displaystyle {\frac {{\rm {d}}\tau }{{\rm {d}}t}}={\frac {\sqrt {1-{\frac {r_{\text{s}}}{r}}}}{\gamma }}}$

donde el numerador es el gravitacional y el denominador es el componente cinemático de la dilatación del tiempo. Para una partícula que cae desde el infinito, el factor izquierdo es igual al factor derecho, ya que la velocidad de caída coincide con la velocidad de escape en este caso. ${\textstyle v}$ ${\textstyle c{\sqrt {\frac {r_{\text{s}}}{r}}}}$

Las dos constantes del momento angular y la energía total de una partícula de prueba con masa están en términos de ${\textstyle L}$ ${\textstyle E}$ ${\textstyle m}$ ${\textstyle v}$

${\displaystyle L=m\ v_{\perp }\ r\ \gamma }$

y

${\displaystyle E=mc^{2}\ {\sqrt {1-{\frac {r_{\text{s}}}{r}}}}\ \gamma }$

dónde

${\displaystyle E=E_{\rm {rest}}+E_{\rm {kin}}+E_{\rm {pot}}}$

y

${\displaystyle E_{\rm {rest}}=mc^{2}\ ,\ \ E_{\rm {kin}}=(\gamma -1)mc^{2}\ ,\ \ E_{\rm {pot}}=\left({\sqrt {1-{\frac {r_{\text{s}}}{r}}}}-1\right)\ \gamma \ mc^{2}}$

Para partículas de prueba masivas, el factor de Lorentz es el tiempo adecuado, mientras que para partículas sin masa, como los fotones, se establece y asume el papel de un parámetro afín. Si la partícula no tiene masa se reemplaza con y con , donde está la constante de Planck y la frecuencia observada localmente. ${\textstyle \gamma }$ ${\textstyle \gamma =1/{\sqrt {1-v^{2}/c^{2}}}}$ ${\textstyle \tau }$ ${\textstyle \gamma }$ ${\textstyle 1}$ ${\textstyle \tau }$ ${\textstyle E_{\rm {rest}}}$ ${\textstyle E_{\rm {kin}}}$ ${\textstyle mc^{2}}$ ${\textstyle hf}$ ${\textstyle h}$ ${\textstyle f}$

Solución exacta usando funciones elípticas.

La ecuación fundamental de la órbita es más fácil de resolver ^{[nota 1]} si se expresa en términos del radio inverso ${\textstyle u={\frac {1}{r}}}$

${\displaystyle \left({\frac {du}{d\varphi }}\right)^{2}={\frac {1}{b^{2}}}-\left(1-ur_{\text{s}}\right)\left({\frac {1}{a^{2}}}+u^{2}\right)}$

El lado derecho de esta ecuación es un polinomio cúbico , que tiene tres raíces , denotadas aquí como , y ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle u_{3}}$

${\displaystyle \left({\frac {du}{d\varphi }}\right)^{2}=r_{\text{s}}\left(u-u_{1}\right)\left(u-u_{2}\right)\left(u-u_{3}\right)}$

La suma de las tres raíces es igual al coeficiente del término. ${\textstyle u^{2}}$

${\displaystyle u_{1}+u_{2}+u_{3}={\frac {1}{r_{\text{s}}}}}$

Un polinomio cúbico con coeficientes reales puede tener tres raíces reales o una raíz real y dos raíces conjugadas complejas . Si las tres raíces son números reales , las raíces se etiquetan de esa manera . Si, en cambio, solo hay una raíz real, entonces se denota como ; las raíces conjugadas complejas están etiquetadas como y . Utilizando la regla de los signos de Descartes , puede haber como máximo una raíz negativa; es negativo si y sólo si . Como se analiza más adelante, las raíces son útiles para determinar los tipos de órbitas posibles. ${\textstyle u_{1}<u_{2}<u_{3}}$ ${\textstyle u_{3}}$ ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle u_{1}}$ ${\textstyle b<a}$

Dado este etiquetado de las raíces, la solución de la ecuación orbital fundamental es

${\displaystyle u=u_{1}+\left(u_{2}-u_{1}\right)\,\mathrm {sn} ^{2}\left({\frac {1}{2}}\varphi {\sqrt {r_{\text{s}}\left(u_{3}-u_{1}\right)}}+\delta \right)}$

donde representa la función amplitudinus sinusal (una de las funciones elípticas de Jacobi ) y es una constante de integración que refleja la posición inicial. El módulo elíptico de esta función elíptica viene dado por la fórmula ${\textstyle \mathrm {sn} }$ ${\textstyle \delta }$ ${\textstyle k}$

${\displaystyle k={\sqrt {\frac {u_{2}-u_{1}}{u_{3}-u_{1}}}}}$

límite newtoniano

Para recuperar la solución newtoniana para las órbitas planetarias, se toma el límite cuando el radio de Schwarzschild llega a cero. En este caso, la tercera raíz se vuelve aproximadamente y mucho más grande que o . Por tanto, el módulo tiende a cero; en ese límite, se convierte en la función seno trigonométrica ${\textstyle r_{\text{s}}}$ ${\textstyle u_{3}}$ ${\textstyle {\frac {1}{r_{\text{s}}}}}$ ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle k}$ ${\textstyle \mathrm {sn} }$

${\displaystyle u=u_{1}+\left(u_{2}-u_{1}\right)\,\sin ^{2}\left({\frac {1}{2}}\varphi +\delta \right)}$

De acuerdo con las soluciones de Newton para los movimientos planetarios, esta fórmula describe una cónica focal de excentricidad. ${\textstyle e}$

${\displaystyle e={\frac {u_{2}-u_{1}}{u_{2}+u_{1}}}}$

Si es un número real positivo, entonces la órbita es una elipse donde y representan las distancias de mayor y mayor aproximación, respectivamente. Si es cero o un número real negativo, la órbita es una parábola o una hipérbola , respectivamente. En estos dos últimos casos, representa la distancia de máxima aproximación; dado que la órbita va al infinito ( ), no hay distancia de máximo acercamiento. ${\textstyle u_{1}}$ ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle u=0}$

Raíces y descripción general de posibles órbitas.

Una raíz representa un punto de la órbita donde la derivada desaparece, es decir, donde . En dicho punto de inflexión, se alcanza un máximo, un mínimo o un punto de inflexión, dependiendo del valor de la segunda derivada, que viene dado por la fórmula ${\textstyle {\frac {du}{d\phi }}=0}$ ${\textstyle u}$

${\displaystyle {\frac {d^{2}u}{d\varphi ^{2}}}={\frac {r_{\text{s}}}{2}}\left[\left(u-u_{2}\right)\left(u-u_{3}\right)+\left(u-u_{1}\right)\left(u-u_{3}\right)+\left(u-u_{1}\right)\left(u-u_{2}\right)\right]}$

Si las tres raíces son números reales distintos, la segunda derivada es positiva, negativa y positiva en u ₁ , u ₂ y u ₃ , respectivamente. De ello se deduce que una gráfica de u versus φ puede oscilar entre u ₁ y u ₂ , o puede alejarse de u ₃ hacia el infinito (lo que corresponde a que r tiende a cero). Si u ₁ es negativo, sólo se producirá realmente una parte de una "oscilación". Esto corresponde a la partícula que viene del infinito, se acerca a la masa central y luego se aleja nuevamente hacia el infinito, como la trayectoria hiperbólica en la solución clásica.

Si la partícula tiene la cantidad justa de energía para su momento angular, u ₂ y u ₃ se fusionarán. Hay tres soluciones en este caso. La órbita puede entrar en espiral hasta , acercándose a ese radio como (asintóticamente) una exponencial decreciente en φ, o . O se puede tener una órbita circular en ese radio. O uno puede tener una órbita que desciende en espiral desde ese radio hasta el punto central. El radio en cuestión se llama radio interior y está entre y 3 veces r _s . También se produce una órbita circular cuando es igual a , y esto se llama radio exterior. Estos diferentes tipos de órbitas se analizan a continuación. ${\textstyle r={\frac {1}{u_{2}}}={\frac {1}{u_{3}}}}$ ${\textstyle \tau }$ ${\textstyle t}$ ${\textstyle {\frac {3}{2}}}$ ${\textstyle u_{2}}$ ${\textstyle u_{1}}$

Si la partícula llega a la masa central con suficiente energía y un momento angular suficientemente bajo, entonces sólo será real. Esto corresponde a la partícula que cae en un agujero negro. La órbita entra en espiral con un cambio finito en φ. ${\textstyle u_{1}}$

Precesión de órbitas

La función sn y su cuadrado sn ² tienen períodos de 4 K y 2 K , respectivamente, donde K está definido por la ecuación ^{[nota 2]}

${\displaystyle K=\int _{0}^{1}{\frac {dy}{\sqrt {\left(1-y^{2}\right)\left(1-k^{2}y^{2}\right)}}}}$

Por lo tanto, el cambio en φ durante una oscilación de (o, equivalentemente, una oscilación de ) es igual a ^[7] ${\textstyle u}$ ${\textstyle r}$

${\displaystyle \Delta \varphi ={\frac {4K}{\sqrt {r_{\text{s}}\left(u_{3}-u_{1}\right)}}}}$

En el límite clásico, u ₃ se aproxima y es mucho mayor que o . Por lo tanto, es aproximadamente ${\textstyle {\frac {1}{r_{\text{s}}}}}$ ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle k^{2}}$

${\displaystyle k^{2}={\frac {u_{2}-u_{1}}{u_{3}-u_{1}}}\approx r_{\text{s}}\left(u_{2}-u_{1}\right)\ll 1}$

Por las mismas razones, el denominador de Δφ es aproximadamente

${\displaystyle {\frac {1}{\sqrt {r_{\text{s}}\left(u_{3}-u_{1}\right)}}}={\frac {1}{\sqrt {1-r_{\text{s}}\left(2u_{1}+u_{2}\right)}}}\approx 1+{\frac {1}{2}}r_{\text{s}}\left(2u_{1}+u_{2}\right)}$

Como el módulo es cercano a cero, el período K se puede expandir en potencias de ; al orden más bajo, esta expansión produce ${\textstyle k}$ ${\textstyle k}$

${\displaystyle K\approx \int _{0}^{1}{\frac {dy}{\sqrt {1-y^{2}}}}\left(1+{\frac {1}{2}}k^{2}y^{2}\right)={\frac {\pi }{2}}\left(1+{\frac {k^{2}}{4}}\right)}$

Sustituyendo estas aproximaciones en la fórmula para Δφ se obtiene una fórmula para el avance angular por oscilación radial.

${\displaystyle \delta \varphi =\Delta \varphi -2\pi \approx {\frac {3}{2}}\pi r_{\text{s}}\left(u_{1}+u_{2}\right)}$

Para una órbita elíptica, representan las inversas de las distancias más larga y más corta, respectivamente. Estos se pueden expresar en términos del semieje mayor de la elipse y su excentricidad orbital . ${\textstyle u_{1}}$ ${\textstyle u_{2}}$ ${\textstyle A}$ ${\textstyle e}$

${\displaystyle {\begin{aligned}r_{\text{max}}&={\frac {1}{u_{1}}}=A(1+e)\\r_{\text{min}}&={\frac {1}{u_{2}}}=A(1-e)\end{aligned}}}$

donación

${\displaystyle u_{1}+u_{2}={\frac {2}{A\left(1-e^{2}\right)}}}$

Sustituyendo la definición de da la ecuación final. ${\textstyle r_{\text{s}}}$

${\displaystyle \delta \varphi \approx {\frac {6\pi GM}{c^{2}A\left(1-e^{2}\right)}}}$

Curvatura de la luz por gravedad.

Diagrama de lentes gravitacionales de un cuerpo compacto.

En el límite, cuando la masa de la partícula m llega a cero (o, de manera equivalente, si la luz se dirige directamente hacia la masa central, cuando la escala de longitud a llega al infinito), la ecuación para la órbita se convierte en

${\displaystyle \varphi =\int {\frac {dr}{r^{2}{\sqrt {{\frac {1}{b^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {1}{r^{2}}}}}}}}$

Expandiéndose en potencias de , el término de orden principal en esta fórmula da la deflexión angular aproximada δ φ para una partícula sin masa que entra desde el infinito y regresa al infinito: ${\textstyle {\frac {r_{\text{s}}}{r}}}$

${\displaystyle \delta \varphi \approx {\frac {2r_{\text{s}}}{b}}={\frac {4GM}{c^{2}b}}.}$

Aquí está el parámetro de impacto , algo mayor que la distancia de máxima aproximación : ^[8] ${\textstyle b}$ ${\textstyle r_{3}}$

${\displaystyle b=r_{3}{\sqrt {\frac {r_{3}}{r_{3}-r_{\text{s}}}}}}$

Aunque esta fórmula es aproximada, es precisa para la mayoría de las mediciones de lentes gravitacionales , debido a la pequeñez de la relación . Para la luz que roza la superficie del sol, la desviación angular aproximada es de aproximadamente 1,75  segundos de arco , aproximadamente una millonésima parte de un círculo. ${\textstyle {\frac {r_{\text{s}}}{r}}}$

De manera más general, las geodésicas de un fotón emitido por una fuente de luz ubicada en una coordenada radial se pueden calcular de la siguiente manera, aplicando la ecuación ${\displaystyle r={1 \over u}\in [r_{s},\infty [}$

${\displaystyle {\left({du \over d\varphi }\right)^{2}}=r_{s}\ u^{3}-u^{2}+{\frac {1}{b^{2}}}}$

Geodésica de un fotón emitido desde una fuente de luz situada en el horizonte de sucesos de un agujero negro y de regreso a él, con parámetro de impacto . ${\displaystyle b>b_{crit}={\frac {3{\sqrt {3}}}{2}}r_{s}}$

Geodésica de un fotón emitido desde una fuente de luz ubicada en el horizonte de sucesos de un agujero negro, con parámetro de impacto y luego desplazándose a la órbita inestable . Si el fotón se libera hacia el infinito. ${\displaystyle b=b_{crit}={\frac {3{\sqrt {3}}}{2}}r_{s}}$ ${\displaystyle r_{crit}={3 \over 2}r_{s}}$ ${\displaystyle b<b_{crit}}$

La ecuación se puede derivar como

${\displaystyle 2\ {\frac {du}{d\varphi }}{\frac {d^{2}u}{d\varphi ^{2}}}=3\ r_{s}u^{2}{\frac {du}{d\varphi }}-2\ u{\frac {du}{d\varphi }}}$

lo que lleva a

${\displaystyle {d^{2}u \over d\varphi ^{2}}={\frac {3}{2}}\ r_{s}\ u^{2}-\ u}$

Esta ecuación con segunda derivada se puede integrar numéricamente de la siguiente manera mediante un método de Runge-Kutta de 4º orden, considerando un tamaño de paso ^y con: ${\displaystyle \Delta \varphi }$

${\displaystyle k_{1}={d^{2}u \over d\varphi ^{2}}(u)}$,

${\displaystyle k_{2}={d^{2}u \over d\varphi ^{2}}{\bigl (}u+{\frac {\Delta \varphi }{2}}{du \over d\varphi }{\bigr )}}$,

${\displaystyle k_{3}={d^{2}u \over d\varphi ^{2}}{\Bigl (}u+{\frac {\Delta \varphi }{2}}{du \over d\varphi }+{\frac {\Delta \varphi ^{2}}{4}}k_{1}{\Bigr )}}$y

${\displaystyle k_{4}={d^{2}u \over d\varphi ^{2}}{\Bigl (}u+\Delta \varphi {du \over d\varphi }+{\frac {\Delta \varphi ^{2}}{2}}k_{2}{\Bigr )}}$.

El valor en el siguiente paso es ${\displaystyle {du \over d\varphi }(\varphi +\Delta \varphi )}$

${\displaystyle {du \over d\varphi }(\varphi )+{\frac {\Delta \varphi }{6}}(k_{1}+2k_{2}+2k_{3}+k_{4})}$

y el valor en el siguiente paso es ${\displaystyle u(\varphi +\Delta \varphi )}$

${\displaystyle u(\varphi )+\Delta \varphi {du \over d\varphi }(\varphi )+{\frac {\Delta \varphi ^{2}}{6}}(k_{1}+k_{2}+k_{3})}$

Se puede elegir que el paso sea constante o adaptativo, dependiendo de la precisión requerida . ${\displaystyle \Delta \varphi }$ ${\displaystyle r={1 \over u}}$

Relación con la física newtoniana

Energía potencial radial efectiva

La ecuación de movimiento de la partícula derivada arriba.

${\displaystyle \left({\frac {dr}{d\tau }}\right)^{2}={\frac {E^{2}}{m^{2}c^{2}}}-c^{2}+{\frac {r_{\text{s}}c^{2}}{r}}-{\frac {L^{2}}{m\mu r^{2}}}+{\frac {r_{\text{s}}L^{2}}{m\mu r^{3}}}}$

se puede reescribir usando la definición del radio de Schwarzschild r _s como

${\displaystyle {\frac {1}{2}}m\left({\frac {dr}{d\tau }}\right)^{2}=\left[{\frac {E^{2}}{2mc^{2}}}-{\frac {1}{2}}mc^{2}\right]+{\frac {GMm}{r}}-{\frac {L^{2}}{2\mu r^{2}}}+{\frac {G(M+m)L^{2}}{c^{2}\mu r^{3}}},}$

que es equivalente a una partícula que se mueve en un potencial efectivo unidimensional

${\displaystyle V(r)=-{\frac {GMm}{r}}+{\frac {L^{2}}{2\mu r^{2}}}-{\frac {G(M+m)L^{2}}{c^{2}\mu r^{3}}}}$

Los dos primeros términos son energías clásicas muy conocidas, siendo el primero la energía potencial gravitacional newtoniana de atracción y el segundo correspondiente a la energía potencial repulsiva "centrífuga" ; sin embargo, el tercer término es una energía de atracción exclusiva de la relatividad general . Como se muestra a continuación y en otros lugares , esta energía cúbica inversa hace que las órbitas elípticas precedan gradualmente en un ángulo δφ por revolución.

${\displaystyle \delta \varphi \approx {\frac {6\pi G(M+m)}{c^{2}A\left(1-e^{2}\right)}}}$

donde es el semieje mayor y es la excentricidad. ${\textstyle A}$ ${\textstyle e}$

El tercer término es atractivo y domina en valores pequeños, dando un radio interno crítico r _interior en el cual una partícula es atraída inexorablemente hacia adentro ; este radio interior es función del momento angular de la partícula por unidad de masa o, de manera equivalente, la escala de longitud definida anteriormente. ${\textstyle r}$ ${\textstyle r=0}$ ${\textstyle a}$

Órbitas circulares y su estabilidad.

Potencial radial efectivo para varios momentos angulares. En radios pequeños, la energía cae precipitadamente, lo que hace que la partícula sea atraída inexorablemente hacia adentro . Sin embargo, cuando el momento angular normalizado es igual a la raíz cuadrada de tres, es posible una órbita circular metaestable en el radio resaltado con un círculo verde. En momentos angulares más altos, hay una barrera centrífuga significativa (curva naranja) y un radio interior inestable, resaltado en rojo. ${\textstyle r=0}$ ${\textstyle {\frac {a}{r_{\text{s}}}}={\frac {L}{mcr_{\text{s}}}}}$

El potencial efectivo se puede reescribir en términos de longitud . ${\textstyle V}$ ${\textstyle a={\frac {h}{c}}}$

${\displaystyle V(r)={\frac {\mu c^{2}}{2}}\left[-{\frac {r_{\text{s}}}{r}}+{\frac {a^{2}}{r^{2}}}-{\frac {r_{\text{s}}a^{2}}{r^{3}}}\right]}$

Las órbitas circulares son posibles cuando la fuerza efectiva es cero.

${\displaystyle F=-{\frac {dV}{dr}}=-{\frac {\mu c^{2}}{2r^{4}}}\left[r_{\text{s}}r^{2}-2a^{2}r+3r_{\text{s}}a^{2}\right]=0}$

es decir, cuando las dos fuerzas de atracción (la gravedad newtoniana (primer término) y la atracción exclusiva de la relatividad general (tercer término)) están exactamente equilibradas por la fuerza centrífuga repulsiva (segundo término). Hay dos radios en los que puede ocurrir este equilibrio, denotados aquí como r _interior y r _exterior.

${\displaystyle {\begin{aligned}r_{\text{outer}}&={\frac {a^{2}}{r_{\text{s}}}}\left(1+{\sqrt {1-{\frac {3r_{\text{s}}^{2}}{a^{2}}}}}\right)\\[3pt]r_{\text{inner}}&={\frac {a^{2}}{r_{\text{s}}}}\left(1-{\sqrt {1-{\frac {3r_{\text{s}}^{2}}{a^{2}}}}}\right)={\frac {3a^{2}}{r_{\text{outer}}}}\end{aligned}}}$

los cuales se obtienen usando la fórmula cuadrática . El radio interior r _interior es inestable, porque la tercera fuerza de atracción se fortalece mucho más rápido que las otras dos fuerzas cuando r se vuelve pequeño; si la partícula se desliza ligeramente hacia adentro desde r _interior (donde las tres fuerzas están en equilibrio), la tercera fuerza domina a las otras dos y atrae a la partícula inexorablemente hacia adentro hasta r = 0. Sin embargo, en el radio exterior, las órbitas circulares son estables; el tercer término es menos importante y el sistema se comporta más como el problema no relativista de Kepler .

Cuando es mucho mayor que (el caso clásico), estas fórmulas se vuelven aproximadamente ${\textstyle a}$ ${\textstyle r_{\text{s}}}$

${\displaystyle {\begin{aligned}r_{\text{outer}}&\approx {\frac {2a^{2}}{r_{\text{s}}}}\\[3pt]r_{\text{inner}}&\approx {\frac {3}{2}}r_{\text{s}}\end{aligned}}}$

Los radios estable e inestable se representan frente al momento angular normalizado en azul y rojo, respectivamente. Estas curvas se encuentran en una órbita circular única (círculo verde) cuando el momento angular normalizado es igual a la raíz cuadrada de tres. A modo de comparación, el radio clásico predicho a partir de la aceleración centrípeta y la ley de gravedad de Newton está representado en negro. ${\textstyle {\frac {a}{r_{\text{s}}}}={\frac {L}{mcr_{\text{s}}}}}$

Sustituyendo las definiciones de y r _s en r _exterior se obtiene la fórmula clásica para una partícula de masa que orbita un cuerpo de masa . ${\textstyle a}$ ${\textstyle m}$ ${\textstyle M}$

${\displaystyle r_{\mathrm {outer} }^{3}={\frac {G(M+m)}{\omega _{\varphi }^{2}}}}$

donde ω _φ es la velocidad angular orbital de la partícula. Esta fórmula se obtiene en mecánica no relativista igualando la fuerza centrífuga a la fuerza gravitacional newtoniana:

${\displaystyle {\frac {GMm}{r^{2}}}=\mu \omega _{\varphi }^{2}r}$

¿Dónde está la masa reducida ? ${\textstyle \mu }$

En nuestra notación, la velocidad angular orbital clásica es igual

${\displaystyle \omega _{\varphi }^{2}\approx {\frac {GM}{r_{\mathrm {outer} }^{3}}}=\left({\frac {r_{\text{s}}c^{2}}{2r_{\mathrm {outer} }^{3}}}\right)=\left({\frac {r_{\text{s}}c^{2}}{2}}\right)\left({\frac {r_{\text{s}}^{3}}{8a^{6}}}\right)={\frac {c^{2}r_{\text{s}}^{4}}{16a^{6}}}}$

En el otro extremo, cuando a ² se acerca a 3 r _s ² desde arriba, los dos radios convergen a un solo valor

${\displaystyle r_{\mathrm {outer} }\approx r_{\mathrm {inner} }\approx 3r_{\text{s}}}$

Las soluciones cuadráticas anteriores aseguran que r _exterior sea siempre mayor que 3 r _s , mientras que r _interior se encuentre entre 3 ⁄ 2 r _s y 3 r _s . Las órbitas circulares menores de 3 ⁄ 2 r _s no son posibles. Para partículas sin masa, a llega al infinito, lo que implica que hay una órbita circular para los fotones en r _interior = 3 ⁄ 2 r _s . La esfera de este radio a veces se conoce como esfera de fotones .   

Precesión de órbitas elípticas.

En el problema no relativista de Kepler , una partícula sigue eternamente la misma elipse perfecta (órbita roja). La relatividad general introduce una tercera fuerza que atrae la partícula ligeramente más fuerte que la gravedad newtoniana, especialmente en radios pequeños. Esta tercera fuerza hace que la órbita elíptica de la partícula precese (órbita cian) en la dirección de su rotación; este efecto se ha medido en Mercurio , Venus y la Tierra. El punto amarillo dentro de las órbitas representa el centro de atracción, como es el Sol .

La tasa de precesión orbital se puede derivar utilizando este potencial radial efectivo V. Una pequeña desviación radial de una órbita circular de radio r _exterior oscilará establemente con una frecuencia angular.

${\displaystyle \omega _{r}^{2}={\frac {1}{m}}\left[{\frac {d^{2}V}{dr^{2}}}\right]_{r=r_{\mathrm {outer} }}}$

que es igual

${\displaystyle \omega _{r}^{2}=\left({\frac {c^{2}r_{\text{s}}}{2r_{\mathrm {outer} }^{4}}}\right)\left(r_{\mathrm {outer} }-r_{\mathrm {inner} }\right)=\omega _{\varphi }^{2}{\sqrt {1-{\frac {3r_{\text{s}}^{2}}{a^{2}}}}}}$

Tomando la raíz cuadrada de ambos lados y realizando una expansión en serie de Taylor se obtiene

${\displaystyle \omega _{r}=\omega _{\varphi }\left[1-{\frac {3r_{\text{s}}^{2}}{4a^{2}}}+{\mathcal {O}}\left({\frac {r_{\text{s}}^{4}}{a^{4}}}\right)\right]}$

Multiplicando por el período T de una revolución se obtiene la precesión de la órbita por revolución.

${\displaystyle \delta \varphi =T\left(\omega _{\varphi }-\omega _{r}\right)\approx 2\pi \left({\frac {3r_{\text{s}}^{2}}{4a^{2}}}\right)={\frac {3\pi m^{2}c^{2}}{2L^{2}}}r_{\text{s}}^{2}}$

donde hemos utilizado ω _φ T = 2 п y la definición de la escala de longitud a . Sustituyendo la definición del radio de Schwarzschild r _s se obtiene

${\displaystyle \delta \varphi \approx {\frac {3\pi m^{2}c^{2}}{2L^{2}}}\left({\frac {4G^{2}M^{2}}{c^{4}}}\right)={\frac {6\pi G^{2}M^{2}m^{2}}{c^{2}L^{2}}}}$

Esto se puede simplificar utilizando el semieje A de la órbita elíptica y la excentricidad e relacionada con la fórmula

${\displaystyle {\frac {h^{2}}{G(M+m)}}=A\left(1-e^{2}\right)}$

para dar el ángulo de precesión

${\displaystyle \delta \varphi \approx {\frac {6\pi G(M+m)}{c^{2}A\left(1-e^{2}\right)}}}$

Derivaciones matemáticas de la ecuación orbital.

Símbolos de Christoffel

Los símbolos de Christoffel que no desaparecen para la métrica de Schwarzschild son: ^[9]

${\displaystyle {\begin{aligned}\Gamma _{rt}^{t}=-\Gamma _{rr}^{r}&={\frac {r_{\text{s}}}{2r(r-r_{\text{s}})}}\\[3pt]\Gamma _{tt}^{r}&={\frac {r_{\text{s}}(r-r_{\text{s}})}{2r^{3}}}\\[3pt]\Gamma _{\phi \phi }^{r}&=(r_{\text{s}}-r)\sin ^{2}(\theta )\\[3pt]\Gamma _{\theta \theta }^{r}&=r_{\text{s}}-r\\[3pt]\Gamma _{r\theta }^{\theta }=\Gamma _{r\phi }^{\phi }&={\frac {1}{r}}\\[3pt]\Gamma _{\phi \phi }^{\theta }&=-\sin(\theta )\cos(\theta )\\[3pt]\Gamma _{\theta \phi }^{\phi }&=\cot(\theta )\end{aligned}}}$

ecuación geodésica

Según la teoría de la relatividad general de Einstein, partículas de masa insignificante viajan a lo largo de geodésicas en el espacio-tiempo. En el espacio-tiempo plano, lejos de una fuente de gravedad, estas geodésicas corresponden a líneas rectas; sin embargo, pueden desviarse de las líneas rectas cuando el espacio-tiempo es curvo. La ecuación de las líneas geodésicas es ^[10]

${\displaystyle {\frac {d^{2}x^{\lambda }}{dq^{2}}}+\Gamma _{\mu \nu }^{\lambda }{\frac {dx^{\mu }}{dq}}{\frac {dx^{\nu }}{dq}}=0}$

donde Γ representa el símbolo de Christoffel y la variable parametriza el camino de la partícula a través del espacio-tiempo , su llamada línea mundial . El símbolo de Christoffel depende únicamente del tensor métrico , o más bien de cómo cambia con la posición. La variable es un múltiplo constante del tiempo adecuado para órbitas temporales (que son recorridas por partículas masivas) y generalmente se considera igual a él. Para órbitas luminosas (o nulas) (que son recorridas por partículas sin masa como el fotón ), el tiempo adecuado es cero y, estrictamente hablando, no puede usarse como variable . Sin embargo, las órbitas similares a la luz se pueden derivar como el límite ultrarelativista de las órbitas temporales, es decir, el límite cuando la masa de la partícula m llega a cero mientras se mantiene fija su energía total . ${\textstyle q}$ ${\textstyle g_{\mu \nu }}$ ${\textstyle q}$ ${\textstyle \tau }$ ${\textstyle q}$

Por lo tanto, para resolver el movimiento de una partícula, la forma más sencilla es resolver la ecuación geodésica, un enfoque adoptado por Einstein ^[11] y otros. ^[12] La métrica de Schwarzschild puede escribirse como

${\displaystyle c^{2}d\tau ^{2}=w(r)c^{2}dt^{2}-v(r)dr^{2}-r^{2}d\theta ^{2}-r^{2}\sin ^{2}\theta d\phi ^{2}\,}$

donde las dos funciones y su recíproco se definen por brevedad. A partir de esta métrica, se pueden calcular los símbolos de Christoffel y sustituir los resultados en las ecuaciones geodésicas. ${\textstyle w(r)=1-{\frac {r_{\text{s}}}{r}}}$ ${\textstyle v(r)={\frac {1}{w(r)}}}$ ${\textstyle \Gamma _{\mu \nu }^{\lambda }}$

${\displaystyle {\begin{aligned}0&={\frac {d^{2}\theta }{dq^{2}}}+{\frac {2}{r}}{\frac {d\theta }{dq}}{\frac {dr}{dq}}-\sin \theta \cos \theta \left({\frac {d\phi }{dq}}\right)^{2}\\[3pt]0&={\frac {d^{2}\phi }{dq^{2}}}+{\frac {2}{r}}{\frac {d\phi }{dq}}{\frac {dr}{dq}}+2\cot \theta {\frac {d\phi }{dq}}{\frac {d\theta }{dq}}\\[3pt]0&={\frac {d^{2}t}{dq^{2}}}+{\frac {1}{w}}{\frac {dw}{dr}}{\frac {dt}{dq}}{\frac {dr}{dq}}\\[3pt]0&={\frac {d^{2}r}{dq^{2}}}-{\frac {1}{v}}{\frac {dv}{dr}}\left({\frac {dr}{dq}}\right)^{2}-{\frac {r}{v}}\left({\frac {d\theta }{dq}}\right)^{2}-{\frac {r\sin ^{2}\theta }{v}}\left({\frac {d\phi }{dq}}\right)^{2}+{\frac {c^{2}}{2v}}{\frac {dw}{dr}}\left({\frac {dt}{dq}}\right)^{2}\end{aligned}}}$

Se puede verificar que es una solución válida sustituyendo la primera de estas cuatro ecuaciones. Por simetría, la órbita debe ser plana y somos libres de disponer el sistema de coordenadas de manera que el plano ecuatorial sea el plano de la órbita. Esta solución simplifica la segunda y cuarta ecuaciones. ${\textstyle \theta ={\frac {\pi }{2}}}$ ${\textstyle \theta }$

Para resolver la segunda y tercera ecuaciones basta con dividirlas entre y , respectivamente. ${\textstyle {\frac {d\phi }{dq}}}$ ${\textstyle {\frac {dt}{dq}}}$

${\displaystyle {\begin{aligned}0&={\frac {d}{dq}}\left[\ln {\frac {d\phi }{dq}}+\ln r^{2}\right]\\[3pt]0&={\frac {d}{dq}}\left[\ln {\frac {dt}{dq}}+\ln w\right],\end{aligned}}}$

lo que produce dos constantes de movimiento.

Enfoque lagrangiano

Debido a que las partículas de prueba siguen las geodésicas en una métrica fija, las órbitas de esas partículas pueden determinarse utilizando el cálculo de variaciones, también llamado enfoque lagrangiano. ^[13] Las geodésicas en el espacio-tiempo se definen como curvas para las cuales pequeñas variaciones locales en sus coordenadas (mientras se mantienen fijos los eventos de sus puntos finales) no producen cambios significativos en su longitud total s . Esto se puede expresar matemáticamente usando el cálculo de variaciones.

${\displaystyle 0=\delta s=\delta \int ds=\delta \int {\sqrt {g_{\mu \nu }{\frac {dx^{\mu }}{d\tau }}{\frac {dx^{\nu }}{d\tau }}}}d\tau =\delta \int {\sqrt {2T}}d\tau }$

donde τ es el tiempo propio , s = cτ es la longitud del arco en el espacio-tiempo y T se define como

${\displaystyle 2T=c^{2}=\left({\frac {ds}{d\tau }}\right)^{2}=g_{\mu \nu }{\frac {dx^{\mu }}{d\tau }}{\frac {dx^{\nu }}{d\tau }}=\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}\left({\frac {dt}{d\tau }}\right)^{2}-{\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\left({\frac {dr}{d\tau }}\right)^{2}-r^{2}\left({\frac {d\varphi }{d\tau }}\right)^{2}}$

en analogía con la energía cinética . Si la derivada con respecto al tiempo propio se representa con un punto por brevedad

${\displaystyle {\dot {x}}^{\mu }={\frac {dx^{\mu }}{d\tau }}}$

T puede escribirse como

${\displaystyle 2T=c^{2}=\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}\left({\dot {t}}\right)^{2}-{\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\left({\dot {r}}\right)^{2}-r^{2}\left({\dot {\varphi }}\right)^{2}}$

Los factores constantes (como c o la raíz cuadrada de dos) no afectan la respuesta al problema variacional; por lo tanto, tomando la variación dentro de la integral se obtiene el principio de Hamilton

${\displaystyle 0=\delta \int {\sqrt {2T}}d\tau =\int {\frac {\delta T}{\sqrt {2T}}}d\tau ={\frac {1}{c}}\delta \int Td\tau .}$

La solución del problema variacional viene dada por las ecuaciones de Lagrange.

${\displaystyle {\frac {d}{d\tau }}\left({\frac {\partial T}{\partial {\dot {x}}^{\sigma }}}\right)={\frac {\partial T}{\partial x^{\sigma }}}.}$

Cuando se aplican a t y φ , estas ecuaciones revelan dos constantes de movimiento

${\displaystyle {\begin{aligned}{\frac {d}{d\tau }}\left[r^{2}{\frac {d\varphi }{d\tau }}\right]&=0,\\{\frac {d}{d\tau }}\left[\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}\right]&=0,\end{aligned}}}$

que puede expresarse en términos de dos escalas de longitud constante, y ${\textstyle a}$ ${\textstyle b}$

${\displaystyle {\begin{aligned}r^{2}{\frac {d\varphi }{d\tau }}&=ac,\\\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}&={\frac {a}{b}}.\end{aligned}}}$

Como se muestra arriba, la sustitución de estas ecuaciones en la definición de la métrica de Schwarzschild produce la ecuación de la órbita.

Enfoque hamiltoniano

Una solución lagrangiana se puede transformar en una forma hamiltoniana equivalente. ^[14] En este caso, el hamiltoniano viene dado por ${\displaystyle H}$

${\displaystyle 2H=c^{2}={\frac {p_{t}^{2}}{c^{2}\left(1-{\frac {r_{\text{s}}}{r}}\right)}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)p_{r}^{2}-{\frac {p_{\theta }^{2}}{r^{2}}}-{\frac {p_{\varphi }^{2}}{r^{2}\sin ^{2}\theta }}}$

Una vez más, la órbita puede estar restringida por simetría. Dado que y no aparecen en el hamiltoniano, sus momentos conjugados son constantes; pueden expresarse en términos de la velocidad de la luz y dos escalas de longitud constante y ${\textstyle \theta ={\frac {\pi }{2}}}$ ${\textstyle t}$ ${\textstyle \varphi }$ ${\textstyle c}$ ${\textstyle a}$ ${\textstyle b}$

${\displaystyle {\begin{aligned}p_{\varphi }&=-ac\\p_{\theta }&=0\\p_{t}&={\frac {ac^{2}}{b}}\end{aligned}}}$

Las derivadas con respecto al tiempo propio están dadas por

${\displaystyle {\begin{aligned}{\frac {dr}{d\tau }}&={\frac {\partial H}{\partial p_{r}}}=-\left(1-{\frac {r_{\text{s}}}{r}}\right)p_{r}\\{\frac {d\varphi }{d\tau }}&={\frac {\partial H}{\partial p_{\varphi }}}={\frac {-p_{\varphi }}{r^{2}}}={\frac {ac}{r^{2}}}\\{\frac {dt}{d\tau }}&={\frac {\partial H}{\partial p_{t}}}={\frac {p_{t}}{c^{2}\left(1-{\frac {r_{\text{s}}}{r}}\right)}}={\frac {a}{b\left(1-{\frac {r_{\text{s}}}{r}}\right)}}\end{aligned}}}$

Dividiendo la primera ecuación por la segunda se obtiene la ecuación orbital.

${\displaystyle {\frac {dr}{d\varphi }}=-{\frac {r^{2}}{ac}}\left(1-{\frac {r_{\text{s}}}{r}}\right)p_{r}}$

El momento radial p _r se puede expresar en términos de r usando la constancia del hamiltoniano ; esto produce la ecuación orbital fundamental ${\textstyle H={\frac {c^{2}}{2}}}$

${\displaystyle \left({\frac {dr}{d\varphi }}\right)^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right)}$

Enfoque de Hamilton-Jacobi

Curvatura de ondas en un campo gravitacional. Debido a la gravedad, el tiempo pasa más lentamente en la parte inferior que en la superior, lo que hace que los frentes de onda (que se muestran en negro) se doblen gradualmente hacia abajo. La flecha verde muestra la dirección de la aparente "atracción gravitacional".

La ecuación orbital se puede derivar de la ecuación de Hamilton-Jacobi . ^[15] La ventaja de este enfoque es que equipara el movimiento de la partícula con la propagación de una onda, y conduce claramente a la derivación de la desviación de la luz por la gravedad en la relatividad general , a través del principio de Fermat . La idea básica es que, debido a la desaceleración gravitacional del tiempo, las partes de un frente de onda más cercanas a una masa gravitante se mueven más lentamente que las más alejadas, desviando así la dirección de propagación del frente de onda.

Utilizando la covarianza general, la ecuación de Hamilton-Jacobi para una sola partícula de masa unitaria se puede expresar en coordenadas arbitrarias como

${\displaystyle g^{\mu \nu }{\frac {\partial S}{\partial x^{\mu }}}{\frac {\partial S}{\partial x^{\nu }}}=c^{2}.}$

Esto es equivalente a la formulación hamiltoniana anterior, donde las derivadas parciales de la acción reemplazan a los momentos generalizados. Usando la métrica de Schwarzschild g ^μν , esta ecuación se convierte en

${\displaystyle {\frac {1}{c^{2}\left(1-{\frac {r_{\text{s}}}{r}}\right)}}\left({\frac {\partial S}{\partial t}}\right)^{2}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {\partial S}{\partial r}}\right)^{2}-{\frac {1}{r^{2}}}\left({\frac {\partial S}{\partial \varphi }}\right)^{2}=c^{2}}$

donde nuevamente orientamos el sistema de coordenadas esféricas con el plano de la órbita. El tiempo t y el ángulo azimutal φ son coordenadas cíclicas, por lo que la solución de la función principal S de Hamilton se puede escribir

${\displaystyle S=-p_{t}t+p_{\varphi }\varphi +S_{r}(r)\,}$

donde y son los momentos generalizados constantes. La ecuación de Hamilton-Jacobi da una solución integral para la parte radial ${\displaystyle p_{t}}$ ${\displaystyle p_{\varphi }}$ ${\displaystyle S_{r}(r)}$

${\displaystyle S_{r}(r)=\int ^{r}{\frac {dr}{1-{\frac {r_{\text{s}}}{r}}}}{\sqrt {{\frac {p_{t}^{2}}{c^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left(c^{2}+{\frac {p_{\varphi }^{2}}{r^{2}}}\right)}}.}$

Tomando la derivada de la función principal S de Hamilton con respecto al momento conservado p _φ se obtiene

${\displaystyle {\frac {\partial S}{\partial p_{\varphi }}}=\varphi +{\frac {\partial S_{r}}{\partial p_{\varphi }}}=\mathrm {constant} }$

que es igual

${\displaystyle \varphi -\int ^{r}{\frac {p_{\varphi }dr}{r^{2}{\sqrt {{\frac {p_{t}^{2}}{c^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left(c^{2}+{\frac {p_{\varphi }^{2}}{r^{2}}}\right)}}}}=\mathrm {constant} }$

Tomando una variación infinitesimal en φ y r se obtiene la ecuación orbital fundamental

${\displaystyle \left({\frac {dr}{d\varphi }}\right)^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right).}$

donde las escalas de longitud conservadas a y b están definidas por los momentos conservados por las ecuaciones

${\displaystyle {\begin{aligned}{\frac {\partial S}{\partial \varphi }}=p_{\varphi }&=-ac\\{\frac {\partial S}{\partial t}}=p_{t}&={\frac {ac^{2}}{b}}\end{aligned}}}$

principio de hamilton

La integral de acción para una partícula afectada únicamente por la gravedad es

${\displaystyle S=\int {-mc^{2}d\tau }=-mc\int {c{\frac {d\tau }{dq}}dq}=-mc\int {{\sqrt {-g_{\mu \nu }{\frac {dx^{\mu }}{dq}}{\frac {dx^{\nu }}{dq}}}}dq}}$

¿Dónde está el tiempo adecuado y cualquier parametrización suave de la línea mundial de la partícula? Si a esto se le aplica el cálculo de variaciones , se obtienen nuevamente las ecuaciones de una geodésica. Para simplificar los cálculos, primero se toma la variación del cuadrado del integrando. Para la métrica y las coordenadas de este caso y suponiendo que la partícula se mueve en el plano ecuatorial , ese cuadrado es ${\textstyle \tau }$ ${\textstyle q}$ ${\textstyle \theta ={\frac {\pi }{2}}}$

${\displaystyle \left(c{\frac {d\tau }{dq}}\right)^{2}=-g_{\mu \nu }{\frac {dx^{\mu }}{dq}}{\frac {dx^{\nu }}{dq}}=\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}\left({\frac {dt}{dq}}\right)^{2}-{\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\left({\frac {dr}{dq}}\right)^{2}-r^{2}\left({\frac {d\varphi }{dq}}\right)^{2}\,.}$

Tomar una variación de esto da

${\displaystyle \delta \left(c{\frac {d\tau }{dq}}\right)^{2}=2c^{2}{\frac {d\tau }{dq}}\delta {\frac {d\tau }{dq}}=\delta \left[\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}\left({\frac {dt}{dq}}\right)^{2}-{\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\left({\frac {dr}{dq}}\right)^{2}-r^{2}\left({\frac {d\varphi }{dq}}\right)^{2}\right]\,.}$

Movimiento en longitud

Varíe con respecto a la longitud sólo para obtener ${\textstyle \varphi }$

${\displaystyle 2c^{2}{\frac {d\tau }{dq}}\delta {\frac {d\tau }{dq}}=-2r^{2}{\frac {d\varphi }{dq}}\delta {\frac {d\varphi }{dq}}\,.}$

Dividir por para obtener la variación del integrando mismo. ${\textstyle 2c{\frac {d\tau }{dq}}}$

${\displaystyle c\,\delta {\frac {d\tau }{dq}}=-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\delta {\frac {d\varphi }{dq}}=-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}{\frac {d\delta \varphi }{dq}}\,.}$

De este modo

${\displaystyle 0=\delta \int {c{\frac {d\tau }{dq}}dq}=\int {c\delta {\frac {d\tau }{dq}}dq}=\int {-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}{\frac {d\delta \varphi }{dq}}dq}\,.}$

Integrando por partes da

${\displaystyle 0=-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\delta \varphi -\int {{\frac {d}{dq}}\left[-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\right]\delta \varphi dq}\,.}$

Se supone que la variación de la longitud es cero en los puntos finales, por lo que el primer término desaparece. La integral puede hacerse distinta de cero mediante una elección perversa de a menos que el otro factor interno sea cero en todas partes. Entonces la ecuación del movimiento es ${\textstyle \delta \varphi }$

${\displaystyle {\frac {d}{dq}}\left[-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\right]=0\,.}$

movimiento en el tiempo

Variar con respecto al tiempo sólo para obtener ${\textstyle t}$

${\displaystyle 2c^{2}{\frac {d\tau }{dq}}\delta {\frac {d\tau }{dq}}=2\left(1-{\frac {r_{\text{s}}}{r}}\right)c^{2}{\frac {dt}{dq}}\delta {\frac {dt}{dq}}\,.}$

Dividir por para obtener la variación del integrando mismo. ${\textstyle 2c{\frac {d\tau }{dq}}}$

${\displaystyle c\delta {\frac {d\tau }{dq}}=c\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}\delta {\frac {dt}{dq}}=c\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}{\frac {d\delta t}{dq}}\,.}$

De este modo

${\displaystyle 0=\delta \int {c{\frac {d\tau }{dq}}dq}=\int {c\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}{\frac {d\delta t}{dq}}dq}\,.}$

Integrando por partes da

${\displaystyle 0=c\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}\delta t-\int {{\frac {d}{dq}}\left[c\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}\right]\delta tdq}\,.}$

Entonces la ecuación del movimiento es

${\displaystyle {\frac {d}{dq}}\left[c\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}\right]=0\,.}$

Momento conservado

Integre estas ecuaciones de movimiento para determinar las constantes de integración obteniendo

${\displaystyle {\begin{aligned}L=p_{\phi }&=mr^{2}{\frac {d\varphi }{d\tau }}\,,\\E=-p_{t}&=mc^{2}\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {dt}{d\tau }}\,.\end{aligned}}}$

Estas dos ecuaciones para las constantes de movimiento (momento angular) y (energía) se pueden combinar para formar una ecuación que es cierta incluso para fotones y otras partículas sin masa para las cuales el tiempo adecuado a lo largo de una geodésica es cero. ${\textstyle L}$ ${\textstyle E}$

${\displaystyle {\frac {d\varphi }{dt}}=\left(1-{\frac {r_{\text{s}}}{r}}\right){\frac {L\,c^{2}}{E\,r^{2}}}\,.}$

movimiento radial

Sustituyendo

${\displaystyle {\frac {d\varphi }{d\tau }}={\frac {L}{m\,r^{2}}}\,}$

y

${\displaystyle {\frac {dt}{d\tau }}={\frac {E}{\left(1-{\frac {r_{\text{s}}}{r}}\right)m\,c^{2}}}\,}$

en la ecuación métrica (y usando ) da ${\textstyle \theta ={\frac {\pi }{2}}}$

${\displaystyle c^{2}={\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\,{\frac {E^{2}}{m^{2}c^{2}}}-{\frac {1}{1-{\frac {r_{\text{s}}}{r}}}}\left({\frac {dr}{d\tau }}\right)^{2}-{\frac {1}{r^{2}}}\,{\frac {L^{2}}{m^{2}}}\,,}$

de donde se puede derivar

${\displaystyle {\left({\frac {dr}{d\tau }}\right)}^{2}={\frac {E^{2}}{m^{2}c^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left(c^{2}+{\frac {L^{2}}{m^{2}r^{2}}}\right)\,,}$

cual es la ecuación de movimiento para . La dependencia de on se puede encontrar dividiendo esto por ${\textstyle r}$ ${\textstyle r}$ ${\textstyle \varphi }$

${\displaystyle {\left({\frac {d\varphi }{d\tau }}\right)}^{2}={\frac {L^{2}}{m^{2}r^{4}}}}$

Llegar

${\displaystyle {\left({\frac {dr}{d\varphi }}\right)}^{2}={\frac {E^{2}r^{4}}{L^{2}c^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {m^{2}c^{2}r^{4}}{L^{2}}}+r^{2}\right)\,}$

lo cual es cierto incluso para partículas sin masa. Si las escalas de longitud están definidas por

${\displaystyle a={\frac {L}{m\,c}}}$

y

${\displaystyle b={\frac {L\,c}{E}}\,,}$

entonces la dependencia de se simplifica a ${\textstyle r}$ ${\textstyle \varphi }$

${\displaystyle {\left({\frac {dr}{d\varphi }}\right)}^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{\text{s}}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right)\,.}$

Ver también

• Portal de astronomía
• Portal de física

• Problema clásico de fuerza central
• Campos de marco en relatividad general
• problema de kepler
• Problema de dos cuerpos en la relatividad general

Notas

1. Esta sustitución de for también es común en los problemas clásicos de fuerza central, ya que también hace que esas ecuaciones sean más fáciles de resolver. Para obtener más información, consulte el artículo sobre el problema clásico de la fuerza central . ${\textstyle u}$ ${\textstyle r}$
2. En la literatura matemática, K se conoce como la integral elíptica completa de primer tipo ; Para obtener más información, consulte el artículo sobre integrales elípticas .

Referencias

1. Kozai, Yoshihide (1998). "Desarrollo de la mecánica celeste en Japón". Planeta. Ciencia espacial . 46 (8): 1031–36. Código Bib : 1998P&SS...46.1031K. doi :10.1016/s0032-0633(98)00033-6.
2. Kaplan, Samuil (1949). "Sobre las órbitas circulares en la teoría de la gravitación de Einstein". J. Exp. Teor. Física . 19 (10): 951–952. arXiv : 2201.07971 . Código Bib : 1949ZhETF..19..951K.
3. Landau y Lifshitz, págs. 299–301.
4. Whittaker 1937.
5. Landau y Lifshitz (1975), págs. 306–309.
6. Gibones, GW; Vyska, M. (29 de febrero de 2012). "La aplicación de funciones elípticas de Weierstrass a geodésicas nulas de Schwarzschild". Gravedad clásica y cuántica . 29 (6): 065016. arXiv : 1110.6508 . Código Bib : 2012CQGra..29f5016G. doi :10.1088/0264-9381/29/6/065016. S2CID  119675906.
7. Synge, págs. 294-295.
8. arXiv.org: gr-qc/9907034v1.
9. Sean Carroll : Apuntes de la conferencia sobre la relatividad general , Capítulo 7, Ec. 7.33
10. Weinberg, pág. 122.
11. Einstein, págs. 95-96.
12. Weinberg, págs. 185-188; Wald, págs. 138-139.
13. Synge, págs. 290–292; Adler, Bazin y Schiffer, págs. 179-182; Whittaker, págs. 390–393; Pauli, pág. 167.
14. Lanczos, págs. 331–338.
15. Landau y Lifshitz, págs. 306–307; Misner, Thorne y Wheeler, págs. 636–679.

Bibliografía

• Schwarzschild, K. (1916). Über das Gravitationsfeld eines Massenpunktes nach der Einstein'schen Theorie. Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften 1 , 189-196.
  • escaneo del papel original
  • texto del artículo original, en Wikisource
  • traducción de Antoci y Loinger
  • un comentario sobre el artículo, dando una derivación más simple
• Schwarzschild, K. (1916). Über das Gravitationsfeld einer Kugel aus inkcompressibler Flüssigkeit. Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften 1 , 424-?.
• Flamm, L (1916). "Beiträge zur Einstein'schen Gravitationstheorie". Physikalische Zeitschrift . 17 : 448–?.
• Adler, R; Bazin M; Schiffer M (1965). Introducción a la Relatividad General . Nueva York: McGraw-Hill Book Company. págs. 177-193. ISBN 978-0-07-000420-7.
• Einstein, A (1956). El significado de la relatividad (5ª ed.). Princeton, Nueva Jersey: Princeton University Press. págs. 92–97. ISBN 978-0-691-02352-6.
• Hagihara, Y (1931). "Teoría de las trayectorias relativistas en un campo gravitacional de Schwarzschild". Revista Japonesa de Astronomía y Geofísica . 8 : 67-176. Código bibliográfico : 1931AOTok..31...67H. ISSN  0368-346X.
• Lanczos, C (1986). Los principios variacionales de la mecánica (4ª ed.). Nueva York: Publicaciones de Dover. págs. 330–338. ISBN 978-0-486-65067-8.
• Landau, LD ; Lifshitz, EM (1975). La teoría clásica de los campos . Curso de Física Teórica . vol. 2 (cuarta edición en inglés revisada). Nueva York: Pergamon Press. págs. 299–309. ISBN 978-0-08-018176-9.
• Misner, CW ; Thorne, K y Wheeler, JA (1973). Gravitación . San Francisco: WH Freeman. Capítulo 25 (págs. 636–687), §33.5 (págs. 897–901) y §40.5 (págs. 1110–1116). ISBN 978-0-7167-0344-0.(Ver Gravitación (libro) ).
• País, A. (1982). Sutil es el Señor: la ciencia y la vida de Albert Einstein . Prensa de la Universidad de Oxford. págs. 253–256. ISBN 0-19-520438-7.
• Pauli, W (1958). Teoría de la relatividad . Traducido por G. Campo. Nueva York: Publicaciones de Dover. págs. 40–41, 166–169. ISBN 978-0-486-64152-2.
• Rindler, W (1977). Relatividad esencial: especial, general y cosmológica (2ª ed. revisada). Nueva York: Springer Verlag. págs. 143-149. ISBN 978-0-387-10090-6.
• Roseveare, NT (1982). El perihelio de Mercurio, de Leverrier a Einstein . Prensa de la Universidad de Oxford. ISBN 0-19-858174-2.

• Synge, JL (1960). Relatividad: la teoría general . Ámsterdam: Editorial de Holanda Septentrional. págs. 289–298. ISBN 978-0-7204-0066-3.
• Wald, RM (1984). Relatividad general . Chicago: Prensa de la Universidad de Chicago. págs. 136-146. ISBN 978-0-226-87032-8.
• Walter, S. (2007). "Rompiendo los 4 vectores: el movimiento tetradimensional en la gravitación, 1905-1910". En Renn, J. (ed.). La Génesis de la Relatividad General . vol. 3. Berlín: Springer. págs. 193-252. Archivado desde el original el 30 de enero de 2009 . Consultado el 13 de julio de 2010 .

• Weinberg, S (1972). Gravitación y cosmología . Nueva York: John Wiley and Sons. págs. 185-201. ISBN 978-0-471-92567-5.
• Whittaker, et al. (1937). Tratado sobre la dinámica analítica de partículas y cuerpos rígidos, con una introducción al problema de los tres cuerpos (4ª ed.). Nueva York: Publicaciones de Dover. págs. 389–393. ISBN 978-1-114-28944-4.

enlaces externos

• Extracto de Reflexiones sobre la relatividad de Kevin Brown.