tecktom.com

この章のまとめ

2026年10月15日12分で読めますAI生成教科書『虎の巻』第2部 第15章 Einstein–Hilbert 作用と Einstein 方程式

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第15章 §2.15.7 ★★★ 作用を選ぶ → 計量で変分 → Einstein 方程式、座標の取り替えの不変性 → 保存則、Newton 極限 → 係数、を1周させる

この章では、重力場の作用に Einstein–Hilbert 作用 −c316πG∫(R+2Λ)−g d4x-\frac{c^3}{16\pi G}\int(R+2\Lambda)\sqrt{-g}\,\mathrm{d}^4x を選び、物質の作用と合わせて逆計量 gμνg^{\mu\nu} で変分した。−g\sqrt{-g} の変分から −12gμνR-\frac12g_{\mu\nu}R が出て RμνR_{\mu\nu} と組み合わさり Einstein テンソルになること、δRμν\delta R_{\mu\nu} の寄与は全発散で境界の項になること、物質の作用の計量微分 TμνT_{\mu\nu} が自動的に対称でゲージ不変なエネルギー運動量テンソルになることを確かめた。作用が座標の取り替えで変わらないことの Noether 恒等式から、重力の側では恒等式 ∇μGμν=0\nabla^\mu G_{\mu\nu}=0、物質の側では保存則 ∇μTμν=0\nabla_\mu T^{\mu\nu}=0 が出て、Einstein 方程式の両辺の発散がそろう。係数 8πG/c48\pi G/c^4 と作用の前の符号は、弱い場で Poisson 方程式に戻ること(重力が引力であること)から決まり、宇宙項 Λ>0\Lambda>0 は距離に比例する斥力、あるいは負の圧力を持つ真空のエネルギーとして読むことができる。図 2.15.2 にこの流れをまとめた。

次の章(第2部第16章)では、真空の方程式 Rμν=0R_{\mu\nu}=0 を球対称で解いて Schwarzschild 解を出し、太陽系での検証に進む。重力波(第2部第18章)は、この章の 2.15.6節と同じく作用を計量のゆらぎの2次まで展開して、今度は時間に依存する解を調べる。

図

図 2.15.2 この章の流れ。作用の変分で Einstein 方程式を出し(上の段)、座標の取り替えの不変性で両辺の発散が 00 になることを確かめ(下の段)、Newton 極限で係数と符号を決める(破線)

この章の道具

  • 変分の公式:δgμν=−gμαgνβδgαβ\delta g_{\mu\nu}=-g_{\mu\alpha}g_{\nu\beta}\delta g^{\alpha\beta}、δ−g=−12−g gμνδgμν\delta\sqrt{-g}=-\frac12\sqrt{-g}\,g_{\mu\nu}\delta g^{\mu\nu}、gμνδRμν=∇σWσg^{\mu\nu}\delta R_{\mu\nu}=\nabla_\sigma W^\sigma、−g ∇σVσ=∂σ(−g Vσ)\sqrt{-g}\,\nabla_\sigma V^\sigma=\partial_\sigma(\sqrt{-g}\,V^\sigma)。
  • Einstein 方程式:Gμν−Λgμν=8πGc4TμνG_{\mu\nu}-\Lambda g_{\mu\nu}=\frac{8\pi G}{c^4}T_{\mu\nu}、同値な形 Rμν=κ(Tμν−12gμνT)−ΛgμνR_{\mu\nu}=\kappa(T_{\mu\nu}-\frac12g_{\mu\nu}T)-\Lambda g_{\mu\nu}。
  • 座標の取り替えの Noether 恒等式:∇μGμν=0\nabla^\mu G_{\mu\nu}=0、∇μTμν=0\nabla_\mu T^{\mu\nu}=0。
  • 弱い場:R00=∇2Φ/c2R_{00}=\nabla^2\Phi/c^2 から κ=8πG/c4\kappa=8\pi G/c^4。Λ>0\Lambda>0 は加速度 +Λc2r/3+\Lambda c^2r/3 の斥力、de Sitter で Λ=3H2/c2\Lambda=3H^2/c^2、R=−4ΛR=-4\Lambda。

総合問題 物質と宇宙項を含む平坦な宇宙 ★★★

