tecktom.com

選ぶもの:平坦な計量のまわりの2次の作用

2026年10月15日16分で読めますAI生成教科書『虎の巻』第2部 第18章 重力波
目次
  1. 展開する作用
  2. 2次の作用
  3. 座標の取り替えがゲージ変換になる
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第18章 §2.18.1 ★★★ g=η+hg=\eta+h と分けて Einstein–Hilbert 作用を hh の2次まで展開すると、hμνh_{\mu\nu} を平坦な時空の上の場とする作用(Fierz–Pauli の形)が得られる

第2部第15章で、重力場そのものの作用として Einstein–Hilbert 作用を選び、その変分から Einstein 方程式を出した。Einstein 方程式は計量 gμνg_{\mu\nu} について非線形で、一般には解けない。ところが太陽系の外を伝わってくる重力の波は、地上に届いたときには計量を平坦な ημν\eta_{\mu\nu} から 10−2110^{-21} 程度しかずらさない。そこで gμν=ημν+hμνg_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu} と分け、作用を hμνh_{\mu\nu} の2次まで展開すると、hμνh_{\mu\nu} を「平坦な時空の上の場」とする線形な場の理論が得られる。この章では、その2次の作用を選ぶ。

当てる道具は、第2部第5章の場の Euler–Lagrange 方程式、第2部第8章のゲージ不変性と拘束による自由度の数え方、第2部第7章のエネルギー運動量テンソル、第2部第14章の測地線偏差の方程式である。2次の作用は座標の取り替えに対応するゲージ変換 hμν→hμν−∂μξν−∂νξμh_{\mu\nu}\to h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu で変わらず、これが電磁場の Aμ→Aμ−∂μΛA_\mu\to A_\mu-\partial_\mu\Lambda とまったく同じ役を果たす。その結果、hμνh_{\mu\nu} の10成分のうち波として伝わる自由度は2つ(h+h_+ と h×h_\times の2つの偏極)しか残らない。

出てくるものは3つある。①横波でトレースが 00 のゲージ(TT ゲージ)での2つの偏極と、それが自由落下する質点の輪を変形させるしかた。②重力波が運ぶエネルギーの流れと、源の質量の四重極モーメントの3階微分で放射の強さが決まる四重極公式。③四重極公式を連星に当てた結果として、公転周期の減少(Hulse–Taylor の連星パルサー)と、合体の直前に周波数が上がるチャープ(GW150914)。最後に、電磁気学の双極子放射・Lorenz ゲージ・Newton の重力のポテンシャルという授業での書き方を、この章の式と並べる。

単位は SI、計量は ημν=diag(+1,−1,−1,−1)\eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1)、xμ=(ct,x)x^\mu=(ct,\bm{x})、□=∂μ∂μ=1c2∂t2−∇2\Box=\partial_\mu\partial^\mu=\frac{1}{c^2}\partial_t^2-\nabla^2 で、上下に1回ずつ現れるギリシャ添字は和を取る(記法の正本どおり)。添字の上げ下げは ημν\eta_{\mu\nu} で行う。観測で測る「ひずみ」は、線素を ds2=c2dt2−(δij+h^ij)dxidxj\mathrm{d} s^2=c^2\mathrm{d} t^2-(\delta_{ij}+\hat h_{ij})\mathrm{d} x^i\mathrm{d} x^j と書いたときの h^ij\hat h_{ij} で表す。この本の流儀では hij=−h^ijh_{ij}=-\hat h_{ij} である(gij=−δij+hijg_{ij}=-\delta_{ij}+h_{ij} のため)。添字の付いた hμνh_{\mu\nu}・h+h_+・h×h_\times と Planck 定数 ℏ\hbar は字が違う。

展開する作用

第2部第15章で選んだ作用は、重力場の部分と物質の部分の和

S=Sg+Sm,Sg=−c316πG∫R−g d4x,Tμν=2c−gδSmδgμν(2.18.1)S=S_{\mathrm g}+S_{\mathrm m},\qquad S_{\mathrm g}=-\frac{c^3}{16\pi G}\int R\sqrt{-g}\,\mathrm{d}^4x,\qquad T_{\mu\nu}=\frac{2c}{\sqrt{-g}}\frac{\delta S_{\mathrm m}}{\delta g^{\mu\nu}}\tag{2.18.1}

