3+1 formalism and bases of numerical relativity - LUTh ...

3+1 formalism and bases of numerical relativity - LUTh ... 3+1 formalism and bases of numerical relativity - LUTh ...

luth.obspm.fr
from luth.obspm.fr More from this publisher
20.07.2013 Views

34 Geometry of hypersurfaces 2.5 Gauss-Codazzi relations We derive here equations that will constitute the basis of the 3+1 formalism for general relativity. They are decompositions of the spacetime Riemann tensor, 4 Riem [Eq. (2.12)], in terms of quantities relative to the spacelike hypersurface Σ, namely the Riemann tensor associated with the induced metric γ, Riem [Eq. (2.34)] and the extrinsic curvature tensor of Σ, K. 2.5.1 Gauss relation Let us consider the Ricci identity (2.34) defining the (three-dimensional) Riemann tensor Riem as measuring the lack of commutation of two successive covariant derivatives with respect to the connection D associated with Σ’s metric γ. The four-dimensional version of this identity is DαDβv γ − DβDαv γ = R γ µαβ vµ , (2.84) where v is a generic vector field tangent to Σ. Let us use formula (2.79) which relates the D-derivative to the ∇-derivative, to write DαDβv γ = Dα(Dβv γ ) = γ µ αγ ν β γγ ρ∇µ(Dνv ρ ). (2.85) Using again formula (2.79) to express Dνv ρ yields DαDβv γ = γ µ α γν β γγ ρ ∇µ γ σ νγρ λ∇σv λ . (2.86) Let us expand this formula by making use of Eq. (2.64) to write ∇µγ σ ν = ∇µ (δ σ ν + nσ nν) = ∇µn σ nν + n σ ∇µnν. Since γ ν β nν = 0, we get DαDβv γ = γ µ α γν β γγ ρ n σ ∇µnν γ ρ λ∇σv λ + γ σ ν∇µn ρ nλ∇σv λ =−vλ +γ ∇σnλ σ νγρ λ λ∇µ∇σv = γ µ αγνβγγ λ∇µnν n σ ∇σv λ − γ µ αγσβγγ ρvλ∇µn ρ ∇σnλ + γ µ αγσβγγ λ λ∇µ∇σv = −Kαβ γ γ λ nσ ∇σv λ − K γ αKβλ v λ + γ µ αγ σ β γγ λ ∇µ∇σv λ , (2.87) where we have used the idempotence of the projection operator γ, i.e. γ γ ργ ρ λ = γγ λ to get the second line and γ µ αγν β∇µnν = −Kβα [Eq. (2.76)] to get the third one. When we permute the indices α and β and substract from Eq. (2.87) to form DαDβvγ − DβDγv γ , the first term vanishes since Kαβ is symmetric in (α,β). There remains DαDβv γ − DβDγv γ = v µ + γ ρ αγσβγγ λ ∇ρ∇σv λ − ∇σ∇ρv λ . (2.88) KαµK γ β − KβµK γ α Now the Ricci identity (2.13) for the connection ∇ gives ∇ρ∇σvλ −∇σ∇ρv λ = 4Rλ µρσvµ . Therefore DαDβv γ − DβDγv γ = v µ + γ ρ αγ σ βγγ 4 λ λ R µρσv µ . (2.89) KαµK γ β − KβµK γ α Substituting this relation for the left-hand side of Eq. (2.84) results in KαµK γ β − KβµK γ α v µ + γ ρ αγ σ βγγ 4 λ λ R µρσv µ = R γ µαβ vµ , (2.90)

or equivalently, since v µ = γ µ σv σ , 2.5 Gauss-Codazzi relations 35 γ µ αγ ν β γγ ργ σ λ 4 R ρ σµνv λ = R γ λαβ vλ + K γ αKλβ − K γ βKαλ v λ . (2.91) In this identity, v can be replaced by any vector of T (M) without changing the results, thanks to the presence of the projector operator γ and to the fact that both K and Riem are tangent to Σ. Therefore we conclude that γ µ αγ ν β γγ ργ σ δ 4 R ρ σµν = R γ δαβ + Kγ αKδβ − K γ β Kαδ . (2.92) This is the Gauss relation. If we contract the Gauss relation on the indices γ and α and use γ µ αγ α ρ = γ µ ρ = δ µ ρ +n µ nρ, we obtain an expression that lets appear the Ricci tensors 4 R and R associated with g and γ respectively: γ µ αγ ν β 4 Rµν + γαµn ν γ ρ βnσ 4 R µ νρσ = Rαβ + KKαβ − KαµK µ β . (2.93) We call this equation the contracted Gauss relation. Let us take its trace with respect to γ, taking into account that K µ µ = K i i = K, KµνK µν = KijK ij and γ αβ γαµn ν γ ρ β nσ 4 R µ νρσ = γ ρ µn ν n σ4 R µ νρσ = 4 R µ νµσ We obtain = 4 Rνσ n ν n σ + 4 R µ νρσn ρ nµn ν n σ = =0 4 Rµνn µ n ν . (2.94) 4 R + 2 4 Rµνn µ n ν = R + K 2 − KijK ij . (2.95) Let us call this equation the scalar Gauss relation. It constitutes a generalization of Gauss’ famous Theorema Egregium (remarkable theorem) [52, 53]. It relates the intrinsic curvature of Σ, represented by the Ricci scalar R, to its extrinsic curvature, represented by K 2 − KijK ij . Actually, the original version of Gauss’ theorem was for two-dimensional surfaces embedded in the Euclidean space R 3 . Since the curvature of the latter is zero, the left-hand side of Eq. (2.95) vanishes identically in this case. Moreover, the metric g of the Euclidean space R 3 is Riemannian, not Lorentzian. Consequently the term K 2 − KijK ij has the opposite sign, so that Eq. (2.95) becomes R − K 2 + KijK ij = 0 (g Euclidean). (2.96) This change of sign stems from the fact that for a Riemannian ambient metric, the unit normal vector n is spacelike and the orthogonal projector is γ α β = δα β −nα nβ instead of γ α β = δα β +nα nβ [the latter form has been used explicitly in the calculation leading to Eq. (2.87)]. Moreover, in dimension 2, formula (2.96) can be simplified by letting appear the principal curvatures κ1 and κ2 of Σ (cf. Sec. 2.3.4). Indeed, K can be diagonalized in an orthonormal basis (with respect to γ) so that Kij = diag(κ1,κ2) and Kij = diag(κ1,κ2). Consequently, K = κ1 + κ2 and KijK ij = κ2 1 + κ22 and Eq. (2.96) becomes R = 2κ1κ2 (g Euclidean, Σ dimension 2). (2.97)

