フリードマン方程式の導出(2)

第2章 FLRW計量に対する重力場の方程式

この記事の目的は、一様・等方な時空に対する重力場の方程式を求めることだった。第1章では一様・等方な時空の計量の一般的な形を求めた。それを重力場の方程式（アインシュタイン方程式）に代入したらどうなるかやってみよう。

重力場の方程式はいつものように混合テンソルによる方程式 $𝐺^{𝜇}_{𝜈} + 𝛬 𝛿^{𝜇}_{𝜈} = \frac{8 𝜋 𝐺}{𝑐^{4}} 𝑇^{𝜇}_{𝜈} (27)$ を使うことにする。 𝐺^𝜇_𝜈 はアインシュタインテンソル、 𝑇^𝜇_𝜈 はエネルギー運動量テンソル、 𝛿^𝜇_𝜈 はクロネッカーのデルタ、 𝛬 は宇宙定数、 𝐺 は万有引力定数、 𝑐 は光速である。

クリストッフェル記号

クリストッフェル記号の0でない成分はすでに第1章の最後の方で(25)式 $\begin{array}{l} 𝛤^{0}_{1}_{1} = \frac{𝑎 \dot{𝑎}}{1 - 𝑘 𝑟^{2}} \\ 𝛤^{0}_{2}_{2} = 𝑎 \dot{𝑎} 𝑟^{2} \\ 𝛤^{0}_{3}_{3} = 𝑎 \dot{𝑎} 𝑟^{2} \sin^{2} 𝜃 \end{array}$ $\begin{array}{l} 𝛤^{1}_{1}_{1} = \frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} \\ 𝛤^{1}_{2}_{2} = - (1 - 𝑘 𝑟^{2}) 𝑟 \\ 𝛤^{1}_{3}_{3} = - (1 - 𝑘 𝑟^{2}) 𝑟 \sin^{2} 𝜃 \\ 𝛤^{2}_{3}_{3} = - \sin 𝜃 \cos 𝜃 \end{array}$ $\begin{array}{l} 𝛤^{1}_{0}_{1} = 𝛤^{1}_{1}_{0} = 𝛤^{2}_{0}_{2} = 𝛤^{2}_{2}_{0} = 𝛤^{3}_{0}_{3} = 𝛤^{3}_{3}_{0} = \frac{\dot{𝑎}}{𝑎} \\ 𝛤^{2}_{1}_{2} = 𝛤^{2}_{2}_{1} = 𝛤^{3}_{1}_{3} = 𝛤^{3}_{3}_{1} = \frac{1}{𝑟} \\ 𝛤^{3}_{2}_{3} = 𝛤^{3}_{3}_{2} = \cot 𝜃 \\ ↑ 10個まとめて(25)式 \end{array}$ で求めてあるので、これを使えばよい。

この後でリッチテンソルを求める際にクリストッフェル記号の微分が必要になるので計算しておく。ただし $\frac{{∂𝛤}^{𝜆}_{𝜇}_{𝜈}}{{∂𝑥}^{𝜌}}$ のすべての成分がいるわけではなく、 𝜆 = 𝜌 または 𝜆 = 𝜈 である成分が求まれば十分であるから、それらだけを計算する。(25)式をただ微分するだけなので、結果だけを書くと、0でない成分は以下である。 $\begin{array}{l} \frac{{∂𝛤}^{0}_{1}_{1}}{{∂𝑥}^{0}} & = & \frac{𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}}{1 - 𝑘 𝑟^{2}} \\ \frac{{∂𝛤}^{0}_{2}_{2}}{{∂𝑥}^{0}} & = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} \\ \frac{{∂𝛤}^{0}_{3}_{3}}{{∂𝑥}^{0}} & = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} \sin^{2} 𝜃 \\ \frac{{∂𝛤}^{1}_{0}_{1}}{{∂𝑥}^{0}} = \frac{{∂𝛤}^{2}_{0}_{2}}{{∂𝑥}^{0}} = \frac{{∂𝛤}^{3}_{0}_{3}}{{∂𝑥}^{0}} & = & \frac{\ddot{𝑎}}{𝑎} - \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} \end{array}$ $\begin{array}{l} \frac{{∂𝛤}^{1}_{1}_{1}}{{∂𝑥}^{1}} & = & \frac{𝑘 (1 + 𝑘 𝑟^{2})}{{(1 - 𝑘 𝑟^{2})}^{2}} \\ \frac{{∂𝛤}^{1}_{2}_{2}}{{∂𝑥}^{1}} & = & - 1 + 3 𝑘 𝑟^{2} \\ \frac{{∂𝛤}^{1}_{3}_{3}}{{∂𝑥}^{1}} & = & (- 1 + 3 𝑘 𝑟^{2}) \sin^{2} 𝜃 \\ \frac{{∂𝛤}^{2}_{1}_{2}}{{∂𝑥}^{1}} = \frac{{∂𝛤}^{3}_{1}_{3}}{{∂𝑥}^{1}} & = & - \frac{1}{𝑟^{2}} \end{array}$ $\begin{array}{l} \frac{{∂𝛤}^{2}_{3}_{3}}{{∂𝑥}^{2}} & = & - \cos^{2} 𝜃 + \sin^{2} 𝜃 \\ \frac{{∂𝛤}^{3}_{2}_{3}}{{∂𝑥}^{2}} & = & - \frac{1}{\sin^{2} 𝜃} \\ …10個まとめて(28) \end{array}$ 実は $\frac{{∂𝛤}^{1}_{1}_{1}}{{∂𝑥}^{1}}$ はこの後で使わないので計算する必要はなかったのだがついでに書いておいた。

