摘要

描述Newton versus Schwarzschild trajectories.gif	English: Comparison of a testparticle's trajectory in Newtonian and Schwarzschild spacetime in the strong gravitational field (r0=10rs=20GM/c²). The initial velocity in both cases is 126% of the circular orbital velocity. φ0 is the launching angle (0° is a horizontal shot, and 90° a radially upward shot). Since the metric is spherically symmetric the frame of reference can be rotated so that Φ is constant and the motion of the test-particle is confined to the r,θ-plane (or vice versa).
日期	2016年五月21日
來源	自己的作品 - Mathematica Code
作者	Yukterez (Simon Tyran, Vienna)
其他版本

Equations of motion

Newton

In spherical coordinates and natural units of ${\displaystyle {\text{G=M=c=1}}}$, where lengths are measured in ${\displaystyle {\text{GM/c}}^{2}}$ and times in ${\displaystyle {\text{GM/c}}^{3}}$, the motion of a testparticle in the presence of a dominant mass is defined by

${\displaystyle {\ddot {\text{r}}}=-{\frac {1}{{\text{r}}^{2}}}+{\text{r}}\ {\dot {\phi }}^{2}\ ,\ \ {\ddot {\phi }}=-{\frac {2\ {\dot {\text{r}}}\ {\dot {\phi }}}{\text{r}}}}$

The initial conditions are

${\displaystyle {\dot {\text{r}}}_{0}=v_{0}\sin(\varphi _{0})\ ,\ \ {\dot {\phi }}_{0}={\frac {v_{0}}{{\text{r}}_{0}}}\cos(\varphi _{0})}$

The overdot stands for the time-derivative. ${\displaystyle \phi }$ is the angular coordinate, ${\displaystyle \varphi }$ the local elevation angle of the test particle, and ${\displaystyle v}$ it's velocity.

${\displaystyle {\rm {E}}}$ and ${\displaystyle {\rm {L}}}$, where the kinetic ${\displaystyle {\text{E}}_{\rm {kin}}=v^{2}/2}$ and potential ${\displaystyle {\text{E}}_{\rm {pot}}=-1/{\rm {r}}}$ component (all in units of ${\displaystyle {\rm {mc}}^{2}}$) give the total energy ${\displaystyle {\rm {E}}={\text{E}}_{\rm {kin}}+{\text{E}}_{\rm {pot}}}$, and the angular momentum, which is given by ${\displaystyle {\rm {L}}=v_{\perp }\ {\rm {r}}}$ (in units of ${\displaystyle {\rm {m}}}$) where ${\displaystyle v_{\perp }={\dot {\phi }}\ {\text{r}}}$ is the transverse and ${\displaystyle v_{\parallel }={\dot {\text{r}}}}$ the radial velocity component, are conserved quantities.

Schwarzschild

The equations of motion ^[1] in Schwarzschild-coordinates are

${\displaystyle {\ddot {\text{r}}}=-{\frac {1}{{\text{r}}^{2}}}+{\text{r}}\ {\dot {\phi }}^{2}-3\ {\dot {\phi }}^{2}\ ,\ \ {\ddot {\phi }}=-{\frac {2\ {\dot {\text{r}}}\ {\dot {\phi }}}{\text{r}}}}$

which is except for the ${\displaystyle -3\ {\dot {\phi }}^{2}}$ term identical with Newton, although the radial coordinate has a different meaning (see farther below). The time dilation is

${\displaystyle {\dot {\text{t}}}={\frac {\sqrt {1+{\dot {\text{r}}}^{2}\ {\text{r}}/({\text{r}}-2)+{\text{r}}^{2}\ {\dot {\phi }}^{2}}}{\sqrt {1-2/{\text{r}}}}}={\frac {1}{{\sqrt {1-2/{\text{r}}}}{\sqrt {1-v^{2}}}}}}$

The coordinates are differentiated by the test particle's proper time ${\displaystyle \tau }$, while ${\displaystyle {\text{t}}}$ is the coordinate time of the bookkeeper at infinity. So the total coordinate time ellapsed between the proper time interval

${\displaystyle \tau _{0}..\tau _{1}}$ is ${\displaystyle {\text{t}}=\int _{\tau _{0}}^{\tau _{1}}{\dot {\text{t}}}\,\mathrm {d} \tau }$

The local velocity ${\displaystyle v}$ (relative to the main mass) and the coordinate celerity are related by

${\displaystyle {\dot {\phi }}={\frac {v_{\perp }}{{\text{r}}\ {\sqrt {1-v^{2}}}}}}$ for the input and ${\displaystyle v_{\perp }={\frac {{\dot {\text{r}}}\ {\dot {\phi }}}{{\dot {\text{t}}}\ {\sqrt {1-2/{\text{r}}}}}}}$ for the output of the transverse ${\displaystyle (\perp )}$ and

${\displaystyle {\dot {r}}=v_{\parallel }{\sqrt {\frac {1-2/{\text{r}}}{1-v^{2}}}}}$ or the other way around ${\displaystyle v_{\parallel }={\frac {\dot {\text{r}}}{{\dot {\text{t}}}\ (1-2/{\text{r}})}}}$ for the radial ${\displaystyle (\parallel )}$ component of motion.

The shapiro-delayed velocity ${\displaystyle {\text{v}}}$ in the bookeeper's frame of reference is

${\displaystyle {\text{v}}_{\perp }=v_{\perp }{\sqrt {1-2/{\text{r}}}}}$ and ${\displaystyle {\text{v}}_{\parallel }=v_{\parallel }(1-2/{\text{r}})}$

The initial conditions in terms of the local physical velocity ${\displaystyle v}$ are therefore