である。この変分から Einstein 方程式 Gμν=8πGc4TμνG_{\mu\nu}=\frac{8\pi G}{c^4}T_{\mu\nu} が出る。ここで計量を

gμν=ημν+hμν,∣hμν∣≪1(2.18.2)g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu},\qquad \lvert h_{\mu\nu}\rvert \ll1\tag{2.18.2}

と分ける。逆計量は1次までで gμν=ημν−hμνg^{\mu\nu}=\eta^{\mu\nu}-h^{\mu\nu}(hμν=ημαηνβhαβh^{\mu\nu}=\eta^{\mu\alpha}\eta^{\nu\beta}h_{\alpha\beta})で、実際 gμαgαν=δμν+O(h2)g^{\mu\alpha}g_{\alpha\nu}=\delta^\mu{}_\nu+\mathcal{O}(h^2) になる。hμνh_{\mu\nu} の跡を h=ημνhμνh=\eta^{\mu\nu}h_{\mu\nu} と書く。

Christoffel 記号は hh の1次で

Γρμν=12ηρσ(∂μhσν+∂νhσμ−∂σhμν)(2.18.3)\Gamma^\rho{}_{\mu\nu}=\tfrac12\eta^{\rho\sigma}(\partial_\mu h_{\sigma\nu}+\partial_\nu h_{\sigma\mu}-\partial_\sigma h_{\mu\nu})\tag{2.18.3}

になり、Riemann テンソルの定義(記法の正本 §4)の ΓΓ\Gamma\Gamma の項は2次なので落ちる。1次の Riemann テンソルを縮約して、1次の Einstein テンソルは

Gμν(1)=12(∂ρ∂μhνρ+∂ρ∂νhμρ−□hμν−∂μ∂νh−ημν∂ρ∂σhρσ+ημν□h)(2.18.4)G^{(1)}_{\mu\nu}=\tfrac12\bigl(\partial^\rho\partial_\mu h_{\nu\rho}+\partial^\rho\partial_\nu h_{\mu\rho}-\Box h_{\mu\nu}-\partial_\mu\partial_\nu h -\eta_{\mu\nu}\partial^\rho\partial^\sigma h_{\rho\sigma}+\eta_{\mu\nu}\Box h\bigr)\tag{2.18.4}

になる。この式は (−,+,+,+) の教科書の式と文字の上でまったく同じ形をしている。hμνh_{\mu\nu}・ημν\eta_{\mu\nu}・□\Box の符号がそろって逆になり、各項で打ち消し合うためである(第2部第11章の対応表で GμνG_{\mu\nu} が両方の流儀で同じであることと整合する)。

2次の作用

SgS_{\mathrm g} を hh の2次まで展開する計算は長いので、結果を書き、それが正しいことを「変分すると Gμν(1)G^{(1)}_{\mu\nu} が出る」ことで確かめる。時間あたり・体積あたりの Lagrangian 密度で書くと(S=∫L dt d3xS=\int\mathcal L\,\mathrm{d} t\,\mathrm{d}^3x)、

L2=c464πG(∂λhμν ∂λhμν−2 ∂μhμλ ∂νhνλ+2 ∂μhμν ∂νh−∂λh ∂λh)(2.18.5)\mathcal L_2=\frac{c^4}{64\pi G}\Bigl(\partial_\lambda h_{\mu\nu}\,\partial^\lambda h^{\mu\nu}-2\,\partial_\mu h^{\mu\lambda}\,\partial^\nu h_{\nu\lambda} +2\,\partial_\mu h^{\mu\nu}\,\partial_\nu h-\partial_\lambda h\,\partial^\lambda h\Bigr)\tag{2.18.5}

である。この形は Fierz と Pauli が質量の無いスピン2の場の作用として書いたものと同じ組み合わせで、ここではそれが Einstein–Hilbert 作用の2次の部分として出てくる。1次の項は全微分になって運動方程式に寄与せず、0次の項は η\eta の曲率が 00 なので消える。

物質の作用の側は、δgμν=−δhμν\delta g^{\mu\nu}=-\delta h^{\mu\nu} と TμνT_{\mu\nu} の定義 (2.18.1) から、hh の1次で

Lint=−12Tμνhμν(2.18.6)\mathcal L_{\mathrm{int}}=-\tfrac12T_{\mu\nu}h^{\mu\nu}\tag{2.18.6}

