General Relativity Fall 2018
Lecture 10: Integration, Einstein-Hilbert action
Yacine Ali-Ha¨ımoud
(Dated: October 3, 2018)
σ
HW comment: Tσ antisym in µν does NOT imply that Tµν is antisym in µν.
µν
- Volu0me element – Consider a LICS with primed coordinates. The 4-volume element is dV = d4x0
- =
- 0
- 0
- 0
dx0 dx1 dx2 dx3 . If we change coordinates, we have
ꢀꢀꢀ
ꢀꢀꢀ
- !
0
∂xµ ∂xµ
d4x0 = det
d4x.
(1)
ꢀꢀ
ꢀꢀ
Now, the metric components change as
- 0
- 0
- 0
- 0
- ∂xµ ∂xν
- ∂xµ ∂xν
0
gµ ν
- 0
- 0
ηµ ν
0
gµν
- =
- =
,
(2)
- ∂xµ ∂xν
- ∂xµ ∂xν
- 0
- 0
since the metric components are ηµ ν in the LICS. Seing this as a matrix operation and taking the determinant, we find
- !
2
0
∂xµ ∂xµ
det(gµν) = − det
.
(3)
Hence, we find that the 4-volume element is
qp
√
dV = − det(gµν)d4x ≡ −g d4x ≡ |g| d4x.
(4) (5)
The integral of a scalar function f is well defined: given any coordinate system (even if only defined locally):
- Z
- Z
p
- f dV =
- f
- |g| d4x.
We can only define the integral of a scalar function. The integral of a vector or tensor field is mean-
ingless in curved spacetime. Think of the integral as a sum. To sum vectors, you need them to belong to the same vector space. There is no common vector space in curved spacetime.
Only in flat spacetime can we define such integrals. First parallel-transport the vector field to a single point of spacetime (it doesn’t matter which one). This operation is independent of the path in flat spacetime, hence well defined. Then sum these components to perform the integral.
Covariant divergence– The covariant divergence of a vector field is
∇µV µ = ∂µV µ + ΓµµνV ν.
(6)
The contracted Christoffel symbol is
- 1
- 1
- 1
2
Γµµν
=
gµλ (∂νgµλ + ∂µgνλ − ∂λgµν) = gµλ∂νgµλ + gµλ∂[λgµ]ν
=
gµλ∂νgµλ
.
(7)
- 2
- 2
Let us now recall that the determinant and inverse of a matrix (here, the metric gµν) can be expressed in terms of the comatrix Cµν, which is the matrix made of the determinants of the (n − 1) × (n − 1) submatrices obtained by deleting row µ and column ν (and then multiplying by (−1)µ+ν)):
X
g ≡ det(gµν) =
gµνCµν
,
- for any fixed ν (not summed over!),
- (8)
µ
1
gµν
=
Cµν
,
(9)
g
2
Therefore, we get
So, we find
p
∂
|g|
- 1 1 ∂g
- 1 1
- 1
2
ν
=
∂νgµλ
=
Cµλ∂νgµλ
=
gµλ∂νgµλ = Γµµν
.
(10)
p
2 g ∂gµλ
2 g
|g|
pp
∂
|g|
V ν
1
ν
∇µV µ = ∂µV µ
- +
- =
∂µ( |g|V µ).
(11)
- p
- p
- |g|
- |g|
Stokes’ theorem – First consider a 4-volume V covered by a single coordinate system. Since ∇µV µ is a scalar field, we may compute its integral over V:
p
- Z
- Z
- Z
∂µ( |g|V µ)
- p
- p
∇µV µ dV =
- |g|d4x =
- ∂µ( |g|V µ)d4x.
(12)
p
|g|
- V
- V
- V
This is just a normal integral in a volume of R4, for which we can use Stokes’ theorem. We denote by ∂V the boundary of V. We denote by nµ the vector normal to the surface, with squared norm ±1, pointing outwards if it is spacelike, and inwards if it is timelike. You can check that this is is required to get back Stokes’ theorem – to convince yourself, use coordinates in which the normal to the surface is a coordinate basis vector ∂(µ) and use the usual Stokes theorem.
By definition, this vector is orthogonal to the vectors spanning the tangent space of ∂V, so the metric can be
- (3)
- (3)
µν
rewritten as gµν = ±n n + g , where g nµ = 0: g(3) is the induced metric on the 3-dimensional manifold ∂V. We
µν
- µ
- ν
(3)
µν
then have det(gµν) = ± det(g ), i.e. |g| = |g(3)|. We then get
- Z
- Z
- Z
qp
∇µV µ dV =
- nµ |g|V µd3x =
- nµV µ |g(3)|d3x.
(13)
- V
- ∂V
- ∂V
We recognize the last two factors as the frame-invariant volume element on the boundary. This last integral is therefore all coordinate-independent, as it should.
Now, if we need to use multiple charts to cover V, we can subdivide it into sub-volumes each with a single coordinate system, and apply Stoke’s theorem there. The integrals on interior surfaces cancel out, and we are left with the boundary of the full volume.
Action formulation of general relativity – The Einstein field equation Gµν = 8πTµν (from now we use
geometric units where G = c = 1) can be derived from an action principle. We want to build an overall action of the form
- Z
- Z
- p
- p
S = Sgravity + SM
≡
d4x |g|Lgravity
+
- d4x |g|LM
- ,
(14) where SM is the matter action, and Sgravity only involves the metric tensor, and L are the corresponding frameinvariant Lagrangian densities. The most natural scalar function of curvature that comes to mind is the Ricci scalar, so we may anticipate that Lgravity ∝ R. The idea then, is to recover the Eintein field equation by extremizing the action with respect to variations in the metric tensor and in the matter fields (whatever they are).
To compute these variations, pick a coordinate system. We want to compute the variation of the action upon
- variation of the metric field by δgµν
- :
- !
p
- ꢁ
- ꢂ
- p
- p
- p
- p
√
- δ( |g|)
- 1
2
gµνδgµν + δR
p
δ( −gR) = Rδ( |g|) + |g| δR = |g|
- R
- + δR
=
|g|
|g|
- ꢁ
- ꢂ
p
1
=
|g| − RgµνδgµνR + δR
,
(15)
2
where we used δ(gµνgµν) = 0 = gµνδgµν + gµνδgµν. Now we have
R = gµνRµν ⇒ δR = Rµνδgµν + gµνδRµν
.
(16) (17)
This implies
p
√
δ( −gR) = |g| (Gµνδgµν + gµνδRµν) .
3
Let us now compute the variation of the Riemann tensor under variation of the metric. To do so, recall that the Riemann tensor is defined in terms of the commutator of the metric-compatible covariant derivative,
RλσµνV σ ≡ ∇µ∇νV λ − ∇ν∇µV λ.
(18)
˜
- Let us denote by ∇ the unique metric-compatible (and torsion-free) covariant derivative associated with gµν + δgµν
- ,
we then have
- σ
- λ
δRλσµνV = ∇µ∇νV − ∇µ∇νV λ − (µ ↔ ν).
˜ ˜
(19)
λ
- ˜
- ˜
We denote by Γµν the associated Christoffel symbols, and δΓ = Γ − Γ the difference relative to those associated with
gµν. We then have, for any vector,
- λ
- λ
∇νV = ∇νV λ + δΓνσV σ,
(20)
˜where the partial derivatives cancel out. Because the left-hand side is a tensor and so is V σ, we conclude that δΓ is a rank (1, 2) tensor (even if each coefficient individually is not!), symmetric in its lower two indices. Similarly, we have for any rank (1, 1) tensor,
- λ
- σ
- σ
∇µX = ∇µXλ + δΓλµσX − δΓµνXλ
(21)
˜
- ν
- ν
- ν
- σ
We then have
- σ
- λ
- λ
- σ
δRλσµνV = ∇µ∇νV − ∇µ∇νV λ + ∇µ(δΓνσV ) − (µ ↔ ν)
- ˜
- ˜
= δΓλ ∇νV σ − δΓµσν∇σV λ + ∇µ(δΓλνσV σ) + O(δΓ2) − (µ ↔ ν)
µσ µσ
= δΓλ ∇νV σ − δΓµσν∇σV λ + ∇µ(δΓλνσ)V σ + δΓνλσ∇µV σ + O(δΓ2) − (µ ↔ ν)
(22)
The second term is symmetric in (µ, nu) so cancels upon antisymmetrization. The sum of the first and fourth terms is also symmetric in (µ, ν), and cancels. We are then left with
δRλ
- = ∇µ(δΓλ ) − ∇ν(δΓλµσ).
- (23)
(24) (25)
- σµν
- νσ
This implies that and so
δRσν = δRλ
= ∇λ(δΓλ ) − ∇ν(δΓλλσ),
- σλν
- νσ
- ꢃ
- ꢄ
gσνδRσν = ∇λ(gσνδΓλ ) − ∇ν(gσνδΓλ ) = ∇µ gσνδΓµνσ − gσµδΓλ
.
- νσ
- λσ
- λσ
Hence, we find that, upon small variations of the metric, we have
- Z
- Z
- Z
- Z
- p
- p
- p
- p
- ꢃ
- ꢄ
δd4x |g|R = d4x |g|Gµνδgµν
+
d4x |g|∇µ gσνδΓµνσ − gσµδΓλ
=
d4x |g|Gµνδgµν
,
(26)
λσ
since the last term is a total divergence and hence integrgates to zero provided metric perturbations (and their partial derivatives) vanish “at infinity”.
The convention is to define the Einstein-Hilbert action with a prefactor of 1/(16π):
Zp
1
S = d4x |g|R + SM
,
(27)
16π
and to define the stress-energy tensor of matter as
Tµν = −
2
δSM
δgµν
p
(28)
|g|
By varying the total action with respect to the metric tensor, we see that we recover the Einstein Field Equation Gµν = 8πTµν. See e.g. Carroll for a couple of worked out examples (scalar field, electromagnetic field).