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

122 Asymptotic flatness and global quantities to the hypersurfaces Σt. Then n · φ = 0 and the integrand in the definition (7.100) is Accordingly Eq. (7.100) becomes ∇ µ φ ν dSµν = ∇µφν(s µ n ν − n µ s ν ) √ q d 2 y = 2∇µφν s µ n ν√ q d 2 y = −2s µ φν∇µn ν√ q d 2 y = 2Kijs i φ j√ q d 2 y. (7.103) JK = 1 Kijs 8π St i φ j√ q d 2 y . (7.104) Remark : Contrary to the 3+1 expression of the Komar mass which turned out to be very different from the expression of the ADM mass, the 3+1 expression of the Komar angular momentum as given by Eq. (7.104) is very similar to the expression of the angular momentum deduced from the Hamiltonian formalism, i.e. Eq. (7.63). The only differences are that it is no longer necessary to take the limit St → ∞ and that there is no trace term Kγijs i φ j in Eq. (7.104). Moreover, if one evaluates the Hamiltonian expression in the asymptotically maximal gauge (7.65) then K = O(r −3 ) and thanks to the asymptotic orthogonality of s and φ, γijs i φ j = O(1), so that Kγijs i φ j does not contribute to the integral and expressions (7.104) and (7.63) are then identical. Example : A trivial example is provided by Schwarzschild spacetime, which among other things is axisymmetric. For the 3+1 foliation associated with the Schwarzschild coordinates (7.16), the extrinsic curvature tensor K vanishes identically, so that Eq. (7.104) yields immediately JK = 0. For other foliations, like that associated with Eddington-Finkelstein coordinates, K is no longer zero but is such that Kijs i φ j = 0, yielding again JK = 0 (as it should be since the Komar angular momentum is independent of the foliation). Explicitely for Eddington-Finkelstein coordinates, Kijs i = − 2m r 2 1 + m r 1 + 2m r , 0, 0 , (7.105) (see e.g. Eq. (D.25) in Ref. [146]) and φ j = (0,0,1), so that obviously Kijs i φ j = 0. Example : The most natural non trivial example is certainly that of Kerr spacetime. Let us use the 3+1 foliation associated with the standard Boyer-Lindquist coordinates (t,r,θ,ϕ) and evaluate the integral (7.104) by choosing for St a sphere r = const. Then y a = (θ,ϕ). The Boyer-Lindquist components of φ are φ i = (0,0,1) and those of s are s i = (s r ,0,0) since γij is diagonal is these coordinates. The formula (7.104) then reduces to JK = 1 Krϕs 8π r=const r√ q dθ dϕ. (7.106)

7.6 Komar mass and angular momentum 123 The extrinsic curvature component Krϕ can be evaluated via formula (4.63), which reduces to 2NKij = Lβ γij since ∂γij/∂t = 0. From the Boyer-Lindquist line element (see e.g. Eq. (5.29) in Ref. [156]), we read the components of the shift: (β r ,β θ ,β ϕ 2amr ) = 0, 0, − (r2 + a2 )(r2 + a2 cos2 θ) + 2a2mr sin2 , (7.107) θ where m and a are the two parameters of the Kerr solution. Then, using Eq. (A.6), Krϕ = 1 2N Lβ γrϕ = 1 ϕ ∂γrϕ ∂β β +γϕϕ 2N ∂ϕ =0 ϕ ∂β + γrϕ ∂r ϕ = ∂ϕ =0 1 2N γϕϕ ∂βϕ . ∂r (7.108) Hence JK = 1 s 16π r=const r N γϕϕ ∂βϕ √ q dθ dϕ. (7.109) ∂r The values of s r , N, γϕϕ and √ q can all be read on the Boyer-Lindquist line element. However this is a bit tedious. To simplify things, let us evaluate JK only in the limit r → ∞. Then s r ∼ 1, N ∼ 1, γϕϕ ∼ r 2 sin 2 θ, √ q ∼ r 2 sin θ and, from Eq. (7.107), βϕ ∼ −2am/r3 , so that JK = 1 16π r 2 2 6am sin r4 r2 sin θ dθ dϕ = 3am π × 2π × sin 8π 0 3 θ dθ. (7.110) Hence, as expected, r=const JK = am. (7.111) Let us now find the 3+1 expression of the volume version (7.101) of the Komar angular momentum. We have n · φ = 0 and, from the 3+1 decomposition (4.10) of T: Hence formula (7.101) becomes with JK = T(n,φ) = −〈p,φ〉. (7.112) Σt 〈p,φ〉 √ γ d 3 x + J H K , (7.113) J H K = 1 Kijs 8π Ht i φ j√ q d 2 y . (7.114) Example : Let us consider a perfect fluid. Then p = (E + P)U [Eq. (5.61)], so that JK = (E + P)U · φ √ γ d 3 x + J H K . (7.115) Σt Taking φ = −y∂x +x∂y (symmetry axis = z-axis), the Newtonian limit of this expression is then JK = ρ(−yU x + xU y )dxdy dz, (7.116) Σt i.e. we recognize the standard expression for the angular momentum around the z-axis.