リッチテンソル 𝑅_𝜇𝜈

リーマンテンソル 𝑅^𝜌_𝜇𝜆𝜈 およびリッチテンソル 𝑅_𝜇𝜈 の定義は $\begin{array}{l} 𝑅^{𝜌}_{𝜇}_{𝜆}_{𝜈} = \frac{{∂𝛤}^{𝜌}_{𝜇}_{𝜈}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜌}_{𝜇}_{𝜆}}{{∂𝑥}^{𝜈}} + 𝛤^{𝜎}_{𝜇}_{𝜈} 𝛤^{𝜌}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{𝜇}_{𝜆} 𝛤^{𝜌}_{𝜎}_{𝜈} & (29) \\ 𝑅_{𝜇 𝜈} = 𝑅^{𝜆}_{𝜇}_{𝜆}_{𝜈} & (30) \end{array}$ である（ただし符号を逆に定義する流儀もある）。そこで(29)式を(30)式に代入すれば $𝑅_{𝜇 𝜈} = \frac{{∂𝛤}^{𝜆}_{𝜇}_{𝜈}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{𝜇}_{𝜆}}{{∂𝑥}^{𝜈}} + 𝛤^{𝜎}_{𝜇}_{𝜈} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{𝜇}_{𝜆} 𝛤^{𝜆}_{𝜎}_{𝜈} (31)$ のようになる。

ではリッチテンソルの各成分の表式を求めよう。(31)式を使って、(25)・(28)式と見比べながら0でない成分を代入していくだけである。この下の式変形では、添え字に具体的な数字（0〜3）を代入する段階で、項の値が0でないものだけを残すようにしている。クリストッフェル記号の縮約 𝛤^𝜆_𝜇𝜆 を求める公式を知っている人は使いたくなるかもしれないが、今あれを使うとかえって繁雑になってめんどくさいのでやめた方がよいと思う。

