Riemann Curvature Tensor of $S^7$¶

A round 7-sphere of radius $r$ has constant sectional curvature $1/r^2$. This notebook builds its metric, Levi-Civita connection, Riemann tensor, Ricci tensor, and scalar curvature in indexed Egison notation.

We use the convention

$$ R^i{}_{jkl} = \partial_k\Gamma^i{}_{jl} - \partial_l\Gamma^i{}_{jk} + \Gamma^m{}_{jl}\Gamma^i{}_{mk} - \Gamma^m{}_{jk}\Gamma^i{}_{ml}. $$

Hyperspherical chart¶

The coordinates are $x=(\alpha, \beta, \gamma, \delta, \xi, \zeta, \eta)$. The standard embedding $X:S^7\hookrightarrow\mathbb{R}^8$ is built by successively multiplying by sines: $X_1=r\cos \alpha$, $X_2=r\sin \alpha\cos \beta$, and so on. It makes $X\mathbin{\cdot}X=r^2$ manifest.

The chart excludes the usual coordinate poles; those singularities are features of hyperspherical coordinates, not of the round geometry.

declare symbol r, α, β, γ, δ, ξ, ζ, η: MathValue

def x : Vector MathValue := [| α, β, γ, δ, ξ, ζ, η |]

def X : Vector MathValue := [| r * cos α, r * sin α * cos β, r * sin α * sin β * cos γ, r * sin α * sin β * sin γ * cos δ, r * sin α * sin β * sin γ * sin δ * cos ξ, r * sin α * sin β * sin γ * sin δ * sin ξ * cos ζ, r * sin α * sin β * sin γ * sin δ * sin ξ * sin ζ * cos η, r * sin α * sin β * sin γ * sin δ * sin ξ * sin ζ * sin η |]

$\begin{pmatrix} \cos(α) r \\ \cos(β) \sin(α) r \\ \cos(γ) \sin(α) \sin(β) r \\ \cos(δ) \sin(α) \sin(β) \sin(γ) r \\ \cos(ξ) \sin(α) \sin(β) \sin(γ) \sin(δ) r \\ \cos(ζ) \sin(α) \sin(β) \sin(γ) \sin(δ) \sin(ξ) r \\ \cos(η) \sin(α) \sin(β) \sin(γ) \sin(δ) \sin(ζ) \sin(ξ) r \\ \sin(α) \sin(β) \sin(γ) \sin(δ) \sin(ζ) \sin(η) \sin(ξ) r\\ \end{pmatrix}$

Metric and inverse metric¶

Differentiating the embedding gives an orthogonal coordinate basis. Its line element is

$$ ds^2=r^2\left(d\alpha^2 + \sin^2 \alpha d\beta^2 + \sin^2 \alpha \sin^2 \beta d\gamma^2 + \sin^2 \alpha \sin^2 \beta \sin^2 \gamma d\delta^2 + \sin^2 \alpha \sin^2 \beta \sin^2 \gamma \sin^2 \delta d\xi^2 + \sin^2 \alpha \sin^2 \beta \sin^2 \gamma \sin^2 \delta \sin^2 \xi d\zeta^2 + \sin^2 \alpha \sin^2 \beta \sin^2 \gamma \sin^2 \delta \sin^2 \xi \sin^2 \zeta d\eta^2\right). $$

Egison differentiates the embedding to obtain the coordinate tangent vectors and constructs every component as their dot product, $g_{ij}=\partial_iX\mathbin{\cdot}\partial_jX$. The inverse metric is then computed from that induced metric.

def e_i_j : Matrix MathValue := ∂/∂ X_j x~i

def g_i_j : Matrix MathValue :=
  generateTensor (\[a, b] -> V.* e_a_# e_b_#) [7, 7]

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

g_1_#

$\begin{pmatrix} r^{2} \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0\\ \end{pmatrix}_{\#}^{\;}$

g_7_#

$\begin{pmatrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ \sin(α)^{2} \sin(β)^{2} \sin(γ)^{2} \sin(δ)^{2} \sin(ζ)^{2} \sin(ξ)^{2} r^{2}\\ \end{pmatrix}_{\#}^{\;}$

g~1~#

$\begin{pmatrix} r^{-2} \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0\\ \end{pmatrix}_{\;}^{\#}$

Levi-Civita connection¶

The Christoffel symbols of the first kind and second kind are

$$ \Gamma_{ijk}=\frac12 (\partial_jg_{ik}+\partial_kg_{ij}-\partial_ig_{jk}), \qquad \Gamma^i{}_{jk}=g^{im}\Gamma_{mjk}. $$

Repeated symbolic indices are contracted by .. The withSymbols block makes the dummy index local to the definition.

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

$-\cos(α) \sin(α)$

Γ~2_1_2

$\cos(α) \sin(α)^{-1}$

Riemann tensor¶

The following definition is a direct transcription of the stated convention. The two output cells sample the same coordinate two-plane with the first two tensor slots exchanged; their different coordinate factors are exactly what one expects in a non-orthonormal coordinate basis.

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

R~1_2_1_2

$\sin(α)^{2}$

R~2_1_1_2

$-1$

Ricci and scalar curvature¶

Contracting the first and third Riemann indices gives

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

For a round $S^7$ the coordinate-free prediction is

$$ \operatorname{Ric}=\frac{6}{r^2}g, \qquad \mathcal{R}=\frac{42}{r^2}. $$

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

def scalarCurvature : MathValue := withSymbols [i, j]
  g~i~j . Ric_i_j

Interpretation¶

The sampled components are coordinate dependent, but their contractions recover the invariant statement: every tangent two-plane has sectional curvature $1/r^2$, the metric is Einstein, and the scalar curvature is $42/r^2$. Factors such as $\sin \alpha$ vanish at chart poles because the coordinate frame degenerates there; the curvature itself remains smooth. The Ricci and scalar-curvature definitions are present, but their full contractions are intentionally not output cells: the selected Riemann components demonstrate the constant-curvature pattern without forcing a large symbolic expansion.