
188 The Einstein field equations
t
The constant is called the cosmological constant, and was originally not present
in Einstein’s equations; he inserted it many years later in order to obtain static
cosmological solutions – solutions for the large-scale behavior of the universe – that
he felt at the time were desirable. Observations of the expansion of the universe sub-
sequently made him reject the term and regret he had ever invented it. However,
recent astronomical observations strongly suggest that it is small but not zero. We
shall return to the discussion of in Ch. (12), but for the moment we shall set
= 0. The justification for doing this, and the possible danger of it, are discussed in
Exer. 18, § 8.6.
The value k = 8π is obtained by demanding that Einstein’s equations predict the correct
behavior of planets in the solar system. This is the Newtonian limit, in which we must
demand that the predictions of GR agree with those of Newton’s theory when the latter
are well tested by observation. We saw in the last chapter that the Newtonian motions are
produced when the metric has the form Eq. (7.8). One of our tasks in this chapter is to
show that Einstein’s equations, Eq. (8.10), do indeed have Eq. (7.8) as a solution when
we assume that gravity is weak (see Exer. 3, § 8.6). We could, of course, keep k arbitrary
until then, adjusting its value to whatever is required to obtain the solution, Eq. (7.8). It
is more convenient, however, for our subsequent use of the equations of this chapter if
we simply set k to 8π at the outset and verify at the appropriate time that this value is
correct.
Eq. (8.10) should be regarded as a system of ten coupled differential equations (not
16, since T
αβ
and G
αβ
are symmetric). They are to be solved for the ten components
g
αβ
when the source T
αβ
is given. The equations are nonlinear, but they have a well-
posed initial-value structure – that is, they determine future values of g
αβ
from given initial
data. However, one point must be made: since {g
αβ
} are the components of a tensor in
some coordinate system, a change in coordinates induces a change in them. In particular,
there are four coordinates, so there are four arbitrary functional degrees of freedom among
the ten g
αβ
. It should be impossible, therefore, to determine all ten g
αβ
from any initial
data, since the coordinates to the future of the initial moment can be changed arbitrarily. In
fact, Einstein’s equations have exactly this property: the Bianchi identities
G
αβ
;β
= 0 (8.11)
mean that there are four differential identities (one for each value of α above) among the
ten G
αβ
. These ten, then, are not independent, and the ten Einstein equations are really
only six independent differential equations for the six functions among the ten g
αβ
that
characterize the geometry independently of the coordinates.
These considerations are of key importance if we want to solve Einstein’s equations to
watch systems evolve in time from some initial state. In this book we will do this in a lim-
ited way for weak gravitational waves in Ch. (9). Because of the complexity of Einstein’s
equations, dynamical situations are usually studied numerically. The field of numerical
relativity has evolved a well-defined approach to the problem of separating the coordinate
freedom in g
αβ
from the true geometric and dynamical freedom. This is described in more
advanced texts, for instance Misner et al. (1973), or Hawking and Ellis (1973), see also
Choquet-Bruhat and York (1980) or the more recent review by Cook (2000). It will suffice