${\displaystyle {\dot {\text{r}}}_{0}={\frac {v_{0}{\sqrt {1-2/{\text{r}}_{0}}}\sin(\varphi _{0})}{\sqrt {1-v_{0}^{2}}}}\ ,\ \ {\dot {\phi }}_{0}={\frac {v_{0}\cos(\varphi _{0})}{{\text{r}}_{0}{\sqrt {1-v_{0}^{2}}}}}}$

The horizontal and vertical components differ by a factor of ${\displaystyle {\sqrt {1-2/{\text{r}}}}}$

because additional to the gravitational time dilation there is also a radial length contraction of the same factor, which means that the physical distance between

${\displaystyle {\text{r}}_{1}}$ and ${\displaystyle {\text{r}}_{2}}$ is not ${\displaystyle {\text{r}}_{2}-{\text{r}}_{1}}$ but ${\displaystyle \int _{{\text{r}}_{1}}^{{\text{r}}_{2}}{\frac {\mathrm {d} {\text{r}}}{\sqrt {1-2/{\text{r}}}}}}$

due to the fact that space around a mass is not euclidean, and a shell of a given diameter contains more volume when a central mass is present than in the absence of a such.

The angular momentum

${\displaystyle {\text{L}}={\text{r}}^{2}\ {\dot {\phi }}=v_{\perp }\ {\rm {r}}/{\sqrt {1-v^{2}}}}$

in units of ${\displaystyle {\text{m}}}$ and the total energy as the sum of rest-, kinetic- and potential energy

${\displaystyle {\text{E}}=1+{\text{E}}_{\rm {kin}}+{\text{E}}_{\rm {pot}}={\dot {\text{t}}}\ \left(1-{\frac {2}{\text{r}}}\right)={\frac {\sqrt {1-2/{\text{r}}}}{\sqrt {1-v^{2}}}}}$

in units of ${\displaystyle {\text{mc}}^{2}}$, where ${\displaystyle {\text{m}}}$ is the test particle's restmass, are the constants of motion. The components of the total energy are

${\displaystyle {\text{E}}_{\rm {kin}}={\frac {1}{\sqrt {1-v^{2}}}}-1}$ for the kinetic plus ${\displaystyle {\text{E}}_{\rm {pot}}={\frac {{\sqrt {1-2/{\text{r}}}}-1}{\sqrt {1-v^{2}}}}}$ for the potential energy plus ${\displaystyle {\text{m}}=1}$, the test particle's invariant rest mass.

The equations of motion in terms of ${\displaystyle {\text{E}}}$ and ${\displaystyle {\text{L}}}$ are

${\displaystyle {\dot {\rm {r}}}=\xi \ ,\ \ {\dot {\phi }}={\frac {\text{L}}{{\text{m}}\ {\text{r}}^{2}}}\ ,\ \ {\dot {\text{t}}}={\frac {\text{E}}{{\text{m}}\ (1-2/{\text{r}})}}}$

or, differentiated by the coordinate time ${\displaystyle {\text{t}}}$

${\displaystyle {\bar {\text{r}}}=\xi \ ,\ \ {\bar {\phi }}={\frac {(1-2/{\text{r}})\ {\text{L}}}{{\text{E}}\ {\text{r}}^{2}}}\ ,\ \ {\bar {\tau }}=1/{\dot {\text{t}}}}$

with

${\displaystyle \xi =\pm {\sqrt {{\frac {{\text{E}}^{2}\ {\text{r}}^{4}}{{\text{L}}^{2}}}-\left(1-{\frac {2}{\text{r}}}\right)\left({\frac {{\text{m}}^{2}\ {\text{r}}^{4}}{{\text{L}}^{2}}}+{\text{r}}^{2}\right)}}}$

where in contrast to the overdot, which stands for ${\displaystyle {\dot {x}}={\rm {d}}x/{\rm {d}}{\tau }}$, the overbar denotes ${\displaystyle {\bar {x}}={\rm {d}}x/{\rm {d}}{\text{t}}}$.

For massless particles like photons ${\displaystyle {\rm {m}}/{\sqrt {1-v^{2}}}}$ in the formula for ${\displaystyle {\rm {E}}}$ and ${\displaystyle {\rm {L}}}$ is replaced with ${\displaystyle {\rm {h\ f}}}$ and the ${\displaystyle {\rm {m}}}$ in the equations of motion set to ${\displaystyle 1}$, with ${\displaystyle {\rm {h}}}$ as Planck's constant and ${\displaystyle {\rm {f}}}$ for the photon's frequency.

授權條款

我，本作品的著作權持有者，決定用以下授權條款發佈本作品：

此檔案採用共享創意 姓名標示-相同方式分享 4.0 國際授權條款。

您可以自由：

• 分享 – 複製、發佈和傳播本作品
• 重新修改 – 創作演繹作品

惟需遵照下列條件：

• 姓名標示 – 您必須指名出正確的製作者，和提供授權條款的連結，以及表示是否有對內容上做出變更。您可以用任何合理的方式來行動，但不得以任何方式表明授權條款是對您許可或是由您所使用。
• 相同方式分享 – 如果您利用本素材進行再混合、轉換或創作，您必須基於如同原先的相同或兼容的條款，來分布您的貢獻成品。

https://creativecommons.org/licenses/by-sa/4.0CC BY-SA 4.0 Creative Commons Attribution-Share Alike 4.0 truetrue

References

Usage in Wikipedia Articles

en.wikipedia.org

• Schwarzschild metric
• Schwarzschild geodesics
• Two-body problem in general relativity

de.wikipedia.org

• Apsidendrehung
• Schwarzschild-Metrik

ru.wikipedia.org

• Метрика Шварцшильда
• Задача Кеплера в относительности

es.wikipedia.org

• La relatividad general

zh.wikipedia.org

• 史瓦西度規