슈바르츠실트 계량

슈바르츠실트 계량(Schwarzschild metric) 혹은 슈바르츠실트 해(Schwarzschild solution)는 아인슈타인 방정식의 엄밀해(exact solution) 중 하나로, 구형 대칭이고 전하와 각운동량이 없다고 가정한 천체의 바깥 시공간을 나타낸다. 이 수학적 모델은 일반 상대성 이론에서 상대적으로 느리게 회전하는, 그리고 대전되지 않은 천체의 중력 작용을 파악할 때 유용한 근사를 제공한다.

슈바르츠실트는 1915년 말 이 해를 처음 발견한 독일 천문학자 카를 슈바르츠실트(Karl Schwarzschild)의 이름을 딴 것이다. 다음은 슈바르츠실트 계량을 지칭하는 세 가지 용어이다.

이 세 가지는 사실상 같은 의미이지만, 이와 같은 뉘앙스의 차이가 있다.

슈바르츠실트 해는 일반 상대성 이론을 기반으로 중심 천체가 만드는 정적 중력장에서의 일체 문제를 해결하기 위한 과정에서 도입되었다. 알베르트 아인슈타인이 1915년 11월 18일에 오랜 난제였던 수성의 근일점 세차 운동 문제를 일반 상대성 이론을 통해 해결한 것[Einstein(1915)]이 일반 상대론 기반 일체 문제의 시초라고 할 수 있지만, 근거로 삼은 중력장 해가 직교 좌표계 기반이었으며 근사적 풀이였다.

일체 문제를 다루는 중력장에 대한 최초의 엄밀한 해는 1915년 12월 독일 물리학자 카를 슈바르츠실트(1873~1916)가 제시하였다.[Schwarzschild(1916)] 러시아 전선에서 복무중이었던 그는 자신의 풀이와 논문을 아인슈타인에게 전달하고 아인슈타인이 1916년 1월 대신 프로이센 학회에 제출해주었다. 이것이 지금의 슈바르츠실트 해이며, 해의 명칭도 그의 이름을 딴 것이다. 다만, 슈바르츠실트는 아인슈타인이 초기에 제시했던 좌표 제한에 따른 방정식을 풀어서 지금의 해와 약간 다르며, 특히 슈바르츠실트 반지름을 좌표 설정으로 $r=0$에 두었다는 차이가 있다.

이후 완성된 (좌표 조건이 제거된) 방정식을 기반으로 수학자 요하네스 드로스테(Johannes Droste, 1886~1963)[Droste(1917)][4], 다비드 힐베르트(David Hilbert, 1862~1943)[Hilbert(1916)][6], 헤르만 바일(Hermann Weyl, 1885~1955)[Weyl(1919)] 등이 차례로 여러 좌표 설정을 통해 동등한 해를 제시하였다. 아인슈타인은 헤르만 바일의 풀이를 선호하여 1921년 프린스턴 강의에서 인용하였다.[8][9] 1922년 발견된 버코프 정리(Birkhoff Theorem)에 의해, 이들이 발견한 해는 모두 동치이며 좌표 변환으로 서로 유도된다는 것이 증명되었다.

슈바르츠실트 계량은 다음 조건으로 구성되어 있다.

(1), (2)는 슈바르츠실트 좌표계



를 도입할 수 있도록 해준다. (이 좌표계에 대한 자세한 내용은 슈바르츠실트 좌표계 문서 참고.) 다음으로 천체가 좌표계에 대하여 정지한 완전 유체(perfect fluid)로 이루어져 있다고 가정하여 스트레스 에너지 텐서를 구한 다음, 좌표계에 대한 아인슈타인 방정식을 풀면 된다. (좁은 의미에서 슈바르츠실트 계량이라 불리는) 외부해, 즉 진공해는 스트레스 에너지 텐서가 0일 때 얻는다.

외부해와 달리 내부해는 별의 상태에 관한 몇 가지 조건(밀도, 압력)에 있어서 보다 신중한 접근을 통해 물리적으로 가능성 있는 모델을 세워나가야 한다. 가령 슈바르츠실트는 유체의 밀도가 항상 일정하다는 조건을 걸어서 내부해를 완성하였다.