まず対角成分は、 $\begin{array}{l} 𝑅_{00} & = & \frac{{∂𝛤}^{𝜆}_{0}_{0}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{0}_{𝜆}}{{∂𝑥}^{0}} + 𝛤^{𝜎}_{0}_{0} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{0}_{𝜆} 𝛤^{𝜆}_{𝜎}_{0} \\ = & 0 - (\frac{{∂𝛤}^{1}_{0}_{1}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{2}_{0}_{2}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{3}_{0}_{3}}{{∂𝑥}^{0}}) + 0 - (𝛤^{1}_{0}_{1} 𝛤^{1}_{1}_{0} + 𝛤^{2}_{0}_{2} 𝛤^{2}_{2}_{0} + 𝛤^{3}_{0}_{3} 𝛤^{3}_{3}_{0}) \\ = & - 3 (\frac{\ddot{𝑎}}{𝑎} - \frac{{\dot{𝑎}}^{2}}{𝑎^{2}}) - 3 {(\frac{\dot{𝑎}}{𝑎})}^{2} \\ = & - 3 \frac{\ddot{𝑎}}{𝑎} + 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} \\ = & - 3 \frac{\ddot{𝑎}}{𝑎} & (32) \\ 𝑅_{11} & = & \frac{{∂𝛤}^{𝜆}_{1}_{1}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{1}_{𝜆}}{{∂𝑥}^{1}} + 𝛤^{𝜎}_{1}_{1} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{1}_{𝜆} 𝛤^{𝜆}_{𝜎}_{1} \\ = & (\frac{{∂𝛤}^{0}_{1}_{1}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{1}_{1}_{1}}{{∂𝑥}^{1}}) - (\frac{{∂𝛤}^{1}_{1}_{1}}{{∂𝑥}^{1}} + \frac{{∂𝛤}^{2}_{1}_{2}}{{∂𝑥}^{1}} + \frac{{∂𝛤}^{3}_{1}_{3}}{{∂𝑥}^{1}}) \\ + {𝛤^{0}_{1}_{1} (𝛤^{1}_{0}_{1} + 𝛤^{2}_{0}_{2} + 𝛤^{3}_{0}_{3}) + 𝛤^{1}_{1}_{1} (𝛤^{1}_{1}_{1} + 𝛤^{2}_{1}_{2} + 𝛤^{3}_{1}_{3})} \\ - (𝛤^{0}_{1}_{1} 𝛤^{1}_{0}_{1} + 𝛤^{1}_{1}_{0} 𝛤^{0}_{1}_{1} + 𝛤^{1}_{1}_{1} 𝛤^{1}_{1}_{1} + 𝛤^{2}_{1}_{2} 𝛤^{2}_{2}_{1} + 𝛤^{3}_{1}_{3} 𝛤^{3}_{3}_{1}) \\ = & \frac{{∂𝛤}^{0}_{1}_{1}}{{∂𝑥}^{0}} - \frac{{∂𝛤}^{2}_{1}_{2}}{{∂𝑥}^{1}} - \frac{{∂𝛤}^{3}_{1}_{3}}{{∂𝑥}^{1}} \\ + 𝛤^{0}_{1}_{1} (𝛤^{2}_{0}_{2} + 𝛤^{3}_{0}_{3} - 𝛤^{1}_{1}_{0}) + 𝛤^{1}_{1}_{1} (𝛤^{2}_{1}_{2} + 𝛤^{3}_{1}_{3}) - {(𝛤^{2}_{1}_{2})}^{2} - {(𝛤^{3}_{1}_{3})}^{2} \\ = & \frac{𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}}{1 - 𝑘 𝑟^{2}} - (- \frac{1}{𝑟^{2}}) - (- \frac{1}{𝑟^{2}}) \\ + \frac{𝑎 \dot{𝑎}}{1 - 𝑘 𝑟^{2}} (\frac{\dot{𝑎}}{𝑎} + \frac{\dot{𝑎}}{𝑎} - \frac{\dot{𝑎}}{𝑎}) + \frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} (\frac{1}{𝑟} + \frac{1}{𝑟}) - {(\frac{1}{𝑟})}^{2} - {(\frac{1}{𝑟})}^{2} \\ = & \frac{𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}}{1 - 𝑘 𝑟^{2}} + \frac{1}{𝑟^{2}} + \frac{1}{𝑟^{2}} + \frac{𝑎 \dot{𝑎}}{1 - 𝑘 𝑟^{2}} \cdot \frac{\dot{𝑎}}{𝑎} + \frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} \cdot \frac{2}{𝑟} - \frac{1}{𝑟^{2}} - \frac{1}{𝑟^{2}} \\ = & \frac{𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}}{1 - 𝑘 𝑟^{2}} + \frac{{\dot{𝑎}}^{2}}{1 - 𝑘 𝑟^{2}} + \frac{2 𝑘}{1 - 𝑘 𝑟^{2}} \\ = & \frac{𝑎 \ddot{𝑎} + 2 {\dot{𝑎}}^{2} + 2 𝑘}{1 - 𝑘 𝑟^{2}} & (33) \\ 𝑅_{22} & = & \frac{{∂𝛤}^{𝜆}_{2}_{2}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{2}_{𝜆}}{{∂𝑥}^{2}} + 𝛤^{𝜎}_{2}_{2} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{2}_{𝜆} 𝛤^{𝜆}_{𝜎}_{2} \\ = & (\frac{{∂𝛤}^{0}_{2}_{2}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{1}_{2}_{2}}{{∂𝑥}^{1}}) - \frac{{∂𝛤}^{3}_{2}_{3}}{{∂𝑥}^{2}} \\ + {𝛤^{0}_{2}_{2} (𝛤^{1}_{0}_{1} + 𝛤^{2}_{0}_{2} + 𝛤^{3}_{0}_{3}) + 𝛤^{1}_{2}_{2} (𝛤^{1}_{1}_{1} + 𝛤^{2}_{1}_{2} + 𝛤^{3}_{1}_{3})} \\ - (𝛤^{0}_{2}_{2} 𝛤^{2}_{0}_{2} + 𝛤^{1}_{2}_{2} 𝛤^{2}_{1}_{2} + 𝛤^{2}_{2}_{0} 𝛤^{0}_{2}_{2} + 𝛤^{2}_{2}_{1} 𝛤^{1}_{2}_{2} + 𝛤^{3}_{2}_{3} 𝛤^{3}_{3}_{2}) \\ = & \frac{{∂𝛤}^{0}_{2}_{2}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{1}_{2}_{2}}{{∂𝑥}^{1}} - \frac{{∂𝛤}^{3}_{2}_{3}}{{∂𝑥}^{2}} \\ + 𝛤^{0}_{2}_{2} (𝛤^{1}_{0}_{1} + 𝛤^{3}_{0}_{3} - 𝛤^{2}_{2}_{0}) + 𝛤^{1}_{2}_{2} (𝛤^{1}_{1}_{1} + 𝛤^{3}_{1}_{3} - 𝛤^{2}_{2}_{1}) - {(𝛤^{3}_{2}_{3})}^{2} \\ = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} + (- 1 + 3 𝑘 𝑟^{2}) - (- \frac{1}{\sin^{2} 𝜃}) \\ + 𝑎 \dot{𝑎} 𝑟^{2} (\frac{\dot{𝑎}}{𝑎} + \frac{\dot{𝑎}}{𝑎} - \frac{\dot{𝑎}}{𝑎}) + {- (1 - 𝑘 𝑟^{2}) 𝑟} (\frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} + \frac{1}{𝑟} - \frac{1}{𝑟}) - \cot^{2} 𝜃 \\ = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} - 1 + 3 𝑘 𝑟^{2} + \frac{1}{\sin^{2} 𝜃} + 𝑎 \dot{𝑎} 𝑟^{2} \frac{\dot{𝑎}}{𝑎} - (1 - 𝑘 𝑟^{2}) 𝑟 \frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} - \frac{\cos^{2} 𝜃}{\sin^{2} 𝜃} \\ = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} + 3 𝑘 𝑟^{2} + {\dot{𝑎}}^{2} 𝑟^{2} - 𝑘 𝑟^{2} \\ = & (𝑎 \ddot{𝑎} + 2 {\dot{𝑎}}^{2} + 2 𝑘) 𝑟^{2} & (34) \\ 𝑅_{33} & = & \frac{{∂𝛤}^{𝜆}_{3}_{3}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{3}_{𝜆}}{{∂𝑥}^{3}} + 𝛤^{𝜎}_{3}_{3} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{3}_{𝜆} 𝛤^{𝜆}_{𝜎}_{3} \\ = & (\frac{{∂𝛤}^{0}_{3}_{3}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{1}_{3}_{3}}{{∂𝑥}^{1}} + \frac{{∂𝛤}^{2}_{3}_{3}}{{∂𝑥}^{2}}) - 0 \\ + {𝛤^{0}_{3}_{3} (𝛤^{1}_{0}_{1} + 𝛤^{2}_{0}_{2} + 𝛤^{3}_{0}_{3}) + 𝛤^{1}_{3}_{3} (𝛤^{1}_{1}_{1} + 𝛤^{2}_{1}_{2} + 𝛤^{3}_{1}_{3}) + 𝛤^{2}_{3}_{3} 𝛤^{3}_{2}_{3}} \\ - (𝛤^{0}_{3}_{3} 𝛤^{3}_{0}_{3} + 𝛤^{1}_{3}_{3} 𝛤^{3}_{1}_{3} + 𝛤^{2}_{3}_{3} 𝛤^{3}_{2}_{3} + 𝛤^{3}_{3}_{0} 𝛤^{0}_{3}_{3} + 𝛤^{3}_{3}_{1} 𝛤^{1}_{3}_{3} + 𝛤^{3}_{3}_{2} 𝛤^{2}_{3}_{3}) \\ = & \frac{{∂𝛤}^{0}_{3}_{3}}{{∂𝑥}^{0}} + \frac{{∂𝛤}^{1}_{3}_{3}}{{∂𝑥}^{1}} + \frac{{∂𝛤}^{2}_{3}_{3}}{{∂𝑥}^{2}} \\ + 𝛤^{0}_{3}_{3} (𝛤^{1}_{0}_{1} + 𝛤^{2}_{0}_{2} - 𝛤^{3}_{3}_{0}) + 𝛤^{1}_{3}_{3} (𝛤^{1}_{1}_{1} + 𝛤^{2}_{1}_{2} - 𝛤^{3}_{3}_{1}) - 𝛤^{3}_{3}_{2} 𝛤^{2}_{3}_{3} \\ = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} \sin^{2} 𝜃 + (- 1 + 3 𝑘 𝑟^{2}) \sin^{2} 𝜃 + (- \cos^{2} 𝜃 + \sin^{2} 𝜃) \\ + 𝑎 \dot{𝑎} 𝑟^{2} \sin^{2} 𝜃 (\frac{\dot{𝑎}}{𝑎} + \frac{\dot{𝑎}}{𝑎} - \frac{\dot{𝑎}}{𝑎}) + {- (1 - 𝑘 𝑟^{2}) 𝑟 \sin^{2} 𝜃} (\frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} + \frac{1}{𝑟} - \frac{1}{𝑟}) \\ - \cot 𝜃 (- \sin 𝜃 \cos 𝜃) \\ = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} \sin^{2} 𝜃 - \sin^{2} 𝜃 + 3 𝑘 𝑟^{2} \sin^{2} 𝜃 - \cos^{2} 𝜃 + \sin^{2} 𝜃 \\ + 𝑎 \dot{𝑎} 𝑟^{2} \sin^{2} 𝜃 \frac{\dot{𝑎}}{𝑎} - (1 - 𝑘 𝑟^{2}) 𝑟 \sin^{2} 𝜃 \frac{𝑘 𝑟}{1 - 𝑘 𝑟^{2}} + \cos^{2} 𝜃 \\ = & (𝑎 \ddot{𝑎} + {\dot{𝑎}}^{2}) 𝑟^{2} \sin^{2} 𝜃 + 3 𝑘 𝑟^{2} \sin^{2} 𝜃 + {\dot{𝑎}}^{2} 𝑟^{2} \sin^{2} 𝜃 - 𝑘 𝑟^{2} \sin^{2} 𝜃 \\ = & (𝑎 \ddot{𝑎} + 2 {\dot{𝑎}}^{2} + 2 𝑘) 𝑟^{2} \sin^{2} 𝜃 & (35) \end{array}$ のようになる。続いて非対角成分は、 $\begin{array}{l} 𝑅_{01} = 𝑅_{10} & = & \frac{{∂𝛤}^{𝜆}_{0}_{1}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{0}_{𝜆}}{{∂𝑥}^{1}} + 𝛤^{𝜎}_{0}_{1} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{0}_{𝜆} 𝛤^{𝜆}_{𝜎}_{1} \\ = & 0 - 0 + 𝛤^{1}_{0}_{1} (𝛤^{1}_{1}_{1} + 𝛤^{2}_{1}_{2} + 𝛤^{3}_{1}_{3}) \\ - (𝛤^{1}_{0}_{1} 𝛤^{1}_{1}_{1} + 𝛤^{2}_{0}_{2} 𝛤^{2}_{2}_{1} + 𝛤^{3}_{0}_{3} 𝛤^{3}_{3}_{1}) \\ = & 𝛤^{1}_{0}_{1} 𝛤^{2}_{1}_{2} + 𝛤^{1}_{0}_{1} 𝛤^{3}_{1}_{3} - 𝛤^{2}_{0}_{2} 𝛤^{2}_{2}_{1} - 𝛤^{3}_{0}_{3} 𝛤^{3}_{3}_{1} \\ = & \frac{\dot{𝑎}}{𝑎} \cdot \frac{1}{𝑟} + \frac{\dot{𝑎}}{𝑎} \cdot \frac{1}{𝑟} - \frac{\dot{𝑎}}{𝑎} \cdot \frac{1}{𝑟} - \frac{\dot{𝑎}}{𝑎} \cdot \frac{1}{𝑟} \\ = & 0 & (36) \\ 𝑅_{02} = 𝑅_{20} & = & \frac{{∂𝛤}^{𝜆}_{0}_{2}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{0}_{𝜆}}{{∂𝑥}^{2}} + 𝛤^{𝜎}_{0}_{2} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{0}_{𝜆} 𝛤^{𝜆}_{𝜎}_{2} \\ = & 0 - 0 + 𝛤^{2}_{0}_{2} 𝛤^{3}_{2}_{3} - 𝛤^{3}_{0}_{3} 𝛤^{3}_{3}_{2} \\ = & \frac{\dot{𝑎}}{𝑎} \cot 𝜃 - \frac{\dot{𝑎}}{𝑎} \cot 𝜃 \\ = & 0 & (37) \\ 𝑅_{03} = 𝑅_{30} & = & \frac{{∂𝛤}^{𝜆}_{0}_{3}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{0}_{𝜆}}{{∂𝑥}^{3}} + 𝛤^{𝜎}_{0}_{3} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{0}_{𝜆} 𝛤^{𝜆}_{𝜎}_{3} \\ = & 0 - 0 + 0 - 0 \\ = & 0 & (38) \\ 𝑅_{12} = 𝑅_{21} & = & \frac{{∂𝛤}^{𝜆}_{1}_{2}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{1}_{𝜆}}{{∂𝑥}^{2}} + 𝛤^{𝜎}_{1}_{2} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{1}_{𝜆} 𝛤^{𝜆}_{𝜎}_{2} \\ = & 0 - 0 + 𝛤^{2}_{1}_{2} 𝛤^{3}_{2}_{3} - 𝛤^{3}_{1}_{3} 𝛤^{3}_{3}_{2} \\ = & \frac{1}{𝑟} \cot 𝜃 - \frac{1}{𝑟} \cot 𝜃 \\ = & 0 & (39) \\ 𝑅_{13} = 𝑅_{31} & = & \frac{{∂𝛤}^{𝜆}_{1}_{3}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{1}_{𝜆}}{{∂𝑥}^{3}} + 𝛤^{𝜎}_{1}_{3} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{1}_{𝜆} 𝛤^{𝜆}_{𝜎}_{3} \\ = & 0 - 0 + 0 - 0 \\ = & 0 & (40) \\ 𝑅_{23} = 𝑅_{32} & = & \frac{{∂𝛤}^{𝜆}_{2}_{3}}{{∂𝑥}^{𝜆}} - \frac{{∂𝛤}^{𝜆}_{2}_{𝜆}}{{∂𝑥}^{3}} + 𝛤^{𝜎}_{2}_{3} 𝛤^{𝜆}_{𝜎}_{𝜆} - 𝛤^{𝜎}_{2}_{𝜆} 𝛤^{𝜆}_{𝜎}_{3} \\ = & 0 - 0 + 0 - 0 \\ = & 0 & (41) \end{array}$ のようにすべて0である。

