フリードマン・ルメートル・ロバートソン・ウォーカー計量

Dec 23 2016

以下のコードは、フリードマン・ルメートル・ロバートソン・ウォーカー計量のリーマン曲率テンソルを計算しています。

-- Parameters
def x := [| w, r, θ, φ |]

-- Scale factor function
def W $r := 1 / '(1 - K * r^2)

-- Metric tensor
def g_i_j :=
  [| [| -1, 0, 0, 0 |]
   , [| 0, `a w^2 * W r, 0, 0 |]
   , [| 0, 0, `a w^2 * r^2, 0 |]
   , [| 0, 0, 0, `a w^2 * r^2 * (sin θ)^2 |]
   |]

def g~i~j := M.inverse g_#_#

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

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

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

-- Ricci curvature
def Ric_i_j := withSymbols [m]
  sum (contract R~m_i_m_j)

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

-- Note: The expected scalar curvature is:
-- (6 * (`a|1|1 w) * (`a w) + 6 * (`a|1 w)^2 + 6 * K) / (`a w)^2

リンク

Egison 数学ノート目次に戻る