이 문서에서는 특별한 언급이 없는 한 기하단위계 $G = c = 1$을 사용한다. 먼저, 별의 내부 구조를 결정하는 몇 가지 매개변수가 필요하다. 별은 완전 유체(Perfect fluid)로 구성되어 있다고 가정하자. 이로부터 다음과 같이 관련 매개변수를 도입할 수 있다. 이들은 모두 유체의 국소적 정지계(rest frame)에서 정의되었으며, $r$에만 의존한다고 가정할 수 있다. (유체에 관한 자세한 설명은 상대론적 유체 참고.)

첫번째로 스트레스 에너지 텐서를 계산해보자. 이를 위해서는 유체의 4-속도를 먼저 결정해야 한다. 별은 (좌표계에 대해) 정지한 상태이므로, $U^r = U^{\theta} = U^{\phi} = 0$으로 설정할 수 있다. 이 때



이므로,



를 얻는다. 따라서, $T_{\mu\nu}$의 각 성분은 다음과 같이 계산된다. 비대각 성분은 모두 $0$이다.



이제, 논의를 단순하게 만들기 위해 유체의 압력 $p$와 밀도 $\rho$가 개수 밀도 $n$에만 의존한다고 가정할 수 있다. 일반적으로는 본래 온도 $T$나 엔트로피 $s$를 추가로 알아야 하나, 상대론이 중요해지는 천체들은 모두 엔트로피를 $r$에 의존하지 않는 상수로 근사시킬 수 있다. (자세한 논의는 MTW gravitation (1973) p.599를 참고할 수 있다.) 이렇게 정의된 다음 방정식을 별의 상태 방정식(equations of state)이라 부른다.

별의 상태에 의해 결정되고, 외부의 시공간 구조를 결정하는 $\Phi(r)$과 $\Lambda(r)$를 일반적으로 구조 방정식(equations of structure)이라 부른다. 아인슈타인 방정식에 본격적으로 대입하기 전에 에너지-운동량 보존법칙



을 살펴볼 수 있다. (일반적으로, 10개의 아인슈타인 방정식을 바로 푸는 것보다 에너지 보존법칙으로 4개를 먼저 해결하는 것이 훨씬 수월하다.) 대칭성에 의해, 이 항등식은 다음 $\mu = r$ 성분만이 남는다.



이제, 아인슈타인 방정식을 풀 차례이다. 아인슈타인 텐서의 $00$ 성분을 살펴보면 (각종 수치는 슈바르츠실트 좌표계 문서 참고.)



에서, 편의상



라 설정한다. 그러면



이 되며, 아인슈타인 방정식 $G_{00} = 8\pi T_{00}$ 및 $G_{rr} = 8\pi T_{rr}$은 다음과 같다.

으로 정리된다. 두번째 식은 $(\rho + p)d\Phi/dr = -dp/dr$로부터



으로 바꿀 수 있는데, 이를 톨만-오펜하이머-볼코프(Tolman–Oppenheimer–Volkoff) 방정식 혹은 TOV 방정식이라 한다. 이 방정식에 대한 자세한 내용은 문서를 참고한다.

두번째 식은 다음과 같이 $r$을 고유 거리로 계산해야 정확히 비교할 수 있다.

이로부터 알 수 있는 가장 중요한 사실은 기본적으로 일반 상대론에서 예측하는 압력의 증가율이 더 크다는 것이다. 이 압력은 기본적으로 별의 자체 중력을 견디기 위한(혹은 좌표계에 대하여 가만히 있기 위해 바깥으로 가속하는[14]) 것이므로, 일반 상대성 이론이 예측하는 중력은 고전 역학이 예측하는 중력보다 더 크다고 압축 설명할 수 있다. 뿐만 아니라, 우변에 $p$가 들어가 있기 때문에 압력을 더 크게 증가시킬수록 더 깊은 곳에서는 훨씬 빠르게 압력을 증가시켜야 한다(self-regenerating 또는 self-defeated).

별의 바깥으로 나가면, $\rho = p = 0$이 된다. 따라서,



이 된다. 즉,



을 얻는다. 여기에서 $M$은 상수이며 두번째 식은 $r \rightarrow \infty$일 때 $\Phi \rightarrow 0$임을 이용한 것이다. (이들은 사실 진공 방정식 $G_{00} = 0, G_{rr} = 0$을 푼 것이다.) 이로부터, 다음 슈바르츠실트 계량(Schwarzschild metric)



을 얻는다. $r$이 충분히 클 때, 슈바르츠실트 계량의 $M$은 뉴턴 역학에서 별의 총 질량과 대응되며 더 나아가 이 시공간은 민코프스키 공간에 수렴한다. 여기서 흔히 $r_{s} = 2M$이라 두어,



