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 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
7.2.2 Asymptotic coordinate freedom<br />
7.3 ADM mass 105<br />
Obviously the above definition <strong>of</strong> asymptotic flatness depends both on the foliation (Σt)t∈R <strong>and</strong><br />
on the coordinates (x i ) chosen on each leaf Σt. It is <strong>of</strong> course important to assess whether this<br />
dependence is strong or not. In other words, we would like to determine the class <strong>of</strong> coordinate<br />
changes (x α ) = (t,x i ) → (x ′α ) = (t ′ ,x ′i ) which preserve the asymptotic properties (7.1)-(7.4).<br />
The answer is that the coordinates (x ′α ) must be related to the coordinates (x α ) by [157]<br />
x ′α = Λ α µx µ + c α (θ,ϕ) + O(r −1 ) (7.7)<br />
where Λ α β is a Lorentz matrix <strong>and</strong> the cα ’s are four functions <strong>of</strong> the angles (θ,ϕ) related to the<br />
coordinates (x i ) = (x,y,z) by the st<strong>and</strong>ard formulæ:<br />
x = r sin θ cos ϕ, y = r sin θ sin ϕ, z = r cos θ. (7.8)<br />
The group <strong>of</strong> transformations generated by (7.7) is related to the Spi group (for Spatial infinity)<br />
introduced by Ashtekar <strong>and</strong> Hansen [25, 24]. However the precise relation is not clear because<br />
the definition <strong>of</strong> asymptotic flatness used by these authors is not expressed as decay conditions<br />
for γij <strong>and</strong> Kij, as in Eqs. (7.1)-(7.4).<br />
Notice that Poincaré transformations are contained in transformation group defined<br />
by (7.7): they simply correspond to the case cα (θ,ϕ) = const. The transformations with<br />
cα (θ,ϕ) = const <strong>and</strong> Λα β = δα β constitute “angle-dependent translations” <strong>and</strong> are called supertranslations.<br />
involves a boost, the transformation (7.7) implies a<br />
Note that if the Lorentz matrix Λα β<br />
change <strong>of</strong> the <strong>3+1</strong> foliation (Σt)t∈R, whereas if Λα β corresponds only to some spatial rotation<br />
<strong>and</strong> the cα ’s are constant, the transformation (7.7) describes some change <strong>of</strong> Cartesian-type<br />
coordinates (xi ) (rotation + translation) within the same hypersurface Σt.<br />
7.3 ADM mass<br />
7.3.1 Definition from the Hamiltonian formulation <strong>of</strong> GR<br />
In the short introduction to the Hamiltonian formulation <strong>of</strong> general <strong>relativity</strong> given in Sec. 4.5,<br />
we have for simplicity discarded any boundary term in the action. However, because the gravitational<br />
Lagrangian density (the scalar curvature 4R) contains second order derivatives <strong>of</strong> the<br />
metric tensor (<strong>and</strong> not only first order ones, which is a particularity <strong>of</strong> general <strong>relativity</strong> with<br />
respect to other field theories), the precise action should be [209, 205, 265, 157]<br />
<br />
<br />
4 √ 4<br />
S = R −g d x + 2 (Y − Y0) √ h d 3 y, (7.9)<br />
V<br />
where ∂V is the boundary <strong>of</strong> the domain V (∂V is assumed to be a timelike hypersurface), Y the<br />
trace <strong>of</strong> the extrinsic curvature (i.e. three times the mean curvature) <strong>of</strong> ∂V embedded in (M,g)<br />
<strong>and</strong> Y0 the trace <strong>of</strong> the extrinsic curvature <strong>of</strong> ∂V embedded in (M,η), where η is a Lorentzian<br />
metric on M which is flat in the region <strong>of</strong> ∂V. Finally √ hd 3 y is the volume element induced<br />
by g on the hypersurface ∂V, h being the induced metric on ∂V <strong>and</strong> h its determinant with<br />
∂V