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

58 3+1 decomposition of Einstein equation Besides, since ∂i and ∂j are tangent to Σt, gij = g(∂i,∂j) = γ(∂i,∂j) = γij. (4.46) Collecting Eqs. (4.43), (4.45) and (4.46), we get the following expression of the metric components in terms of 3+1 quantities: gαβ = g00 g0j gi0 gij or, in terms of line elements [using Eq. (4.40)], −N2 + βkβ = k βj βi γij , (4.47) gµν dx µ dx ν = −N 2 dt 2 + γij(dx i + β i dt)(dx j + β j dt) . (4.48) The components of the inverse metric are given by the matrix inverse of (4.47): g αβ g00 g = 0j = g i0 g ij − 1 N2 βi N2 βj N2 γij − βiβj N2 . (4.49) Indeed, it is easily checked that the matrix product g αµ gµβ is equal to the identity matrix δ α β . Remark : Notice that gij = γij but that in general g ij = γ ij . One can deduce from the above formulæ a simple relation between the determinants of g and γ. Let us first define the latter ones by g := det(gαβ) , (4.50) γ := det(γij) . (4.51) Notice that g and γ depend upon the choice of the coordinates (x α ). They are not scalar quantities, but scalar densities. Using Cramer’s rule for expressing the inverse (g αβ ) of the matrix (gαβ), we have g 00 = C00 det(gαβ) C00 = , (4.52) g where C00 is the element (0,0) of the cofactor matrix associated with (gαβ). It is given by C00 = (−1) 0 M00 = M00, where M00 is the minor (0,0) of the matrix (gαβ), i.e. the determinant of the 3 × 3 matrix deduced from (gαβ) by suppressing the first line and the first column. From Eq. (4.47), we read M00 = det(γij) = γ. (4.53) Hence Eq. (4.52) becomes g 00 = γ . (4.54) g Expressing g 00 from Eq. (4.49) yields then g = −N 2 γ, or equivalently, √ −g = N √ γ . (4.55)

4.3 3+1 Einstein equation as a PDE system 59 4.2.4 Choice of coordinates via the lapse and the shift We have seen above that giving a coordinate system (x α ) on M such that the hypersurfaces x 0 = const. are spacelike determines uniquely a lapse function N and a shift vector β. The converse is true in the following sense: setting on some hypersurface Σ0 a scalar field N, a vector field β and a coordinate system (x i ) uniquely specifies a coordinate system (x α ) in some neighbourhood of Σ0, such that the hypersurface x 0 = 0 is Σ0. Indeed, the knowledge of the lapse function a each point of Σ0 determines a unique vector m = Nn and consequently the location of the “next” hypersurface Σδt by Lie transport along m (cf. Sec. 3.3.2). Graphically, we may also say that for each point of Σ0 the lapse function specifies how far is the point of Σδt located “above” it (“above” meaning perpendicularly to Σ0, cf. Fig. 3.2). Then the shift vector tells how to propagate the coordinate system (x i ) from Σ0 to Σδt (cf. Fig. 4.1). This way of choosing coordinates via the lapse function and the shift vector is one of the main topics in 3+1 numerical relativity and will be discussed in detail in Chap. 9. 4.3 3+1 Einstein equation as a PDE system 4.3.1 Lie derivatives along m as partial derivatives Let us consider the term LmK which occurs in the 3+1 Einstein equation (4.15). Thanks to Eq. (4.28), we can write LmK = L∂t K − Lβ K. (4.56) This implies that L∂t K is a tensor field tangent to Σt, since both LmK and Lβ K are tangent to Σt, the former by the property (3.32) and the latter because β and K are tangent to Σt. Moreover, if one uses tensor components with respect to a coordinate system (x α ) = (t,x i ) adapted to the foliation, the Lie derivative along ∂t reduces simply to the partial derivative with respect to t [cf. Eq. (A.3)]: L∂tKij = ∂Kij . (4.57) ∂t By means of formula (A.6), one can also express Lβ K in terms of partial derivatives: with k ∂Kij ∂β Lβ Kij = β + Kkj ∂xk k ∂β + Kik ∂xi k . (4.58) ∂xj Similarly, the relation (3.24) between Lmγ and K becomes L∂t γ − Lβ γ = −2NK, (4.59) L∂tγij = ∂γij . (4.60) ∂t and, evaluating the Lie derivative with the connection D instead of partial derivatives [cf. Eq. (A.8)]: i.e. Lβ γij = β k Dkγij +γkjDiβ k + γikDjβ k , (4.61) =0 Lβ γij = Diβj + Djβi. (4.62)

58 <strong>3+1</strong> decomposition <strong>of</strong> Einstein equation<br />

Besides, since ∂i <strong>and</strong> ∂j are tangent to Σt,<br />

gij = g(∂i,∂j) = γ(∂i,∂j) = γij. (4.46)<br />

Collecting Eqs. (4.43), (4.45) <strong>and</strong> (4.46), we get the following expression <strong>of</strong> the metric components<br />

in terms <strong>of</strong> <strong>3+1</strong> quantities:<br />

gαβ =<br />

g00 g0j<br />

gi0 gij<br />

or, in terms <strong>of</strong> line elements [using Eq. (4.40)],<br />

<br />

−N2 + βkβ<br />

=<br />

k βj<br />

βi<br />

γij<br />

<br />

, (4.47)<br />

gµν dx µ dx ν = −N 2 dt 2 + γij(dx i + β i dt)(dx j + β j dt) . (4.48)<br />

The components <strong>of</strong> the inverse metric are given by the matrix inverse <strong>of</strong> (4.47):<br />

g αβ <br />

g00 g<br />

=<br />

0j<br />

<br />

=<br />

g i0 g ij<br />

− 1<br />

N2 βi N2 βj N2 γij − βiβj N2 <br />

. (4.49)<br />

Indeed, it is easily checked that the matrix product g αµ gµβ is equal to the identity matrix δ α β .<br />

Remark : Notice that gij = γij but that in general g ij = γ ij .<br />

One can deduce from the above formulæ a simple relation between the determinants <strong>of</strong> g<br />

<strong>and</strong> γ. Let us first define the latter ones by<br />

g := det(gαβ) , (4.50)<br />

γ := det(γij) . (4.51)<br />

Notice that g <strong>and</strong> γ depend upon the choice <strong>of</strong> the coordinates (x α ). They are not scalar<br />

quantities, but scalar densities. Using Cramer’s rule for expressing the inverse (g αβ ) <strong>of</strong> the<br />

matrix (gαβ), we have<br />

g 00 = C00<br />

det(gαβ)<br />

C00<br />

= , (4.52)<br />

g<br />

where C00 is the element (0,0) <strong>of</strong> the c<strong>of</strong>actor matrix associated with (gαβ). It is given by<br />

C00 = (−1) 0 M00 = M00, where M00 is the minor (0,0) <strong>of</strong> the matrix (gαβ), i.e. the determinant<br />

<strong>of</strong> the 3 × 3 matrix deduced from (gαβ) by suppressing the first line <strong>and</strong> the first column. From<br />

Eq. (4.47), we read<br />

M00 = det(γij) = γ. (4.53)<br />

Hence Eq. (4.52) becomes<br />

g 00 = γ<br />

. (4.54)<br />

g<br />

Expressing g 00 from Eq. (4.49) yields then g = −N 2 γ, or equivalently,<br />

√ −g = N √ γ . (4.55)

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

Saved successfully!

Ooh no, something went wrong!