라고 표기한다. 이 $r_{s}$를 슈바르츠실트 반경(Schwarzschild radius)이라고 한다. 여기서는 $G=c=1$이라 두었으므로 단위를 살리면 $\displaystyle r_s = \frac{2GM}{c^2}$이다.

을 각각 $r$에 대하여 적분함으로써 $m(r)$, $p(r)$, $\Phi(r)$을 구할 수 있다. 이들 적분값을 하나로 정하기 위해서는 공통적으로 $r = 0$에서의 초기조건이 필요하다. 먼저, 일반적인 경우 $m(0) = 0$이라 둘 수 있다. 한편, $p_c := p(0)$은 별의 중심 압력(central pressure)으로써 TOV 방정식의 초기 조건이 되며, 압력은 바깥으로 갈수록 단조 감소하다 별의 경계에서 $0$이 된다. 내부해와 외부해의 연속성을 보장하기 위해, $p(R) = 0$이 되는 $r=R$을 별의 반지름으로 설정한다. 이 때, 별의 총 질량은



이 된다. 마지막으로



을 적분 상수로 두면 세 가지 방정식에 대한 해가 모두 완성된다. 하지만 내부해를 결정하기 위해서는 별의 상태 방정식 $\rho(r), \,\, p(r)$이 먼저 주어져야 하는데, 이에 대한 모델이 다음에 제시되어 있다.

처음 진공해를 제시한 슈바르츠실트는 얼마 지나지 않아 내부해 또한 제시하였는데, 그는 별의 에너지 밀도 $\rho$가 언제나 일정하다고 가정하였다. 이는 지나치게 단순화된 모델이긴 하나, 중성자별의 경우 실제로 내부 밀도는 거의 일정하며, 다른 모델을 살펴볼 때에도 참고가 될 수 있다.

이 경우, $m(r)$은 즉각 다음과 같이 결정된다. ($R$은 별의 반지름이다.)



$r \geq R$에서 $\rho = 0$이므로



이 된다. 이것은 이 모델에서 별의 총 질량을 나타내며, 슈바르츠실트 질량(Schwarzschild Mass)이라고도 한다.

한편, TOV 방정식은



이므로, 중심 압력 $p_c$로부터 다음과 같이 적분할 수 있다.



이로부터, 다음 관계식들을 얻을 수 있다.



만약 $M/R \rightarrow 4/9$이면 $p_c$는 무한대로 발산함을 알 수 있다. 마지막으로, $g_{00}(R) = -(1-2M/R)$이라 두면



를 얻는다.

측지선은 아인슈타인 방정식의 각각의 특수해를 깊이 이해하기 위한 가장 핵심적인 도구이다. 측지선에 속하는 것은 빛(광자), 입자, 작은 행성 등이 있으며, 이들은 다양한 조건의 천체 주위에서 특징적인 궤도를 그린다. 슈바르츠실트 시공간은 매우 느리게 회전하는 구형 대칭 천체 주위의 공간을 나타내며, 적색 편이, 근일점 세차운동, 빛의 굴절 등 고전 역학이 예측하지 못하는 다양한 현상을 가장 단순하고도 선명하게 보여준다. 이러한 비고전적 효과는 천문학계의 다양한 노력으로 검증이 이루어졌고, 반대로 이러한 현상을 활용하여 천체(펄사 등)의 보다 정확한 질량을 예측하는 등의 활용도 동시에 이루어지고 있다.

일반 상대성 이론에서는 아무런 힘을 받지 않는 측지선도 다음 항등식



즉 $g_{\mu\nu}$의 변동에 의해 4-운동량 $p^{\mu}$의 성분이 보존되지 않는다. 이는 고전 역학의 중력 작용에 대응된다. 비슷하게 어떤 벡터장을 따라 시공간의 점들을 옮겼을 때 시공간이 정확히 똑같은 형태이면, 즉 $g_{\mu\nu}$가 벡터장을 따라 보존되면 입자는 이 벡터장을 따라서는 힘(중력)을 받지 않으며 4-운동량의 해당 성분이 보존된다. 이는 그 자체로 측지선을 따라 보존되는 양, 즉 운동 상수(constants of motion)를 제공한다. 이러한 $g_{\mu\nu}$의 대칭성이 많을수록 운동 상수는 많아지고 4-운동량의 표현은 간단해진다. 따라서, 효율적으로 측지선을 분석하는 데 있어서 운동 상수를 찾는 것은 최우선순위가 된다. 일반적인 경우에 이렇게 얻는 운동 상수는 사실 고전 물리학의 에너지, 운동량, 각운동량과 대응된다.