オリジナル空間が平坦で一様・等方な宇宙 ds2=c2dt2−a(t)2(dx2+dy2+dz2)\mathrm{d} s^2=c^2\mathrm{d} t^2-a(t)^2(\mathrm{d} x^2+\mathrm{d} y^2+\mathrm{d} z^2)(x0=ctx^0=ct)に、宇宙とともに動く塵(質量密度 ρ(t)\rho(t)、圧力 00)と宇宙項 Λ>0\Lambda>0 がある。H=a˙/aH=\dot a/a とする。

  1. (1) (使う道具:第2部第14章の曲率の計算)この計量で R=−6(aa¨+a˙2)c2a2R=-\dfrac{6(a\ddot a+\dot a^2)}{c^2a^2}、G00=3a˙2c2a2G_{00}=\dfrac{3\dot a^2}{c^2a^2}、Gij=−2aa¨+a˙2c2δijG_{ij}=-\dfrac{2a\ddot a+\dot a^2}{c^2}\delta_{ij} となることを、Christoffel 記号から確かめよ。
  2. (2) (使う道具:2.15.3節)Einstein 方程式 (2.15.11) の 0000 成分と ijij 成分から、H2=8πG3ρ+Λc23H^2=\frac{8\pi G}{3}\rho+\frac{\Lambda c^2}{3} と 2a¨a+H2=Λc2\frac{2\ddot a}{a}+H^2=\Lambda c^2 を出せ。
  3. (3) (使う道具:2.15.4節・例題 2.15.4)∇μTμν=0\nabla_\mu T^{\mu\nu}=0 から ρa3\rho a^3 が一定であることを示し、(2) の2式がこれと矛盾しないこと(第1式を時間微分すると第2式が出ること)を確かめよ。
  4. (4) (使う道具:2.15.5節)加速の式 a¨a=−4πG3ρ+Λc23\frac{\ddot a}{a}=-\frac{4\pi G}{3}\rho+\frac{\Lambda c^2}{3} を示せ。今の値 Ωm=8πGρ03H02=0.315\Omega_{\mathrm m}=\frac{8\pi G\rho_0}{3H_0^2}=0.315、ΩΛ=Λc23H02=0.685\Omega_\Lambda=\frac{\Lambda c^2}{3H_0^2}=0.685 を使って、減速パラメータ q0=−a¨aa˙2∣0q_0=-\frac{\ddot a a}{\dot a^2}\big|_0 と、膨張が加速に転じた時の赤方偏移 zz(1+z=a0/a1+z=a_0/a)を求めよ。
  5. (5) (使う道具:2.15.1節・2.15.5節)a(t)=(ΩmΩΛ)1/3sinh⁡2/3(32ΩΛH0t)a(t)=\Bigl(\frac{\Omega_{\mathrm m}}{\Omega_\Lambda}\Bigr)^{1/3}\sinh^{2/3}\Bigl(\frac32\sqrt{\Omega_\Lambda}H_0t\Bigr)(a0=1a_0=1)が (2) の第1式を満たすことを示せ。宇宙の年齢 t0t_0(a=1a=1 になる時刻)を H0=67.4 km s−1 Mpc−1H_0=67.4\,\mathrm{km\,s^{-1}\,Mpc^{-1}} で求め、遠い未来の膨張率 H∞H_\infty とスカラー曲率 R∞R_\infty を求めよ(放射は無視する)。
解答

