解析力学から様々な理論分野を組み立てる 第2部 第11章 §2.11.1 ★★★ 慣性質量と重力質量が等しいので、重力は m m m に依らない計量 g μ ν g_{\mu\nu} g μν に入れる。弱い場では g 00 = 1 + 2 Φ / c 2 g_{00}=1+2\Phi/c^2 g 00 = 1 + 2Φ/ c 2 で Newton の − m Φ -m\Phi − m Φ に戻る
第2部第3章では、自由粒子の作用を世界線の固有時に − m c 2 -mc^2 − m c 2 を掛けたもの S = − m c ∫ d s S=-mc\int\mathrm{d} s S = − m c ∫ d s と選んだ。そこでは世界間隔をd s 2 = η μ ν d x μ d x ν \mathrm{d} s^2=\eta_{\mu\nu}\mathrm{d} x^\mu\mathrm{d} x^\nu d s 2 = η μν d x μ d x ν と平坦な計量で測った。この章では、重力の中の粒子の作用を
S = − m c ∫ g μ ν ( x ) d x μ d x ν S=-mc\int\sqrt{g_{\mu\nu}(x)\,\mathrm{d} x^\mu\mathrm{d} x^\nu} S = − m c ∫ g μν ( x ) d x μ d x ν
と選ぶ。違いは、定数の行列 η μ ν \eta_{\mu\nu} η μν を場所によって変わる行列 g μ ν ( x ) g_{\mu\nu}(x) g μν ( x ) に替えたことだけである。この1行を選ぶ理由は、重力の中ではどんな物体も同じように落ちる、という実験事実にある。Newton 力学の Lagrangian
L = 1 2 m v 2 − m Φ L=\frac12mv^2-m\Phi L = 2 1 m v 2 − m Φ では、慣性を表す m m m と重力を受ける強さを表す m m m がたまたま同じ字で、運動方程式から m m m が消える。消える量で軌道を区別できないなら、重力は物体の性質でなく時空の性質として書くのが自然である。g μ ν g_{\mu\nu} g μν の中に重力を入れれば、m m m は作用全体に掛かる係数になり、軌道に入らないことが最初から保証される。
この作用に第1部の道具を当てると、次のものが順に出てくる。座標の取り替えと「Lagrangian は全微分だけの違いを許す」という性質から、落下するエレベーターの中では一様な重力が消えること(等価原理)と、潮汐力だけは消せないことが出る。第2部第3章の双曲運動を観測者の族に広げると、一様に加速する観測者の座標(Rindler 座標)で計量が g 00 = ( 1 + a ξ / c 2 ) 2 g_{00}=(1+a\xi/c^2)^2 g 00 = ( 1 + a ξ / c 2 ) 2 になり、平坦な時空の中にも「重力」が座標の取り方として現れる。時刻をずらす変換に対する Noether の保存量からは、高い所の時計が速く進むこと(重力赤方偏移)が出る。最後に、この効果を GPS の衛星の時計と、1 cm の高さの違いを読み取る光格子時計の数値で確かめる。
この章では計量の弱い場の形 g 00 ≈ 1 + 2 Φ / c 2 g_{00}\approx1+2\Phi/c^2 g 00 ≈ 1 + 2Φ/ c 2 を、等価原理から決まる部分だけ使う。空間の成分 g i j g_{ij} g ij の補正は Einstein 方程式から決まり、第2部第15・16章で扱う。測地線方程式と Christoffel 記号は第2部第12章で、この章の作用から Euler–Lagrange 方程式として出す。計量は η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) (この本の流儀)とし、上下に1回ずつ現れるギリシャ添字だけ和をとる。この章では Newton のポテンシャルを Φ \Phi Φ (Φ < 0 \Phi<0 Φ < 0 )、地表の重力加速度を g g g 、ロケットの固有加速度を a a a と書く。g g g は添字の付いた計量 g μ ν g_{\mu\nu} g μν とは添字の有無で区別する。
慣性質量と重力質量
Newton 力学で重力の中の粒子の Lagrangian を書くと、
L = 1 2 m i v 2 − m g Φ ( x ) (2.11.1) L=\frac12m_{\mathrm i}v^2-m_{\mathrm g}\Phi(\bm{x})\tag{2.11.1} L = 2 1 m i v 2 − m g Φ ( x ) ( 2.11.1 )
となる。m i m_{\mathrm i} m i は運動エネルギーに付く質量(慣性質量)、m g m_{\mathrm g} m g は重力のポテンシャルエネルギーに付く質量(重力質量)で、定義の上では別の量である。電磁気で L = 1 2 m v 2 − q ϕ L=\frac12mv^2-q\phi L = 2 1 m v 2 − q ϕ と書くとき、m m m と q q q が別の量であるのと同じ立場である。Euler–Lagrange 方程式は
x ¨ = − m g m i ∇ Φ (2.11.2) \ddot{\bm{x}}=-\frac{m_{\mathrm g}}{m_{\mathrm i}}\nabla\Phi\tag{2.11.2} x ¨ = − m i m g ∇Φ ( 2.11.2 )
で、加速度は比 m g / m i m_{\mathrm g}/m_{\mathrm i} m g / m i によって物体ごとに違ってよいはずである。電場の中の粒子の加速度が q / m q/m q / m で物体ごとに違うのと同じである。
ところが実験では、この比は物質によらない。2種類の物体の加速度の差を平均で割った量を η \eta η と書くと、人工衛星の中で自由落下する白金合金とチタン合金の試験体を比べた実験は ∣ η ∣ \lvert \eta\rvert ∣ η ∣ が 10 − 14 10^{-14} 1 0 − 14 より小さいことを示した1 1 P. Touboul et al. (2022), Phys. Rev. Lett. 129 , 121102(MICROSCOPE 衛星の最終結果。η = [ − 1.5 ± 2.3 ( 統計 ) ± 1.5 ( 系統 ) ] × 10 − 15 \eta=[-1.5\pm2.3(\text{統計})\pm1.5(\text{系統})]\times10^{-15} η = [ − 1.5 ± 2.3 ( 統計 ) ± 1.5 ( 系統 )] × 1 0 − 15 )。 。比が物体によらないなら単位の取り方で m g = m i ≡ m m_{\mathrm g}=m_{\mathrm i}\equiv m m g = m i ≡ m とおけて、式 (2.11.2) から m m m が消える。軌道は初めの位置と速度だけで決まり、どんな物体でも同じになる。これを自由落下の普遍性 と呼ぶ。
選ぶ作用
軌道が物体によらないなら、重力は物体の側の量(電荷のような結合定数)を通してでなく、粒子が動く舞台の側に書くのが自然である。第2部第3章の作用の中で、舞台を表しているのは世界間隔を測る計量 η μ ν \eta_{\mu\nu} η μν だけなので、そこに重力を入れる。
定義 2.11.1(重力の中の粒子の作用)
質量 m > 0 m>0 m > 0 の粒子の作用を、場所に依存する計量 g μ ν ( x ) g_{\mu\nu}(x) g μν ( x ) (対称、符号は ( + , − , − , − ) (+,-,-,-) ( + , − , − , − ) )を使って
S [ x ] = − m c ∫ d s , d s 2 = g μ ν ( x ) d x μ d x ν (2.11.3) S[x]=-mc\int\mathrm{d} s,\qquad \mathrm{d} s^2=g_{\mu\nu}(x)\,\mathrm{d} x^\mu\mathrm{d} x^\nu\tag{2.11.3} S [ x ] = − m c ∫ d s , d s 2 = g μν ( x ) d x μ d x ν ( 2.11.3 )
と選ぶ。時刻 t t t を径数にとると x 0 = c t x^0=ct x 0 = c t で、Lagrangian は
L = − m c g μ ν ( x ) x ˙ μ x ˙ ν , x ˙ μ = d x μ d t = ( c , v ) (2.11.4) L=-mc\sqrt{g_{\mu\nu}(x)\,\dot x^\mu\dot x^\nu},\qquad \dot x^\mu=\frac{\mathrm{d} x^\mu}{\mathrm{d} t}=(c,\bm{v})\tag{2.11.4} L = − m c g μν ( x ) x ˙ μ x ˙ ν , x ˙ μ = d t d x μ = ( c , v ) ( 2.11.4 )
となる。重力は g μ ν g_{\mu\nu} g μν だけに入り、粒子の性質 m m m は作用全体に掛かる係数である。
m m m が作用全体の係数なので、Euler–Lagrange 方程式の両辺を m m m で割ることができ、軌道は m m m によらない。自由落下の普遍性は、この作用を選んだ時点で式の形として保証される。電磁場との結合 − q ∫ A μ d x μ -q\int A_\mu\mathrm{d} x^\mu − q ∫ A μ d x μ (第2部第3章)は m m m と別の係数 q q q を持つので、q / m q/m q / m の違う粒子は違う軌道を描く。この違いが、重力だけを計量に入れる理由である。
弱い場で Newton に戻ること
時間によらない弱い重力の中で、計量を
g 00 = 1 + 2 Φ ( x ) c 2 , g 0 i = 0 , g i j = − δ i j (2.11.5) g_{00}=1+\frac{2\Phi(\bm{x})}{c^2},\qquad g_{0i}=0,\qquad g_{ij}=-\delta_{ij}\tag{2.11.5} g 00 = 1 + c 2 2Φ ( x ) , g 0 i = 0 , g ij = − δ ij ( 2.11.5 )
とおいてみる(∣ Φ ∣ / c 2 ≪ 1 \lvert \Phi\rvert /c^2\ll1 ∣ Φ ∣ / c 2 ≪ 1 )。g i j g_{ij} g ij にも Φ / c 2 \Phi/c^2 Φ/ c 2 の大きさの補正が付くが、遅い粒子の運動には g 00 g_{00} g 00 の補正だけが現れることをすぐ下で確かめる。式 (2.11.4) に入れると L = − m c ( 1 + 2 Φ / c 2 ) c 2 − v 2 = − m c 2 1 + ( 2 Φ − v 2 ) / c 2 L=-mc\sqrt{(1+2\Phi/c^2)c^2-v^2}=-mc^2\sqrt{1+(2\Phi-v^2)/c^2} L = − m c ( 1 + 2Φ/ c 2 ) c 2 − v 2 = − m c 2 1 + ( 2Φ − v 2 ) / c 2 で、1 + x = 1 + x / 2 − x 2 / 8 + ⋯ \sqrt{1+x}=1+x/2-x^2/8+\cdots 1 + x = 1 + x /2 − x 2 /8 + ⋯ を使って展開すると
L = − m c 2 + 1 2 m v 2 − m Φ + m ( 2 Φ − v 2 ) 2 8 c 2 + O ( c − 4 ) (2.11.6) L=-mc^2+\frac12mv^2-m\Phi+\frac{m(2\Phi-v^2)^2}{8c^2}+\mathcal{O}(c^{-4})\tag{2.11.6} L = − m c 2 + 2 1 m v 2 − m Φ + 8 c 2 m ( 2Φ − v 2 ) 2 + O ( c − 4 ) ( 2.11.6 )
を得る。定数 − m c 2 -mc^2 − m c 2 は運動に関係しないので、c → ∞ c\to\infty c → ∞ で Newton の Lagrangian (2.11.1) (m g = m i = m m_{\mathrm g}=m_{\mathrm i}=m m g = m i = m )に戻る。逆に言えば、Newton 力学に戻るという要請が g 00 g_{00} g 00 の 1 / c 2 1/c^2 1/ c 2 の係数を 2 Φ 2\Phi 2Φ に決めた。g i j g_{ij} g ij を − ( 1 + 2 Ψ / c 2 ) δ i j -(1+2\Psi/c^2)\delta_{ij} − ( 1 + 2Ψ/ c 2 ) δ ij と補正しても、L L L には − m c 2 ⋅ ( − Ψ v 2 / c 4 ) -mc^2\cdot(-\Psi v^2/c^4) − m c 2 ⋅ ( − Ψ v 2 / c 4 ) 、つまり m v 2 Ψ / c 2 mv^2\Psi/c^2 m v 2 Ψ/ c 2 の大きさの項しか加わらず、これは式 (2.11.6) の最後の項と同じく 1 / c 2 1/c^2 1/ c 2 だけ小さい。したがって、Newton 極限で決まるのは g 00 g_{00} g 00 だけである。Ψ \Psi Ψ は光の曲がりのように v ∼ c v\sim c v ∼ c の運動で現れ、Einstein 方程式から Ψ = − Φ \Psi=-\Phi Ψ = − Φ と決まる(第2部第16章)。
近似の裏付け
型 B(極限) | 小さいもの ∣ Φ ∣ / c 2 → 0 \lvert \Phi\rvert /c^2\to0 ∣ Φ ∣ / c 2 → 0 かつ v / c → 0 v/c\to0 v / c → 0 (重力のポテンシャルと速さが光速の2乗・光速に比べて小さい) | 証明 証明済み(有限の時間区間で ∇ Φ \nabla\Phi ∇Φ が Lipschitz 連続なとき、運動方程式の解のパラメータへの連続依存。g i j g_{ij} g ij の補正は Einstein 方程式から第2部第16章で決める) | 行き先 近似の数理 第2章
例題 2.11.1 地球の公転で 1 / c 2 1/c^2 1/ c 2 の項はどれだけ小さいか ★★
オリジナル 太陽のまわりを半径 r r r の円軌道で回る地球を考える。Φ = − G M ⊙ / r \Phi=-GM_\odot/r Φ = − G M ⊙ / r 、v 2 = G M ⊙ / r v^2=GM_\odot/r v 2 = G M ⊙ / r である。
(1) 式 (2.11.6) の最後の項と運動エネルギー 1 2 m v 2 \frac12mv^2 2 1 m v 2 の比を G M ⊙ / ( r c 2 ) GM_\odot/(rc^2) G M ⊙ / ( r c 2 ) で表せ。
(2) G M ⊙ = 1.327 × 10 20 m 3 s − 2 GM_\odot=1.327\times10^{20}\,\mathrm{m^3\,s^{-2}} G M ⊙ = 1.327 × 1 0 20 m 3 s − 2 、r = 1.496 × 10 11 m r=1.496\times10^{11}\,\mathrm m r = 1.496 × 1 0 11 m として数値を求めよ。
解答 (1) 2 Φ − v 2 = − 2 G M ⊙ / r − G M ⊙ / r = − 3 G M ⊙ / r 2\Phi-v^2=-2GM_\odot/r-GM_\odot/r=-3GM_\odot/r 2Φ − v 2 = − 2 G M ⊙ / r − G M ⊙ / r = − 3 G M ⊙ / r なので、最後の項は m ⋅ 9 ( G M ⊙ / r ) 2 / ( 8 c 2 ) m\cdot9(GM_\odot/r)^2/(8c^2) m ⋅ 9 ( G M ⊙ / r ) 2 / ( 8 c 2 ) である。1 2 m v 2 = 1 2 m G M ⊙ / r \frac12mv^2=\frac12mGM_\odot/r 2 1 m v 2 = 2 1 m G M ⊙ / r で割ると
9 ( G M ⊙ / r ) 2 / ( 8 c 2 ) G M ⊙ / ( 2 r ) = 9 4 G M ⊙ r c 2 . \frac{9(GM_\odot/r)^2/(8c^2)}{GM_\odot/(2r)}=\frac94\,\frac{GM_\odot}{rc^2}. G M ⊙ / ( 2 r ) 9 ( G M ⊙ / r ) 2 / ( 8 c 2 ) = 4 9 r c 2 G M ⊙ .
(2) G M ⊙ / ( r c 2 ) = 1.327 × 10 20 / ( 1.496 × 10 11 × 8.988 × 10 16 ) = 9.87 × 10 − 9 GM_\odot/(rc^2)=1.327\times10^{20}/(1.496\times10^{11}\times8.988\times10^{16})=9.87\times10^{-9} G M ⊙ / ( r c 2 ) = 1.327 × 1 0 20 / ( 1.496 × 1 0 11 × 8.988 × 1 0 16 ) = 9.87 × 1 0 − 9 で、比は 2.22 × 10 − 8 2.22\times10^{-8} 2.22 × 1 0 − 8 である。地球の公転では Newton 力学からのずれは 10 − 8 10^{-8} 1 0 − 8 の大きさで、軌道の形を何周にもわたって積み重ねたときに初めて見える。水星の近日点の移動(第2部第16章)は、この大きさの効果が1世紀にわたって積み重なったものである。
計量の流儀:ほかの教科書との対応
一般相対論の教科書は計量の符号の流儀が2つに分かれる。この本は特殊相対論の章と同じ ( + , − , − , − ) (+,-,-,-) ( + , − , − , − ) で、Landau–Lifshitz 第2巻と同じである。Carroll・Wald・Misner–Thorne–Wheeler(MTW)は ( − , + , + , + ) (-,+,+,+) ( − , + , + , + ) を使う。曲率の定義はこの本と MTW で同じ式を使うので、Christoffel 記号・Riemann テンソル・Ricci テンソルは成分まで一致し、符号が変わるのはスカラー曲率 R R R と、それを含む Einstein–Hilbert 作用の前の係数である。これらの量はまだ定義していないが(第2部第13〜15章)、ほかの教科書と読み比べるときの目印として、表 2.11.1 を先に置く。
表 2.11.1 計量の符号の流儀による違い。「同じ」は成分の値まで一致することを表す 量 この本・Landau–Lifshitz ( + , − , − , − ) (+,-,-,-) ( + , − , − , − ) Carroll・Wald・MTW ( − , + , + , + ) (-,+,+,+) ( − , + , + , + ) 平坦な計量 d i a g ( + 1 , − 1 , − 1 , − 1 ) \mathrm{diag}(+1,-1,-1,-1) diag ( + 1 , − 1 , − 1 , − 1 ) d i a g ( − 1 , + 1 , + 1 , + 1 ) \mathrm{diag}(-1,+1,+1,+1) diag ( − 1 , + 1 , + 1 , + 1 ) (成分は全部 − 1 -1 − 1 倍)u μ u μ u_\mu u^\mu u μ u μ + c 2 +c^2 + c 2 − c 2 -c^2 − c 2 粒子の Lagrangian − m c g μ ν x ˙ μ x ˙ ν -mc\sqrt{g_{\mu\nu}\dot x^\mu\dot x^\nu} − m c g μν x ˙ μ x ˙ ν − m c − g μ ν x ˙ μ x ˙ ν -mc\sqrt{-g_{\mu\nu}\dot x^\mu\dot x^\nu} − m c − g μν x ˙ μ x ˙ ν Γ ρ μ ν \Gamma^\rho{}_{\mu\nu} Γ ρ μν 同じ 同じ R ρ σ μ ν R^\rho{}_{\sigma\mu\nu} R ρ σ μν ・R μ ν R_{\mu\nu} R μν 同じ 同じ(Wald は添字の並びが違うが値は同じ流儀) スカラー曲率 R R R 符号が逆 符号が逆(正定値の球面ではどちらも + 2 / a 2 +2/a^2 + 2/ a 2 ) G μ ν G_{\mu\nu} G μν ・Einstein 方程式G μ ν = + 8 π G T μ ν / c 4 G_{\mu\nu}=+8\pi GT_{\mu\nu}/c^4 G μν = + 8 π G T μν / c 4 同じ 宇宙項 G μ ν − Λ g μ ν G_{\mu\nu}-\Lambda g_{\mu\nu} G μν − Λ g μν 、作用は R + 2 Λ R+2\Lambda R + 2Λ G μ ν + Λ g μ ν G_{\mu\nu}+\Lambda g_{\mu\nu} G μν + Λ g μν 、作用は R − 2 Λ R-2\Lambda R − 2Λ Einstein–Hilbert 作用の前 − c 3 / ( 16 π G ) -c^3/(16\pi G) − c 3 / ( 16 π G ) + c 3 / ( 16 π G ) +c^3/(16\pi G) + c 3 / ( 16 π G ) T μ ν T_{\mu\nu} T μν の定義+ ( 2 c / − g ) δ S m / δ g μ ν +(2c/\sqrt{-g})\,\delta S_{\mathrm m}/\delta g^{\mu\nu} + ( 2 c / − g ) δ S m / δ g μν − ( 2 / − g ) δ S m / δ g μ ν -(2/\sqrt{-g})\,\delta S_{\mathrm m}/\delta g^{\mu\nu} − ( 2/ − g ) δ S m / δ g μν
粒子の Lagrangian の行は、この章ですぐ使う。( − , + , + , + ) (-,+,+,+) ( − , + , + , + ) の流儀では時間的な世界線で g μ ν x ˙ μ x ˙ ν < 0 g_{\mu\nu}\dot x^\mu\dot x^\nu<0 g μν x ˙ μ x ˙ ν < 0 なので、根号の中に負号を付けて正にする。どちらの流儀でも作用は「− m c 2 × -mc^2\times − m c 2 × 固有時」で同じ値をとり、物理の結果は変わらない。弱い場の式 (2.11.5) は ( − , + , + , + ) (-,+,+,+) ( − , + , + , + ) では g 00 = − ( 1 + 2 Φ / c 2 ) g_{00}=-(1+2\Phi/c^2) g 00 = − ( 1 + 2Φ/ c 2 ) と書かれる。Weinberg の教科書は Riemann テンソル全体の符号がこの本と逆なので、式を写すときは注意する。
1 P. Touboul et al. (2022), Phys. Rev. Lett. 129 , 121102(MICROSCOPE 衛星の最終結果。η = [ − 1.5 ± 2.3 ( 統計 ) ± 1.5 ( 系統 ) ] × 10 − 15 \eta=[-1.5\pm2.3(\text{統計})\pm1.5(\text{系統})]\times10^{-15} η = [ − 1.5 ± 2.3 ( 統計 ) ± 1.5 ( 系統 )] × 1 0 − 15 )。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com