tecktom.com

この章のまとめ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第2部 第16章 Schwarzschild 解と太陽系の検証

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第16章 §2.16.7 ★★★ Rμν=0R_{\mu\nu}=0 を解いて計量を決め、その上で第1部の Noether と有効ポテンシャルを当てると、太陽系の検証がすべて rs/rr_{\mathrm s}/r の1次の量として出る

この章で選んだものは、物質の無い球対称な領域の真空の Einstein 方程式 Rμν=0R_{\mu\nu}=0 である。球対称と無限遠での平坦さを課すと、Ricci テンソルの3成分から α+β=0\alpha+\beta=0 と (re−2β)′=1(r\mathrm{e}^{-2\beta})'=1 の2つの式が出て、積分定数を Newton 極限で決めると Schwarzschild 計量 (2.16.10) が得られた。Birkhoff の定理により、球対称なら静的という仮定は要らない。

この計量の上の運動には、第2部第12章の Lagrangian L=12gμνx˙μx˙νL=\frac12g_{\mu\nu}\dot x^\mu\dot x^\nu を当てた。計量が tt と φ\varphi によらないので、Noether の定理でエネルギー E~\tilde E と角運動量 ℓ\ell が保存し、動径の運動は有効ポテンシャル Veff=−GM/r+ℓ2/(2r2)−GMℓ2/(c2r3)V_{\mathrm{eff}}=-GM/r+\ell^2/(2r^2)-GM\ell^2/(c^2r^3) の1次元運動になる。Newton との違いは最後の1項だけで、軌道の式では 3GMu2/c23GMu^2/c^2 になる。この項を摂動として扱うと、質量のある粒子では共鳴によって永年項が出て、それが近日点の回転 6πGM/[c2a(1−e2)]6\pi GM/[c^2a(1-e^2)] を表す(水星で 1 世紀に 42.98′′42.98''、観測 42.56′′±0.94′′42.56''\pm0.94'')。光では共鳴が無く、無限遠の向きがずれるだけで、そのずれが曲がり角 4GM/(c2b)4GM/(c^2b) になる(太陽の縁で 1.75′′1.75'')。光の曲がりは、屈折率 n=1+2GM/(c2r)n=1+2GM/(c^2r) の Fermat の原理からも出て、その半分は g00g_{00}、半分は空間の曲がり gijg_{ij} による。同じ屈折率から Shapiro 遅延と重力レンズのレンズ方程式も出る。

r=rsr=r_{\mathrm s} と r=3GM/c2r=3GM/c^2 では何が起きるか、つまり Schwarzschild 半径が座標の特異点にすぎないこと、光子球、円軌道が不安定になる 6GM/c26GM/c^2 は、第2部第17章で扱う。連星パルサーの近星点の移動は同じ (2.16.21) の強い場の版で、公転周期の減少とともに第2部第18章で扱う。

総合問題 地球の軌道の時計と近日点を、円に近い軌道の振動から読む ★★★

オリジナル太陽(GM⊙=1.32712440018×1020 m3 s−2GM_\odot=1.32712440018\times10^{20}\,\mathrm{m^3\,s^{-2}})のまわりの Schwarzschild 計量 (2.16.10) の赤道面で、半径 rr の円軌道とそれに近い軌道を考える。x≔GM/(c2r)x\coloneqq GM/(c^2r) とおく。

  1. (1) (使う道具:第2部第12章・2.16.2 節)L=12gμνx˙μx˙νL=\frac12g_{\mu\nu}\dot x^\mu\dot x^\nu から E~\tilde E と ℓ\ell が保存することを示し、円軌道の ℓ2\ell^2 と、公転の角振動数 ωφ≔φ˙\omega_\varphi\coloneqq\dot\varphi(固有時あたり)を xx で表せ。
  2. (2) (使う道具:第2部第11章・例題 2.16.2)円軌道を回る時計の進み dτdt\frac{\mathrm{d} \tau}{\mathrm{d} t} が 1−3x\sqrt{1-3x} であることを示せ。地球の軌道(r=1 au=1.495978707×1011 mr=1\,\mathrm{au}=1.495978707\times10^{11}\,\mathrm m、地球自身の重力は無視)の時計は、遠方の時計に比べて1年(3.15576×107 s3.15576\times10^7\,\mathrm s)で何秒遅れるか。
  3. (3) (使う道具:第2部第1章・第1部第8章)円軌道のまわりの小さな動径の振動の角振動数 κ\kappa(固有時あたり)を κ2=Veff′′(r)\kappa^2=V_{\mathrm{eff}}''(r) から求め、κ2/ωφ2=1−6x\kappa^2/\omega_\varphi^2=1-6x が厳密に成り立つことを示せ。
  4. (4) (使う道具:2.16.3 節)(3) から、円に近い軌道の近日点が1周に 2π[(1−6x)−1/2−1]2\pi[(1-6x)^{-1/2}-1] 進むことを導き、x→0x\to0 で (2.16.21) の e→0e\to0 の値と一致することを確かめよ。地球の軌道(a=1.00000261 aua=1.00000261\,\mathrm{au}、e=0.0167086e=0.0167086、公転周期 365.256363365.256363 日)で、1 世紀(3652536525 日)あたりの近日点の進みを (2.16.21) で求めよ。
  5. (5) (使う道具:第2部第17章への予告)κ2\kappa^2 が負になる半径を求め、そのとき円軌道に何が起きるか述べよ。1.41.4 倍の太陽質量の中性子星では、その半径は何 km か。
