ブラックホールのシミュレーション① 時間発展の式
「シミュレーションで学ぶ 相対論入門」を参考にブラックホールのシミュレーションを作っていきます。
目次
測地線方程式
シミュレーションのためにシュワルツシルトブラックホール周りの光子や質点の運動を求めたいので、それを記述する測地線方程式(時間発展方程式)を求めます。 シュワルツシルト計量は以下の通り。
はシュワルツシルト半径です。 ラグランジアンは
ただし、ドット()は固有時による微分を表しています。 は質点の場合1、光子の場合は0を表しています。
オイラー・ラグランジュ方程式から測地線方程式を求めます。 とは循環座標なので運動の定数(保存量)が存在し、それぞれと表すことにします。 まずはそれらから考えると良いです。 はひとまず式(1)を用いて次のように表します。
ニュートン力学的には角運動量が保存しているので運動面を (平面)に限定してもよいです。 実際、の微分方程式について、,を初期条件として解くと、 (一定)の解を得ます。 しかし、シミュレーションでは一つの平面に固定しない3次元的な動きを表現したいのでここではを固定しないで求めます。
また、固有時間と観測者にとっての時間での時間発展では運動のふるまいが異なることはみなさんご存知だと思います。 ここで固有時間は注目する物体に伴って移動する座標系で計測した時間を表し、 は遠方で静止している観測者の時間を表しています。 簡便のため以下では出来るだけを粒子の時間、を観測者の時間と呼ぶことにします。 他にもを座標時間と読んだりもします。 シミュレーションでは両方の立場で運動を表現したいので両方の時間発展方程式を求めます。
光子の計算では固有時間は0になってしまうので、ドット()は、ではなく適当なパラメータによる微分とします。 このパラメータを削除して積分すれば光子の軌跡を求めることができるが、時間発展を求めたいので固有時間による光子の運動は考えないことにします。 これらの式を使ってルンゲクッタ法により運動を求めることもできますが、 プログラムの仕様上の2階微分も求めます。
粒子の時間による時間発展方程式
質点の固有時間の時間発展は
初期条件は、を与えればよい。 保存量はから求めます。
観測者の時間による時間発展方程式
続いて遠方の静止系による表式を求めます。 遠方の静止形での時間として、のによる2階微分を求めます。 次の関係を利用して上記の式を変形していきます。
による微分をプライム()で表すことにすると以下のようになります。
ただし
としています。 初期条件は書籍による座標の設定では以下の図の通りで、
と表されます。 運動の定数は初期条件を使って、
と表せます。 ただし、光子の場合はです。 初期条件は、を与えればよいです。