Table of Contents  Back to Lecture 13: Three Basic Premises of General Relativity  On to Lecture 15: The Schwarzschild Metric and Event Horizons 
14. The Einstein Field Equations and Derivation of Newton's Law
Einstein's field equations show how the sources of gravitational fields alter the metric. They can actually be motivated by Newton's law for gravitational potential , with which we begin this discussion.
First, Newton's law postulates the existence of a certain scalar field , called gravitational potential which exerts a force on a unit mass given by
Further, satisfies
where _{} is the mass density and G is a constant. (The divergence theorem then gives the more familiar F = = GM/r^{2} for a spherical source of mass M  see the exercise set.) In relativity, we need an invariant analogue of (I). First, we generalize the mass density to energy density (recall that energy and mass are interchangeable according to relativity), which in turn is only one of the components of the stressenergy tensor T. Thus we had better use the whole of T.
Question What about the mysterious gravitational potential ?
Answer That is a more subtle issue. Since the second principle of general relativity tells us that particles move along geodesics, we should interpret the gravitational potential as somehow effecting the geodesics. But the most fundamental determinant of geodesics is the underlying metric g. Thus we will generalize to g. In other words, Einstein replaced a mysterious "force" by a purely geometric quantity. Put another way, gravity is nothing but a distortion of the local geometry in spacetime. But we are getting ahead of ourselves...
Finally, we generalize the (second order differential) operator to some yettobedetermined second order differential operator . This allows us to generalize (I) to
where k is some constant. In an MCRF, (g) is some linear combination of g^{ab}_{,ij}, g^{ab}_{,i} and g^{ab}, and must also be symmetric (since T is). Examples of such a tensors are the Ricci tensors R^{ab}, g^{ab}R, as well as g^{ab}. Let us take a linear combination as our candidate:
We now apply the conservation laws T^{ab}_{b} = 0, giving
since g^{ab}_{b} = 0 already (Exercise Set 8 #4). But in Lecture 9 we also saw that
(R^{ab}    1 2 
g^{ab}R)_{b}  =  0,  ... (b) 
where the term in parentheses is the Einstein tensor G^{ab}. Calculating (a)  (b), using the product rule for differentiation and the fact that g^{ab}_{b} = 0, we find
(  +  1 2 
)g^{ab}R_{b}  =  0, 
giving (upon multiplication by g^{**})
(  +  1 2 
)R_{j}  =  0, 
Finally, the requirement that these equations reduce to Newton's for v/c << 1 tells us that k = 8 (discussed below) so that we have
G^{ab} + g^{ab} = 8T^{ab} 
The constant is called the cosmological constant. Einstein at first put = 0, but later changed his mind when looking at the large scale behavior of the universe. Later still, he changed his mind again, and expressed regret that he had ever come up with it in the first place. The cosmological constant remains a problem child to this day . We shall set it equal to zero in what follows.
Solution of Einstein's Equations for Static Spherically Symmetric Stars
In the case of spherical symmetry, we use polar coordinates (r, , , t) with origin thought of as at the center of the star as our coordinate system (note it is singular there, so in fact this coordinate system does not include the origin) and restrict attention to g of the form
g_{**}  = 

, 
or
where each of the coordinates is a function of r and t only. In other words, at any fixed time t, the surfaces = const, = const and r = const are all orthogonal. (This causes the zeros to be in the positions shown.)
Question Explain why the nonzeros terms have the above form.
Answer For motivation, let us first look at the standard metric on a 2sphere of radius r: (see Example 5.2(d))
g_{**}  = 

. 
If we throw r in as the third coordinate, we could calculate
g_{**}  = 

. 
Moving into Minkowski space, we have
ds^{2}  =  dx^{2} + dy^{2} + dz^{2}  dt^{2} 
=  dr^{2} + r^{2}(d^{2} + sin^{2} d^{2})  dt^{2}, 
Minkowski Space Metrtic in Polar Coordinates

For the general spherically symmetric stellar medium, we can still define the radial coordinate to make g_{} = r^{2} (through adjustment by scaling if necessary). Further, we take as the definition of spherical symmetry, that the geometry of the surfaces r = t = const. are spherical, thus foring us to have the central 22 block.
For static spherical symmetry, we also require, among other things, (a) that the geometry be unchanged under timereversal, and (b) that g be independent of time t. For (a), if we change coordinates using
then the metric remains unchanged; that is, = g. But changing coordinates in this way amounts to multiplying on the left and right (we have an order 2 tensor here) by the changeofcoordinates matrix diag (1, 1, 1, 1), giving
_{**}  = 

. 
Setting = g gives g_{rt} = 0. Combining this with (b) results in g of the form
g_{**}  = 

, 
where we have introduced the exponentials to fix the signs, and where = (r), and = (r). Using this version of g, we can calculate the Einstein tensor to be (see the exercise set!)

We also need to calculate the stress energy tensor,
In the static case, there is assumed to be no flow of star material in our frame, so that u^{1} = u^{2} = u^{3} = 0. Further, the normal condition for four velocity, "u, u' = 1, gives
[0, 0, 0, u^{4}] 


=  1, 
whence
so that T^{44} = (_{}+p)e^{2} + p(e^{2}) (note that we are using g^{**} here). Hence,
T^{**}  =  (_{} + p)u^{*}u^{*} + pg^{**}  
= 


= 

(a) Equations of Motion T^{ab}_{b} = 0
To solve these, we first notice that we are not in an inertial frame (the metric g is not nice at the origin; in fact, nothing is even defined there!) so we need the Christoffel symbols, and use
T^{ab}_{b}  =  x^{b} 
+  _{k}^{a}_{b}T^{kb}  +  _{b}^{b}_{k}T^{ak}, 
where
_{h}^{p}_{k}  =  1 2 
g^{lp}  x^{h} 
+  x^{k} 
  x^{l} 
. 
Now, lots of the terms in T^{ab}_{b} vanish by symmetry, and the restricted nature of the functions. We shall focus on a = 1, the rcoordinate. We have:
and we calculate these terms oneatatime.
a = 1, b = 1:  T^{11}_{1}  =  x^{1} 
+  _{1}^{1}_{1}T^{11} + _{1}^{1}_{1}T^{11}. 
To evaluate this, first look at the term _{1}^{1}_{1}:
_{1}^{1}_{1 }  =  1 2 
g^{l1}(g_{1l,1} + g_{l1,1}  g_{11,l})  
=  1 2 
g^{11}(g_{11,1} + g_{11,1}  g_{11,1})  (because g is diagonal, whence l = 1)  
=  1 2 
g^{11}(g_{11,1})  
=  1 2 
e^{2}e^{2}.2'(r) = '(r) 
Hence,
T^{11}_{1}  =  dr 
e^{2} + ( 2p'(r)e^{2}) + 2'(r)pe^{2}  =  dr 
e^{2}. 
Now for the next term:
a = 1, b = 2:  T^{12}_{2}  = 


= 


= 


= 


= 

Similarly (exercise set)
Finally,
a = 1, b = 4:  T^{14}_{4}  = 


= 


= 


= 

Hence, the conservation equation becomes
T^{1a}_{a} = 0 



This gives the pressure gradient required to keep the plasma static in a star.
Note In classical mechanics, the term on the right has _{} rather than _{}+p. Thus, the pressure gradient is larger in relativistic theory than in classical theory. This increased pressure gradient corresponds to greater values for p, and hence bigger values for all the components of T. By Einstein's field equations, this now leads to even greater values of (manifested as gravitational force) thereby causing even larger values of the pressure gradient. If p is large to begin with (big stars) this vicious cycle diverges, ending in the gravitational collapse of a star, leading to neutron stars or, in extreme cases, black holes. You can go directly to Lecture 15 on stellar collapse to find out more.
(b) Einstein Field Equations G^{ab} = 8T^{ab}
Looking at the (4, 4) component first, and substituting from the expressions for G and T, we find
r^{2} 
e^{2}  dr 
[r(1e^{2}]  =  8_{}e^{2}. 
If we define
1 2 
r(1e^{2})  =  m(r), 
then the equation becomes
r^{2} 
e^{2}  dr 
=  4_{}e^{2}, 
or
dr 
=  4r^{2}_{}  ... (I) 
This looks like an equation for classical mass, since classically,
M(R)  =  R 0 
4r^{2}_{}(r) dr  , 
where the integrand is the mass of a shell whose thickness is dr. Thus,
dr 
=  4^{2}_{}(r). 
Here, _{} is energy density, and by our choice of units, energy is equal to rest mass, so we interpret m(r) as the total mass of the star enclosed by a sphere of radius r.
Now look at the (1, 1) component:
2 r 
'e^{4}    r^{2} 
(1e^{2})  =  8pe^{2} 
2 r 
'    r^{2} 
(1e^{2})  =  8pe^{2} 
2r'  e^{2} (1e^{2}) = 8r^{2}pe^{2} 
' = e^{2}  2r 
In the expression for m, solve for e^{2} to get
e^{2}  =  12m/r 
, 
giving
dr 
=  2r(12m/r) 
or
dr 
=  r(r2m) 
... (II) 
It can be checked using the Bianchi identities that we in fact get no additional information from the (2,2) and (3,3) components, so we ignore them.
Consequences of the Field Equations: Outside the Star
Outside the star we take p = _{} 0, and m(r) = M, the total stellar mass, getting
(I):  dr  =  0  (nothing new, since m = M = constant) 
(II):  dr 
=  r(r2M) 
, 
which is a separable first order differential equation with solution
e^{2}  =  1   r 
, 
if we impose the boundary condition 0 as r+. (See the exercise Set).
Recalling from the definition of m that
e^{2}  =  12M/r 
, 
we can now express the metric outside a star as follows:
Schwarzschild Metric

In the exercise set, you will see how this leads to Newton's Law of Gravity.
Exercise Set 13
1. Use ^{2} = 4G_{} and the divergence theorem to deduce Newton's law = GM/r^{2} for a spherical mass of uniform density _{}.
2. Calculate the Einstein tensor for the metric g = diag(e^{2}, r^{2}, r^{2}sin, e^{2}), and verify that it agrees with that in the notes.
3. Referring to the notes above, show that T^{13}_{3} = 0.
4. Show that T^{i4}_{4} = 0 for i = 2, 3, 4.
5. If we impose the condition that, far from the star, spacetime is flat, show that this is equivalent to saying that _{r}_{}_{+}(r) = _{r}_{}_{+}(r) = 0. Hence obtain the formula e^{2} = 1 2M/r.
6. A Derivation of Newton's Law of Gravity
(a) Show that, at a large distance R from a static stable star, the Schwarzschild metric can be approximated as
g^{**} 

(b) (Schutz, p. 272 #9) Define a new coordinate by R = (1+M/)^{2}, and deduce that, in terms of the new coordinates (ignoring terms of order 1/R^{2})
g^{**} 

(c) Now convert to Cartesian coordinates, (x, y, z, t) to obtain
g^{**} 

(d) Now refer to the last formula in Lecture 13, and obtain Newton's Law of Gravity. To how many kilograms does one unit of M correspond?
Table of Contents  Back to Lecture 13: Three Basic Premises of General Relativity  On to Lecture 15: The Schwarzschild Metric and Event Horizons 