(1) x0=ctx^0=ct で g00=1g_{00}=1、gij=−a2δijg_{ij}=-a^2\delta_{ij}。0でない Christoffel 記号は Γ0ij=aa˙cδij\Gamma^0{}_{ij}=\frac{a\dot a}{c}\delta_{ij}、Γi0j=a˙caδij\Gamma^i{}_{0j}=\frac{\dot a}{ca}\delta^i{}_j である(Γ0ij=−12g00∂0gij=12c∂t(a2)δij\Gamma^0{}_{ij}=-\frac12g^{00}\partial_0g_{ij}=\frac{1}{2c}\partial_t(a^2)\delta_{ij}、Γi0j=12gik∂0gkj=a˙caδij\Gamma^i{}_{0j}=\frac12g^{ik}\partial_0g_{kj}=\frac{\dot a}{ca}\delta^i{}_j)。これらを Ricci テンソルの定義に入れると R00=−3a¨c2aR_{00}=-\frac{3\ddot a}{c^2a}、Rij=aa¨+2a˙2c2δijR_{ij}=\frac{a\ddot a+2\dot a^2}{c^2}\delta_{ij} で、R=R00−1a2∑iRii=−3a¨c2a−3(aa¨+2a˙2)c2a2=−6(aa¨+a˙2)c2a2R=R_{00}-\frac{1}{a^2}\sum_iR_{ii}=-\frac{3\ddot a}{c^2a}-\frac{3(a\ddot a+2\dot a^2)}{c^2a^2}=-\frac{6(a\ddot a+\dot a^2)}{c^2a^2}。G00=R00−12R=−3a¨c2a+3(aa¨+a˙2)c2a2=3a˙2c2a2G_{00}=R_{00}-\frac12R=-\frac{3\ddot a}{c^2a}+\frac{3(a\ddot a+\dot a^2)}{c^2a^2}=\frac{3\dot a^2}{c^2a^2}、Gij=Rij+a22R δij=aa¨+2a˙2−3aa¨−3a˙2c2δij=−2aa¨+a˙2c2δijG_{ij}=R_{ij}+\frac{a^2}{2}R\,\delta_{ij}=\frac{a\ddot a+2\dot a^2-3a\ddot a-3\dot a^2}{c^2}\delta_{ij}=-\frac{2a\ddot a+\dot a^2}{c^2}\delta_{ij}。Λ>0\Lambda>0 の de Sitter(a=eHta=\mathrm{e}^{Ht})では R=−12H2/c2R=-12H^2/c^2 で、正定値の計量の球面の R>0R>0 とは逆の符号になる。(+,−,−,−) で空間の部分を負に取っているためである(記法の正本 §4)。

(2) 塵の T00=ρc2T_{00}=\rho c^2、Tij=0T_{ij}=0。0000 成分:3a˙2c2a2−Λ=8πGc4ρc2\frac{3\dot a^2}{c^2a^2}-\Lambda=\frac{8\pi G}{c^4}\rho c^2 から H2=8πG3ρ+Λc23H^2=\frac{8\pi G}{3}\rho+\frac{\Lambda c^2}{3}。ijij 成分:−2aa¨+a˙2c2δij−Λ(−a2δij)=0-\frac{2a\ddot a+\dot a^2}{c^2}\delta_{ij}-\Lambda(-a^2\delta_{ij})=0 から 2a¨a+H2=Λc2\frac{2\ddot a}{a}+H^2=\Lambda c^2。第1式は例題 2.15.3 で NN の方程式として出したものと同じである。

(3) 例題 2.15.4 で P=0P=0 とすると ρ˙+3Hρ=0\dot\rho+3H\rho=0、つまり ddt(ρa3)=a3(ρ˙+3Hρ)=0\frac{\mathrm{d} }{\mathrm{d} t}(\rho a^3)=a^3(\dot\rho+3H\rho)=0。第1式を時間微分すると 2HH˙=8πG3ρ˙=−8πGHρ2H\dot H=\frac{8\pi G}{3}\dot\rho=-8\pi GH\rho、H˙=a¨a−H2\dot H=\frac{\ddot a}{a}-H^2 なので a¨a=H2−4πGρ\frac{\ddot a}{a}=H^2-4\pi G\rho。第1式で ρ\rho を消すと 4πGρ=32(H2−Λc23)4\pi G\rho=\frac32(H^2-\frac{\Lambda c^2}{3}) だから a¨a=−12H2+Λc22\frac{\ddot a}{a}=-\frac12H^2+\frac{\Lambda c^2}{2} で、第2式と一致する。ijij 成分の方程式が 0000 成分と保存則から出ることは、2.15.4節の恒等式 ∇μGμν=0\nabla^\mu G_{\mu\nu}=0 の表れである。

