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 ...
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)
- Page 1: arXiv:gr-qc/0703035v1 6 Mar 2007 3+
- Page 4 and 5: 4 CONTENTS 3.3.6 Evolution of the o
- Page 6 and 7: 6 CONTENTS 8 The initial data probl
- Page 8 and 9: 8 CONTENTS
- Page 10 and 11: 10 CONTENTS
- Page 12 and 13: 12 Introduction 3D (no symmetry at
- Page 14 and 15: 14 Introduction
- Page 16 and 17: 16 Geometry of hypersurfaces in two
- Page 18 and 19: 18 Geometry of hypersurfaces 2.2.3
- Page 20 and 21: 20 Geometry of hypersurfaces Figure
- Page 22 and 23: 22 Geometry of hypersurfaces •
- Page 24 and 25: 24 Geometry of hypersurfaces The ei
- Page 26 and 27: 26 Geometry of hypersurfaces Figure
- Page 28 and 29: 28 Geometry of hypersurfaces The no
- Page 30 and 31: 30 Geometry of hypersurfaces Since
- Page 32 and 33: 32 Geometry of hypersurfaces 2.4.3
- Page 36 and 37: 36 Geometry of hypersurfaces Exampl
- Page 38 and 39: 38 Geometry of hypersurfaces
- Page 40 and 41: 40 Geometry of foliations Figure 3.
- Page 42 and 43: 42 Geometry of foliations 3.3.2 Nor
- Page 44 and 45: 44 Geometry of foliations means Eq.
- Page 46 and 47: 46 Geometry of foliations Remark :
- Page 48 and 49: 48 Geometry of foliations Note that
- Page 50 and 51: 50 Geometry of foliations
- Page 52 and 53: 52 3+1 decomposition of Einstein eq
- Page 54 and 55: 54 3+1 decomposition of Einstein eq
- Page 56 and 57: 56 3+1 decomposition of Einstein eq
- Page 58 and 59: 58 3+1 decomposition of Einstein eq
- Page 60 and 61: 60 3+1 decomposition of Einstein eq
- Page 62 and 63: 62 3+1 decomposition of Einstein eq
- Page 64 and 65: 64 3+1 decomposition of Einstein eq
- Page 66 and 67: 66 3+1 decomposition of Einstein eq
- Page 68 and 69: 68 3+1 decomposition of Einstein eq
- Page 70 and 71: 70 3+1 decomposition of Einstein eq
- Page 72 and 73: 72 3+1 equations for matter and ele
- Page 74 and 75: 74 3+1 equations for matter and ele
- Page 76 and 77: 76 3+1 equations for matter and ele
- Page 78 and 79: 78 3+1 equations for matter and ele
- Page 80 and 81: 80 3+1 equations for matter and ele
- Page 82 and 83: 82 3+1 equations for matter and ele
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)