以上でリッチテンソル 𝑅_𝜇𝜈 の表式が求まった。

リッチテンソル 𝑅^𝜇_𝜈

次に、1個目の添え字を上にあげた 𝑅^𝜇_𝜈 = 𝑔^𝜇𝜌𝑅_𝜌𝜈 を計算する。(13)式と(32)〜(41)式を代入すればよい。

まず対角成分は次のようになる。 $\begin{array}{l} 𝑅^{0}_{0} & = & 𝑔^{0 𝜌} 𝑅_{𝜌 0} \\ = & 𝑔^{00} 𝑅_{00} \\ = & (- 1) (- 3 \frac{\ddot{𝑎}}{𝑎}) \\ = & 3 \frac{\ddot{𝑎}}{𝑎} & (42) \\ 𝑅^{1}_{1} & = & 𝑔^{1 𝜌} 𝑅_{𝜌 1} \\ = & 𝑔^{11} 𝑅_{11} \\ = & \frac{1 - 𝑘 𝑟^{2}}{𝑎^{2}} \cdot \frac{𝑎 \ddot{𝑎} + 2 {\dot{𝑎}}^{2} + 2 𝑘}{1 - 𝑘 𝑟^{2}} \\ = & \frac{\ddot{𝑎}}{𝑎} + 2 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 2 \frac{𝑘}{𝑎^{2}} & (43) \\ 𝑅^{2}_{2} & = & 𝑔^{2 𝜌} 𝑅_{𝜌 2} \\ = & 𝑔^{22} 𝑅_{22} \\ = & \frac{1}{𝑎^{2} 𝑟^{2}} \cdot (𝑎 \ddot{𝑎} + 2 {\dot{𝑎}}^{2} + 2 𝑘) 𝑟^{2} \\ = & \frac{\ddot{𝑎}}{𝑎} + 2 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 2 \frac{𝑘}{𝑎^{2}} \\ = & 𝑅^{1}_{1} & (44) \\ 𝑅^{3}_{3} & = & 𝑔^{3 𝜌} 𝑅_{𝜌 3} \\ = & 𝑔^{33} 𝑅_{33} \\ = & \frac{1}{𝑎^{2} 𝑟^{2} \sin^{2} 𝜃} \cdot (𝑎 \ddot{𝑎} + 2 {\dot{𝑎}}^{2} + 2 𝑘) 𝑟^{2} \sin^{2} 𝜃 \\ = & \frac{\ddot{𝑎}}{𝑎} + 2 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 2 \frac{𝑘}{𝑎^{2}} \\ = & 𝑅^{1}_{1} & (45) \end{array}$ (44)・(45)式は、 𝑅^²_₂ や 𝑅^³_₃ を計算した結果を 𝑅^¹_₁ と見比べてみたら 𝑅^¹_₁ に等しいことがわかった、という意味である。