(4) 第2式と第1式から a¨a=12(Λc2−H2)=12(Λc2−8πG3ρ−Λc23)=−4πG3ρ+Λc23\frac{\ddot a}{a}=\frac12(\Lambda c^2-H^2)=\frac12\Bigl(\Lambda c^2-\frac{8\pi G}{3}\rho-\frac{\Lambda c^2}{3}\Bigr)=-\frac{4\pi G}{3}\rho+\frac{\Lambda c^2}{3}。第1項は Newton の引力で、第2項は 2.15.5節の距離に比例する斥力 Λc23r\frac{\Lambda c^2}{3}r そのものである。今の値を H02H_0^2 で割ると q0=−a¨aH02=Ωm2−ΩΛ=0.1575−0.685=−0.528q_0=-\frac{\ddot a}{aH_0^2}=\frac{\Omega_{\mathrm m}}{2}-\Omega_\Lambda=0.1575-0.685=-0.528 で、今の膨張は加速している。ρ∝a−3=(1+z)3\rho\propto a^{-3}=(1+z)^3 なので、a¨=0\ddot a=0 は Ωm2(1+z)3=ΩΛ\frac{\Omega_{\mathrm m}}{2}(1+z)^3=\Omega_\Lambda、z=(2×0.6850.315)1/3−1=0.632z=\Bigl(\frac{2\times0.685}{0.315}\Bigr)^{1/3}-1=0.632。これより前(z>0.632z>0.632)は減速していた。

(5) τ≔32ΩΛH0t\tau\coloneqq\frac32\sqrt{\Omega_\Lambda}H_0t、A≔(Ωm/ΩΛ)1/3A\coloneqq(\Omega_{\mathrm m}/\Omega_\Lambda)^{1/3} とおくと a=Asinh⁡2/3τa=A\sinh^{2/3}\tau、H=a˙a=23coth⁡τ⋅32ΩΛH0=ΩΛH0coth⁡τH=\frac{\dot a}{a}=\frac23\coth\tau\cdot\frac32\sqrt{\Omega_\Lambda}H_0=\sqrt{\Omega_\Lambda}H_0\coth\tau。右辺は H02(Ωma−3+ΩΛ)=H02ΩΛ(sinh⁡−2τ+1)=H02ΩΛcoth⁡2τH_0^2\bigl(\Omega_{\mathrm m}a^{-3}+\Omega_\Lambda\bigr)=H_0^2\Omega_\Lambda\bigl(\sinh^{-2}\tau+1\bigr)=H_0^2\Omega_\Lambda\coth^2\tau で、左辺 H2H^2 と一致する(ρ=ρ0a−3\rho=\rho_0a^{-3} を使った)。a=1a=1 は sinh⁡τ0=ΩΛ/Ωm\sinh\tau_0=\sqrt{\Omega_\Lambda/\Omega_{\mathrm m}} のときで、H0t0=23ΩΛarsinh⁡ΩΛ/Ωm=0.951H_0t_0=\frac{2}{3\sqrt{\Omega_\Lambda}}\operatorname{arsinh}\sqrt{\Omega_\Lambda/\Omega_{\mathrm m}}=0.951。1/H0=14.5 Gyr1/H_0=14.5\,\mathrm{Gyr} なので t0=13.8 Gyrt_0=13.8\,\mathrm{Gyr} である。t→∞t\to\infty で coth⁡τ→1\coth\tau\to1、H∞=ΩΛH0=cΛ/3=55.8 km s−1 Mpc−1H_\infty=\sqrt{\Omega_\Lambda}H_0=c\sqrt{\Lambda/3}=55.8\,\mathrm{km\,s^{-1}\,Mpc^{-1}} で、a∝eH∞ta\propto\mathrm{e}^{H_\infty t} の de Sitter 宇宙に近づく。そのとき物質は薄まり、式 (2.15.12) のトレースから R∞=−4Λ=−4.36×10−52 m−2R_\infty=-4\Lambda=-4.36\times10^{-52}\,\mathrm{m^{-2}} になる。図 2.15.3 に a(t)a(t) を物質だけの宇宙と比べて示した。

図

図 2.15.3 同じ H0H_0 を持つ2つの平坦な宇宙のスケール因子。実線は Ωm=0.315\Omega_{\mathrm m}=0.315・ΩΛ=0.685\Omega_\Lambda=0.685(総合問題 (5) の厳密解)、破線は物質だけ(a=(3H0t/2)2/3a=(3H_0t/2)^{2/3}、今は H0t=2/3H_0t=2/3)。宇宙項のある宇宙は z=0.632z=0.632 を境に加速に転じ、遠い未来に指数関数で膨張する

木下の勉強ステータス

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

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

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