일반 상대성 이론에서 운동 상수를 찾는 가장 일반화된 과정은 킬링 벡터장(Killing vector field)이란 것이 제공한다.[15] 간단하게만 언급하자면 킬링 벡터장이란, $g_{\mu\nu}$의 리미분(Lie Derivatives)이 0인 벡터장을 말한다.[16] 다시 말해 킬링 벡터장 $\xi^{\mu}$는 $\mathcal{L}_{\xi} g_{\mu\nu} = 0$을 만족시키며, 이를 풀면 다음 킬링 방정식(killing equation)을 얻는다.



따라서, 측지선 $x^{\mu}(\lambda)$에 접하는 벡터(4-속도 혹은 4-운동량) $\displaystyle p^{\mu} = \frac{dx^{\mu}}{d\lambda}$을 정의하면 이는 다음을 만족시킨다.



첫 항은 $\sigma, \mu$에 관한 대칭성과 킬링 방정식 때문에, 그리고 두번째 항은 측지선 방정식 때문에 사라진다. 따라서 $\xi_{\mu}\,p^{\mu}$, 즉 $p^{\mu}$의 $\xi$ 성분은 측지선을 따라 보존되는 양이다. 이것은 킬링 벡터장의 방향에 따라 고전 역학에 잘 대응되는 운동 상수를 정의할 수 있다. 예를 들어, 시간 방향 킬링 벡터장은 에너지를, 회전 방향 킬링 벡터장은 각 운동량을 정의한다.

만약 킬링 벡터장이 특정 좌표계의 한 기저 벡터장 $\vec{e}_{\alpha}$와 일치할 경우,



이므로 다음과 같이 표현할 수도 있다.



즉, "공변" 4-속도 혹은 4-운동량의 해당 성분이 측지선을 따라 보존된다고 말할 수 있다.[17] ($p_{\mu}$의 보존은 보다 초보적인 방법으로 증명할 수 있다. 다음을 살펴보자.)

[펼치기 · 접기]

$\displaystyle p^{\mu} = \frac{dx^{\mu}}{d\lambda}$가 정의하는 측지선 방정식

$\displaystyle p^{\alpha}p^{\mu}_{\,\,\,;\alpha} = 0$
$\displaystyle p^{\alpha}p_{\mu;\alpha} = 0$

에서

$\displaystyle \begin{aligned}p^{\alpha}p_{\mu,\alpha} &= p^{\alpha}\Gamma^{\beta}_{\mu\alpha}p_{\beta} = \frac{1}{2}p^{\alpha}g^{\beta\sigma}(g_{\sigma\mu, \alpha} + g_{\sigma\alpha, \mu} - g_{\mu\alpha, \sigma})p_{\beta} \\&= \frac{1}{2}(g_{\sigma\mu, \alpha} + g_{\sigma\alpha, \mu} - g_{\mu\alpha, \sigma})p^{\alpha}p^{\sigma} \end{aligned}$

이다. 그런데 $p^{\alpha}p^{\sigma}$는 $\alpha, \sigma$에 대하여 대칭이고, $g_{\sigma\mu, \alpha} - g_{\mu\alpha, \sigma}$는 반대칭이므로 그 곱은 사라진다. 즉,

$\displaystyle p^{\alpha}p_{\mu,\alpha} = \frac{dp_{\mu}}{d\lambda} = \frac{1}{2}g_{\sigma\alpha, \mu}p^{\alpha}p^{\sigma}$

를 얻는다. 즉, 메트릭 텐서의 성분 $g_{\mu\nu}$가 $x^{\alpha}$ 방향에 대하여 변화하지 않으면, $p_{\alpha}$는 측지선 $x^{\mu}(\lambda)$ 위에서 보존된다.

이제 본론으로 돌아오면, 슈바르츠실트 해에는 4개의 킬링 벡터장이 있다. 하나는 시간 대칭을, 세 개는 회전 대칭을 나타낸다. 고전역학에서 시간 대칭성을 나타내는 킬링 벡터장은 에너지(스칼라)의 보존, 회전 대칭성을 띠는 킬링 벡터장은 각운동량(3-벡터)의 보존과 대응된다. 이 때, 각운동량의 보존은 크기의 보존과 방향의 보존으로 나눌 수 있는데, 방향의 보존은 입자 또는 광자가 (원점을 지나는) 정해진 평면에서만 운동한다는 것으로 해석할 수 있다.