続いて非対角成分であるが、計量テンソル 𝑔^𝜇𝜈 もリッチテンソル 𝑅_𝜇𝜈 も非対角成分はすべて0なので、 𝑅^𝜇_𝜈 = 𝑔^𝜇𝜌𝑅_𝜌𝜈 も非対角成分はすべて0である。

以上でリッチテンソル 𝑅^𝜇_𝜈 の表式が求まった。

スカラー曲率

スカラー曲率（リッチスカラー） 𝑅 の定義は $\begin{array}{l} 𝑅 & = & 𝑔^{𝜇 𝜈} 𝑅_{𝜇 𝜈} \\ = & 𝑅^{𝜈}_{𝜈} \end{array}$ であるから、(42)〜(45)式を代入して計算すると次のようになる。 $\begin{array}{l} 𝑅 & = & 𝑅^{𝜈}_{𝜈} \\ = & 𝑅^{0}_{0} + 𝑅^{1}_{1} + 𝑅^{2}_{2} + 𝑅^{3}_{3} \\ = & 𝑅^{0}_{0} + 𝑅^{1}_{1} + 𝑅^{1}_{1} + 𝑅^{1}_{1} \\ = & 𝑅^{0}_{0} + 3 𝑅^{1}_{1} \\ = & 3 \frac{\ddot{𝑎}}{𝑎} + 3 (\frac{\ddot{𝑎}}{𝑎} + 2 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 2 \frac{𝑘}{𝑎^{2}}) \\ = & 3 \frac{\ddot{𝑎}}{𝑎} + 3 \frac{\ddot{𝑎}}{𝑎} + 6 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 6 \frac{𝑘}{𝑎^{2}} \\ = & 6 \frac{\ddot{𝑎}}{𝑎} + 6 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 6 \frac{𝑘}{𝑎^{2}} & (46) \end{array}$

