解析力学から様々な理論分野を組み立てる 第2部 第12章 §2.12.1 ★★★ 作用 −mc∫ds の ημν を gμν(x) に替える。時計の進みが場所と速さで変わることは、この1行に入っている
第2部第3章では、自由粒子の作用を世界線に沿った固有時に −mc2 を掛けたもの
S=−mc∫ds,ds2=ημνdxμdxν
と選んだ(第2部第3章)。第2部第11章では、等価原理から、重力は時空の計量 gμν(x) の場所による違いとして書くしかないことを見た(第2部第11章)。この章で選ぶものは、この2つを合わせた1行である。
S[x]=−mc∫gμν(x)dλdxμdλdxνdλ.
重力の中の粒子は、力を受けて曲がるのではなく、この作用を停留させる世界線、つまり固有時を停留させる世界線を進む。この世界線を測地線と呼ぶ。
この作用に Euler–Lagrange 方程式を当てると、計量の1階微分の組み合わせとして Christoffel 記号 Γρμν が現れ、測地線方程式が出る。平方根の付いた L は第2部第3章で見たとおり径数の取り替えで変わらず、H が恒等的に 0 になる。径数を固有時に固定すると、平方根を外した L=21gμνx˙μx˙ν に置き換えてよいことが分かり、光の世界線もこの形で扱える。計量を Newton のポテンシャル Φ で g00≈1+2Φ/c2 と書くと、測地線方程式の先頭の項は x¨=−∇Φ で、L を展開すると −mc2+21mv2−mΦ が出る。計量が座標の取り替えで変わらないとき(Killing ベクトルがあるとき)には、第1部第10章の Noether の定理からエネルギーと角運動量の保存が出る。最後に、回転する座標系の Christoffel 記号が遠心力と Coriolis 力になることを確かめ、学部の授業で習う式と対応させる。
この章を通して計量の符号は (+,−,−,−)(Landau–Lifshitz と同じ。他の教科書との対応は第2部第11章の表)とし、上下に1回ずつ現れるギリシャ添字だけ和をとる。曲率(Riemann テンソル)はこの章には出てこない。測地線方程式は計量とその1階微分だけで書け、曲率は2階微分を含む量として第2部第14章で導入する。
計量と固有時
座標を xμ=(x0,x1,x2,x3) と書く。近い2つの出来事の間の世界間隔を
ds2=gμν(x)dxμdxν(2.12.1)
で定める。gμν は場所によってよい 4×4 の対称行列で、各点で固有値の符号が (+,−,−,−) になっているものとする。その行列式を g=detgμν<0 と書き、逆行列を gμν(gμλgλν=δμν)と書く。添字の上げ下げは g で行う(uμ=gμνuν)。平坦な時空を慣性系の座標で書いたときが gμν=ημν である。
ds2>0 の向き(時間的な向き)に進む粒子について、固有時を
dτ=cds=c1gμνdxμdxν(2.12.2)
で定める。第2部第3章と同じく、dτ は粒子が持ち歩く時計の進みにあたる。η を g に替えても「世界線の長さを時計が測る」という意味は変わらず、変わるのは長さの測り方が場所によることだけである。4元速度は uμ=dτdxμ で、式 (2.12.2) を dτ2 で割ると
gμνuμuν=c2(2.12.3)
を満たす。
選ぶ作用
定義 2.12.1(曲がった時空の自由粒子の作用)
質量 m>0 の粒子の世界線を任意の径数 λ で xμ(λ) と書き、x˙μ=dλdxμ とする。作用を
S[x]=−mc∫λ1λ2gμν(x)x˙μx˙νdλ=−mc2∫dτ,L=−mcgμν(x)x˙μx˙ν(2.12.4)
と選ぶ。両端の出来事 xμ(λ1), xμ(λ2) を固定して δS=0 となる世界線を時間的な測地線と呼ぶ。
この作用を選ぶ理由は2つある。1つ目は、gμν=ημν のとき第2部第3章の作用に戻ることである。2つ目は、∫dτ がどの座標で計算しても同じ値になることで、そのため停留の条件も座標によらない。座標を x→x′ と取り替えると gμν の成分は変わるが、ds2 が同じ値になるように変わる(成分の変わり方は第2部第13章で扱う)。この章で「重力」と呼んできたものは、すべて gμν(x) に入っている。L には重力のポテンシャルを別に足す項が無い。
静的な弱い場の計量
具体的な計量として、この章では次のものを繰り返し使う。
ds2=(1+c22Φ(x))c2dt2−(1−c22Φ(x))dx2,x0=ct.(2.12.5)
Φ は Newton のポテンシャルで、質量 M の球の外では Φ=−GM/r<0 である。g00=1+2Φ/c2 の部分は、第2部第11章の重力赤方偏移(高い所の時計が速く進む)を計量に書いたものになっている。gij の部分は等価原理だけでは決まらず、Einstein 方程式を弱い場で解いて得る(第2部第15章)。この章ではこの計量を与えられたものとして使い、∣Φ∣/c2 の1次までを残す。地表では GM/(Rc2)=7×10−10 程度で、太陽の表面でも 2×10−6 程度なので、1次で十分よい。
この計量の中を速度 v=dtdx で動く時計の進みは、式 (2.12.2) に式 (2.12.5) を入れて
dtdτ=1+c22Φ−(1−c22Φ)c2v2=1+c2Φ−2c2v2+O(c−4)(2.12.6)
となる。第1項の補正 Φ/c2 は低い所ほど時計が遅れることを、第2項 −v2/(2c2) は特殊相対論の時間の遅れを表す。2つの補正が同じ1行から出る。
円軌道を回る時計と地表の時計の比べ方(同じ速さになるのは r=3R/2)は、第2部第11章の式 (2.11.23) で扱った。ここでは自転する地球の表面に置いた時計どうしを比べる。
例題 2.12.1 ジオイドの上の時計はどこでも同じ速さで進む ★★
オリジナル地球は自転軸のまわりに角速度 Ω で回り、外の計量を式 (2.12.5) とする(Φ は自転しない慣性系で書いた重力のポテンシャルで、時刻によらないとする)。地表に置いた時計は、慣性系から見ると軸からの距離 ρ の円を速さ v=Ωρ で回る。
- (1) 地表に置いた時計の dtdτ を 1/c2 の1次まで求め、遠心力のポテンシャルを含む W=Φ−21Ω2ρ2 だけで決まることを示せ。海面(ジオイド)は W が一定の面なので、その上の時計はどこでも同じ速さで進むことを説明せよ。
- (2) 地球の扁平さを J2 の項まで入れた重力のポテンシャル
Φ(r,β)=−rGM[1−J2(ra)223sin2β−1]
(β は緯度)を使い、赤道(r=a)と極(r=b)で Φ と v2/2 を求めよ。数値は GM=3.986004418×1014m3s−2、a=6378137m、b=6356752.3m、J2=1.08263×10−3、Ω=7.292115×10−5s−1 とする1世界測地系 WGS84 の定数。National Imagery and Mapping Agency (2000), Department of Defense World Geodetic System 1984, Technical Report TR8350.2, 3rd ed.。3. (3) 赤道と極の W の差を、Φ の差・v2/2 の差と比べよ。時計の進みの割合の差は1日あたり何 ns か。
解答
(1) 式 (2.12.6) に v=Ωρ を入れて
dtdτ=1+c21(Φ−21Ω2ρ2)+O(c−4)=1+c2W+O(c−4).
重力による遅れ Φ/c2 と速さによる遅れ −v2/(2c2) が、W の1つにまとまった。W は地球と一緒に回る座標で見た重力と遠心力のポテンシャルの和である(回る座標の計量の g00 が遠心力のポテンシャルを含むことは例題 2.12.7 で確かめる)。回る座標で静止した海水の表面は、重力と遠心力の合力に垂直な面、つまり W が一定の面に落ち着く。したがって海面の上ではどこでも dtdτ が同じで、時計は同じ速さで進む。
(2) 赤道は β=0、r=a、ρ=a なので
Φ赤道=−aGM(1+2J2)=−6.252864×107m2s−2,
v=Ωa=465.1ms−1,2v2=1.0816×105m2s−2.
極は β=π/2、r=b、ρ=0 なので
Φ極=−bGM(1−J2b2a2)=−6.263670×107m2s−2,v=0.
(3) Φ赤道−Φ極=1.0807×105m2s−2 で、赤道の時計は重力のぶんだけ極の時計より 1.0807×105/c2=1.202×10−12 の割合で速く進む。一方、赤道の時計は速さのぶん v2/(2c2)=1.203×10−12 だけ遅れる。2つを合わせた
W赤道−W極=1.0807×105−1.0816×105=−94m2s−2
は、どちらの項に比べても 10−3 倍の大きさである。残りの 94 は、Φ を J2 の1次で打ち切ったことによるずれ(J22GM/a≈73m2s−2 の大きさ)と同じ程度で、この精度で W は一致している。時計の進みの割合の差は −94/c2=−1.05×10−15、1日あたり 0.09ns で、2つの項のそれぞれ(1日あたり約 104ns)よりずっと小さい。地球を半径の違う球とみなして J2 を落とすと、Φ の差は GM/b−GM/a=2.10×105m2s−2 になり v2/2 と釣り合わない。地球が扁平になったのは自転の遠心力のためで、その形が重力の項と速さの項をちょうど釣り合わせている。
1 世界測地系 WGS84 の定数。National Imagery and Mapping Agency (2000), Department of Defense World Geodetic System 1984, Technical Report TR8350.2, 3rd ed.
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com