시간 대칭을 표현하는 킬링 벡터장은 $\vec{\xi}_t = \vec{e}_{t} = (1, 0, 0, 0)$, 회전 대칭을 표현하는 킬링 벡터장은 측지선이 놓이는 평면을 $\displaystyle \theta = \frac{\pi}{2}$로 고정시켰을 때 $\vec{\xi}_\phi = \vec{e}_{\phi} = (0, 0, 0, 1)$로 정할 수 있다. 이처럼 킬링 벡터장이 어떤 좌표 기저 벡터장과 일치하므로, 운동 상수는 공변 4-운동량의 성분으로 표현할 수 있다. 이제, 각각에 대응하는 운동 상수를 다음과 같이 정의한다.

슈바르츠실트 좌표계 $(r, \theta, \phi)$에 대하여 정지한, 즉 순간적인 4-속도 $\vec{U}$가 $\vec{e}_0$과 나란한 관찰자를 생각하자. 관찰자는 4-운동량이 $\vec{p}$ 인 신호(빛)을 방출하거나, 수신한다고 가정했을 때 빛의 에너지는 다음과 같이 계산할 수 있다.



이 때, $\vec{U}\cdot\vec{U} = g_{00}(U^0)^2 = -1$임을 이용하면



가 된다. 한편, $p_0 = g_{00}(\vec{e}_0)^0p^0 = -E$라 두면, 이는 빛의 측지선 상에서 보존된다.

를 얻는다. $r \rightarrow \infty$일 때 $E_{\text{observer}} = E$이므로, $E$는 무한히 멀리 떨어진 곳에 있는 관찰자가 받아들인 빛의 에너지임을 알 수 있다.

$r = r_1$에 놓인 관찰자 1과 $r = r_2$에 놓인 관찰자 2가 동일한 빛을 공유했을 때, 이들이 받아들이는 진동수의 비율은



가 된다. 반대로 적색편이의 값은



이 된다. $r_s$는 매우 작으므로, 이는



과 같다.

일반 상대성 이론에 따른 주변 천체나, 빛의 궤도는 측지선을 통해 파악할 수 있다. 먼저, 슈바르츠실트 해는 정적이고 회전 대칭성을 띠므로 $p_0$과 $p_{\phi}$가 보존된다. 각각에 대하여 입자의 질량 당 에너지 $\tilde E$와 각운동량 $\tilde L$을 다음과 같이 정의한다. 광자는 질량이 없으므로 그냥 에너지($E$)와 각운동량($L$)에 대응된다.



한편, 슈바르츠실트 해는 구형 대칭이므로 입자와 광자는 곡면 $\theta = \pi/2$ 위에서 운동한다고 가정할 수 있다. 이 때, $g_{\phi\phi} = r^2$이다. 이제, 입자의 운동량 성분을 다음과 같이 구할 수 있다. 광자는 에너지와 각운동량의 물결표를 제거하고 $\tau$가 정의되지 않으므로 어떤 아핀 매개변수(affine parameter)인 $\lambda$로 바꾸면 동일한 식을 얻는다.



여기서 결정되지 않은 것은 시간에 대한 $r$의 변화, 즉 $\displaystyle \frac{\mathrm{d}r}{\mathrm{d}\tau}$이다. 운동량의 크기가 항상 $p^{\mu}p_{\mu} = -m^2$(광자는 $m = 0$)임을 이용하면 다음을 얻는다.

이것이 슈바르츠실트의 측지선(운동 방정식)을 분석하는 가장 기본적인 방정식이다. 여기에서, 유효 퍼텐셜(effective potential) $\tilde V^2(r)$ (제곱이 아니다.)을 다음과 같이 정의한다.

유효 펴텐셜은 공통적으로 $\displaystyle \left(\frac{\mathrm{d}r}{\mathrm{d}\tau}\right)^2 = \tilde E^2 - \tilde V^2$을 만족시킨다. 좌변은 입자(광자)의 운동 에너지, 우변은 총 에너지와 퍼텐셜의 차이에 대응된다. 유효 퍼텐셜은 총 에너지보다 클 수 없으며, $r$ 방향의 속도(운동 에너지)는 유효 퍼텐셜이 작아질수록 커지고 유효 퍼텐셜이 커질수록 작아진다. 실제로 이 식의 양변을 미분하여 비교해보면



