Appendix A The Laplacian in a spherical coordinate system
In order to be able to deduce the most important physical consequences from the Poisson equation (12.5), which represents the Newtonian limit of Einstein’s field equations, we must know the form of the Laplacian in a spherical coordinate system. This is because most important applications of any theory of gravitation involve physical systems, for example planets, stars, black holes and the whole universe, that are spherically symmetric. Of course it is possible to find the wanted form of the Laplacian in collections of mathematical formulae. However, in the spirit of the rest of this book, we here offer a detailed deduction. The Laplace operator, or more commonly called the Laplacian, is defined as the divergence of the gradient, i.e.
r2 div grad : (A.1)
We shall find an expression of the Laplacian valid in an arbitrary orthogonal coordinate system, and then specialize to a spherical coordinate system. In order to calculate the Laplacian in an orthogonal coordinate system, we must first find expressions for the gradient and the divergence in such coordinate systems. We start with the gradient. An operator is something which acts upon a function and changes it in a prescribed way. The gradient operator acts upon a scalar function by differentiating it, and gives out a vector called the gradient of the function. In an arbitrary coordinate system the gradient operator is defined as a vector with covariant components @ r : (A.2) i @xi According to the rule (5.77) for raising an index, the contravariant components are
@ ri D gij r D gij : (A.3) j @xj
Ø. Grøn and A. Næss, Einstein’s Theory: A Rigorous Introduction 313 for the Mathematically Untrained, DOI 10.1007/978-1-4614-0706-5, © Springer Science+Business Media, LLC 2011 314 A The Laplacian in a spherical coordinate system
We shall need to know the expression for the gradient in the class of coor- dinate systems with orthogonal coordinate axes. Such coordinate systems have diagonal metric tensors, i.e. gij D 0 for i ¤ j . In this case the only non-vanishing contravariant components of the metric tensor are given in terms of the covariant components by Eq. (5.75).
ii g D 1=gii: Inserting this into Eq. (A.3)gives
1 @ ri D : i gii @x
Next we shall find an expression for the divergence in terms of ordinary partial derivatives and the components of the metric tensor in an orthogonal coordinate system. The divergence is given as a covariant derivative in Eq. (10.3). Inserting the expression (7.15) for the covariant derivative into Eq. (10.3), leads to
@F i div FE D C F k i : (A.4) @xi ki
i Hence, we must find an expression for the Christoffel symbols ki in an orthogonal coordinate system with a diagonal metric tensor. For this purpose we use Eq. (7.30) with D i, D k, D i, which for a diagonal metric tensor reduces to 1 @g @g @g i D ii C ik ik : ki k i i 2gii @x @x @x The last two terms cancel each other, so we are left with
1 @g i D ii : ki k (A.5) 2gii @x
Inserting this into Eq. (A.4)weget
@F i 1 @g FE D C F k ii : div i i (A.6) @x 2gii @x
We now define g g11 g22 g33: (A.7)
Differentiation g by means of the product rule (2.22) we first put u D g11 g22 and v D g33 and get @g @.g g / @g D 11 22 g C g g 33 : (A.8) @xk @xk 33 11 22 @xk A The Laplacian in a spherical coordinate system 315
Using the product rule once more, we have
@.g g / @g @g 11 22 D 11 g C g 22 : @xk @xk 22 11 @xk Inserting this into Eq. (A.8)gives
@g @g @g @g D 11 g g C g 22 g C g g 33 : @xk @xk 22 33 11 @xk 33 11 22 @xk Multiplying the numerators and the denominators of the terms at the right-hand side by g11, g22,andg33, respectively, reordering, and using Eq. (A.7), we get
@g g @g g @g g @g D 11 C 22 C 33 : k k k k @x g11 @x g22 @x g33 @x
Using Einstein’s summation convention this may be written as
@g g @g D ii : k k @x gii @x
Dividing by g and exhanging the left-hand and right-hand sides, we get
1 @g 1 @g ii D : k k gii @x g @x
Inserting this into Eq. (A.6)weget
@F i 1 @g div FE D C F k : (A.9) @xk 2g @xk p In order to simplify this expression, we differentiate g. Using the rule (2.36) with n D 1=2 and the chain rule (2.31), we obtain
p 1=2 @ g @g 1 @g 1 @g D D g 1=2 D p : @xk @xk 2 @xk 2 g @xk p Dividing by g gives p 1 @ g 1 @g p D : g @xk 2g @xk Thus, Eq. (A.9) can be written as p @F i 1 @ g div FE D C F i p : (A.10) @xi g @xi 316 A The Laplacian in a spherical coordinate system
Using, once more, the product rule for differentiation, we get p p 1 @. gFi / 1 p @F i @ g p D p g C F i g @xi g @xi @xi p @F i 1 @ g D C F i p : @xi g @xi
Comparing with Eq. (A.10) we see that this equation can be written as
1 @ p div FE D p gFi : (A.11) g @xi
Inserting F i Dri from Eq. (A.3) and the expression (A.11) for the divergence into Eq. (A.1), we ultimately arrive at the expression for the Laplacian in an arbitrary orthogonal coordinate system p 1 @ g @ r2 : p i (A.12) g @xi gii @x
In the case of a Cartesian coordinate system, the only non-vanishing components of p themetrictensoraregxx D gyy D gzz D 1.Then g D 1, and performing the summation over j in Eq. (A.12), we get
@2 @2 @2 r2 D C C : @x2 @y 2 @z2
Using Einstein’s summation convention this may be written as
@2 r2 D i (A.13) @xi @x valid in Cartesian coordinates. The expression (A.12) shall now be used to calculate the Laplace operator in a spherical coordinate system. Then we need the line element of Euclidean 3-dimen- sional space as expressed in a spherical coordinate system. The basis vectors of this coordinate system are given in terms of the basis vectors of a Cartesian coordinate system in Eq. (6.22). The components of the metric tensor in this coordinate system are given by the scalar products of the basis vectors in Eq. (6.22). Since the vectors are orthogonal to each other only the products of each vector with itself are different from zero. Using Eq. (4.11)wefind
2 2 2 2 2 grr DEer Eer D sin cos ' C sin sin ' C cos D sin2 cos2 ' C sin2 ' C cos2 D sin2 C cos2 D 1; A The Laplacian in a spherical coordinate system 317 2 2 2 2 2 2 g DEe Ee D r cos cos ' C cos sin ' C sin D r2 cos2 cos2 ' C sin2 ' C sin2 D r2 cos2 C sin2 D r2; 2 2 2 2 2 g'' DEe' Ee' D r sin sin ' C sin cos ' D r2 sin2 sin2 ' C cos2 ' D r2 sin2 :
Thus the line element of flat space in spherical coordinates has the form
d`2 D dr2 C r2 d 2 C r2 sin2 d'2:
The components of the metric tensor in a spherical coordinate system are therefore
2 2 2 grr D 1; g D r ; and g'' D r sin : (A.14)
Performing the summation over i in Eq. (A.12), with x1 D r, x2 D ,andx3 D ', we have p p 1 @ g @ 1 @ g @ r2 D p C p g @r grr @r g @ g @ p 1 @ g @ C p : g @' g'' @'
Inserting the expressions (A.14) and using that g D r4 sin2 ,weget 1 @ @ r2 D r2 sin r2 sin @r @r 1 @ r2 sin @ C r2 sin @ r2 @ 1 @ r2 sin @ C : r2 sin @' r2 sin2 @'
In the first term sin is constant during the differentiation with respect to r,and can be moved to the numerator in front of @[email protected]