7.6 Komar mass <strong>and</strong> angular momentum 123<br />

The extrinsic curvature component Krϕ can be evaluated via formula (4.63), which reduces<br />

to 2NKij = Lβ γij since ∂γij/∂t = 0. From the Boyer-Lindquist line element (see e.g.<br />

Eq. (5.29) in Ref. [156]), we read the components <strong>of</strong> the shift:<br />

(β r ,β θ ,β ϕ <br />

2amr<br />

) = 0, 0, −<br />

(r2 + a2 )(r2 + a2 cos2 θ) + 2a2mr sin2 <br />

, (7.107)<br />

θ<br />

where m <strong>and</strong> a are the two parameters <strong>of</strong> the Kerr solution. Then, using Eq. (A.6),<br />

Krϕ = 1<br />

2N Lβ γrϕ = 1<br />

<br />

ϕ ∂γrϕ ∂β<br />

β +γϕϕ<br />

2N ∂ϕ<br />

<br />

=0<br />

ϕ ∂β<br />

+ γrϕ<br />

∂r ϕ <br />

=<br />

∂ϕ<br />

<br />

=0<br />

1<br />

2N γϕϕ<br />

∂βϕ .<br />

∂r<br />

(7.108)<br />

Hence<br />

JK = 1<br />

<br />

s<br />

16π r=const<br />

r<br />

N γϕϕ<br />

∂βϕ √<br />

q dθ dϕ. (7.109)<br />

∂r<br />

The values <strong>of</strong> s r , N, γϕϕ <strong>and</strong> √ q can all be read on the Boyer-Lindquist line element.<br />

However this is a bit tedious. To simplify things, let us evaluate JK only in the limit<br />

r → ∞. Then s r ∼ 1, N ∼ 1, γϕϕ ∼ r 2 sin 2 θ, √ q ∼ r 2 sin θ <strong>and</strong>, from Eq. (7.107),<br />

βϕ ∼ −2am/r3 , so that<br />

JK = 1<br />

<br />

16π<br />

r 2 2 6am<br />

sin<br />

r4 r2 sin θ dθ dϕ = 3am<br />

π<br />

× 2π × sin<br />

8π 0<br />

3 θ dθ. (7.110)<br />

Hence, as expected,<br />

r=const<br />

JK = am. (7.111)<br />

Let us now find the <strong>3+1</strong> expression <strong>of</strong> the volume version (7.101) <strong>of</strong> the Komar angular<br />

momentum. We have n · φ = 0 <strong>and</strong>, from the <strong>3+1</strong> decomposition (4.10) <strong>of</strong> T:<br />

Hence formula (7.101) becomes<br />

with<br />

<br />

JK =<br />

T(n,φ) = −〈p,φ〉. (7.112)<br />

Σt<br />

〈p,φ〉 √ γ d 3 x + J H K , (7.113)<br />

J H K = 1<br />

<br />

Kijs<br />

8π Ht<br />

i φ j√ q d 2 y . (7.114)<br />

Example : Let us consider a perfect fluid. Then p = (E + P)U [Eq. (5.61)], so that<br />

<br />

JK = (E + P)U · φ √ γ d 3 x + J H K . (7.115)<br />

Σt<br />

Taking φ = −y∂x +x∂y (symmetry axis = z-axis), the Newtonian limit <strong>of</strong> this expression<br />

is then<br />

<br />

JK = ρ(−yU x + xU y )dxdy dz, (7.116)<br />

Σt<br />

i.e. we recognize the st<strong>and</strong>ard expression for the angular momentum around the z-axis.

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

Saved successfully!

Ooh no, something went wrong!