으로, $r$ 방향 가속도(힘)은 유효 퍼텐셜의 $r$ 방향 미분의 음수에 대응되며, 이러한 논의는 고전 역학과 잘 대응된다. (현재 정의의 절반을 유효 퍼텐셜로 정의한다면 더 잘 맞는다.)

한편, 유효 퍼텐셜의 각 항을 다음과 같이 분석할 수 있다.



첫번째는 상수(정지 에너지), 두번째는 고전적 중력 퍼텐셜, 세번째는 구심력에 대응된다. 마지막은 일반 상대성 이론에서 추가되는 보정항으로, 뉴턴의 중력에 추가되는 힘(상대론적 효과)이라고 볼 수 있다. 혹은,



로 정리하여 일반 상대성 이론의 중력 퍼텐셜로 보고, 뉴턴 중력은 첫번째 항만을 설명하는 근사 법칙으로 볼 수도 있다.

슈바르츠실트 측지선을 통해 입자 및 광자의 궤도를 분석하는 방법은 고전 역학에서 중심력을 다루는 방법과 크게 다르지 않다. 다만 3차 항이 추가되면서 그래프의 모양 및 극점의 위치가 조금 변화했을 뿐이다. 마찬가지로, 입자의 운동은 주어진 $\tilde L$로 $\tilde V^2(r)$의 그래프가 결정된 후 $\tilde E^2$의 크기에 따라 정해지며, 예를 들어, 그림의 상황에서 입자의 궤도는 다음 세 가지 유형을 살펴볼 수 있다. 출발한 $r$의 값과도 연관이 있으나, 같은 방식으로 분석할 수 있다. 주의할 점은, 이것은 $\phi$ 방향을 생략한 그림이라는 것이다. 실제로 $\tilde L$을 갖는 입자(광자)는 회전하면서 이러한 그래프를 그리게 된다.

일반적으로 입자는 가속도가 $0$인, 즉 $\displaystyle \frac{\mathrm{d}\tilde V^2}{\mathrm{d}r} = 0$인 곳에서 $r$이 고정된 궤도(즉, 원궤도)를 갖게 되는데, 그림에서는 $\tilde V^2 = \tilde E^2_A$와 $\tilde V^2 = \tilde E^2_B$가 그러하다. 그러나 $\tilde V^2 = \tilde E^2_A$의 경우 양쪽으로 조금만 움직여도 가속도 방향이 바깥쪽을 향하므로 불안정한 궤도이다. 반면 $\tilde V^2 = \tilde E^2_B$의 경우 양쪽에서 가속도 방향이 안쪽을 향하므로, 안정적인 궤도가 된다. (따라서, 원궤도로부터의 약간의 오차는 $r$의 진동을 만들어 타원 궤도와 유사한 형태가 나타난다.) 이것을 만족시키는 반지름의 해를 구하려면 다음 방정식을 풀면 된다.



다음은 각 $\tilde L$에 대한 $\tilde V^2(r)$의 그래프를 나타낸 것이며, 아래에서부터 차례대로 $(\tilde L/r_s)^2 = 1, 2, 3, 4, 5$에 해당하는 그래프이다. 점선은 각 그래프 위의 극점을 표시한 것으로 위 방정식의 해를 나타내며, 안정적이거나 불안정한 원궤도의 반지름과 같다. 표와 같이 세 가지 경우가 있다.

그림에서 살펴볼 수 있듯이, $\tilde L$이 클수록 원심력 효과로 가능한 궤도의 선택폭이 넓어진다. $\tilde L < \sqrt{3}r_s$부터는 입자가 궤도를 유지할 수 없고 언제나 중심 별로 빨려 들어간다.



마지막으로, 입자의 각운동량 $\tilde L$과 (불)안정적 궤도 $r$의 관계를 나타낸 그래프를 통해 고전 역학과 일반 상대성 이론을 비교해보자. 고전 역학에서는 어느 각운동량에서나 안정적 궤도가 하나 존재하며, 궤도의 반지름은 임의적이다. 그러나 일반 상대성 이론에서는 어느 이상의 각운동량이 확보되지 않으면 궤도는 존재하지 않으며, 그 이상부터는 안정적 궤도(큰 값), 불안정 궤도(작은 값) 두 개가 존재하게 된다. 그리고 그 경계에서는 불안정한 궤도 $r = 3r_s$가 하나 존재한다. 이로부터, 일반 상대성 이론에서 입자가 가질 수 있는 최소의 안정적 궤도 반지름은 $r = 3r_s$임을 알 수 있다.

