[論文レビュー] Spherically symmetric cosmology: resource paper

A. A. Coley, W C Lim+1 more|ArXiv.org|Mar 6, 2008

Cosmology and Gravitation Theories参考文献 3被引用数 6

ひとこと要約

本稿では、1+3フレーム形式を用いて球対称な宇宙論的モデルを定式化し、時間と径方向座標の2変数における適切に定義された一階偏微分方程式系を導出することで、強固な数値的・定性的解析が可能となる。主な貢献は、完全流体およびダスト（例えば、Lemaître-Tolman-Bondiモデル）の力学を単純化する体系的で、共変的かつゲージ不変な枠組みを提供することであり、相対論的流体に対しては $ \gamma \geq 1 $ の場合に適切に定義された系が示され、非渦度のあるダスト（irrotational dust）の正確な解は、遅刻にミルン型の振る舞いへと等方化することを示している。

ABSTRACT

We use the 1+3 frame formalism to write down the evolution equations for spherically symmetric models as a well-posed system of first order PDEs in two variables, suitable for numerical and qualitative analysis.

研究の動機と目的

• 球対称な宇宙論的モデルを分析するための体系的で、共変的かつ数値的に取り扱いやすい枠組みを構築すること。
• 1+3形式を用いてアインシュタイン場方程式を2変数（時間と径方向座標）における適切に定義された一階偏微分方程式系に再定式化すること。
• 特にLemaître-Tolman-Bondiダストモデルにおける不均一モデルの構造形成の定性的および数値的解析を可能にすること。
• 特に方程式状態パラメータ $ \gamma \geq 1 $ の相対論的完全流体に対して、進化系の適切性を確立すること。
• LTBモデルの遅刻における漸近的振る舞いを示し、等方的ミルン型解への収束を示すこと。

提案手法

• 球対称な時空に1+3フレーム形式を適用し、時間的単位ベクトルと正規直交空間トライアドを用いて計量と接続を時間的および空間的成分に分解する。
• 運動学的変数（膨張、歪み、加速度）および力学的量（エネルギー密度、圧力、熱フラックス、非等方的応力）を正規直交フレームで表現し、対称性により独立成分の数を削減する。
• 進化方程式はアインシュタイン場方程式から導出され、時間と径方向座標に関する一階偏微分方程式系が得られ、ガウスおよびコーディー方程式からの制約が含まれる。
• ガウス曲率 $ {}^{2}K = (e_2^2)^2 $ を導入することで $ e_2^2 $ を置き換え、$ K $ を基本変数として用いることで、角度依存性を排除する。
• 空間微分係数行列の固有値を用いた適切性の分析により、$ \gamma > 1 $ の場合に双曲型（したがって適切に定義された）であることが示され、$ \gamma < 1 $ の場合に楕円型となる。
• LTBモデル（ダスト、$ \gamma = 1 $、$ v = 0 $）では、系が空間微分から分離され、自明に適切に定義されたものとなり、遅刻における振る舞いの正確な解析が可能となる。

実験結果

リサーチクエスチョン

• RQ1 球対称な宇宙論的モデルを、適切に定義された一階偏微分方程式系として体系的に定式化する方法は何か？
• RQ2 球対称性下での相対論的完全流体の進化方程式の適切性を保証する条件は何か？
• RQ3 不均一なダストモデル（例：LTB）は、遅刻に至る過程でどのように等方化に向かうか？
• RQ4 ガウス曲率 $ {}^{2}K $ は、力学的方程式の簡略化に果たす役割は何か？
• RQ5 1+3形式は、不均一宇宙論における構造形成の数値的および定性的解析をどのように容易にするか？

主な発見

• 完全流体に対して $ \gamma > 1 $ の場合、進化方程式系は適切に定義される。空間微分係数行列は実数で相異なる固有値を持ち、対角化可能であり、音速は $ c_s = \sqrt{\gamma - 1} $ である。
• $ \gamma < 1 $ の場合、系は楕円型となり、適切に定義されない。これは因果構造と予測可能性の破綻を示している。
• 非渦度のあるダスト（LTB）モデルでは、進化方程式が空間微分に依存しないため、系は自明に適切に定義され、正確な解を用いた解析が可能である。
• 遅刻において、常に膨張を続けるLTB解は等方化し、ミルン解に収束する。$ \beta \sim e^{-\tau} $、$ \sigma_+ / \beta \sim \tau e^{-\tau} $、$ \mu / (3\beta^2) \sim e^{-\tau} $ であり、均一性および等方性への収束が示唆される。
• LTBモデルにおける密度対比および歪みパラメータは、FLRWモデルのそれらに類似した動的挙動を示し、構造形成の3次元力学系解析が可能となる。
• $ e_2^2 $ の代わりに $ {}^{2}K $ を用いることで、$ \vartheta $-依存性が効果的に隠され、方程式が簡略化されつつも幾何的意味が保持される。