34 Geometry <strong>of</strong> hypersurfaces<br />

2.5 Gauss-Codazzi relations<br />

We derive here equations that will constitute the basis <strong>of</strong> the <strong>3+1</strong> <strong>formalism</strong> for general <strong>relativity</strong>.<br />

They are decompositions <strong>of</strong> the spacetime Riemann tensor, 4 Riem [Eq. (2.12)], in terms <strong>of</strong><br />

quantities relative to the spacelike hypersurface Σ, namely the Riemann tensor associated with<br />

the induced metric γ, Riem [Eq. (2.34)] <strong>and</strong> the extrinsic curvature tensor <strong>of</strong> Σ, K.<br />

2.5.1 Gauss relation<br />

Let us consider the Ricci identity (2.34) defining the (three-dimensional) Riemann tensor Riem<br />

as measuring the lack <strong>of</strong> commutation <strong>of</strong> two successive covariant derivatives with respect to the<br />

connection D associated with Σ’s metric γ. The four-dimensional version <strong>of</strong> this identity is<br />

DαDβv γ − DβDαv γ = R γ<br />

µαβ vµ , (2.84)<br />

where v is a generic vector field tangent to Σ. Let us use formula (2.79) which relates the<br />

D-derivative to the ∇-derivative, to write<br />

DαDβv γ = Dα(Dβv γ ) = γ µ αγ ν β γγ ρ∇µ(Dνv ρ ). (2.85)<br />

Using again formula (2.79) to express Dνv ρ yields<br />

DαDβv γ = γ µ α γν β γγ ρ ∇µ<br />

<br />

γ σ νγρ λ∇σv λ<br />

. (2.86)<br />

Let us exp<strong>and</strong> this formula by making use <strong>of</strong> Eq. (2.64) to write ∇µγ σ ν = ∇µ (δ σ ν + nσ nν) =<br />

∇µn σ nν + n σ ∇µnν. Since γ ν β nν = 0, we get<br />

DαDβv γ = γ µ α γν β γγ ρ<br />

<br />

n σ ∇µnν γ ρ<br />

λ∇σv λ + γ σ ν∇µn ρ nλ∇σv λ<br />

<br />

=−vλ +γ<br />

∇σnλ<br />

σ νγρ <br />

λ<br />

λ∇µ∇σv = γ µ αγνβγγ λ∇µnν n σ ∇σv λ − γ µ αγσβγγ ρvλ∇µn ρ ∇σnλ + γ µ αγσβγγ λ<br />

λ∇µ∇σv = −Kαβ γ γ<br />

λ nσ ∇σv λ − K γ αKβλ v λ + γ µ αγ σ β γγ<br />

λ ∇µ∇σv λ , (2.87)<br />

where we have used the idempotence <strong>of</strong> the projection operator γ, i.e. γ γ ργ ρ<br />

λ = γγ<br />

λ to get<br />

the second line <strong>and</strong> γ µ αγν β∇µnν = −Kβα [Eq. (2.76)] to get the third one. When we permute<br />

the indices α <strong>and</strong> β <strong>and</strong> substract from Eq. (2.87) to form DαDβvγ − DβDγv γ , the first term<br />

vanishes since Kαβ is symmetric in (α,β). There remains<br />

DαDβv γ − DβDγv γ <br />

=<br />

v µ + γ ρ αγσβγγ <br />

λ ∇ρ∇σv λ − ∇σ∇ρv λ<br />

. (2.88)<br />

<br />

KαµK γ<br />

β − KβµK γ α<br />

Now the Ricci identity (2.13) for the connection ∇ gives ∇ρ∇σvλ −∇σ∇ρv λ = 4Rλ µρσvµ . Therefore<br />

DαDβv γ − DβDγv γ =<br />

<br />

v µ + γ ρ αγ σ βγγ 4 λ<br />

λ R µρσv µ . (2.89)<br />

<br />

KαµK γ<br />

β − KβµK γ α<br />

Substituting this relation for the left-h<strong>and</strong> side <strong>of</strong> Eq. (2.84) results in<br />

<br />

KαµK γ<br />

β − KβµK γ α<br />

<br />

v µ + γ ρ αγ σ βγγ 4 λ<br />

λ R µρσv µ = R γ<br />

µαβ vµ , (2.90)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!