アインシュタインテンソル

アインシュタインテンソル 𝐺^𝜇_𝜈 の定義は $𝐺^{𝜇}_{𝜈} = 𝑅^{𝜇}_{𝜈} - \frac{1}{2} 𝛿^{𝜇}_{𝜈} 𝑅$ である。では(42)〜(46)式を代入して各成分の表式を求めよう。

まず対角成分は次のようになる。 $\begin{array}{l} 𝐺^{0}_{0} & = & 𝑅^{0}_{0} - \frac{1}{2} 𝛿^{0}_{0} 𝑅 \\ = & 3 \frac{\ddot{𝑎}}{𝑎} - \frac{1}{2} (6 \frac{\ddot{𝑎}}{𝑎} + 6 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 6 \frac{𝑘}{𝑎^{2}}) \\ = & 3 \frac{\ddot{𝑎}}{𝑎} - 3 \frac{\ddot{𝑎}}{𝑎} - 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - 3 \frac{𝑘}{𝑎^{2}} \\ = & - 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - 3 \frac{𝑘}{𝑎^{2}} \\ 𝐺^{1}_{1} & = & 𝑅^{1}_{1} - \frac{1}{2} 𝛿^{1}_{1} 𝑅 \\ = & (\frac{\ddot{𝑎}}{𝑎} + 2 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 2 \frac{𝑘}{𝑎^{2}}) - \frac{1}{2} (6 \frac{\ddot{𝑎}}{𝑎} + 6 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 6 \frac{𝑘}{𝑎^{2}}) \\ = & \frac{\ddot{𝑎}}{𝑎} + 2 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} + 2 \frac{𝑘}{𝑎^{2}} - 3 \frac{\ddot{𝑎}}{𝑎} - 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - 3 \frac{𝑘}{𝑎^{2}} \\ = & - 2 \frac{\ddot{𝑎}}{𝑎} - \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - \frac{𝑘}{𝑎^{2}} \\ 𝐺^{2}_{2} & = & 𝑅^{2}_{2} - \frac{1}{2} 𝛿^{2}_{2} 𝑅 \\ = & 𝑅^{1}_{1} - \frac{1}{2} 𝛿^{1}_{1} 𝑅 \\ = & 𝐺^{1}_{1} \\ 𝐺^{3}_{3} & = & 𝑅^{3}_{3} - \frac{1}{2} 𝛿^{3}_{3} 𝑅 \\ = & 𝑅^{1}_{1} - \frac{1}{2} 𝛿^{1}_{1} 𝑅 \\ = & 𝐺^{1}_{1} \end{array}$

