光や物質で満たされた宇宙に対するフリードマン方程式を解く(1)

前回の記事は「真空の宇宙に対するフリードマン方程式を解く」であった。今回は一歩進んで、一様に光や物質で満たされているような宇宙に対して、状態方程式の具体的な形を仮定してフリードマン方程式の解を求めてみる。現実の宇宙と合うか合わないかといった検討は後回しにして、計算上出てくる解は全部出してみよう。今の段階では宇宙定数は特に意味を考えずに単なる定数とみなす。とりあえず数学的にただ方程式を解くことに集中し、観測データとの関係とか宇宙項の意味といった話はまだここでは扱わない。

目次

第1章 フリードマン方程式等

フリードマン方程式およびそれに付随する方程式の一般的な形は フリードマン方程式: 𝑎˙2𝑎2= 8𝜋𝐺3𝑐4𝜀 𝑘𝑎2 +𝛬3 (1) 加速度方程式: 𝑎¨𝑎= 4𝜋𝐺3𝑐4 (𝜀+3𝑝) +𝛬3 (2) エネルギー保存則: 𝜀˙ +3𝑎˙𝑎(𝜀+𝑝) =0 (3) 状態方程式: 𝑝=𝑝(𝜀) (4) であった。 𝑎 と 𝜀 と 𝑝 は時間の関数であり、 𝑎 はスケール因子、 𝜀 は宇宙を満たす完全流体のエネルギー密度、 𝑝 は宇宙を満たす完全流体の圧力である。ドット˙は時間座標 𝑐𝑡 による微分 1𝑐dd𝑡 である。その他の文字は定数であり、 𝐺 は万有引力定数、 𝑐 は光速、 𝛬 は宇宙定数、 𝑘 は空間の曲率である。

4個の式を書いたが独立な条件は3個である。解が恒等的に 𝑎˙=0 (𝑎 が定数)でない限り、フリードマン方程式(1)式が成り立っている条件の下で、加速度方程式(2)式とエネルギー保存則(3)式は一方が成り立てば自動的に他方も成り立つのでどちらか一方は無視してよい。ただし解が恒等的に 𝑎˙=0 (𝑎 が定数)であるときは、フリードマン方程式(1)式が成り立てば自動的にエネルギー保存則(3)式も成り立つけれども、加速度方程式(2)式は独立になるので無視してはならない。ここまでは「フリードマン方程式の導出」の記事でやったことだ。

具体的にこの連立方程式を解こうとすると、 𝑎 が定数関数であるか否かによってやることが異なってくる。そこで第2章で 𝑎 が定数関数の場合を扱い、第3章でそうでない場合を扱う。

なお、仮定により 𝑎 ≧ 0 (等号成立は宇宙の始まりと終わりの瞬間だけ)である。

第2章 スケール因子が一定の解

スケール因子 𝑎 が一定の場合、 𝑎˙=𝑎¨=0 であるから解くべき方程式は フリードマン方程式: 0= 8𝜋𝐺3𝑐4𝜀 𝑘𝑎2 +𝛬3 (5) 加速度方程式: 0= 4𝜋𝐺3𝑐4 (𝜀+3𝑝) +𝛬3 (6) 状態方程式: 𝑝=𝑝(𝜀) (7) である。フリードマン方程式(5)式が成り立てばエネルギー保存則の式は自動的に成り立つ。(5)・(6)式はもはや微分方程式でなくただの代数方程式であり、とても簡単な式変形により 𝛬= 4𝜋𝐺𝑐4 (𝜀+3𝑝) 𝑘𝑎2= 4𝜋𝐺𝑐4 (𝜀+𝑝) であることがわかる。 𝑎 の他に 𝜀 と 𝑝 が未知関数だったのでそれらについて解けば 𝜀= 𝑐48𝜋𝐺 (3𝑘𝑎2𝛬) (8) 𝑝= 𝑐48𝜋𝐺 (𝛬𝑘𝑎2) (9) となり、当然ながら定数関数になる。ここで 𝑎 と 𝑘 に関してはその値を単独で見てもあまり意味はなく、2つを組み合わせた 𝑎𝑘 が空間の曲率半径の実際の大きさである。だから式の中に現れる 𝑘𝑎2 が曲率半径の−2乗を表しているということだ。

もしも(8)・(9)式の 𝜀 と 𝑝 が状態方程式(7)式を満たしていれば、これらは解である。満たしていなければ、スケール因子 𝑎 が一定となる解は存在しない。

ここでは状態方程式の例として1つだけ挙げておく。それは流体のエネルギー密度に比べて圧力が無視できる場合、すなわち宇宙にある物質の速度が光速に比べて十分に小さい場合で、そのときは近似的に 𝑝=0 とすることができる(この状態方程式については3.2節の冒頭で再度説明する)。このとき空間の曲率半径は 𝑎𝑘 =1𝛬 =𝑐24𝜋𝐺𝜀 であり、エネルギー密度 𝜀 に着目すれば 𝜀 =𝑐4𝛬4𝜋𝐺 =𝑐44𝜋𝐺𝑘𝑎2 となっていればよい。このようになっている宇宙を「アインシュタインの静止宇宙」とか単に「アインシュタイン宇宙」と呼ぶ。

ほかの状態方程式も試したければ、(7)式を好きな形に仮定してそこに(8)・(9)式を代入すれば良い。特に難しいことはなく 𝑘𝑎2 と 𝛬 が満たすべき関係が出てくるはずだ。 𝑎 が一定の解についてはこれくらいにしておこう。

⛭ 数式の表示設定 (S)