解答

(1) LL は tt と φ\varphi を含まないので、Euler–Lagrange 方程式 ddτ∂L∂t˙=∂L∂t=0\frac{\mathrm{d} }{\mathrm{d} \tau}\frac{\partial L}{\partial \dot t}=\frac{\partial L}{\partial t}=0、ddτ∂L∂φ˙=0\frac{\mathrm{d} }{\mathrm{d} \tau}\frac{\partial L}{\partial \dot\varphi}=0 から∂L∂t˙=(1−rs/r)c2t˙=E~\frac{\partial L}{\partial \dot t}=(1-r_{\mathrm s}/r)c^2\dot t=\tilde E と −∂L∂φ˙=r2φ˙=ℓ-\frac{\partial L}{\partial \dot\varphi}=r^2\dot\varphi=\ell が一定である(第1部第10章の循環座標の Noether の定理)。円軌道は例題 2.16.2 の (1) より ℓ2=GMr/(1−3x)\ell^2=GMr/(1-3x)。ωφ=ℓ/r2\omega_\varphi=\ell/r^2 なので

ωφ2=ℓ2r4=GMr311−3x.\omega_\varphi^2=\frac{\ell^2}{r^4}=\frac{GM}{r^3}\frac{1}{1-3x}.

(2) 規格化 (2.16.11) を t˙2\dot t^2 で割ると、円軌道(r˙=0\dot r=0)では (1−2x)c2−r2Ω2=c2(dτ/dt)2(1-2x)c^2-r^2\Omega^2=c^2(\mathrm{d}\tau/\mathrm{d} t)^2。例題 2.16.2 の (2) の Ω2=GM/r3\Omega^2=GM/r^3、つまり r2Ω2=c2xr^2\Omega^2=c^2x を入れて (dτ/dt)2=1−3x(\mathrm{d}\tau/\mathrm{d} t)^2=1-3x。このうち 2x2x は重力の位置による遅れ(第2部第11章の重力赤方偏移)、xx は速さ v2=GM/rv^2=GM/r による特殊相対論の遅れである。地球の軌道では x=1476.6/1.495979×1011=9.871×10−9x=1476.6/1.495979\times10^{11}=9.871\times10^{-9} で、1−dτ/dt≈32x=1.481×10−81-\mathrm{d}\tau/\mathrm{d} t\approx\frac32x=1.481\times10^{-8}。1年で 1.481×10−8×3.15576×107=0.467 s1.481\times10^{-8}\times3.15576\times10^7=0.467\,\mathrm s 遅れる。

(3) Veff′′(r)=−2GMr3+3ℓ2r4−12GMℓ2c2r5V_{\mathrm{eff}}''(r)=-\frac{2GM}{r^3}+\frac{3\ell^2}{r^4}-\frac{12GM\ell^2}{c^2r^5}。ℓ2/r4=GMr311−3x\ell^2/r^4=\frac{GM}{r^3}\frac{1}{1-3x} を入れると

κ2=Veff′′(r)=GMr3[−2+3−12x1−3x]=GMr31−6x1−3x.\kappa^2=V_{\mathrm{eff}}''(r)=\frac{GM}{r^3}\Bigl[-2+\frac{3-12x}{1-3x}\Bigr]=\frac{GM}{r^3}\frac{1-6x}{1-3x}.

(1) の ωφ2\omega_\varphi^2 で割ると κ2/ωφ2=1−6x\kappa^2/\omega_\varphi^2=1-6x で、近似を使っていない。Newton の場合(x→0x\to0)は κ=ωφ\kappa=\omega_\varphi で、動径の振動1回と公転1回が同じ時間で終わるので軌道が閉じる(第1部第8章の ωr=ωθ\omega_r=\omega_\theta)。

(4) 動径の振動の1周期 2π/κ2\pi/\kappa の間に φ\varphi は 2πωφ/κ=2π(1−6x)−1/22\pi\omega_\varphi/\kappa=2\pi(1-6x)^{-1/2} 進む。近日点から近日点までの角度がこれなので、1周あたりの進みは 2π[(1−6x)−1/2−1]≈6πx2\pi[(1-6x)^{-1/2}-1]\approx6\pi x。e→0e\to0 で a(1−e2)→ra(1-e^2)\to r なので (2.16.21) は 6πGM/(c2r)=6πx6\pi GM/(c^2r)=6\pi x となり一致する。地球の軌道では1周に 0.0384′′0.0384''、1 世紀に 100.0100.0 周するので 3.84′′3.84'' である。

(5) κ2<0\kappa^2<0 となるのは 1−6x<01-6x<0、すなわち r<6GM/c2=3rsr<6GM/c^2=3r_{\mathrm s}。この内側では円軌道のまわりの小さなずれが指数関数的に増え、円軌道は不安定になる。r=6GM/c2r=6GM/c^2 は安定な円軌道の最も内側(最内安定円軌道)で、第2部第17章で扱う。1.4M⊙1.4M_\odot の中性子星では 6×1.4×1476.6 m=12.4 km6\times1.4\times1476.6\,\mathrm m=12.4\,\mathrm{km} で、太陽では 8.9 km8.9\,\mathrm{km} と太陽の内部に埋もれて意味を持たない。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