FLRW計量のリーマン曲率

フリードマン–ルメートル–ロバートソン–ウォーカー（FLRW）時空は、 均質かつ等方的な宇宙を記述します。共動座標$(w,r,\theta,\phi)$と 単位系$c=1$を用いると、

$$ ds^2=-dw^2+a(w)^2\left( \frac{dr^2}{1-Kr^2}+r^2d\theta^2+r^2\sin^2\theta\,d\phi^2 \right). $$

尺度因子$a(w)$は任意の記号関数のまま扱い、$K$は一定の 空間曲率パラメータとします。

計量の設定

$W(r)=(1-Kr^2)^{-1}$とおきます。選んだ接続や曲率の成分を 短時間で計算できるよう、逆計量を明示的に入力します。

declare symbol w, r, θ, φ, K: MathValue

def x : Vector MathValue := [| w, r, θ, φ |]
def a : MathValue := function (w)

def W (r: MathValue) : MathValue := 1 / `(1 - K * r^2)

def g_i_j : Matrix MathValue :=
  [| [| -1, 0, 0, 0 |]
   , [| 0, a^2 * W r, 0, 0 |]
   , [| 0, 0, a^2 * r^2, 0 |]
   , [| 0, 0, 0, a^2 * r^2 * (sin θ)^2 |]
   |]_i_j

def g~i~j : Matrix MathValue :=
  [| [| -1, 0, 0, 0 |]
   , [| 0, 1 / (a^2 * W r), 0, 0 |]
   , [| 0, 0, 1 / (a^2 * r^2), 0 |]
   , [| 0, 0, 0, 1 / (a^2 * r^2 * (sin θ)^2) |]
   |]~i~j

W r
$(-K r^{2} + 1)^{-1}$
g_#_#
$\begin{pmatrix} -1 & 0 & 0 & 0 \\ 0 & (-K r^{2} + 1)^{-1} a^{2} & 0 & 0 \\ 0 & 0 & a^{2} r^{2} & 0 \\ 0 & 0 & 0 & \sin(θ)^{2} a^{2} r^{2} \\ \end{pmatrix}_{\#\#}^{\;\;}$
g~#~#
$\begin{pmatrix} -1 & 0 & 0 & 0 \\ 0 & (-K r^{2} + 1) a^{-2} & 0 & 0 \\ 0 & 0 & a^{-2} r^{-2} & 0 \\ 0 & 0 & 0 & \sin(θ)^{-2} a^{-2} r^{-2} \\ \end{pmatrix}_{\;\;}^{\#\#}$

宇宙膨張と接続

空間計量を時間微分すると、$a'(w)$に比例するクリストッフェル記号が 現れます。特に、動径方向の膨張は$\Gamma^w{}_{rr}$と $\Gamma^r{}_{wr}$の両方に現れます。

def Γ_i_j_k : Tensor MathValue :=
  (1 / 2) * (∂/∂ g_i_k x~j + ∂/∂ g_i_j x~k - ∂/∂ g_j_k x~i)

def Γ~i_j_k : Tensor MathValue := withSymbols [m]
  g~i~m . Γ_m_j_k

Γ~1_2_2
$(-K r^{2} + 1)^{-1} a \frac{\partial a}{\partial 1}$
Γ~2_1_2
$a^{-1} \frac{\partial a}{\partial 1}$

リーマンテンソルと縮約

時間と空間が混ざる成分は尺度因子の加速度を測り、空間成分だけからなる 成分では$(a')^2$と$K$が組み合わさります。完全な縮約は

$$ \operatorname{Ric}_{ij}=R^m{}_{imj},\qquad \mathcal R=g^{ij}\operatorname{Ric}_{ij}. $$

def R~i_j_k_l : Tensor MathValue := withSymbols [m]
  ∂/∂ Γ~i_j_l x~k - ∂/∂ Γ~i_j_k x~l
    + Γ~m_j_l . Γ~i_m_k - Γ~m_j_k . Γ~i_m_l

def Ric_i_j : Matrix MathValue := withSymbols [m]
  sum (contract R~m_i_m_j)

def scalarCurvature : MathValue := withSymbols [i, j]
  expandAll' (g~i~j . Ric_i_j)

R~1_2_1_2
$(-K r^{2} + 1)^{-1} a \frac{\partial^2 a}{\partial 1^2}$
R~2_3_2_3
$\frac{\partial a}{\partial 1}^{2} r^{2} + K r^{2}$

スカラー曲率とその意味

この規約では、期待されるスカラー曲率は

$$ \mathcal R =6\left(\frac{a''(w)}{a(w)}+ \frac{a'(w)^2+K}{a(w)^2}\right) =\frac{6\bigl(a''a+(a')^2+K\bigr)}{a^2}. $$

です。選んだリーマンテンソルの成分から、加速膨張と空間切片の曲率という 2つの幾何学的要素が見て取れます。scalarCurvatureを評価すると、 数式処理系は任意関数$a$について完全な縮約を構成して簡約するため、 数分かかる場合があります。そのため、意図的に自動実行する出力セルには せず、定義だけを示しています。

リンク

Egison 数学ノート目次に戻る