을 풀면, 다음과 같이 단 하나의 경우만 남는다.



이는 불안정한 궤도이나, 빛이 하나의 궤도를 이루는 해를 가진다는 점은 주목할 만하다. 블랙홀의 빛 구 개념은 여기에서 온 것이다.

뉴턴 역학에서 행성의 궤도는 한바퀴, 즉 $\Delta \phi = 2\pi$마다 주기가 완료되며, 닫힌 타원이다. 따라서 $r$도 $2\pi$라는 주기를 갖는다. 일반 상대성 이론에서는 보정항 때문에 $r$의 주기가 약간 길어지며, 근일점(perihelion)[19]도 약간 회전한다. 이 때 $r$이 최소가 되는 점을 근일점이라 하므로, 근일점이 얼마나 회전하는지 알려면 $r(\phi)$의 주기를 구하면 된다. 운동 방정식으로부터 $\displaystyle \frac{\mathrm{d}r}{\mathrm{d}\phi} = \frac{\mathrm{d}r}{\mathrm{d}\tau}/\frac{\mathrm{d}\phi}{\mathrm{d}\tau}$를 계산하고, 이를 정리하면 $1/r$이 $k\phi$에 대한 삼각함수로 정리된다. 따라서 $2 \pi/k - 2\pi$가 근일점의 주기 당 이동량이 된다.

먼저



에서



이고, 여기서 $u=1/r$이라 두면




을 얻는다. 먼저, 뉴턴 중력은 $u^3$ 항을 무시했을 때 얻게 된다. $u$가 상수인 원형 궤도를 가정했을 때 $u = r_s/2\tilde L^2$이므로, 이로부터의 오차를 표현하기 위해 $\displaystyle s = u - \frac{r_s}{2\tilde L^2}$이라 둔다. 그러면,



을 얻는다. 양변을 적분하고, 초기값을 $\phi = -\pi/2$라 두면 다음을 얻게 된다.



이므로



이다. 이는 이심률이 $\displaystyle e = \frac{2(\tilde E^2 - 1)}{r_s} + \frac{r_s}{2\tilde L^2}$인 타원이다. 이 때 근일점을 $r_{\mathrm P} = a(1-e)$, 원일점을 $r_{\mathrm A} = a(1+e)$라 하면[20], 이들은



의 두 해이다. 이로부터 $\displaystyle \frac{1}{r_{\mathrm P}} + \frac{1}{r_{\mathrm A}} = \frac{2}{a(1-e^2)} = \frac{r_s}{\tilde L^2}$을 얻을 수 있다.

$u^3$을 무시하지 않았을 때, 식은 좀 더 복잡해져서 (모두 전개한 뒤 $s^3$을 무시하면)



을 얻는다. 하지만 이 역시, 삼각함수의 주기를 조금 조정하면 적합한 해를 얻을 수 있다.



여기에서 주기 당 근일점의 이동량을 알기 위해 상수 $A, B, C, k$ 중 가장 중요한 건 $k$이다. 이는 $\Delta \phi = 2\pi/k$ 당 $r$의 한 주기가 끝남을 의미하며, 따라서 $r$이 최소가 되는 근일점에서 근일점으로 돌아옴을 의미한다. 따라서, $2\pi/k - 2\pi$가 주기 당 근일점의 회전 각도를 의미하게 되며, 이는 약



이다. 행성의 궤도가 타원과 유사하다면, 근일점 $r_{\mathrm P}$와 원일점 $r_{\mathrm A}$의 공식을 그대로 사용할 수 있다. 이 때, $\displaystyle \tilde L^2 \approx r_sa(1-e^2)/2$이므로



을 얻는다.

역사적으로, 수성의 근일점 세차운동은 1915년 11월 18일 아인슈타인이 처음으로 관련 계산을 해냄으로써 일반 상대성 이론이 뉴턴 이론에 대한 실험적 강점을 증명한 첫번째 증거였다. 태양계의 행성들에 대한 관련 계산 및 관측이 이루어졌으며, 그 결과는 다음과 같다.