運動方程式
アリスの運動を表す方程式を求める。この章ではボブの視点の座標系における計算結果を検証したいので、S(0)系やアリスの視点の座標系のことはいったん忘れて、わざわざ一般相対性理論の手法を使ってS(1)系で計算するのである。
アリスにはS(1)系で見て重力以外の力がかかっていないから、その運動を表す方程式は自由粒子の運動方程式であり測地線の方程式である。それは
のような2階微分方程式である。クリストッフェル記号
は前のページの最後で求めたのでそれを代入すればよい。左辺の見かけ上の第2項は添え字の 𝜇 と 𝜈 がそれぞれ上下にペアで現れるので、規約によりそれぞれに 0,1,2,3 を入れて和をとるので実体は16個の項である。𝜆 は上下にペアで現れないので和をとってはいけない。この式は 𝜆 = 0,1,2,3 の4個の場合の式をまとめて書いたものであり、実体は4個の連立方程式なのである。これを解けば
の4成分
を 𝜏𝐴 によるパラメータ表示で表すことができる。しかしこれだけだとパラメータ 𝜏𝐴 は自由なので、 𝜏𝐴 をアリスの固有時と一致させるために線素の式
も条件に加える。計量テンソル
は前のページの終わりのほうで求めたのでそれを代入すればよい。𝜇 と 𝜈 はそれぞれに 0,1,2,3 を入れて和をとるので右辺の実体は16個の項である。今から(65)・(66)式の実数解を求めるのである。また、これ以降、4元速度を
と書くこととする。(この微分方程式を解くだけの作業に興味がない人は、一般解が求まる表8まで飛ばしてもらって構わない。)
簡単なものから先に解く。まず(65)式の第2成分 (𝜆 = 2) は、
であるが、 𝜇 と 𝜈 が何であろうと はすべて0だから
となる。
次に(65)式の第3成分 (𝜆 = 3) は第2成分とまったく同様に、
となる。
次に(65)式の第0成分 (𝜆 = 0) は、
であり、 のうち0でないものは
と
だから
となる。 が恒等的に0であると仮定すればそれはこの式を満たす。そうでないと仮定すれば
で割って、
であるが、改めて ±𝐾₁₀″ = 𝐾₁₀ とおき、さらに が恒等的に0の場合と合わせれば、積分定数𝐾₁₀はすべての実数として
である。 すなわち
を求めるためにはもう一度積分する必要があるが、まだその計算はできない。その前に
すなわち
を求める。しかし(65)式の第1成分 (𝜆 = 1) を積分することは難しいので、先に(66)式を使う。(66)式に(61)式を代入すると、
のようになるが、 𝐾₁₀ = 0 だと右辺が常に負となり実数解がなくなってしまうので、𝐾₁₀ ≠ 0 である。
ここで(72)式の右辺が恒等的に0であるかないかによって場合分けをする。(72)式の右辺が恒等的に0であると仮定すると、
であるが、今 だから複号が付いている項は正でなければならないので、
である。この
は定数関数であるから(72)式の左辺も恒等的に0となり、(72)式を満たす。一方、(72)式の右辺が恒等的に0でないと仮定すれば、
であるが、今 だから複号は正を採用し、
である。微分すると4元速度の𝑥成分は
である。後の計算に備えてもう一度微分すると4元加速度の𝑥成分は
である。ここまでで(73)式と(74)式の2とおりの が得られたが、まだ(65)式の第1成分 (𝜆 = 1) の条件を使っていないのであった。それが満たされているかどうかを確認する。(65)式の第1成分は、
であり、 のうち0でないものは
だけだから
である。まず(73)式を(77)式の左辺に代入すると、
となって0にならないので、これは(77)式を満たさない。したがって(65)式の第1成分を満たさないので(73)式の は(65)式の解ではない。次に(74)式およびそれから導いた(76)式を(77)式の左辺に代入すると、
であるからこれは(77)式を満たす。したがって(65)式の第1成分を満たすので(74)式の は(65)式の解である。
すなわち
が求まったので、後回しにしていた
すなわち
を求める。(71)式に(74)式を代入すると、
である。最後の等号は微分の公式
を使った。
ここまでで2階微分方程式(65)式(に線素の条件(66)式を加えたもの)の一般解が求まったので表8にまとめておく。
表8. (65)・(66)式の一般解
式番号 | 一般解(0階微分と1階微分) |
(79) |
|
(74) |
|
(68) |
|
(70) |
|
(78)(71) |
|
(75) |
|
(67) |
|
(69) |
|
| (𝐾₁₀, 𝐾₁₁, 𝐾₁₂, 𝐾₁₃, 𝐾₁₄, 𝐾₁₅, 𝐾₁₆ は積分定数) |
なお と
の最後の等号は(74)式を使えば確認できる。
次に初期条件を適用して積分定数を定める。今の出発局面の場合、初期条件はその最初の時刻
のときのアリスの位置と4元速度と固有時である。ボブが加速を始める最初の瞬間まではアリスとボブの運動は共通であるから、
における値は
であり、またアリスの固有時も出発時を基準にしたいので のとき 𝜏𝐴 = 0 とする。
これらの初期条件のうち第2・3成分(𝑦, 𝑧 成分)と 𝜏𝐴 を(67)〜(70)式に代入すれば、簡単な計算により 𝐾₁₃ = 𝐾₁₄ = 𝐾₁₅ = 𝐾₁₆ = 0 のように定まる。それを表8の各式に代入すると次のようになる。
第2・3成分(𝑦, 𝑧 成分)はずっと 0 であることがわかったので、これ以降この節では書かない。 さらに初期条件のうち第0・1成分(𝑐𝑡, 𝑥 成分)と 𝜏𝐴 を(80)〜(83)式に代入すれば、簡単な計算により
, 𝐾₁₁ = 𝐾₁₂ = 0 のように定まる。
積分定数の値を代入する前に、後で使うために(80)〜(83)式を固有時 𝜏𝐴 を使ったパラメータ表示から座標時
の関数としての表示に変形する。(80)式を 𝜏𝐴 について解くと、
である。(81)式より、
である。(82)式より、
である。(83)式より、
である。以上の(85)・(86)・(88)・(89)式が座標時 の関数としての表示である。
先ほど求めた積分定数の値をこれらに代入すれば、S(1)系で見たアリスの固有時と位置と4元速度は次のようになる。
今は
の範囲のみを考えることにしているが、(91)式より確かにその範囲に収まっていることがわかる。(90)式より、 𝜏𝐴 は
で頭打ちになり、S(1)系で測ると が大きくなったときアリスの時間はほぼ止まってしまうことになる。
出発局面の終了時では、 であるからこれを(90)〜(93)式に代入すると、アリスの固有時は
、位置は
、4元速度は
である。