になる。これは電磁場の結合 −JμAμ-J^\mu A_\mu と同じ形で、エネルギー運動量テンソルが重力の「源の流れ」の役をしている。

L2+Lint\mathcal L_2+\mathcal L_{\mathrm{int}} を hμνh^{\mu\nu} について変分すると、場の Euler–Lagrange 方程式(第2部第5章)は

c416πGGμν(1)−12Tμν=0⟺Gμν(1)=8πGc4Tμν(2.18.7)\frac{c^4}{16\pi G}G^{(1)}_{\mu\nu}-\frac12T_{\mu\nu}=0\quad\Longleftrightarrow\quad G^{(1)}_{\mu\nu}=\frac{8\pi G}{c^4}T_{\mu\nu}\tag{2.18.7}

となり、Einstein 方程式の1次の部分(線形化した Einstein 方程式)に一致する(検算は verify で一般の hμνh_{\mu\nu} について行った)。

座標の取り替えがゲージ変換になる

座標を x′μ=xμ+ξμ(x)x'^\mu=x^\mu+\xi^\mu(x) と少しずらすと、計量の成分は gμν′(x′)=∂xα∂x′μ∂xβ∂x′νgαβ(x)g'_{\mu\nu}(x')=\frac{\partial x^\alpha}{\partial x'^\mu}\frac{\partial x^\beta}{\partial x'^\nu}g_{\alpha\beta}(x) で変わる。ξ\xi と hh がともに小さいとして1次を残すと

hμν ⟶ hμν−∂μξν−∂νξμ(ξμ=ημνξν)(2.18.8)h_{\mu\nu}\ \longrightarrow\ h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu\qquad(\xi_\mu=\eta_{\mu\nu}\xi^\nu)\tag{2.18.8}

になる。同じ物理的な時空を別の座標で書いただけなので、測れる量はこの変換で変わってはならない。実際、1次の Riemann テンソルは (2.18.8) で変わらず、したがって Gμν(1)G^{(1)}_{\mu\nu} も2次の作用 (2.18.5) も(全微分を除いて)変わらない。電磁場の Aμ→Aμ−∂μΛA_\mu\to A_\mu-\partial_\mu\Lambda では任意関数が1つだったのに対し、ここでは任意関数が4つ(ξ0,…,ξ3\xi_0,\dots,\xi_3)ある。第2部第8章の Noether の第2定理をそのまま当てると、運動方程式の左辺の発散が恒等的に 00 になる、つまり

∂μGμν(1)=0(すべての hμν について)(2.18.9)\partial^\mu G^{(1)}_{\mu\nu}=0\qquad(\text{すべての }h_{\mu\nu}\text{ について})\tag{2.18.9}

が成り立つ。これは縮約した Bianchi 恒等式の1次の部分で、方程式 (2.18.7) が解を持つには ∂μTμν=0\partial^\mu T_{\mu\nu}=0 が要る。電荷の保存が Maxwell 方程式の整合条件だったのと同じ構造である(第2部第8章)。

静的な弱い場 h00=2Φ/c2h_{00}=2\Phi/c^2、hij=2Ψδij/c2h_{ij}=2\Psi\delta_{ij}/c^2 では、L2\mathcal L_2 が 第2部第15章の式 (2.15.21) になり、Ψ=Φ\Psi=\Phi とおくと Newton の重力の作用 −(∇Φ)2/(8πG)-(\nabla\Phi)^2/(8\pi G) に、結合 (2.18.6) は −ρΦ-\rho\Phi に戻る。これは電磁気のスカラーポテンシャル ϕ\phi との対応である。結合 −12Tμνhμν-\frac12T_{\mu\nu}h^{\mu\nu} には物質の流れ T0iT_{0i} と h0ih_{0i} の項もあり、これがベクトルポテンシャル A\bm{A} との対応になる。次の例題でそれを確かめる。

例題 2.18.1 ゆっくり回る球殻のまわりの重力磁気 ★★

オリジナル質量 MM・半径 RR の薄い球殻が、zz 軸のまわりに一定の角速度 ω\bm{\omega} でゆっくり回っている(ωR≪c\omega R\ll c)。殻の面密度は σ=M/(4πR2)\sigma=M/(4\pi R^2)、殻の上の点 x′\bm{x}' の速度は v=ω×x′\bm{v}=\bm{\omega}\times\bm{x}' で、物質の分布は時刻によらない。2.18.2節で示すとおり、調和ゲージでは線形化した Einstein 方程式が波動方程式 (2.18.14) になる。ここではそれを使う。

  1. (1) 遅い塵の Tμν=ρuμuνT^{\mu\nu}=\rho u^\mu u^\nu、uμ≈(c,v)u^\mu\approx(c,\bm{v}) から T0iT_{0i} を求め、hˉ0i\bar h_{0i} の静的な方程式を書け。静磁場の ∇2A=−μ0J\nabla^2\bm{A}=-\mu_0\bm{J} と比べ、解を積分で書け。
  2. (2) 一様に帯電した球殻(電荷 QQ、面電荷密度 σQ\sigma_Q)が角速度 ω\bm{\omega} で回るときのベクトルポテンシャルは、殻の外で A=μ04πm×xr3\bm{A}=\frac{\mu_0}{4\pi}\frac{\bm{m}\times\bm{x}}{r^3}(m=13QR2ω\bm{m}=\frac13QR^2\bm{\omega})、殻の内で A=μ0σQR3ω×x\bm{A}=\frac{\mu_0\sigma_QR}{3}\bm{\omega}\times\bm{x} である。これを読み替えて hˉ0i\bar h_{0i} を求め、殻の外では殻の角運動量 J=23MR2ω\bm{J}=\frac23MR^2\bm{\omega} だけで書けることを示せ。
  3. (3) 殻の内では h0ih_{0i} が ω×x\bm{\omega}\times\bm{x} に比例する。第2部第12章の例題 2.12.7 の回転座標の計量 g0i=−(Ω×x)i/cg_{0i}=-(\bm{\Omega}\times\bm{x})_i/c と比べ、殻の内の局所慣性系が遠方の座標に対してどちら向きに、どれだけの角速度で回るかを求めよ。地球の質量と半径(GM=3.986×1014 m3 s−2GM=3.986\times10^{14}\,\mathrm{m^3\,s^{-2}}、R=6.371×106 mR=6.371\times10^6\,\mathrm m)を持ち、自転の角速度 ω=7.292×10−5 s−1\omega=7.292\times10^{-5}\,\mathrm{s^{-1}} で回る球殻では、その角速度は1年あたり何秒角か。
解答

(1) T0i=ρu0ui=ρc viT^{0i}=\rho u^0u^i=\rho c\,v^i で、添字を下げると T0i=η00ηijT0j=−ρc viT_{0i}=\eta_{00}\eta_{ij}T^{0j}=-\rho c\,v^i である。静的なら □=∂02−∇2→−∇2\Box=\partial_0^2-\nabla^2\to-\nabla^2 なので、式 (2.18.14) の (0i)(0i) 成分は

∇2hˉ0i=16πGc4T0i=−16πGc3ρ vi\nabla^2\bar h_{0i}=\frac{16\pi G}{c^4}T_{0i}=-\frac{16\pi G}{c^3}\rho\,v^i

になる。静磁場の ∇2Ai=−μ0Ji\nabla^2A^i=-\mu_0J^i と同じ形で、hˉ0i↔Ai\bar h_{0i}\leftrightarrow A^i、μ0Ji↔(16πG/c3)ρvi\mu_0J^i\leftrightarrow(16\pi G/c^3)\rho v^i と読み替えればよい。静磁場の解 A=μ04π∫J(x′)∣x−x′∣d3x′\bm{A}=\frac{\mu_0}{4\pi}\int\frac{\bm{J}(\bm{x}')}{\lvert \bm{x}-\bm{x}'\rvert }\mathrm{d}^3x' から

hˉ0i(x)=4Gc3∫ρ(x′) vi(x′)∣x−x′∣ d3x′\bar h_{0i}(\bm{x})=\frac{4G}{c^3}\int\frac{\rho(\bm{x}')\,v^i(\bm{x}')}{\lvert \bm{x}-\bm{x}'\rvert }\,\mathrm{d}^3x'

である。η0i=0\eta_{0i}=0 なので h0i=hˉ0ih_{0i}=\bar h_{0i} で、これが計量の g0ig_{0i} になる。

(2) 読み替えは μ0σQ→16πGσ/c3\mu_0\sigma_Q\to16\pi G\sigma/c^3 である。殻の外では μ04πm=μ0σQR43ω\frac{\mu_0}{4\pi}\bm{m}=\frac{\mu_0\sigma_QR^4}{3}\bm{\omega} を 16πGσR43c3ω=4GMR23c3ω\frac{16\pi G\sigma R^4}{3c^3}\bm{\omega}=\frac{4GMR^2}{3c^3}\bm{\omega} に替えて

hˉ0i=4GMR23c3(ω×x)ir3=2Gc3(J×x)ir3\bar h_{0i}=\frac{4GMR^2}{3c^3}\frac{(\bm{\omega}\times\bm{x})_i}{r^3}=\frac{2G}{c^3}\frac{(\bm{J}\times\bm{x})_i}{r^3}

になる。磁気双極子のベクトルポテンシャル μ04πm×xr3\frac{\mu_0}{4\pi}\frac{\bm{m}\times\bm{x}}{r^3} と同じ形で、磁気モーメント m\bm{m} の役を角運動量 J\bm{J} がしている(μ0m/(4π)↔2GJ/c3\mu_0\bm{m}/(4\pi)\leftrightarrow2G\bm{J}/c^3)。殻の外の場は、殻の大きさによらず角運動量だけで決まる。殻の内では

hˉ0i=16πGσR3c3(ω×x)i=4GM3Rc3(ω×x)i\bar h_{0i}=\frac{16\pi G\sigma R}{3c^3}(\bm{\omega}\times\bm{x})_i=\frac{4GM}{3Rc^3}(\bm{\omega}\times\bm{x})_i

で、r=Rr=R で外の式とつながる。磁場に当たる ∇×\nabla\times は殻の内で一様になる。v/cv/c の1次で現れる計量のこの部分を重力磁気(gravitomagnetism)と呼ぶ。

(3) 殻の内では Φ\Phi が一定なので、h00h_{00} と hijh_{ij} は定数で、時間と長さの目盛りを定数倍するだけである。残る g0i=h0i=4GM3Rc2(ω×x)icg_{0i}=h_{0i}=\frac{4GM}{3Rc^2}\frac{(\bm{\omega}\times\bm{x})_i}{c} を g0i=−(Ω×x)i/cg_{0i}=-(\bm{\Omega}\times\bm{x})_i/c と比べると

Ω=−4GM3Rc2 ω\bm{\Omega}=-\frac{4GM}{3Rc^2}\,\bm{\omega}

である。遠方で慣性系に一致する座標が、殻の内では局所慣性系に対して −4GM3Rc2ω-\frac{4GM}{3Rc^2}\bm{\omega} で回る座標に見える。言い換えると、殻の内の局所慣性系は遠方の座標に対して殻と同じ向きに角速度

Ωd=4GM3Rc2 ω\Omega_{\mathrm d}=\frac{4GM}{3Rc^2}\,\omega

で回る。殻の中の自由なこまの軸は、遠方の星に対してこの角速度でゆっくり向きを変える。回る殻が慣性系を引きずるこの効果は Thirring が線形化した方程式から求めた1H. Thirring (1918), Physikalische Zeitschrift 19, 33。下の脚注へ。数値は GM/(Rc2)=6.96×10−10GM/(Rc^2)=6.96\times10^{-10} から Ωd=43×6.96×10−10×7.292×10−5=6.77×10−14 s−1\Omega_{\mathrm d}=\frac43\times6.96\times10^{-10}\times7.292\times10^{-5}=6.77\times10^{-14}\,\mathrm{s^{-1}} で、1年(3.156×107 s3.156\times10^7\,\mathrm s)で 2.14×10−6 rad=0.44′′2.14\times10^{-6}\,\mathrm{rad}=0.44'' である。

近似の裏付け

型 B(極限) | 小さいもの ∣hμν∣→0\lvert h_{\mu\nu}\rvert \to0(平坦な計量からのずれが小さい。地上の重力波で 10−2110^{-21}) | 証明 要確認(線形化した解が非線形の解の近似になる条件を原論文で確かめていない) | 行き先 近似の数理 第1章

1 H. Thirring (1918), Physikalische Zeitschrift 19, 33。

木下の勉強ステータス

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

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

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