続いて非対角成分であるが、リッチテンソル 𝑅^𝜇_𝜈 もクロネッカーのデルタ 𝛿^𝜇_𝜈 も非対角成分はすべて0なので、アインシュタインテンソル 𝐺^𝜇_𝜈 も非対角成分はすべて0である。

方程式の完成

ここまででアインシュタインテンソル 𝐺^𝜇_𝜈 の各成分の表式が求まったので、それらを(27)式に代入すれば、 $\begin{array}{l} 第(0, 0)成分: & - 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - 3 \frac{𝑘}{𝑎^{2}} + 𝛬 = \frac{8 𝜋 𝐺}{𝑐^{4}} 𝑇^{0}_{0} & (47) \\ 第(1, 1), (2, 2), (3, 3)成分: & - 2 \frac{\ddot{𝑎}}{𝑎} - \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - \frac{𝑘}{𝑎^{2}} + 𝛬 = \frac{8 𝜋 𝐺}{𝑐^{4}} 𝑇^{𝜇}_{𝜈} & (𝜇 = 𝜈 = 1, 2, 3) & (48) \\ 非対角成分: & 0 = \frac{8 𝜋 𝐺}{𝑐^{4}} 𝑇^{𝜇}_{𝜈} & (𝜇 \neq 𝜈) & (49) \end{array}$ のようになる。これらが一様・等方な時空が満たすべき重力場の方程式である。

右辺についてはまだ何の条件も考えていなかったが、(47)〜(49)式からエネルギー運動量テンソル 𝑇^𝜇_𝜈 に条件が課されることがわかる。その条件とは非対角成分が0であることと、 𝑇^¹_₁ = 𝑇^²_₂ = 𝑇^³_₃ である。それと左辺は時間座標 𝑤 の関数であるから、右辺の 𝑇^𝜇_𝜈 も 𝑤 の関数である。したがって $(𝑇^{𝜇}_{𝜈}) = (\begin{matrix} 𝑈 (𝑤) & 0 & 0 & 0 \\ 0 & 𝑉 (𝑤) & 0 & 0 \\ 0 & 0 & 𝑉 (𝑤) & 0 \\ 0 & 0 & 0 & 𝑉 (𝑤) \end{matrix}) (50)$ のような形になっているべきである。 𝑈(𝑤) および 𝑉(𝑤) は今のところ未定の 𝑤 の関数である。

ところで話は変わるが、一般に完全流体のエネルギー運動量テンソルは $𝑇^{𝜇 𝜈} = (𝜌 + \frac{𝑝}{𝑐^{2}}) 𝑢^{𝜇} 𝑢^{𝜈} + 𝑝 𝑔^{𝜇 𝜈} (51)$ と表されるのだった（ただしどれかの項の符号が逆になる流儀もある）。ここで 𝑢^𝜇 は流体要素の4元速度、 𝜌 は質量密度、 𝑝 は圧力である。 𝜌 と 𝑝 は流体要素が瞬間的に静止する局所慣性系で測った値であるからスカラーである。(51)式の片方の添え字をおろして混合テンソルにすれば $𝑇^{𝜇}_{𝜈} = (𝜌 + \frac{𝑝}{𝑐^{2}}) 𝑢^{𝜇} 𝑢_{𝜈} + 𝑝 𝛿^{𝜇}_{𝜈} (52)$ のようになる。もし今考えている座標系で流体要素が静止していれば、 𝑢^¹ = 𝑢^² = 𝑢^³ = 𝑢_₁ = 𝑢_₂ = 𝑢_₃ = 0, 𝑢^⁰ = 𝑐, 𝑢_₀ = −𝑐 であるからこれを(52)式に代入すると、 $(𝑇^{𝜇}_{𝜈}) = (\begin{matrix} - 𝑐^{2} 𝜌 & 0 & 0 & 0 \\ 0 & 𝑝 & 0 & 0 \\ 0 & 0 & 𝑝 & 0 \\ 0 & 0 & 0 & 𝑝 \end{matrix}) (53)$ となる。ここで、 𝜌 と 𝑝 が 𝑤 のみの関数になっていることにすれば(53)式は(50)式の条件に合致している。空間は一様・等方だから 𝜌 と 𝑝 が空間座標に依存しないことは理にかなっている。よって、静止している完全流体のエネルギー運動量テンソルである(53)式を、(47)〜(49)式の 𝑇^𝜇_𝜈 として採用することにする。そういえば第1章の最後で、この座標系で静止している粒子の世界線は測地線であることを確認したのだった。そのことも流体要素が静止し続けることと整合している。「静止している完全流体」以外で(50)式の形になることはあり得ないのかどうか知らないが、そういう話が載っている本を見たことがないので、たぶん気にしなくてよいのだろう。

ただし今後の説明の都合上、ここでは質量密度 𝜌 の代わりにエネルギー密度 𝜀 を使うことにする。質量とエネルギーは等価であり 𝜀 = 𝜌𝑐^² だから値が 𝑐^² 倍違うだけで式の意味は同じであり、(53)式は $(𝑇^{𝜇}_{𝜈}) = (\begin{matrix} - 𝜀 & 0 & 0 & 0 \\ 0 & 𝑝 & 0 & 0 \\ 0 & 0 & 𝑝 & 0 \\ 0 & 0 & 0 & 𝑝 \end{matrix}) (54)$ となる。そういうわけで(54)式を(47)〜(49)式に代入すると、 $\begin{array}{l} 第(0, 0)成分: & - 3 \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - 3 \frac{𝑘}{𝑎^{2}} + 𝛬 = - \frac{8 𝜋 𝐺}{𝑐^{4}} 𝜀 & (55) \\ 第(1, 1), (2, 2), (3, 3)成分: & - 2 \frac{\ddot{𝑎}}{𝑎} - \frac{{\dot{𝑎}}^{2}}{𝑎^{2}} - \frac{𝑘}{𝑎^{2}} + 𝛬 = \frac{8 𝜋 𝐺}{𝑐^{4}} 𝑝 & (56) \end{array}$ となり、非対角成分は 0 = 0 となって消える。 𝑎 と 𝜀 と 𝑝 は時間座標 𝑤 の関数であり、その他の文字は定数である。

このページでやったことを改めて振り返っておこう。重力場の方程式の一般的な形である(27)式の左辺に、一様・等方な時空の計量であるFLRW計量を代入したものが、(47)〜(49)式である。そしてこれらに矛盾しないようなエネルギー運動量テンソルとして(54)式を採用して右辺に代入したものが、(55)・(56)式である。一様・等方な時空は(55)・(56)式を満たすはずだ。