/

/
/

Z
Z
Z

.
.
.
.
.
`
``
`
`
`
`
`
`
`
`
`

~
~
~
-
-
-

.
.
. .
.
.
.
IFT
Instituto de Fsica Teorica
Universidade Estadual Paulista
An Introduction to
GENERAL RELATIVITY
R. Aldrovandi and J. G. Pereira
March-April/2004
A Preliminary Note
These notes are intended for a two-month, graduate-level course. Ad-
dressed to future researchers in a Centre mainly devoted to Field Theory,
they avoid the ex cathedra style frequently assumed by teachers of the sub-
ject. Mainly, General Relativity is not presented as a nished theory.
Emphasis is laid on the basic tenets and on comparison of gravitation
with the other fundamental interactions of Nature. Thus, a little more space
than would be expected in such a short text is devoted to the equivalence
principle.
The equivalence principle leads to universality, a distinguishing feature of
the gravitational eld. The other fundamental interactions of Naturethe
electromagnetic, the weak and the strong interactions, which are described
in terms of gauge theoriesare not universal.
These notes, are intended as a short guide to the main aspects of the
subject. The reader is urged to refer to the basic texts we have used, each
one excellent in its own approach:
L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Perg-
amon, Oxford, 1971)
C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation (Freeman,
New York, 1973)
S. Weinberg, Gravitation and Cosmology (Wiley, New York, 1972)
R. M. Wald, General Relativity (The University of Chicago Press,
Chicago, 1984)
J. L. Synge, Relativity: The General Theory (North-Holland, Amster-
dam, 1960)
i
Contents
1 Introduction 1
1.1 General Concepts . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Some Basic Notions . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 The Equivalence Principle . . . . . . . . . . . . . . . . . . . . 3
1.3.1 Inertial Forces . . . . . . . . . . . . . . . . . . . . . . . 5
1.3.2 The Wake of Non-Trivial Metric . . . . . . . . . . . . . 10
1.3.3 Towards Geometry . . . . . . . . . . . . . . . . . . . . 13
2 Geometry 18
2.1 Dierential Geometry . . . . . . . . . . . . . . . . . . . . . . . 18
2.1.1 Spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.1.2 Vector and Tensor Fields . . . . . . . . . . . . . . . . . 29
2.1.3 Dierential Forms . . . . . . . . . . . . . . . . . . . . . 35
2.1.4 Metrics . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.2 Pseudo-Riemannian Metric . . . . . . . . . . . . . . . . . . . . 44
2.3 The Notion of Connection . . . . . . . . . . . . . . . . . . . . 46
2.4 The LeviCivita Connection . . . . . . . . . . . . . . . . . . . 50
2.5 Curvature Tensor . . . . . . . . . . . . . . . . . . . . . . . . . 53
2.6 Bianchi Identities . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.6.1 Examples . . . . . . . . . . . . . . . . . . . . . . . . . 57
3 Dynamics 63
3.1 Geodesics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3.2 The Minimal Coupling Prescription . . . . . . . . . . . . . . . 71
3.3 Einsteins Field Equations . . . . . . . . . . . . . . . . . . . . 76
3.4 Action of the Gravitational Field . . . . . . . . . . . . . . . . 79
3.5 Non-Relativistic Limit . . . . . . . . . . . . . . . . . . . . . . 82
3.6 About Time, and Space . . . . . . . . . . . . . . . . . . . . . 85
3.6.1 Time Recovered . . . . . . . . . . . . . . . . . . . . . . 85
3.6.2 Space . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
ii
3.7 Equivalence, Once Again . . . . . . . . . . . . . . . . . . . . . 90
3.8 More About Curves . . . . . . . . . . . . . . . . . . . . . . . . 92
3.8.1 Geodesic Deviation . . . . . . . . . . . . . . . . . . . . 92
3.8.2 General Observers . . . . . . . . . . . . . . . . . . . . 93
3.8.3 Transversality . . . . . . . . . . . . . . . . . . . . . . . 95
3.8.4 Fundamental Observers . . . . . . . . . . . . . . . . . . 96
3.9 An Aside: Hamilton-Jacobi . . . . . . . . . . . . . . . . . . . 99
4 Solutions 107
4.1 Transformations . . . . . . . . . . . . . . . . . . . . . . . . . . 107
4.2 Small Scale Solutions . . . . . . . . . . . . . . . . . . . . . . . 111
4.2.1 The Schwarzschild Solution . . . . . . . . . . . . . . . 111
4.3 Large Scale Solutions . . . . . . . . . . . . . . . . . . . . . . . 128
4.3.1 The Friedmann Solutions . . . . . . . . . . . . . . . . . 128
4.3.2 de Sitter Solutions . . . . . . . . . . . . . . . . . . . . 135
5 Tetrad Fields 141
5.1 Tetrads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
5.2 Linear Connections . . . . . . . . . . . . . . . . . . . . . . . . 146
5.2.1 Linear Transformations . . . . . . . . . . . . . . . . . . 146
5.2.2 Orthogonal Transformations . . . . . . . . . . . . . . . 148
5.2.3 Connections, Revisited . . . . . . . . . . . . . . . . . . 150
5.2.4 Back to Equivalence . . . . . . . . . . . . . . . . . . . 154
5.2.5 Two Gates into Gravitation . . . . . . . . . . . . . . . 159
6 Gravitational Interaction of the Fundamental Fields 161
6.1 Minimal Coupling Prescription . . . . . . . . . . . . . . . . . 161
6.2 General Relativity Spin Connection . . . . . . . . . . . . . . . 162
6.3 Application to the Fundamental Fields . . . . . . . . . . . . . 164
6.3.1 Scalar Field . . . . . . . . . . . . . . . . . . . . . . . . 164
6.3.2 Dirac Spinor Field . . . . . . . . . . . . . . . . . . . . 165
6.3.3 Electromagnetic Field . . . . . . . . . . . . . . . . . . 166
7 General Relativity with Matter Fields 170
7.1 Global Noether Theorem . . . . . . . . . . . . . . . . . . . . . 170
7.2 EnergyMomentum as Source of Curvature . . . . . . . . . . . 171
7.3 EnergyMomentum Conservation . . . . . . . . . . . . . . . . 173
7.4 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
7.4.1 Scalar Field . . . . . . . . . . . . . . . . . . . . . . . . 175
7.4.2 Dirac Spinor Field . . . . . . . . . . . . . . . . . . . . 176
iii
7.4.3 Electromagnetic Field . . . . . . . . . . . . . . . . . . 177
8 Closing Remarks 179
Bibliography 180
iv
Chapter 1
Introduction
1.1 General Concepts
1.1 All elementary particles feel gravitation the same. More specically,
particles with dierent masses experience a dierent gravitational force, but
in such a way that all of them acquire the same acceleration and, given the
same initial conditions, follow the same path. Such universality of response
is the most fundamental characteristic of the gravitational interaction. It is a
unique property, peculiar to gravitation: no other basic interaction of Nature
has it.
Due to universality, the gravitational interaction admits a description
which makes no use of the concept of force. In this description, instead of
acting through a force, the presence of a gravitational eld is represented
by a deformation of the spacetime structure. This deformation, however,
preserves the pseudo-riemannian character of the at Minkowski spacetime
of Special Relativity, the non-deformed spacetime that represents absence of
gravitation. In other words, the presence of a gravitational eld is supposed
to produce curvature, but no other kind of spacetime deformation.
A free particle in at space follows a straight line, that is, a curve keeping
a constant direction. A geodesic is a curve keeping a constant direction on
a curved space. As the only eect of the gravitational interaction is to bend
spacetime so as to endow it with curvature, a particle submitted exclusively
to gravity will follow a geodesic of the deformed spacetime.
1
This is the approach of Einsteins General Relativity, according to which
the gravitational interaction is described by a geometrization of spacetime.
It is important to remark that only an interaction presenting the property of
universality can be described by such a geometrization.
1.2 Some Basic Notions
1.2 Before going further, let us recall some general notions taken from
classical physics. They will need renements later on, but are here put in a
language loose enough to make them valid both in the relativistic and the
non-relativistic cases.
Frame: a reference frame is a coordinate system for space positions, to which
a clock is bound.
Inertia: a reference frame such that free (unsubmitted to any forces) mo-
tion takes place with constant velocity is an inertial frame; in classical
physics, the force law in an inertial frame is m
dv
k
dt
= F
k
; in Special
Relativity, the force law in an inertial frame is
m
d
ds
U
a
= F
a
, (1.1)
where U is the four-velocity U = (, v/c), with = 1/
_
1 v
2
/c
2
(as
U is dimensionless, F above has not the mechanical dimension of a force
only Fc
2
has). Incidentally, we are stuck to cartesian coordinates to
discuss accelerations: the second time derivative of a coordinate is an
acceleration only if that coordinate is cartesian.
Transitivity: a reference frame moving with constant velocity with respect
to an inertial frame is also an inertial frame;
Relativity: all the laws of nature are the same in all inertial frames; or,
alternatively, the equations describing them are invariant under the
transformations (of space coordinates and time) taking one inertial
frame into the other; or still, the equations describing the laws of Nature
in terms of space coordinates and time keep their forms in dierent
inertial frames; this principle can be seen as an experimental fact; in
non-relativistic classical physics, the transformations referred to belong
to the Galilei group; in Special Relativity, to the Poincare group.
2
Causality: in non-relativistic classical physics the interactions are given by
the potential energy, which usually depends only on the space coordi-
nates; forces on a given particle, caused by all the others, depend only
on their position at a given instant; a change in position changes the
force instantaneously; this instantaneous propagation eect or ac-
tion at a distance is a typicallly classical, non-relativistic feature; it
violates special-relativistic causality; Special Relativity takes into ac-
count the experimental fact that light has a nite velocity in vacuum
and says that no eect can propagate faster than that velocity.
Fields: there have been tentatives to preserve action at a distance in a
relativistic context, but a simpler way to consider interactions while
respecting Special Relativity is of common use in eld theory: interac-
tions are mediated by a eld, which has a well-dened behaviour under
transformations; disturbances propagate, as said above, with nite ve-
locities.
1.3 The Equivalence Principle
Equivalence is a guiding principle, which inspired Einstein in his construction
of General Relativity. It is rmly rooted on experience.

In its most usual form, the Principle includes three subprinciples: the
weak, the strong and that which is called Einsteins equivalence principle.
We shall come back and forth to them along these notes. Let us shortly list
them with a few comments.
1.3 The weak equivalence principle: universality of free fall, or inertial
mass = gravitational mass.
In a gravitational eld, all pointlike structureless particles fol-
low one same path; that path is xed once given (i) an initial
position x(t
0
) and (ii) the correspondent velocity x(t
0
).
This leads to a force equation which is a second order ordinary dierential
equation. No characteristic of any special particle, no particular property

Those interested in the experimental status will nd a recent appraisal in C. M. Will,

The Confrontation between General Relativity and Experiment, arXiv:gr-qc/0103036 12
Mar 2001. Theoretical issues are discussed by B. Mashhoon, Measurement Theory and
General Relativity, gr-qc/0003014, and Relativity and Nonlocality, gr-qc/0011013 v2.
3
appears in the equation. Gravitation is consequently universal. Being uni-
versal, it can be seen as a property of space itself. It determines geometrical
properties which are common to all particles. The weak equivalence princi-
ple goes back to Galileo. It raises to the status of fundamental principle a
deep experimental fact: the equality of inertial and gravitational masses of
all bodies.
The strong equivalence principle: (Einsteins lift) says that
Gravitation can be made to vanish locally through an appro-
priate choice of frame.
It requires that, for any and every particle and at each point x
0
, there exists
a frame in which x

= 0.
Einsteins equivalence principle requires, besides the weak principle,
the local validity of Poincare invariance that is, of Special Relativity. This
invariance is, in Minkowski space, summed up in the Lorentz metric. The
requirement suggests that the above deformation caused by gravitation is a
change in that metric.
In its complete form, the equivalence principle
1. provides an operational denition of the gravitational interaction;
2. geometrizes it;
3. xes the equation of motion of the test particles.
1.4 Use has been made above of some undened concepts, such as path,
and local. A more precise formulation requires more mathematics, and will
be left to later sections. We shall, for example, rephrase the Principle as a
prescription saying how an expression valid in Special Relativity is changed
once in the presence of a gravitational eld. What changes is the notion of
derivative, and that change requires the concept of connection. The prescrip-
tion (of minimal coupling) will be seen after that notion is introduced.
4
1.5 Now, forces equally felt by all bodies were known since long. They are
the inertial forces, whose name comes from their turning up in non-inertial
frames. Examples on Earth (not an inertial system !) are the centrifugal
force and the Coriolis force. We shall begin by recalling what such forces
are in Classical Mechanics, in particular how they appear related to changes
of coordinates. We shall then show how a metric appears in an non-inertial
frame, and how that metric changes the law of force in a very special way.
1.3.1 Inertial Forces
1.6 In a frame attached to Earth (that is, rotating with a certain angular
velocity ), a body of mass m moving with velocity

X on which an external
force F
ext
acts will actually experience a strange total force. Let us recall
in rough brushstrokes how that happens.
A simplied model for the motion of a particle in a system attached to
Earth is taken from the classical formalism of rigid body motion.

It runs as
follows:
The rotating
Earth
Start with an inertial cartesian system, the space system (inertial means
we insist devoid of proper acceleration). A point particle will
have coordinates {x
i
}, collectively written as a column vector x = (x
i
).
Under the action of a force f , its velocity and acceleration will be, with
respect to that system, x and x. If the particle has mass m, the force
will be f = m x.
Consider now another coordinate system (the body system) which rotates
around the origin of the rst. The point particle will have coordinates
X in this system. The relation between the coordinates will be given
by a rotation matrix R,
X = R x.
The forces acting on the particle in both systems are related by the same

The standard approach is given in H. Goldstein,Classical Mechanics, AddisonWesley,

Reading, Mass., 1982. A modern description can be found in J. L. McCauley,Classical
Mechanics, Cambridge University Press, Cambridge, 1997.
5
relation,
F = R f .
We are using symbols with capitals (X, F, , . . . ) for quantities re-
ferred to the body system, and the corresponding small letters (x, f ,
, . . . ) for the same quantities as seen from the space system.
Now comes the crucial point: as Earth is rotating with respect to the space
system, a dierent rotation is necessary at each time to pass from that
system to the body system; this is to say that the rotation matrix R
is time-dependent. In consequence, the velocity and the acceleration
seen from Earths system are given by

X =

R x + R x

X =

R x + 2

R x + R x. (1.2)
Introduce the matrix = R
1

R. It is an antisymmetric 3 3 matrix,
consequently equivalent to a vector. That vector, with components

k
=
1
2

k
ij

ij
(1.3)
(which is the same as
ij
=
ijk

k
), is Earths angular velocity seen
from the space system. is, thus, a matrix version of the angular
velocity. It will correspond, in the body system, to
= RR
1
=

R R
1
.
Comment 1.1 Just in case,
ijk
is the 3-dimensional Kronecker symbol in 3-
dimensional space:
123
= 1; any odd exchange of indices changes the sign;
ijk
= 0
if there are repeated indices. Indices are raised and lowered with the Kronecker
delta
ij
, dened by
ii
= 1 and
ij
= 0 if i = j. In consequence,
ijk
=
ijk
=

i
jk
, etc. The usual vector product has components given by (v u)
i
= (v u)
i
=
ijk
u
j
v
k
. An antisymmetric matrix like , acting on a vector will give
ij
v
j
=

ijk

k
v
j
= ( v)
i
.
A few relations turn out without much ado:
2
= R
2
R
1
,

= R R
1
and
= R
1

R ,
6
or

R = R [ ] .
Substitutions put then Eq. (1.2) into the form

X+ 2

X+ [

+
2
] X = R x
The above relationship between 3 3 matrices and vectors takes matrix
action on vectors into vector products: x = x, etc. Transcribing
into vector products and multiplying by the mass, the above equation
acquires its standard form in terms of forces,
m

X = m X
. .
2m

X
. .
m

X
. .
+ F
ext
.
centrifugal Coriolis fluctuation
We have indicated the usual names of the contributions. A few words
on each of them
uctuation force: in most cases can be neglected for Earth, whose angular
velocity is very nearly constant.
centrifugal force: opposite to Earths attraction, it is already taken into
account by any balance (you are fatter than you think, your mass is
larger than suggested by your your weight by a few grams ! the ratio
is 3/1000 at the equator).
Coriolis force: responsible for trade winds, rivers one-sided overows, as-
symmetric wear of rails by trains, and the eect shown by the Foucault
pendulum.
1.7 Inertial forces have once been called cticious, because they disap-
pear when seen from an inertial system at rest. We have met them when
we started from such a frame and transformed to coordinates attached to
Earth. We have listed the measurable eects to emphasize that they are
actually very real forces, though frame-dependent.
1.8 The remarkable fact is that each body feels them the same. Think of
the examples given for the Coriolis force: air, water and iron feel them, and
7
in the same way. Inertial forces are universal, just like gravitation. This
has led Einstein to his formidable stroke of genius, to conceive gravitation as
an inertial force.
1.9 Nevertheless, if gravitation were an inertial eect, it should be ob-
tained by changing to a non-inertial frame. And here comes a problem. In
Classical Mechanics, time is a parameter, external to the coordinate system.
In Special Relativity, with Minkowskis invention of spacetime, time under-
went a violent conceptual change: no more a parameter, it became the fourth
coordinate (in our notation, the zeroth one).
Classical non-inertial frames are obtained from inertial frames by trans-
formations which depend on time. Relativistic non-inertial frames should be
obtained by transformations which depend on spacetime. Timedependent
coordinate changes ought to be special cases of more general transforma-
tions, dependent on all the spacetime coordinates. In order to be put into
a position closer to inertial forces, and concomitantly respect Special Rela-
tivity, gravitation should be related to the dependence of frames on all the
coordinates.
1.10 Universality of inertial forces has been the rst hint towards General
Relativity. A second ingredient is the notion of eld. The concept allows the
best approach to interactions coherent with Special Relativity. All known
forces are mediated by elds on spacetime. Now, if gravitation is to be
represented by a eld, it should, by the considerations above, be a universal
eld, equally felt by every particle. It should change spacetime itself. And,
of all the elds present in a space the metric the rst fundamental form,
as it is also called seemed to be the basic one. The simplest way to
change spacetime would be to change its metric. Furthermore, the metric
does change when looked at from a non-inertial frame.
1.11 The Lorentz metric of Special Relativity is rather trivial. There
is a coordinate system (the cartesian system) in which the line element of
Lorentz
metric
Minkowski space takes the form
ds
2
=
ab
dx
a
dx
b
= dx
0
dx
0
dx
1
dx
1
dx
2
dx
2
dx
3
dx
3
8
= c
2
dt
2
dx
2
dy
2
dz
2
. (1.4)
Take two points P and Q in Minkowski spacetime, and consider the in-
tegral
_
Q
P
ds =
_
Q
P
_

ab
dx
a
dx
b
.
Its value depends on the path chosen. In consequence, it is actually a func-
tional on the space of paths between P and Q,
S[
PQ
] =
_

PQ
ds. (1.5)
An extremal of this functional would be a curve such that S[] =
_
ds
= 0. Now,
ds
2
= 2 ds ds = 2
ab
dx
a
dx
b
,
so that
ds =
ab
dx
a
ds
dx
b
=
ab
U
a
dx
b
.
Thus, commuting d and and integrating by parts,
S[] =
_
Q
P

ab
dx
a
ds
dx
b
ds
ds =
_
Q
P

ab
d
ds
dx
a
ds
x
b
ds
=
_
Q
P

ab
d
ds
U
a
x
b
ds.
The variations x
b
are arbitrary. If we want to have S[] = 0, the integrand
must vanish. Thus, an extremal of S[] will satisfy
d
ds
U
a
= 0. (1.6)
This is the equation of a straight line, the force law (1.1) when F
a
= 0.
The solution of this dierential equation is xed once initial conditions are
given. We learn here that a vanishing acceleration is related to an extremal
of S[
PQ
].
1.12 Let us see through an example what happens when a force is present.
For that it is better to notice beforehand that, when considering elds, it is
9
in general the action which is extremal. Simple dimensional analysis shows
that, in order to have a real physical action, we must take
S = mc
_
ds (1.7)
instead of the length. Consider the case of a charged test particle. The
coupling of a particle of charge e to an electromagnetic potential A is given
by A
a
j
a
= e A
a
U
a
, so that the action along a curve is
S
em
[] =
e
c
_

A
a
U
a
ds =
e
c
_

A
a
dx
a
.
The variation is
S
em
[] =
e
c
_

A
a
dx
a

e
c
_

A
a
dx
a
=
e
c
_

A
a
dx
a
+
e
c
_

dA
b
x
b
=
e
c
_

b
A
a
x
b
dx
a
+
e
c
_

a
A
b
x
b
dx
a
=
e
c
_

[
b
A
a

a
A
b
]x
b
dx
a
ds
ds
=
e
c
_

F
ba
U
a
x
b
ds .
Combining the two pieces, the variation of the total action
S = mc
_
Q
P
ds
e
c
_
Q
P
A
a
dx
a
(1.8)
is
S =
_
Q
P
_

ab
mc
d
ds
U
a

e
c
F
ba
U
a
_
x
b
ds.
The extremal satises
Lorentz
force law
mc
d
ds
U
a
=
e
c
F
a
b
U
b
, (1.9)
which is the Lorentz force law and has the form of the general case (1.1).
1.3.2 The Wake of Non-Trivial Metric
Let us see now in another example that the metric changes when
viewed from a non-inertial system. This fact suggests that, if gravitation is
to be related to non-inertial systems, a gravitational eld is to be related to
a non-trivial metric.
10
1.13 Consider a rotating disc (details can be seen in Mllers book

), seen
as a system performing a uniform rotation with angular velocity on the x,
y plane:
x = r cos( + t) ; y = r sin( + t) ; Z = z ;
X = Rcos ; Y = Rsin .
This is the same as
x = X cos t Y sin t ; y = Y cos t + X sin t .
As there is no contraction along the radius (the motion being orthogonal
to it), R = r. Both systems coincide at t = 0. Now, given the standard
Minkowski line element
ds
2
= c
2
dt
2
dx
2
dy
2
dz
2
in cartesian (space, inertial) coordinates (x
0
, x
1
, x
2
, x
3
) = (ct, x, y, z), how
will a body observer on the disk see it ?
It is immediate that
dx = dr cos( + t) r sin( + t)[d + dt]
dy = dr sin( + t) + r cos( + t)[d + dt]
dx
2
= dr
2
cos
2
( + t) + r
2
sin
2
( + t)[d + dt]
2
2rdr cos( + t) sin( + t)[d + dt] ;
dy
2
= dr
2
sin
2
( + t) + r
2
cos
2
( + t)[d + dt]
2
+2rdr sin( + t) cos( + t)[d + dt]
dx
2
+ dy
2
= dR
2
+ R
2
(d
2
+
2
dt
2
+ 2ddt).
It follows from
dX
2
+ dY
2
= dR
2
+ R
2
d
2
,
that
dx
2
+ dy
2
= dX
2
+ dY
2
+ R
2

2
dt
2
+ 2R
2
ddt.

C. Mller, The Theory of Relativity, Oxford at Clarendon Press, Oxford, 1966, mainly
in 8.9.
11
A simple check shows that
XdY Y dX = R
2
d,
so that
dx
2
+ dy
2
= dX
2
+ dY
2
+ R
2

2
dt
2
+ 2XdY dt 2Y dXdt.
Thus,
ds
2
= (1

2
R
2
c
2
) c
2
dt
2
dX
2
dY
2
2XdY dt + 2Y dXdt dZ
2
.
In the moving body system, with coordinates (X
0
, X
1
, X
2
, X
3
) = (ct, X, Y, Z =
z) the metric will be
ds
2
= g

dX

dX

,
where the only non-vanishing components of the modied metric g are:
g
11
= g
22
= g
33
= 1; g
01
= g
10
= Y/c; g
02
= g
20
= X/c;
g
00
= 1

2
R
2
c
2
.
This is better visualized as the matrix
g = (g

) =
_
_
_
_
_
1

2
R
2
c
2
Y/c X/c 0
Y/c 1 0 0
X/c 0 1 0
0 0 0 1
_
_
_
_
_
. (1.10)
We can go one step further an dene the body time coordinate T to be
such that dT =
_
1
2
R
2
/c
2
dt, that is,
T =
_
1
2
R
2
/c
2
t .
This expression is physically appealing, as it is the same as T =
_
1 v
2
/c
2
t,
the time-contraction of Special Relativity, if we take into account the fact
that a point with coordinates (R, ) will have squared velocity v
2
=
2
R
2
.
We see that, anyhow, the body coordinate system can be used only for points
12
satisfying the condition R < c. In the body coordinates (cT, X, Y, Z), the
line element becomes
ds
2
= c
2
dT
2
dX
2
dY
2
dZ
2
+ 2[Y dX XdY ]
dT
_
1
2
R
2
/c
2
.
(1.11)
Time, as measured by the accelerated frame, diers from that measured in
the inertial frame. And, anyhow, the metric has changed. This is the point
we wanted to make: when we change to a non-inertial system the metric
undergoes a signicant transformation, even in Special Relativity.
Comment 1.2 Put = R/c. Matrix (1.10) and its inverse are
g = (g

) =
_
_
1
2 Y
R

X
R
0
Y
R
1 0 0

X
R
0 1 0
0 0 0 1
_
_
; g
1
= (g

) =
_
_
_
1
Y
R

X
R
0
Y
R

2 Y
2
R
2
1

2
XY
R
2
0

X
R


2
XY
R
2

2 X
2
R
2
1 0
0 0 0 1
_
_
_ .
1.3.3 Towards Geometry
1.14 We have said that the only eect of a gravitational eld is to bend
spacetime, so that straight lines become geodesics. Now, there are two quite
distinct denitions of a straight line, which coincide on at spaces but not
on spaces endowed with more sophisticated geometries. A straight line going
from a point P to a point Q is
1. among all the lines linking P to Q, that with the shortest length;
2. among all the lines linking P to Q, that which keeps the same direction
all along.
There is a clear problem with the rst denition: length presupposes a
metric a real, positive-denite metric. The Lorentz metric does not dene
lengths, but pseudo-lengths. There is always a zero-length path between
any two points in Minkowski space. In Minkowski space,
_
ds is actually
maximal for a straight line. Curved lines, or broken ones, give a smaller
pseudo-length. We have introduced a minus sign in Eq.(1.7) in order to
conform to the current notion of minimal action.
The second denition can be carried over to spacetime of any kind, but
at a price. Keeping the same direction means keeping the tangent velocity
13
vector constant. The derivative of that vector along the line should vanish.
Now, derivatives of vectors on non-at spaces require an extra concept, that
of connection which, will, anyhow, turn up when the rst denition is
used. We shall consequently feel forced to talk a lot about connections in
what follows.
1.15 Consider an arbitrary metric g, dening the interval by
general
metric
ds
2
= g

dx

dx

.
What happens now to the integral of Eq.(1.7) with a point-dependent metric?
Consider again a charged test particle, but now in the presence of a non-trivial
metric. We shall retrace the steps leading to the Lorentz force law, with the
action
S = mc
_

PQ
ds
e
c
_

PQ
A

dx

, (1.12)
but now with ds =
_
g

dx

dx

.
1. Take rst the variation
ds
2
= 2dsds = [g

dx

dx

] = dx

dx

+ 2g

dx

dx

ds =
1
2
dx

ds
dx

ds

ds + g

dx

ds
dx

ds
ds
We have conveniently divided and multiplied by ds.
2. We now insert this in the rst piece of the action and integrate by parts
the last term, getting
S = mc
_

PQ
_
1
2
dx

ds
dx

ds

d
ds
(g

dx

ds
)

ds

e
c
_

PQ
[A

dx

+ A

dx

]. (1.13)
3. The derivative
d
ds
(g

dx

ds
) is
d
ds
(g

dx

ds
) =
dx

ds
d
ds
g

+ g

d
ds
U

= U

+ g

d
ds
U

= g

d
ds
U

+ U

= g

d
ds
U

+
1
2
U

].
14
4. Collecting terms in the metric sector, and integrating by parts in the
electromagnetic sector,
S = mc
_

PQ
_
g

d
ds
U

1
2
U

)
_
x

ds

e
c
_

PQ
[

dx

dx

] = (1.14)
mc
_

PQ
g

d
ds
U

_
1
2
g

)
_
_
x

ds

e
c
_

PQ
[

dx

dx

]. (1.15)
5. We meet here an important character of all metric theories. The ex-
pression between curly brackets is the Christoel symbol, which will be
Christoel
symbol
indicated by the notation

=
1
2
g

) . (1.16)
6. After arranging the terms, we get
S =
_

PQ
_
mc g

_
d
ds
U

e
c
(

)U

_
x

ds.
(1.17)
7. The variations x

, except at the xed endpoints, is quite arbitrary. To

have S = 0, the integrand must vanish. Which gives, after contracting
with g

,
mc
_
d
ds
U

_
=
e
c
F

. (1.18)
8. This is the Lorentz law of force in the presence of a non-trivial metric.
We see that what appears as acceleration is now

=
d
ds
U

. (1.19)
15
The Christoel symbol is a non-tensorial quantity, a connection. We
shall see later that a reference frame can be always chosen in which it
vanishes at a point. The law of force
mc
_
d
ds
U

_
= F

(1.20)
will, in that frame and at that point, reduce to that holding for a trivial
metric, Eq. (1.1).
9. In the absence of forces, the resulting expression,
geodesic
equation
d
ds
U

= 0, (1.21)
is the geodesic equation, dening the straightest possible line on a
space in which the metric is non-trivial.
Comment 1.3 An accelerated frame creates the illusion of a force. Suppose a point P is
at rest. It may represent a vessel in space, far from any other body. An astronaut in
the spacecraft can use gyros and accelerometers to check its state of motion. It will never
be able to say that it is actually at rest, only that it has some constant velocity. Its own
reference frame will be inertial. Assume another craft approaches at a velocity which is
constant relative to P, and observes P. It will measure the distance from P, see that the
velocity x is constant. That observer will also be inertial.
Suppose now that the second vessel accelerates towards P. It will then see x = 0, and
will interpret this result in the normal way: there is a force pulling P. That force is clearly
an illusion: it would have opposite sign if the accelerated observer moved away from P.
No force acts on P, the force is due to the observers own acceleration. It comes from the
observer, not from P.
Comment 1.4 Curvature creates the illusion of a force. Two old travellers (say, Hero-
dotus and Pausanias) move northwards on Earth, starting from two distinct points on the
equator. Suppose they somehow communicate, and have a means to evaluate their relative
distance. They will notice that that distance decreases with their progress until, near the
pole, they will see it dwindle to nothing. Suppose further they have ancient notions, and
think the Earth is at. How would they explain it ? They would think there were some
force, some attractive force between them. And what is the real explanation ? It is simply
that Earths surface is a curved space. The force is an illusion, born from the atness
prejudice.
16
1.16 Gravitation is very weak. To present time, no gravitational bending
in the trajectory of an elementary particle has been experimentally observed.
Only large agglomerates of fermions have been seen to experience it. Never-
theless, an eect on the phase of the wave-function has been detected, both
for neutrons and atoms.

1.17 Suppose that, of all elementary particles, one single existed which did
not feel gravitation. That would be enough to change all the picture. The
underlying spacetime would remain Minkowskis, and the metric responsible
for gravitation would be a eld g

on that, by itself at, spacetime.

Spacetime is a geometric construct. Gravitation should change the geom-
etry of spacetime. This comes from what has been said above: coordinates,
metric, connection and frames are part of the dierential-geometrical struc-
ture of spacetime. We shall need to examine that structure. The next chapter
is devoted to the main aspects of dierential geometry.

The so-called COW experiment with neutrons is described in R. Colella, A. W.

Overhauser and S. A. Werner, Phys. Rev. Lett. 34, 1472 (1975). See also U. Bonse
and T. Wroblevski, Phys. Rev. Lett. 51, 1401 (1983). Experiments with atoms are
reviewed in C. J. Borde, Matter wave interferometers: a synthetic approach, and in B.
Young, M. Kasevich and S. Chu, Precision atom interferometry with light pulses, in Atom
Interferometry, P. R. Bergman (editor) (Academic Press, San Diego, 1997).
17
Chapter 2
Geometry
The basic equations of Physics are dierential equations. Now, not every
space accepts dierentials and derivatives. Every time a derivative is written
in some space, a lot of underlying structure is assumed, taken for granted. It
is supposed that that space is a dierentiable (or smooth) manifold. We shall
give in what follows a short survey of the steps leading to that concept. That
will include many other notions taken for granted, as that of coordinate,
parameter, curve, continuous, and the very idea of space.
2.1 Dierential Geometry
Physicists work with sets of numbers, provided by experiments, which they
must somehow organize. They make always implicitly a large number
of assumptions when conceiving and preparing their experiments and a few
more when interpreting them. For example, they suppose that the use of
coordinates is justied: every time they have to face a continuum set of
values, it is through coordinates that they distinguish two points from each
other. Now, not every kind of point-set accept coordinates. Those which do
accept coordinates are specically structured sets called manifolds. Roughly
speaking, manifolds are sets on which, at least around each point, everything
looks usual, that is, looks Euclidean.
2.1 Let us recall that a distance function is a function d taking any pair
(p, q) of points of a set X into the real line R and satisfying the following four
distance
function
conditions : (i) d(p, q) 0 for all pairs (p, q); (ii) d(p, q) = 0 if and only if
p = q; (iii) d(p, q) = d(q, p) for all pairs (p, q); (iv) d(p, r) + d(r, q) d(p, q)
for any three points p, q, r. It is thus a mapping d: XX R
+
. A space on
18
which a distance function is dened is a metric space. For historical reasons,
a distance function is here (and frequently) called a metric, though it would
be better to separate the two concepts (see in section 2.1.4, page 40, how a
positive-denite metric, which is a tensor eld, can dene a distance).
2.2 The Euclidean spaces are the basic spaces we shall start with. The 3-
dimensional space E
3
consists of the set R
3
of ordered triples of real numbers p
= (p
1
, p
2
, p
3
), q = (q
1
, q
2
, q
3
), etc, endowed with the distance function d(p, q)
=
_
3
i=1
(p
i
q
i
)
2

1/2
. A r-ball around p is the set of points q such that
d(p, q) < r, for r a positive number. These open balls dene a topology, that
is, a family of subsets of E
3
leading to a well-dened concept of continuity.
It was thought for much time that a topology was necessarily an ospring of
a distance function. This is not true. The modern concept, presented below
(2.7), is more abstract and does without any distance function.
euclidean
spaces
Non-relativistic elds live on space E
3
or, if we prefer, on the direct
product spacetime E
3
E
1
, with the extra E
1
accounting for time. In non-
relativistic physics space and time are independent of each other, and this is
encoded in the directproduct character: there is one distance function for
space, another for time. In relativistic theories, space and time are blended
together in an inseparable way, constituting a real spacetime. The notion of
spacetime was introduced by Poincare and Minkowski in Special Relativity.
Minkowski spacetime, to be described later, is the paradigm of every other
spacetime.
For the n-dimensional Euclidean space E
n
, the point set is the set R
n
of
ordered n-uples p = (p
1
, p
2
, ..., p
n
) of real numbers and the topology is the
ball-topology of the distance function d(p, q) = [

n
i=1
(p
i
q
i
)
2
]
1/2
. E
n
is the
basic, initially assumed space, as even dierential manifolds will be presently
dened so as to generalize it. The introduction of coordinates on a general
space S will require that S resemble some E
n
around each point.
2.3 When we say around each point, mathematicians say locally. For
example, manifolds are locally Euclidean sets. But not every set of points
can resemble, even locally, an Euclidean space. In order to do so, a point set
must have very special properties. To begin with, it must have a topology. A
set with such an underlying structure is a topological space. Manifolds are
19
topological spaces with some particular properties which make them locally
Euclidean spaces.
The procedure then runs as follows:
it is supposed that we know everything on usual Analysis, that is,
Analysis on Euclidean spaces. Structures are then progressively
added up to the point at which it becomes possible to transfer
notions from the Euclidean to general spaces. This is, as a rule,
only possible locally, in a neighborhood around each point.
2.4 We shall later on represent physical systems by elds. Such elds are
present somewhere in space and time, which are put together in a unied
spacetime. We should say what we mean by that. But there is more. Fields
are idealized objects, which we represent mathematically as members of some
other spaces. We talk about vectors, matrices, functions, etc. There will be
spaces of vectors, of matrices, of functions. And still more: we operate with
these elds. We add and multiply them, sometimes integrate them, or take
their derivatives. Each one of these operations requires, in order to have
a meaning, that the objects they act upon belong to spaces with specic
properties.
2.1.1 Spaces
2.5 Thus, rst task, it will be necessary to say what we understand by
spaces in general. Mathematicians have built up a systematic theory of
spaces, which describes and classies them in a progressive order of complex-
ity. This theory uses two primitive notions - sets, and functions from one
set to another. The elements belonging to a space may be vectors, matrices,
functions, other sets, etc, but the standard language calls simply points
the members of a generic space.
A space S is an organized set of points, a point set plus a structure.
This structure is a division of S, a convenient family of subsets. Dierent
general
notion
purposes require dierent kinds of subset families. For example, in order to
arrive at a well-dened notion of integration, a measure space is necessary,
which demands a special type of sub-division called -algebra. To make of
20
S a topological space, we decompose it in another peculiar way. The latter
will be our main interest because most spaces used in Physics are, to start
with, topological spaces.
2.6 That this is so is not evident at every moment. The customary ap-
proach is just the contrary. The physicist will implant the object he needs
without asking beforehand about the possibilities of the underlying space.
He can do that because Physics is an experimental science. He is justi-
ed in introducing an object if he obtains results conrmed by experiment.
A well-succeeded experiment brings forth evidence favoring all the assump-
tions made, explicit or not. Summing up: the additional objects (say, elds)
dened on a certain space (say, spacetime) may serve to probe into the un-
derlying structure of that space.
2.7 Topological spaces are, thus, the primary spaces. Let us begin with
them.
Given a point set S, a topology is a family T of subsets of S
topology
to which belong: (a) the whole set S and the empty set ; (b)
the intersection

k
U
k
of any nite sub-family of members U
k
of
T; (c) the union

k
U
k
of any sub-family (nite or innite) of
members.
A topological space (S, T) is a set of points S on which a
topology T is dened.
The members of the family T are, by denition, the open sets of (S, T).
Notice that a topological space is indicated by the pair (S, T). There are, in
general, many dierent possible topologies on a given point set S, and each
one will make of S a dierent topological space. Two extreme topologies
are always possible on any S. The discrete space is the topological space
(S, P(S)), with the power set P(S) the set of all subsets of S as the
topology. For each point p, the set {p} containing only p is open. The other
extreme case is the indiscrete (or trivial) topology T = {, S}.
Any subset of S containing a point p is a neighborhood of p. The comple-
ment of an open set is (by denition) a closed set. A set which is open in a
21
topology may be closed in another. It follows that and S are closed (and
open!) sets in all topologies.
Comment 2.1 The space (S, T) is connected if and S are the only sets which are
simultaneously open and closed. In this case S cannot be decomposed into the union of
two disjoint open sets (this is dierent from path-connectedness). In the discrete topology
all open sets are also closed, so that unconnectedness is extreme.
2.8 Let f : A B be a function between two topological spaces A (the
domain) and B (the target). The inverse image of a subset X of B by f is the
set f
<1>
(X) = {a A such that f(a) X}. The function f is continuous
if the inverse images of all the open sets of the target space B are open sets
of the domain space A. It is necessary to specify the topology whenever one
continuity
speaks of a continuous function. A function dened on a discrete space is
automatically continuous. On an indiscrete space, a function is hard put to
be continuous.
2.9 A topology is a metric topology when its open sets are the open balls
B
r
(p) = {q S such that d(q, p) < r} of some distance function. The
simplest example of such a ball-topology is the discrete topology P(S): it
can be obtained from the so-called discrete metric: d(p, q) = 1 if p = q, and
d(p, q) = 0 if p = q. In general, however, topologies are independent of any
distance function: the trivial topology cannot be given by any metric.
2.10 A caveat is in order here. When we say metric we mean a positive-
denite distance function as above. Physicists use the word metrics for
some invertible bilinear forms which are not positive-denite, and this prac-
tice is progressively infecting mathematicians. We shall follow this seemingly
inevitable trend, though it should be clear that only positive-denite metrics
can dene a topology. The fundamental bilinear form of relativistic Physics,
the Lorentz metric on Minkowski space-time, does not dene true distances
between points.
2.11 We have introduced Euclidean spaces E
n
in 2.2. These spaces, and
Euclidean half-spaces (or upper-spaces) E
n
+
are, at least for Physics, the
most important of all topological spaces. This is so because Physics deals
22
mostly with manifolds, and a manifold (dierentiable or not) will be a space
which can be approximated by some E
n
or E
n
+
in some neighborhood of
each point (that is, locally). The half-space E
n
+
has for point set R
n
+
=
{p = (p
1
, p
2
, ..., p
n
) R
n
such that p
n
0}. Its topology is that induced
by the ball-topology of E
n
(the open sets are the intersections of R
n
+
with
the balls of E
n
). This space is essential to the denition of manifolds-with-
boundary.
2.12 A bijective function f : A B will be a homeomorphism if it is
continuous and has a continuous inverse. It will take open sets into open sets
homeo
morphism
and its inverse will do the same. Two spaces are homeomorphic when there
exists a homeomorphism between them. A homeomorphism is an equiva-
lence relation: it establishes a complete equivalence between two topological
spaces, as it preserves all the purely topological properties. Under a home-
omorphism, images and pre-images of open sets are open, and images and
pre-images of closed sets are closed. Two homeomorphic spaces are just the
same topological space. A straight line and one branch of a hyperbola are
the same topological space. The same is true of the circle and the ellipse.
A 2-dimensional sphere S
2
can be stretched in a continuous way to become
an ellipsoid or a tetrahedron. From a purely topological point of view, these
three surfaces are indistinguishable. There is no homeomorphism, on the
other hand, between S
2
and a torus T
2
, which is a quite distinct topological
space.
Take again the Euclidean space E
n
. Any isometry (distancepreserving
mapping) will be a homeomorphism, in particular any translation. Also
homothecies with reason = 0 are homeomorphisms. From these two prop-
erties it follows that each open ball of E
n
is homeomorphic to the whole E
n
.
Suppose a space S has some open set U which is homeomorphic to an open
set (a ball) in some E
n
: there is a homeomorphic mapping : U ball,
f(p U) = x = (x
1
, x
2
, ..., x
n
). Such a local homeomorphism , with E
n
as
target space, is called a coordinate mapping and the values x
k
are coordinates coordinates
of p.
2.13 S is locally Euclidean if, for every point p S, there exists an open
set U to which p belongs, which is homeomorphic to either an open set in
23
some E
s
or an open set in some E
s
+
. The number s is the dimension of S at
the point p.
2.14 We arrive in this way at one of the concepts announced at the begin-
ning of this chapter: a (topological) manifold is a connected space on which
coordinates make sense.
A manifold is a topological space S which is manifold
(i) locally Euclidean;
(ii) has the same dimension s at all points, which is then the
dimension of S, s = dim S.
Points whose neighborhoods are homeomorphic to open sets of E
s
+
and not
to open sets of E
s
constitute the boundary S of S. Manifolds including
points of this kind are manifoldswithboundary.
The local-Euclidean character will allow the denition of coordinates and
will have the role of a complementarity principle: in the local limit, a
dierentiable manifold will look still more Euclidean than the topological
manifolds. Notice that we are indicating dimensions by m, n, s, etc, and
manifolds by the corresponding capitals: dim M = m; dim N = n, dim S =
s, etc.
2.15 Each point p on a manifold has a neighborhood U homeomorphic to
an open set in some E
n
, and so to E
n
itself. The corresponding homeomor-
phism
: U open set in E
n
will give local coordinates around p. The neighborhood U is called a co-
ordinate neighborhood of p. The pair (U, ) is a chart, or local system of
coordinates (LSC) around p.
We must be more specic. Take E
n
itself: an open neighborhood V of a
point q E
n
is homeomorphic to another open set of E
n
. Each homeomor-
phism u: V V

included in E
n
denes a system of coordinate functions
(what we usually call coordinate systems: Cartesian, polar, spherical, ellip-
tic, stereographic, etc.). Take the composite homeomorphism x: S E
n
,
x(p) = (x
1
, x
2
, ..., x
n
) = (u
1
(p), u
2
(p), ..., u
n
(p)). The functions
24
x
i
= u
i
: U E
1
will be the local coordinates around p. We shall use
the simplied notation (U, x) for the chart. Dierent systems of coordinate
functions require dierent number of charts to plot the space S. For E
2
itself, coordinates
one Cartesian system is enough to chart the whole space: V = E
2
, u = the
identity mapping. The polar system, however, requires at least two charts.
For the sphere S
2
, stereographic coordinates require only two charts, while
the cartesian system requires four.
Comment 2.2 Suppose the polar system with only one chart: E
2
R
1
+
(0, 2). Intu-
itively, close points (r, 0+) and (r, 2 ), for small, are represented by faraway points.
Technically, due to the necessity of using open sets, the whole half-line (r, 0) is absent, not
represented. Besides the chart above, it is necessary to use E
2
R
1
+
(, +2 ), with
arbitrary in the interval (0, 2 ).
Comment 2.3 Classical Physics needs coordinates to distinguish points. We see that
the method of coordinates can only work on locally Euclidean spaces.
2.16 As we have said, every time we write a derivative, a dierential, a
Laplacian we are assuming an additional underlying structure for the space
we are working on: it must be a dierentiable (or smooth) manifold. And
manifolds and smooth manifolds can be introduced by imposing progres-
sively restrictive conditions on the decomposition which has led to topologi-
cal spaces. Just as not every space accepts coordinates (that is, not not every
space is a manifold), there are spaces on which to dierentiate is impossible.
We arrive nally at the crucial notion by which knowledge on dierentiability
on Euclidean spaces is translated into knowledge on dierentiability on more
general spaces. We insist that knowledge of Analysis on Euclidean spaces is
taken for granted.
A given point p S can in principle have many dierent coordinate neigh-
borhoods and charts. Given any two charts (U, x) and (V, y) with U

V = ,
to a given point p in their intersection, p U

V , will correspond coordi-

nates x = x(p) and y = y(p). These coordinates will be related by a homeo-
morphism between open sets of E
n
,
y x
<1>
: E
n
E
n
which is a coordinate transformation, usually written y
i
= y
i
(x
1
, x
2
, . . . , x
n
).
Its inverse is x y
<1>
, written x
j
= x
j
(y
1
, y
2
, ..., y
n
). Both the coordinate
25
transformation and its inverse are functions between Euclidean spaces. If
both are C

(dierentiable to any order) as functions from E

n
into E
n
, the
two local systems of coordinates are said to be dierentially related. An atlas
on the manifold is a collection of charts {(U
a
, y
a
)} such that

a
U
a
= S.
If all the charts are dierentially related in their intersections, it will be a
dierentiable atlas.

The chain rule

i
k
=
y
i
x
j
x
j
y
k
says that both Jacobians are = 0.

An extra chart (W, x), not belonging to a dierentiable atlas A, is said

to be admissible to A if, on the intersections of W with all the coordinate-
neighborhoods of A, all the coordinate transformations from the atlas LSCs
to (W, x) are C

. If we add to a dierentiable atlas all its admissible charts,

we get a complete atlas, or maximal atlas, or C

structure. The extension

of a dierentiable atlas, obtained in this way, is unique (this is a theorem).
A topological manifold with a complete dierentiable atlas is
a dierentiable manifold.
dierentiable
manifold
2.17 A function f between two smooth manifolds is a dierentiable func-
tion (or smooth function) when, given the two atlases, there are coordinates
systems in which y f x
<1>
is dierentiable as a function between
Euclidean spaces.
2.18 A curve on a space S is a function a : I S, a : t a(t), taking the
interval I = [0, 1] E
1
into S. The variable t I is the curve parameter.
If the function a is continuous, then a is a path. If the function a is also
curves
dierentiable, we have a smooth curve.

When a(0) = a(1), a is a closed

This requirement of innite dierentiability can be reduced to k-dierentiability (to

give a C
k
atlas).

If some atlas exists on S whose Jacobians are all positive, S is orientable. When 2
dimensional, an orientable manifold has two faces. The M obius strip and the Klein bottle
are non-orientable manifolds.

The trajectory in a brownian motion is continuous (thus, a path) but is not dieren-
tiable (not smooth) at the turning points.
26
curve, or a loop, which can be alternatively dened as a function from the
circle S
1
into S. Some topological properties of a space can be grasped by
studying its possible paths.
Comment 2.4 This is the subject matter of homotopy theory. We shall need one concept
contractibility for which the notion of homotopy is an indispensable preliminary.
Let f, g : X Y be two continuous functions between the topological spaces X
and Y. They are homotopic to each other (f g) if there exists a continuous function
F : X I Y such that F(p, 0) = f(p) and F(p, 1) = g(p) for every p X. The
function F(p, t) is a one-parameter family of continuous functions interpolating between
f and g, a homotopy between f and g. Homotopy is an equivalence relation between
continuous functions and establishes also a certain equivalence between spaces. Given any
space Z, let id
Z
: Z Z be the identity mapping on Z, id
Z
(p) = p for every p Z. A
continuous function f : X Y is a homotopic equivalence between X and Y if there exists
a continuous function g : Y X such that g f id
X
and f g id
Y
. The function
g is a kind of homotopic inverse to f. When such a homotopic equivalence exists, X
and Y are homotopic. Every homeomorphism is a homotopic equivalence but not every
homotopic equivalence is a homeomorphism.
Comment 2.5 A space X is contractible if it is homotopically equivalent to a point. More
precisely, there must be a continuous function h : X I X and a constant function
f : X X, f(p) = c (a xed point) for all p X, such that h(p, 0) = p = id
X
(p) and
h(p, 1) = f(p) = c. Contractibility has important consequences in standard, 3-dimensional
vector analysis. For example, the statements that divergenceless uxes are rotational
(div v = 0 v = rot w) and irrotational uxes are potential (rot v = 0 v = grad )
are valid only on contractible spaces. These properties generalize to dierential forms (see
page 38).
2.19 We have seen that two spaces are equivalent from a purely topologi-
cal point of view when related by a homeomorphism, a topology-preserving
transformation. A similar role is played, for spaces endowed with a dier-
entiable structure, by a dieomorphism: a dieomorphism is a dierentiable
dieo
morphism
homeomorphism whose inverse is also smooth. When some dieomorphism
exists between two smooth manifolds, they are said to be dieomorphic. In
this case, besides being topologically the same, they have equivalent dier-
entiable structures. They are the same dierentiable manifold.
2.20 Linear spaces (or vector spaces) are spaces allowing for addition and
rescaling of their members. This means that we know how to add two vectors
vector
space
27
so that the result remains in the same space, and also to multiply a vector by
some number to obtain another vector, also a member of the same space. In
the cases we shall be interested in, that number will be a complex number.
In that case, we have a vector space V over the eld C of complex numbers.
Every vector space V has a dual V

, another linear space formed by all the

linear mappings taking V into C. If we indicate a vector V by the ket
|v >, a member of the dual can be indicated by the bra < u|. The latter will
be a linear mapping taking, for example, |v > into a complex number, which
we indicate by < u|v >. Being linear means that a vector a|v > + b|w > will
be taken by < u| into the complex number a < u|v > + b < u|w >. Two
linear spaces with the same nite dimension (= maximal number of linearly
independent vectors) are isomorphic. If the dimension of V is nite, V and
V

have the same dimension and are, consequently, isomorphic.

Comment 2.6 Every vector space is contractible. Many of the most remarkable proper-
ties of E
n
come from its being, besides a topological space, a vector space. E
n
itself and
any open ball of E
n
are contractible. This means that any coordinate open set, which is
homeomorphic to some such ball, is also contractible.
Comment 2.7 A vector space V can have a norm, which is a distance function and
denes consequently a certain topology called the norm topology. In this case, V is a
metric space. For instance, a norm may come from an inner product, a mapping from
the Cartesian set product V V into C, V V C, (v, u) < v, u > with suitable
properties. The number v = (| < v, v > |)
1/2
will be the norm of v V induced by
the inner product. This is a special norm, as norms can be dened independently of inner
products. When the norm comes from an inner space, we have a Hilbert space. When
not, a Banach space. When the operations (multiplication by a scalar and addition) keep
a certain coherence with the topology, we have a topological vector space.
Once in possession of the means to dene coordinates, we can proceed to
transfer to manifolds all the (supposedly wellknown) results of usual vector
and tensor analysis on Euclidean spaces. Because a manifold is equivalent
to an Euclidean space only locally, this will be possible only in a certain
neighborhood of each point. This is the basic dierence between Euclidean
spaces and general manifolds: properties which are global on the rst hold
only locally on the latter.
28
2.1.2 Vector and Tensor Fields
2.21 The best means to transfer the concepts of vectors and tensors from
Euclidean spaces to general dierentiable manifolds is through the mediation
of spaces of functions. We have talked on function spaces, such as Hilbert
spaces separable or not, and Banach spaces. It is possible to dene many
distinct spaces of functions on a given manifold M, diering from each other
by some characteristics imposed in their denitions: squareintegrability for
example, or dierent kinds of norms. By a suitable choice of conditions we
can actually arrive at a space of functions containing every information on
M. We shall not deal with such involved subjects. At least for the time
being, we shall need only spaces with poorly dened structures, such as the
space of real functions on M, which we shall indicate by R(M).
2.22 Of the many equivalent notions of a vector on E
n
, the directional vectors
derivative is the easiest to adapt to dierentiable manifolds. Consider the
set R(E
n
) of real functions on E
n
. A vector V = (v
1
, v
2
, . . . , v
n
) is a linear
operator on R(E
n
): take a point p E
n
and let f R(E
n
) be dierentiable
in a neighborhood of p. The vector V will take f into the real number
V (f) = v
1
_
f
x
1
_
p
+ v
2
_
f
x
2
_
p
+ + v
n
_
f
x
n
_
p
.
This is the directional derivative of f along the vector V at p. This action
of V on functions respects two conditions:
1. linearity: V (af + bg) = aV (f) + bV (g), a, b E
1
and f, g R(E
n
);
2. Leibniz rule: V (f g) = f V (g) + g V (f).
2.23 This conception of vector an operator acting on functions can
be dened on a dierential manifold N as follows. First, introduce a curve
through a point p N as a dierentiable curve a : (1, 1) N such that
a(0) = p (see page 27). It will be denoted by a(t), with t (1, 1). When t
varies in this interval, a 1-dimensional continuum of points is obtained on N.
In a chart (U, x) around p, these points will have coordinates a
i
(t) = x
i
(t).
29
Consider now a function f R(N). The vector V
p
tangent to the curve a(t)
at p is given by
V
p
(f) =
_
d
dt
(f a)(t)
_
t=0
=
_
dx
i
dt
_
t=0

x
i
f .
V
p
is independent of f, which is arbitrary. It is an operator V
p
: R(N) E
1
.
Now, any vector V
p
, tangent at p to some curve on N, is a tangent vector
to N at p. In the particular chart used above,
dx
k
dt
is the k-th component of
V
p
. The components are chart-dependent, but V
p
itself is not. From its very
denition, V
p
satises the conditions (1) and (2) above. A tangent vector on
N at p is just that, a mapping V
p
: R(N) E
1
which is linear and satises
the Leibniz rule.
2.24 The vectors tangent to N at p constitute a linear space, the tan-
gent space T
p
N to the manifold N at p. Given some coordinates x(p) =
(x
1
, x
2
, . . . , x
n
) around the point p, the operators {

x
i
} satisfy conditions (1)
tangent
space
and (2) above. More than that, they are linearly independent and conse-
quently constitute a basis for the linear space: any vector can be written in
the form
V
p
= V
i
p

x
i
.
The V
i
p
s are the components of V
p
in this basis. Notice that each coordinate
x
j
belongs to R(N). The basis {

x
i
} is the natural, holonomic, or coordinate
basis associated to the coordinate system {x
j
}. Any other set of n vectors
{e
i
} which are linearly independent will provide a base for T
p
N. If there is
no coordinate system {y
k
} such that e
k
=

y
k
, the base {e
i
} is anholonomic
or non-coordinate.
2.25 T
p
N and E
n
are nite vector spaces of the same dimension and are
consequently isomorphic. The tangent space to E
n
at some point will be
itself an E
n
. Euclidean spaces are dieomorphic to their own tangent spaces,
and that explains in part their simplicity in equations written on such
spaces, one can treat indices related to the space itself and to the tangent
spaces on the same footing. This cannot be done on general manifolds.
These tangent vectors are called simply vectors, or contravariant vectors.
The members of the dual cotangent space T

p
N, the linear mappings
p
: T
p
N
E
n
, are covectors, or covariant vectors.
30
2.26 Given an arbitrary basis {e
i
} of T
p
N, there exists a unique basis
{
j
} of T

p
N, its dual basis, with the property
j
(e
i
) =
j
i
. Any
p
T

p
N,
is written
p
=
p
(e
i
)
i
. Applying V
p
to the coordinates x
i
, we nd V
i
p
= V
p
(x
i
), so that V
p
= V
p
(x
i
)

x
i
= (V
p
)e
i
. The members of the basis
dual to the natural basis {

x
i
} are indicated by {dx
i
}, with dx
j
(

x
i
) =

j
i
. This notation is justied in the usual cases, and extended to general
manifolds (when f is a function between general dierentiable manifolds, df
takes vectors into vectors). The notation leads also to the reinterpretation of
the usual expression for the dierential of a function, df =
f
x
i
dx
i
, as a linear
operator:
df(V
p
) =
f
x
i
dx
i
(V
p
).
In a natural basis,

p
=
p
(

x
i
)dx
i
.
2.27 The same order of ideas can be applied to tensors in general: a tensors
tensor at a point p on a dierentiable manifold M is dened as a tensor on
T
p
M. The usual procedure to dene tensors covariant and contravariant
on Euclidean vector spaces can be applied also here. A covariant tensor of
order s, for example, is a multilinear mapping taking the Cartesian product
T
s
p
M = T
p
M T
p
M T
p
M of T
p
M by itself s-times into the set of real
numbers. A contravariant tensor of order r will be a multilinear mapping
taking the the Cartesian product T
r
p
M = T

p
M T

p
M T

p
M of T

p
M
by itself r-times into E
1
. A mixed tensor, s-times covariant and r-times
contravariant, will take the Cartesian product T
s
p
M T
r
p
M multilinearly
into E
1
. Basis for these spaces are built as the direct product of basis for the
corresponding vector and covector spaces. The whole lore of tensor algebra is
in this way transmitted to a point on a manifold. For example, a symmetric
covariant tensor of order s applies to s vectors to give a real number, and is
indierent to the exchange of any two arguments:
T(v
1
, v
2
, . . . , v
k
, . . . , v
j
, . . . , v
s
) = T(v
1
, v
2
, . . . , v
j
, . . . , v
k
, . . . , v
s
).
An antisymmetric covariant tensor of order s applies to s vectors to give a
real number, and change sign at each exchange of two arguments:
T(v
1
, v
2
, . . . , v
k
, . . . , v
j
, . . . , v
s
) = T(v
1
, v
2
, . . . , v
j
, . . . , v
k
, . . . , v
s
).
31
2.28 Because they will be of special importance, let us say a little more on
such antisymmetric covariant tensors. At each xed order, they constitute
a vector space. But the tensor product of two antisymmetric tensors
and of orders p and q is a (p + q)-tensor which is not antisymmetric,
so that the antisymmetric tensors do not constitute a subalgebra with the
tensor product.
2.29 The wedge product is introduced to recover a closed algebra. First
we dene the alternation Alt(T) of a covariant tensor T, which is an anti-
symmetric tensor given by
Alt(T)(v
1
, v
2
, . . . , v
s
) =
1
s!

(P)
(sign P)T(v
p
1
, v
p
2
, . . . , v
p
s
),
the summation taking place on all the permutations P = (p
1
, p
2
, . . . , p
s
) of
the numbers (1,2,. . . , s) and (sign P) being the parity of P. Given two
antisymmetric tensors, of order p and of order q, their exterior product,
or wedge product, indicated by , is the (p+q)-antisymmetric tensor
=
(p + q)!
p! q!
Alt( ).
With this operation, the set of antisymmetric tensors constitutes the exte-
rior algebra, or Grassmann algebra, encompassing all the vector spaces of
Grassmann
algebra
antisymmetric tensors. The following properties come from the denition:
( + ) = + ; (2.1)
( + ) = + ; (2.2)
a( ) = (a) = (a), a R; (2.3)
( ) = ( ); (2.4)
= ()

. (2.5)
In the last property,

and

are the respective orders of and . If

{
i
} is a basis for the covectors, the space of s-order antisymmetric tensors
has a basis
{
i
1

i
2

i
s
}, 1 i
1
, i
2
, . . . , i
s
dim T
p
M, (2.6)
32
in which an antisymmetric covariant s-tensor will be written
=
1
s!

i
1
i
2
...i
s

i
1

i
2

i
s
.
In a natural basis {dx
j
},
=
1
s!

i
1
i
2
...i
s
dx
i
1
dx
i
2
dx
i
s
.
2.30 Thus, a tensor at a point p M is a tensor dened on the tangent
space T
p
M. One can choose a chart around p and use for T
p
M and T

p
M the
natural bases {

x
i
} and {dx
j
}. A general tensor will be written
T = T
i
1
i
2
...i
r
j
1
j
2
...j
s

x
i
1

x
i
2


x
i
r
dx
j
1
dx
j
2
dx
j
s
.
In another chart, with natural bases {

x
i

} and (dx
j

), the same tensor will

be written
T = T
i

1
i

2
...i

r
j

1
j

2
...j

x
i

x
i

2


x
i

r
dx
j

1
dx
j

2
dx
j

s
= T
i

1
i

2
...i

r
j

1
j

2
...j

s
x
i
1
x
i

x
i
2
x
i

2

x
i
r
x
i

x
j

1
x
j
1

x
j

2
x
j
2

x
j

s
x
j
s

x
i
1

x
i
2


x
i
r
dx
j
1
dx
j
2
dx
j
s
, (2.7)
which gives the transformation of the components under changes of coordi-
nates in the charts intersection. We nd frequently tensors dened as entities
whose components transform in this way, with one Lame coecient
x
j

r
x
j
r
for
each index. It should be understood that a tensor is always a tensor with re-
spect to a given group. Just above, the group of coordinate transformations
was involved. General base transformations constitute another group.
2.31 Vectors and tensors have been dened at a xed point p of a dier-
entiable manifold M. The natural basis we have used is actually {
_

x
i

p
}. A
vector at p M has been dened as the tangent to a curve a(t) on M, with
a(0) = p. We can associate a vector to each point of the curve by allowing
the variation of the parameter t: X
a(t)
(f) =
d
dt
(f a)(t). X
a(t)
is then the
vector
elds
tangent eld to a(t), and a(t) is the integral curve of X through p. In general,
this only makes sense locally, in a neighborhood of p. When X is tangent to
a curve globally, X is a complete eld.
33
2.32 Let us, for the sake of simplicity take a neighborhood U of p and
suppose a(t) U, with coordinates (a
1
(t), a
2
(t), , a
m
(t)). Then, X
a(t)
=
da
i
dt

a
i
, and
da
i
dt
is the component X
i
a(t)
. In this sense, the eld whose integral
curve is a(t) is given by the velocity
da
dt
. Conversely, if a eld is given
by its components X
k
(x
1
(t), x
2
(t), . . . , x
m
(t)) in some natural basis, its
integral curve x(t) is obtained by solving the system of dierential equations
X
k
=
dx
k
dt
. Existence and uniqueness of solutions for such systems hold in
general only locally, as most elds exhibit singularities and are not complete.
Most manifolds accept no complete vector elds at all. Those which do are
called parallelizable. Toruses are parallelizable, but, of all the spheres S
n
,
only S
1
, S
3
and S
7
are parallelizable. S
2
is not.

2.33 At a point p, V
p
takes a function belonging to R(M) into some real
number, V
p
: R(M) R. When we allow p to vary in a coordinate neigh-
borhood, the image point will change as a function of p. By using successive
cordinate transformations and as long as singularities can be surounded, V
can be extended to M. Thus, a vector eld is a mapping V : R(M) R(M).
In this way we arrive at the formal denition of a eld:
a vector eld V on a smooth manifold M is a linear mapping V : R(M)
R(M) obeying the Leibniz rule:
X(f g) = f X(g) + g X(f), f, g R(M).
We can say that a vector eld is a dierentiable choice of a member of T
p
M
at each p of M. An analogous reasoning can be applied to arrive at tensors
elds of any kind and order.
2.34 Take now a eld X, given as X = X
i
x
i
. As X(f) R(M), another
eld as Y = Y
i
x
i
can act on X(f). The result,
Y Xf = Y
j
X
i
x
j
f
x
i
+ Y
j
X
i

2
f
x
j
x
i
,
does not belong to the tangent space because of the last term, but the com-
mutator
[X, Y ] := (XY Y X) =
_
X
i
Y
j
x
i
Y
i
X
j
x
i
_

x
j

This is the hedgehog theorem: you cannot comb a hedgehog so that all its prickles
stay at; there will be always at least one singular point, like the head crown.
34
does, and is another vector eld. The operation of commutation denes a
Lie
algebra
linear algebra. It is also easy to check that
[X, X] = 0, (2.8)
[[X, Y ], Z] + [[Z, X], Y ] + [[Y, Z], X] = 0, (2.9)
the latter being the Jacobi identity. An algebra satisfying these two condi-
tions is a Lie algebra. Thus, the vector elds on a manifold constitute, with
the operation of commutation, a Lie algebra.
2.1.3 Dierential Forms
2.35 Dierential forms

are antisymmetric covariant tensor elds on dif-

ferentiable manifolds. They are of extreme interest because of their good
behavior under mappings. A smooth mapping between M and N take dif-
ferential forms on N into dierential forms on M (yes, in that inverse order)
while preserving the operations of exterior product and exterior dierenti-
ation (to be dened below). In Physics they have acquired the status of
a new vector calculus: they allow to write most equations in an invariant
(coordinate- and frame-independent) way. The covector elds, or Pfaan
forms, or still 1-forms, provide basis for higher-order forms, obtained by ex-
terior product [see eq. (2.6)]. The exterior product, whose properties have
been given in eqs.(2.1)-(2.5), generalizes the vector product of E
3
to spaces of
any dimension and thus, through their tangent spaces, to general manifolds.
2.36 The exterior product of two members of a basis {
i
} is a 2-form,
typical member of a basis {
i

j
} for the space of 2-forms. In this basis,
a 2-form F, for instance, will be written F =
1
2
F
ij

i

j
. The basis for the
m-forms on an m-dimensional manifold has a unique member,
1

2

m
.
The nonvanishing m-forms are called volume elements of M, or volume forms.

On the subject, a beginner should start with H. Flanders, Dierential Forms, Aca-
demic Press, New York, l963; and then proceed with C. Westenholz, Dierential Forms in
Mathematical Physics, North-Holland, Amsterdam, l978; or W. L. Burke, Applied Dier-
ential Geometry, Cambridge University Press, Cambridge, l985; or still with R. Aldrovandi
and J. G. Pereira, Geometrical Physics, World Scientic, Singapore, l995.
35
2.37 The name dierential forms is misleading: most of them are not
dierentials of anything. Perhaps the most elementary form in Physics is
the mechanical work, a Pfaan form in E
3
. In a natural basis, it is written
W = F
k
dx
k
, with the components F
k
representing the force. The total work
realized in taking a particle from a point a to point b along a line is
W
ab
[] =
_

W =
_

F
k
dx
k
,
and in general depends on the chosen line. It will be path-independent only
when the force comes from a potential U as a gradient, F
k
= (grad U)
k
.
In this case W = dU, truly the dierential of a function, and W
ab
=
U(a) U(b). An integrability criterion is: W
ab
[] = 0 for any closed curve.
Work related to displacements in a non-potential force eld is a typical non-
dierential 1-form. Another well-known example is heat exchange.
2.38 In a more geometric mood, the form appearing in the integrand of the
arc length
_
x
a
ds is not the dierential of a function, as the integral obviously
depends on the trajectory from a to x, and is a multi-valued function of x.
The elementary length ds is a prototype form which is not a dierential,
despite its conventional appearance. A 1-form is exact if it is a gradient, like
= dU. Being exact is not the same as being integrable. Exact forms are
integrable, but non-exact forms may also be integrable if they are of the form
fdU.
2.39 The 0-form f has the dierential df =
f
x
i
dx
i
=
f
x
i
dx
i
, which is a
1-form. The generalization of this dierential of a function to forms of any
order is the dierential operator d with the following properties:
1. when applied to a k-form, d gives a (k+1)-form;
2. d( + ) = d + d ;
3. d( ) = (d) + ()

d(), where

is the order of ;
4. d
2
= dd 0 for any form .
36
2.40 The invariant, basis-independent denition of the dierential of a
k-form is given in terms of vector elds:
d(X
0
, X
1
, . . . , X
k
) =
k

i=0
()
i
X
i
_
(X
0
, X
1
, . . . , X
i1
,

X
i
, X
i+1
. . . , X
k
)
_
+

i<j
()
i+j
([X
i
, X
j
], X
0
, X
1
, . . . ,

X
i
, . . . ,

X
j
. . . , X
k
).(2.10)
Wherever it appears, the notation

X
n
means that X
n
is absent. From this
denition, or from the systematic use of the dening conditions, we can
obtain the rst examples of derivatives:
if f is a function (0-form), df =
i
f dx
i
(gradient) ;
if A = A
i
dx
i
is a covector (1-form), then dA =
1
2
(
i
A
j

j
A
i
) dx
i
dx
j
(rotational)
Comment 2.8 To grasp something about the meaning of
d
2
0,
which is usually called the Poincare lemma, let us examine the simplest case, a 1-form in
a natural basis: =
i
dx
i
. Its dierential is d = (d
i
) dx
i
+
i
d(dx
i
) =

i
x
j
dx
i
dx
j
=
1
2
(

i
x
j


j
x
i
) dx
i
dx
j
. If is exact, = df (in components,
i
=
i
f) then
d = d
2
f =
1
2
_

2
f
x
i
x
j


2
f
x
j
x
i
_
dx
i
dx
j
and the property d
2
f 0 is just the symmetry of the mixed second derivatives of a
function. Along the same lines, if is not exact, we can consider
d
2
=

2

i
x
j
x
k
dx
j
dx
k
dx
i
=
1
2
_

2

i
x
j
x
k


2

i
x
k
x
j
_
dx
j
dx
k
dx
i
= 0.
Thus, the condition d
2
0 comes from the equality of mixed second derivatives of the
functions
i
, related to integrability conditions.
2.41 A form such that d = 0 is said to be closed. A form which can
be written as = d for some is said to be exact.
37
Comment 2.9 It is natural to ask whether every closed form is exact. The answer, given
by the inverse Poincare lemma, is: yes, but only locally. It is yes in Euclidean spaces,
and dierentiable manifolds are locally Euclidean. Every closed form is locally exact. The
precise meaning of locally is the following: if d = 0 at the point p M, then there
exists a contractible (see below) neighborhood of p in which there exists a form (the
local integral of ) such that = d. But attention: if is another form of the same
order of and satisfying d = 0, then also = d( +). There are innite forms of which
an exact form is the dierential.
The inverse Poincare lemma gives an expression for the local integral of = d. In
order to state it, we have to introduce still another operation on forms. Given in a natural
basis the p-form
(x) =
i
1
i
2
i
3
...i
p
(x)dx
i
1
dx
i
2
dx
i
3
dx
i
p
the transgression of is the (p-1)-form
T =
p

j=1
()
j1
_
1
0
dtt
p1
x
i
j

i
1
i
2
i
3
...i
p
(tx)
dx
i
1
dx
i
2
. . . dx
i
j1
dx
i
j+1
dx
i
p
. (2.11)
Notice that, in the x-dependence of , x is replaced by (tx) in the argument. As t ranges
from 0 to 1, the variables are taken from the origin up to x. This expression is frequently
referred to as the homotopy formula.
The operation T is meaningful only in a star-shaped region, as x is linked to the origin
by the straight line tx, but can be generalized to a contractible region. Contractibility
has been dened in Comment 2.5. Consider the interval I = [0, 1]. A space or domain X
is contractible if there exists a continuous function h : XI X and a constant function
f : X X, f(p) = c (a xed point) for all p X, such that h(p, 0) = p = id
X
(p) and
h(p, 1) = f(p) = c. Intuitively, X can be continuously contracted to one of its points.
E
n
is contractible (and, consequently, any coordinate neighborhood), but spheres S
n
and
toruses T
n
are not. The limitation to the result given below comes from this strictly local
property. Well, the lemma then says that, in a contractible region, any form can be
written in the form
= dT +Td. (2.12)
When
d = 0 , (2.13)
= dT, (2.14)
so that is indeed exact and the integral looked for is just = T, always up to s such
that d = 0. Of course, the formulae above hold globally on Euclidean spaces, which are
38
contractible. The condition for a closed form to be exact on the open set V is that V
be contractible (say, a coordinate neighborhood). On a smooth manifold, every point has
an Euclidean (consequently contractible) neighborhood and the property holds at least
locally. The sphere S
2
requires at least two neighborhoods to be charted, and the lemma
holds only on each of them. The expression stating the closedness of , d = 0 becomes,
when written in components, a system of dierential equations whose integrability (i.e.,
the existence of a unique integral ) is granted locally. In vector analysis on E
3
, this
includes the already mentioned fact that an irrotational ux (dv = rot v = 0) is potential
(v = grad U = dU). If one tries to extend this from one of the S
2
neighborhoods, a
singularity inevitably turns up.
2.42 Let us nally comment on the mappings between dierential mani-
folds and the announced goodbehavior of forms. A C

function f : M N
between dierentiable manifolds M and N induces a mapping between the
tangent spaces:
f

: T
p
M T
f(p)
N.
If g is an arbitrary real function on N, g R(N), this mapping is dened by
[f

(X
p
)](g) = X
p
(g f) (2.15)
for every X
p
T
p
M and all g R(N). When M = E
m
and N = E
n
, f

is the
jacobian matrix. In the general case, f

is a homomorphism (a mapping which

preserves the algebraic structure) of vector spaces, called the dierential of f.
It is also frequently written df. When f and g are dieomorphisms, then
(f g)

X = f

X. Still more important, a dieomorphism f preserves

the commutator:
f

[X, Y ] = [f

X, f

Y ]. (2.16)
Consider now an antisymmetric s-tensor w
f(p)
on the vector space T
f(p)
N.
Then f determines a tensor on T
p
M by
(f

)
p
(v
1
, v
2
, . . . , v
s
) =
f(p)
(f

v
1
, f

v
2
, . . . , f

v
s
). (2.17)
Thus, the mapping f induces a mapping f

between the tensor spaces, work-

ing however in the inverse sense. f

is called a pull-back and f

, by extension,
pull
back
push-forward. The pull-back has some wonderful properties:
39
f

is linear;
(f g)

= g

.
f

preserves the exterior product: f

( ) = f

;
f

preserves the exterior derivative: f

(d) = d(f

).
Mappings between dierential manifolds preserve, consequently, the most
important aspects of dierential exterior algebra. The well-dened behavior
when mapped between dierent manifolds makes of the dierential forms
the most interesting of all tensors. Notice that all these properties apply also
when f is simply some dierentiable transformation mapping M into itself.
2.1.4 Metrics
2.43 The Euclidean space E
3
consists of the set of triples R
3
with the ball-
topology. The balls come from the Euclidean metric, a symmetric second-
order positive-denite tensor g whose components are, in global Cartesian
coordinates, given by g
ij
=
ij
. Thus, E
3
is R
3
plus the Euclidean metric.
We use this metric to measure lengths in our everyday life, but it happens
frequently that another metric is simultaneously at work on the same R
3
.
Suppose, for example, that the space is permeated by a medium endowed
with a point-dependent isotropic refractive index n(p). Light rays will feel
the metric g

ij
= n
2
(p)
ij
. To feel means that they will bend, acquire a
curved aspect. Fermats principle says simply that light rays will become
geodesics of the new metric, the straightest possible curve if measurements
are made using g

ij
instead of g
ij
. As long as we proceed to measurements
using only light rays, distances optical lengths will be dierent from those
given by the Euclidean metric. Suppose further that the medium is some
compressible uid, with temperature gradients and all which is necessary to
render point-dependent the derivative of the pressure with respect to the
uid density at xed entropy, c
s
=
_
p

_
S
. In that case, sound propagation
will be governed by still another metric, g

ij
=
1
c
s

ij
. Nevertheless, in both
cases we use also the Euclidean metric to make measurements, and much
of geometric optics and acoustics comes from comparing the results in both
40
metrics involved. These words are only to call attention to the fact that there
is no such a thing like the metric of a space. It happens frequently that more
than one is important in a given situation.
2.44 Bilinear forms are covariant tensors of second order. The tensor
product of two linear forms w and z is dened by (wz) (X,Y) = w(X)z(Y ).
The most fundamental bilinear form appearing in Physics is the Lorentz
metric on R
4
(see the nal of this section).
Given a basis {
j
} for the space of 1-forms, the products w
i
w
j
, with
i, j = 1, 2, . . . , m, constitute a basis for the space of covariant 2-tensors, in
terms of which a bilinear form g is written g = g
ij
w
i
w
j
. In a natural
basis, g = g
ij
dx
i
dx
j
.
A metric on a smooth manifold is a bilinear form, denoted g(X, Y ), X Y
or < X, Y >, satisfying the following conditions:
1. it is indeed bilinear:
X (Y + Z) = X Y + X Z
(X + Y ) Z = X Z + Y Z;
2. it is symmetric:
X Y = Y X;
3. it is non-singular: if X Y = 0 for every eld Y, then X = 0.
In a basis, g(X
i
, X
j
) = X
i
X
j
= g
mn

m
(X
i
)
n
(X
j
), so that
g
ij
= g
ji
= g(X
i
, X
j
) = X
i
X
j
.
It is standard notation to write simply w
i
w
j
= w
(i
w
j)
=
1
2
(w
i
w
j
+w
j
w
i
)
for the symmetric part of the bilinear basis, so that
g = g
ij
w
i
w
j
or, in a natural basis,
g = g
ij
dx
i
dx
j
.
41
2.45 Given a eld Y = Y
i
e
i
and a form z = z
j
w
j
in the dual basis,
z(Y ) = < z, Y > = z
j
Y
j
. A metric establishes a relation between vector
and covector elds: Y is said to be the contravariant image of a form z
if, for every X, g(X, Y ) = z(X). Then g
ij
Y
j
= z
i
. In this case, we write
simply z
j
= Y
j
. This is the usual role of the covariant metric, to lower
indices, taking a vector into the corresponding covector. If the mapping
Y z so dened is onto, the metric is non-degenerate. This is equivalent
to saying that the matrix (g
ij
) is invertible. A contravariant metric g can
then be introduced whose components (denoted by g
rs
) are the elements
of the matrix inverse to (g
ij
). If w and z are the covariant images of X
and Y, dened in a way inverse to the image given above, then g(w, z) =
g(X, Y ). All this denes on the spaces of vector and covector elds an internal
product (X, Y ) := (w, z) := g(X, Y ) = g(w, z). Invertible metrics are called
semi-Riemannian. Although physicists usually call them just Riemannian,
mathematicians more frequently reserve this denomination to non-degenerate
positive-denite metrics, with values in the positive real line R
+
.
As the Lorentz metric is not positive denite, it does not dene balls and
is consequently unable to provide for a topology on Minkowski space-time
(whose topology is, by the way, unknown). A Riemannian manifold is a
smooth manifold on which a Riemannian metric is dened. A theorem due
to Whitney states that it is always possible to dene at least one Riemannian
metric on an arbitrary dierentiable manifold.
A positive denite metric is presupposed in any measurement: lengths,
angles, volumes, etc. The length of a vector X is introduced as
X = (X, X)
1/2
.
A metric is indenite when X = 0 does not imply X = 0. It is the case of
Lorentz metric, which attributes zero length to vectors on the light cone.
The length of a curve : (a, b) M is then dened as
L

=
_
b
a

d
dt
dt.
Given two points p, q on a Riemannian manifold M, consider all the piecewise
dierentiable curves with (a) = p and (b) = q. The distance between
42
p and q is dened as the inmum of the lengths of all such curves between
them:
d(p, q) = inf
(t)
_
b
a

dC
dt
dt. (2.18)
In this way a metric tensor denes a distance function on M.
2.46 The metrics referred to in the introduction of this section, concerned
with simplied models for the behavior of light rays and sound waves, are
both obtained by multiplying all the components of the Euclidean metric
by a given function. A transformation like g
ij
g

ij
= f(p) g
ij
is called a
conformal transformation. In angle measurements, the metric appears in a
numerator and in a denominator and in consequence two metrics diering
by a conformal transformation will give the same angles. Conformal trans-
formations preserve the angles, or the cones.
Comment 2.10 To nd the angle made by two vector elds U and V at each point,
calculate ||U V ||
2
= ||U||
2
+ ||V ||
2
- 2 U V = U U - V V - 2 ||U||||V || cos
UV
, that is,
g

(U V )

(U V )

= g

+ g

-
_
g

_
g

cos
UV
. Then,

UV
= arccos
g

+g

(U V )

(U V )

_
g

_
g

,
which does not change if g

is replaced by g

= f(p) g

.
2.47 Geometry has had a very strong historical bond to metric. Geome-
tries have been synonymous of kinds of metric manifolds. This comes from
the impression that we measure something (say, distance from the origin)
when attributing coordinates to a point. We do not. Only homeomorphisms
are needed in the attribution, and they have nothing to do with metrics.
We hope to have made it clear that a metric on a dierentiable manifold is
chosen at convenience.
2.48 Minkowski spacetime is a 4-dimensional connected manifold on which
a certain indenite metric (the Lorentz metric) is dened. Points on space-
Minkowski
spacetime
time are called events. Being indenite, the metric denes only a pseudo-
distance function for any pair of points. Given two events x and y with Carte-
sian coordinates (x
0
, x
1
, x
2
, x
3
) and (y
0
, y
1
, y
2
, y
3
), their pseudo-distance will
43
be
s
2
=

= (x
0
y
0
)
2
(x
1
y
1
)
2
(x
2
y
2
)
2
(x
3
y
3
)
2
. (2.19)
This pseudo-distance is called the interval between x and y. Notice the
usual practice of attributing the rst place, with index zero, to the time
related coordinate. The Lorentz metric does not dene a topology, but es-
tablishes a partial ordering of the events: causality. A general spacetime will
be any dierentiable manifold S such that, at each point p, the tangent space
T
p
S is a Minkowski spacetime. This will induce on S another metric, with
the same set of signs (+,-,-,-) in the diagonalized form. A metric with that
set of signs is said to be Lorentzian.
In General Relativity, Einsteins equations determine a Lorentzian metric
which will be felt by any particle or wave travelling in spacetime. They are
non-linear second-order dierential equations for the metric, with as source
an energy-momentum density of the other elds in presence. Thus, the metric
depends on the source and on the assumed boundary conditions.
2.2 Pseudo-Riemannian Metric
2.49 Each spacetime is a 4dimensional pseudoRiemannian manifold. Its
main character is the fundamental form, or metric
g(x) = g

dx

dx

. (2.20)
This metric has signature 2. Being symmetric, the matrix g(x) = (g

) can
be diagonalized. Signature concerns the signs of the eigenvalues: it is the
numbers of eigenvalues with one sign minus the number of eigenvalues with
the opposite sign. It is important because it is an invariant under changes of
coordinates and vector bases. In the convention we shall adopt this means
that, at any selected point P, it is possible to choose coordinates {x

} in
terms of which g

takes the form

g(P) =
_
+|g
00
| 0 0 0
0 |g
11
| 0 0
0 |g
22
| 0
0 0 0 |g
33
|
_
. (2.21)
44
2.50 The rst example has been given above: it is the Lorentz metric of
Minkowski space, for which we shall use the notation
(x) =
ab
dx
a
dx
b
. (2.22)
We are using indices , , , . . . for Riemannian spacetime, and a, b, c, . . . for
Minkowski spacetime.
Minkowski space is the simplest, standard spacetime. Up to the signa-
ture, is an Euclidean space, and as such can be covered by a single, global
coordinate system. This system the cartesian system is the father of
all coordinate systems and just puts in the diagonal form
g(P) =
_
+1 0 0 0
0 1 0 0
0 1 0
0 0 0 1
_
. (2.23)
Comment 2.11 A metric is a real symmetric nonsingular bilinear form. A bilinear form
takes two vectors into a real number: g(V, U) = g

R. When symmetric, the

numbers g(V, U) and g(U, V ) are the same. A metric denes orthogonality. Two vectors
U and V are orthogonal to each other by g if g(V, U) = 0. The metric components can be
disposed in a matrix (g

). It will have, consequently, ten components on a 4dimensional

spacetime.
2.51 Given the metric, a vector eld V is timelike, spacelike or a null
vector, depending on the sign of g

:
g

_
> 0 timelike
< 0 spacelike
=0 null
(2.24)
V

are the components of the vector V in the coordinate system{x

}. Higher
indices indicate a contravariant vector. The metric can be used to lower
indices. From V , obtain a covariant vector, or covector, whose components
are V

= g

. Einsteins summation convention is being used: repeated

higherlower indices are summed over everytime they appear, without the
summation symbol.
2.52 As presented above, with lower indices, the metric is sometimes called
the contravariant metric. The elements of its inverse matrix are indicated by
g

. Thus, with the Kronecker delta

(which is = 1 when = and zero

otherwise), we have
g

= g

.
45
The set {g

} is then called the covariant metric, and can be used to raise

indices, or to get a contravariant vector from a covariant vector: V

= g

.
The same holds for indices of general tensors. Frequently used notations are
g = |g| = det(g

). Of course, det(g

) = g
1
.
2.53 The norm ds of the innitesimal displacement dx

, whose square is
ds
2
= g

dx

dx

(2.25)
is the innitesimal interval, or simply interval. Given a curve with extreme interval
points a and b, the integral
L[] =
_

ds =
_
b
a
ds (2.26)
along is a function(al) of , called its length.
Comment 2.12 Again a name used by extension of the strictly Riemannian case, in
which this integrals is a true length. A strictly Riemannian metric determines a true
distance between two points. In the case above, the distance between a and b would be
the inmum of L[], all curves considered.
Comment 2.13 The determinant of a metric matrix like (2.23) is always negative. In
cartesian coordinates, an integration over 4-space has the form
_
V
4
d
4
x.
In another coordinate system, a Jacobian turns up. We recall that what appears in
integration measures is the Jacobian up to the sign. Integration using a coordinate system
in which the innitesimal length takes the form (2.25), and in which the metric has a
negative determinant, is given by the expression
_
V
4

g d
4
x,
which holds in every case.
2.3 The Notion of Connection
Besides wellbehaved entities like tensors (including metrics, vectors and
functions), a manifold contains other, not so well-behaved objects. The most
important are connections, essential to the notion of parallelism. We proceed
now to present the physicits approach to connections.
46
2.54 What we understand by good behavior is: covariance under change
of coordinates. A scalar eld is invariant under change of coordinates. Take
the next case in complexity, a vector eld V . It will have components V

in
a coordinate system {x

} and components V

in a coordinate system {y

}.
The two sets are related by
V

=
y

. (2.27)
This is the standard behavior. It denes a vector by the group of coordinate
transformations. Tensors of any order reproduce it, index by index. The
metric, for example, has its components changed according to
g

=
x

. (2.28)
Notice by the way that, contracting (2.27) with the gradient operator

y

,
V

= V

= V

.
Thus, the expression
V = V

(2.29)
is invariant under change of coordinates. The vector eld V , once conceived
as such a directional derivative, is an invariant concept. This is the notion of
vector eld used by the mathematicians: a directional derivative acting on
the functions dened on the manifold.
Take now the derivative of (2.27):

=
y

+

y

_
y

_
V

=
y

+
x

_
y

_
V

=
y

+
x

2
y

=
x

_
y

+
y

2
y

_
.
47

=
x

_

x

+
x

2
y

_
. (2.30)
If alone, the rst term in the righthand side would confer to the derivative a
good tensor status. The second term breaks that behavior: the derivative of
a vector is not a tensor. In other words, the derivative is not covariant. On
a manifold, it is impossible to tell, for example, whether a vector is constant
or not.
The solution is to change the very denition of derivative by adding an-
other structure. We add to each derivative an extra term involving a new
object :
covariant
derivative

D
Dy

=

y

D
Dx

=

x

(2.31)
and impose good behavior of the modied derivative:
D
Dy

=
x

D
Dx

,
or

=
x

_

x

_
. (2.32)
We then compare with (2.30) and look for conditions on the object . These
conditions x the behavior of under coordinate transformations. must
transform according to

=
x

+
x

2
y

_
,
or

=
y

+
x

2
y

. (2.33)
This noncovariant behavior of the connection makes of (2.31) a well
behaved, covariant derivative. We have used a vector eld to nd how
should behave but, once is known, covariant derivatives can be dened on
general tensors.
48
There are actually innite objects satisfying conditions (2.33), that is,
there are innite connections. Take another one,

. It is immediate that
the dierence

is a tensor under coordinate transformations.

The covariant derivative of a function (tensor of zero degree) is the usual
derivative, which in the case is automatically covariant. Take a third order
mixed tensor T. Its covariant derivative will be given by
D

. (2.34)
The rules to calculate, involving terms with contractions for each original
index, are fairly illustrated in this example. Notice the signs: positive for
upper indices, negative for lower indices. The metric tensor, in particular,
will have the covariant derivative
D

. (2.35)
When the covariant derivative of T is zero on a domain, T is selfparallel
parallel
transport
on the domain, or paralleltransported. An intuitive view of this notion will
be given soon (see below, Fig.2.1, page 51). It exactly translates to curved
space the idea of a straight line as a curve with maintains its direction along
all its length.
If the metric is paralleltransported, the equation above gives the metric-
ity condition

= 2
()
, (2.36)
where the symbol with lowered index is dened by

= g

and the
compact notation for the symmetrized part

()
=
1
2
{

}, (2.37)
has been introduced. The analogous notation for the antisymmetrized part

[]
=
1
2
{

} (2.38)
is also very useful.
Another convenient notational device: it it usual to indicate the common
derivative by a comma. We shall adopt also one of the two current notations
49
for the covariant derivative, the semi-colon. The metricity condition, for
example, will have the expressions
g
;
= g
,
2
()
= 0. (2.39)
The bar notation for the covariant derivative, g
|
= g
;
, is also found in
the literature.
2.4 The LeviCivita Connection
There are, actually, innite connections on a manifold, innite objects behav-
ing according to (2.33). And, given a metric, there are innite connections
satisfying the metricity condition. One of them, however, is special. It is
given by

=
1
2
g

} . (2.40)
It is the single connection satisfying metricity and which is symmetric in
the last two indices. This symmetry has a deep meaning. The torsion of a
connection of components

is a tensor T of components
T

= 2

[]
. (2.41)
Connection (2.40) is called the LeviCivita connection. Its components are
just the Christoel symbols we have met in 1.15. It has, as said, a special
relationship to the metric and is the only metricpreserving connection with
zero torsion. Standard General Relativity works only with such connections.
We can now give a clear image of parallel-transport. If the connection is
related to a denite-positive metric so that angles can be measured a
parallel-transported vector keeps the angle with the curve the same all along
(Fig.2.1).
Notice that, with the notations introduced above, a general connection
will have components of the form

()
2 T

. (2.42)
Connections with zero torsion,

()
are usually called symmetric
connections.
We have said that a curve on a manifold M is a mapping from the real
line on M. The mapping will be continuous, dierentiable, etc, depending
50
X
V

X
X
V
V

Figure 2.1: If the connection is related to a denite-positive metric, a parallel-

transported vector keeps the angle with the curve all along.
on the point of interest. We shall be mostly interested in curves with xed
curves
again
endpoints. A curve with initial endpoint p
0
and nal endpoint p
1
is better
dened as a mapping from the closed interval I = [0, 1] into M:
: [0, 1] M ; u (u) ; (0) = p
0
, (1) = p
1
. (2.43)
The variable u [0, 1] is the curve parameter. A loop, or closed curve, with
p
0
= p
1
, can be alternatively dened as a mapping of the circle S
1
into M. We
shall be primarily concerned with dierentiable, or smooth curves, dened
as those for which the above mapping is a dierentiable function. A smooth
curve has always a tangent vector at each point. Suppose the tangent vectors
at all the points is time-like (see Eq. (2.24)). In that case the curve itself is
said to be time-like. In the same way, the curve is said to be space-like or
light-like when the tangent vectors are all respectively space-like and null.
We can now consider derivatives along a curve , whose points have co-
ordinates x

(u). The usual derivative along is dened as

d
du
=
dx

du

The covariant derivative along also called absolute derivative is dened

by
D
Du
=
dx

du
D

, (2.44)
and apply to tensors just in the same way as covariant derivatives. Thus, for
example,
51
DV

Du
=
dV

du
+

dx

du
. (2.45)
Let us go back to the invariant expression of a vector eld, Eq.(2.29).
It is an operator acting on the functions dened on the manifold. If the
manifold is a dierentiable manifold, and {x

} is a coordinate system, the

set of operators {

x

} is linearly independent, and

x

. We say then
that the set of derivative operators {

x

} constitute a base for vector elds.

Such a base is called natural, or coordinate base, as it is closely related to the
coordinate system {x

}. Any vector eld can be expressed as in Eq.(2.29),

with components V

. Under a change of coordinate system,

x

=
y

.
The base member e

=

x

is in this way expressed in terms of another base,

the set {

y

}. The latter constitutes another coordinate base, naturally

related to the coordinate system {y

}. The components V

transform in the
converse way, so that V is, as already said, an invariant object. Vector bases
on 4dimensional spacetime are usually called tetrads (also vierbeine, and
fourlegs).
Tetrad elds are actually much more general they need not be related
to any coordinate system. Put h

=
y

, so that e

=

x

= h

. Of
course, the members of the base {e

} commute with each other, [e

, e

] = 0.
This is also the sucient the condition for a base to be natural in some
coordinate system: if [e

, e

] = 0, then there exists a coordinate system {x

}
such that e

=

x

for = 0, 1, 2, 3. These things are simpler in the socalled

dual formalism, which uses dierential forms.
As said in 2.20, to every real vector space V corresponds another vector
space V

, which is its dual. This V

is dened as the set of linear mappings

dual
space
again
from V into the real line R. Given a base {e

} of V, there exists always

a base {e

} of V

which is its dual, by which we mean that e

(e

) =

.
The base dual to {

x

} is the set of dierential forms {dx

}, each dx

being
understood as a linear mapping satisfying dx

) =

.
Take the dierential of a function f: df =
f
x

dx

. Applying the above

rule, it follows that df(

x

) =
f
x

. An arbitrary dierential 1form (also

called a covector eld) is written =

dx

in base {dx

}. The above dual

base will have members e

= dx

= h

dy

, with h

=
x

. As it always
happen that d
2
0, the condition for a covector base to be naturally related
to a coordinate system is that de

= 0. Each dx

transforms according to
dx

=
x

dy

, and the above expression for a covector is invariant. We

shall later examine general tetrads in detail.
52
2.5 Curvature Tensor
2.55 A connection denes covariant derivatives of general tensorial ob-
jects. It goes actually a little beyond tensors. A connection denes a
covariant derivative of itself. This gives, rather surprisingly, a tensor, the
Riemann curvature tensor of the connection:
R

. (2.46)
It is important to notice the position of the indices in this denition. Au-
thors dier in that point, and these dierences can lead to dierences in the
signs (for example, in the scalar curvature dened below). We are using
all along notations consistent with the dierential forms. There is a clear
antisymmetry in the last two indices,
R

= R

.
2.56 Notice that what exists is the curvature of a connection. Many con-
nections are dened on a given space, each one with its curvature. It is
common language to speak of the curvature of space, but this only makes
sense if a certain connection is assumed to be included in the very denition
of that space.
2.57 The above formulas hold on spaces of any dimension. The meaning
of the curvature tensor can be understood from the diagram of Figure 2.2.
First, build an innitesimal parallelogram formed by pieces of geodesics,
indicated by dx

, dx

, dx

and dx

.
Take a vector eld X with components X

at a corner point P. Parallel-

transport X along dx

. At its extremity, X will have the value X

dx

Parallel-transport that value along dx

. It will lead to a value which

we denote X

.
Back at the starting point, parallel-transport X now rst along dx

and then along dx

. This will lead to a eld value which we denote

X

.
53
In a at case, X

= X

. On a curved case, the dierence between

them is non-vanishing, and given by
X

= X

= R

dx

dx

. (2.47)
This is the innitesimal case, in the limit of vanishing parallelogram
and with the value of R

at the left-lower corner.

dx

dx

dx

dx'

dx'
Figure 2.2: The meaning of the Riemann tensor.
2.58 Other tensors can be obtained from the Riemann curvature tensor
by contraction. The most important is the Ricci tensor
R

= R

. (2.48)
This tensor is symmetric in the case of the Levi-Civita connection:

R
=

R
. (2.49)
In this case, which has a special relation to the metric, the contraction with
it gives the scalar curvature

R
= g

R
. (2.50)
54
2.6 Bianchi Identities
2.59 Take the denition (2.46) in the case of the Levi-Civita connection,

. (2.51)
The metric can be used to lower the index . Calculation shows that

R
= g

, (2.52)
where

= g

. The curvature of a Levi-Civita connection has some

special symmetries in the indices, which can be obtained from the detailed
expression in terms of the metric:

R
=

R
=

R
=

R
; (2.53)

R
=

R
. (2.54)
In consequence of these symmetries, the Ricci tensor (2.48) is essentially the
only contracted second-order tensor obtained from the Riemann tensor, and
the scalar curvature (2.50) is essentially the only scalar.
2.60 A detailed calculation gives the simplest way to exhibit curvature.
Consider a vector eld U with components U

and take twice the covari-

ant derivative, getting U

;;
. Reverse then the order to obtain U

;;
and
compare. The result is
U

;;
U

;;
=

. (2.55)
Curvature turns up in the commutator of two covariant derivatives.
2.61 Detailed calculations lead also to some identities. Of of them is

= 0 . (2.56)
Another is

R;
+

R;
+

R;
= 0 (2.57)
55
(notice, in both cases, the cyclic rotation of the three last indices). The last
expression is called the Bianchi identity. As the metric has zero covariant
derivative, it can be inserted in this identity to contract indices in a conve-
nient way. Contracting with g

, it comes out

R;

R;
+

;
= 0.
Further contraction with g

yields

R;

;
= 0,
which is the same as

R;
2

;
= 0
or
_

1
2

R
_
;
= 0. (2.58)
This expression is the contracted Bianchi identity. The tensor thus co-
variantly conserved will have an important role. Its totally covariant form,
G

R

1
2
g

R
, (2.59)
is called the Einstein tensor. Its contraction with the metric gives the scalar
curvature (up to a sign).
g

R
. (2.60)
2.62 When the Ricci tensor is related to the metric tensor by

R
= g

, (2.61)
where is a constant, It is usual to say that we have an Einstein space. In
that case,

R
= 4 and G

= g

. Spaces in which

R
is a constant are
said to be spaces of constant curvature. This is the standard language. We
insist that there is no such a thing as the curvature of space. Curvature is a
characteristic of a connection, and many connections are dened on a given
space.
Some formulas above hold only in a four-dimensional space. We shall in
the following give some lower-dimensional examples. It should be kept in
mind that, in a d-dimensional space, g

= d. When d = 2 for example,

g

0. On a two-dimensional Einstein space, G

0.
56
2.6.1 Examples
Calculations of the above objects will be necessary to write the dynamical
equations of General Relativity. Such calculations are rather tiresome. In
order to get a feeling, let us examine some low-dimensional examples. We
shall look at the 2-dimensional sphere and at two 2-dimensional hyperboloids.
In eect, two surfaces of revolution can be got from the hyperbola shown in
Figure 2.3: one by rotating around the vertical axis, the other by rotating
around the horizontal axis (Figure 2.4). The rst will have two separate
sheets, the other only one. They are obtained from each other by exchanging
the squared vertical coordinate (below, z
2
) by the squared coordinates of the
resulting surface (below, x
2
+ y
2
). The spacetimes of de Sitter are higher-
dimensional versions of these hyperbolic surfaces.
Figure 2.3: Hyperbola. Two distinct surfaces of revolution can be obtained,
by rotating either around the vertical or the horizontal axis.
2.63 The sphere S
2
is dened as the set of points of E
3
satisfying x
2
+
y
2
+z
2
= a
2
in cartesian coordinates. This means z =
_
a
2
x
2
y
2
, and
consequently
dz =
dxx
_
a
2
x
2
y
2

dy y
_
a
2
x
2
y
2
.
57
H
H
2
(1,1)
H
Figure 2.4: The two hyperbolic surfaces. The second is the complement of
the rst, rotated of 90
o
.
The interval dl
2
= dx
2
+ dy
2
+ dz
2
becomes then
dl
2
=
a
2
dx
2
+ a
2
dy
2
dy
2
x
2
+ 2 dxdy xy dx
2
y
2
a
2
x
2
y
2
It is convenient to change to spherical coordinates
x = a sin cos ; y = a sin sin ; z = a cos .
The interval becomes
dl
2
= a
2
_
d
2
+ sin
2
d
2
_
.
The corresponding metric is given by
g = (g

) =
_
a
2
0
0 a
2
sin
2

_
,
with obvious inverse. The only non-vanishing Christoel symbols are

1
22
=
- sin cos and

2
12
=

2
21
= cot . The only non-vanishing components of
58
the Riemann tensor are (up to symmetries in the indices)

R
1
212
= sin
2
and

R
2
112
= 1. The Ricci tensor has

R11
= 1 and

R22
= sin
2
, so that it can
be represented by the matrix
(

R
) =
_
1 0
0 sin
2

_
=
1
a
2
(g

) .
Consequently, the sphere is an Einstein space, with =
1
a
2
. Finally, the
scalar curvature is

R
=
2
a
2
.
Not surprisingly, a sphere is a space of constant curvature. As previously said,
the Einstein tensor is of scarce interest for two-dimensional spaces. Though
we shall not examine the geodesics, we note that their equations are

= sin cos

2
;

= 2 cot

,
and that constant and are obvious particular solutions. Actually, the
solutions are the great arcs.
2.64 The two-sheeted hyperboloid H
2
is dened as the locus of those
points of E
3
satisfying x
2
+y
2
z
2
= a
2
in cartesian coordinates. The line
element has the form
dx
2
+ dy
2
dz
2
=
a
2
dx
2
+ a
2
dy
2
dy
2
x
2
+ 2 dxdy xy dx
2
y
2
a
2
x
2
y
2
.
Changing to coordinates
x = a sinh cos ; y = a sinh sin ; z = a cosh
the interval becomes
dl
2
= a
2
_
d
2
+ sinh
2
d
2
_
.
The metric will be
g = (g

) =
_
a
2
0
0 a
2
sinh
2

_
.
59
The only non-vanishing Christoel symbols are

1
22
= - sinh cosh and

2
12
=

2
21
= coth . The only non-vanishing components of the Riemann tensor
are (up to symmetries in the indices)

R
1
212
= sinh
2
and

R
2
112
= 1. The
Ricci tensor constitutes the matrix
(

R
) =
_
1 0
0 sinh
2

_
.
We see that H
2
is an Einstein space, with =
1
a
2
. The scalar curvature is

R
=
2
a
2
.
H
2
is a space of constant, but negative, curvature. It is similar to the sphere,
but with a imaginary angle. An old name for it is pseudo-sphere.
2.65 The one-sheeted hyperboloid H
(1,1)
is dened as the locus of those
points of E
3
satisfying x
2
+ y
2
z
2
= a
2
in cartesian coordinates. The line
element has the form
dz
2
dx
2
dy
2
=
a
2
dx
2
+ a
2
dy
2
dy
2
x
2
+ 2 dxdy xy dx
2
y
2
a
2
x
2
y
2
.
Changing to coordinates
x = a cosh cos ; y = a cosh sin ; z = a sinh
the interval becomes
dl
2
= a
2
_
d
2
+ sinh
2
d
2
_
.
The metric will be
g = (g

) =
_
a
2
0
0 a
2
cosh
2

_
.
The only non-vanishing Christoel symbols are

1
22
= sinh cosh and

2
12
=

2
21
= tanh . The only non-vanishing components of the Riemann tensor
are (up to symmetries in the indices)

R
1
212
= cosh
2
and

R
2
112
= 1. The
Ricci tensor constitutes the matrix
(

R
) =
_
1 0
0 cosh
2

_
.
60
We see that H
(1,1)
is an Einstein space, with =
1
a
2
. The scalar curvature is

R
=
2
a
2
.
H
(1,1)
is a space of constant positive curvature, like the sphere.
2.66 The rotating disk We shall now consider the rotating disk of 1.13
from two points of view: in (2 + 1) dimensions (two for the rotation plane,
one for time); and only the 3-dimensional space section. Curiously, the two
cases will have quite dierent results: the (2, 1) case gives zero curvature,
while the 3-dimensional space is curved.
Take rst the (2, 1) case, with coordinates (ct, R, ). The metric will be
g = (g

) =
_
_
_
1

2
R
2
c
2
0
R
2
c
0 1 0

R
2
c
0 R
2
_
_
_
.
The only non-vanishing Christoel symbols are

2
11
=

2
R
c
2
,

2
13
=

2
31
=
R
c
,

2
33
= R,

3
12
=

3
21
=

cR
, and

3
23
=

3
32
=
1
R
. All the com-
ponents of the Riemman tensor vanish, so that the 3-dimensional spacetime
is at.
Take now the 3-dimensional space, with coordinates (x, y, z). The metric
will be
g = (g

) =
_
_
_
1 + fy
2
fxy 0
fxy 1 + fx
2
0
0 0 1
_
_
_
,
with f =

2
c
2
c
2

2
(x
2
+y
2
)
2
. All the Christoel symbols of type

3
ij
are zero, but
the other have some rather lengthy expressions. The same is true of the Ricci
tensor. We shall only quote the scalar curvature, which is (with r
2
= x
2
+y
2
)
R =
2 c
4
r
4

6
2 c
8

2
(3 + 2 r
2

2
) + c
6
(4 r
2

4
2 r
4

6
)
(c
2
r
2

2
)
2
(r
2

2
c
2
(1 + r
2

2
))
2
.
This means that the space is curved. A simpler expression turns up for the
particular values = 1, c = 1:
R =
6
1 r
2
.
61
We see that there is a singularity when r c. Exactly the same results
come out if we consider the pure 2-dimensional case, the plane rotating disk.
This corresponds to dropping the last column and row of the metric above.
2.67 In all these examples,
1. a set point S is rst dened by a constraint on the points of an ambient
space E;
2. the metric is dened by the restriction, on the subset, of the metric of
the ambient space;
3. such a metric is said to be induced by the imbedding of S in E.
62
Chapter 3
Dynamics
3.1 Geodesics
3.1 Curves dened on a manifold provide tests for many of its properties.
In Physics, they do still more: any observer will be ultimately represented
by some special curve on spacetime. We have obtained the geodesic equation
before. We had then used implicitly the assumption that ds = 0. The ap-
proach which follows

though more involved, has the advantage of including

the null geodesics, for which ds = 0.
Our aim is to obtain certain priviledged curves as extremals of certain
functionals, or functions dened on the space of curves. Such functional
involve integrals along each curve, like
F[] =
_

f .
A simple stratagem allows to manipulate such functionals as if they were
functions. To do it, it is necessary to give a label to each curve, so that
varying the curve becomes a simple change of label. This can be done by
considering, beyond the initial family of possible, curves, another family of
transversal curves.
3.2 We shall be interested in families of curves, like the curves , and
in Figure 3.1. We have indicated by a
0
and a
1
the initial and the nal
congru
ences

See J. L. Synge, Relativity: The general theory, NorthHolland, Amsterdam, 1960.

63

u
v
u
1
u
0
a
1
b
1
c
1
a
0
b
0
c
0
Figure 3.1: The family of curves , , is parametrized by u. The crossed
family of curves , is parametrized by v. The second family is a variation
of the rst.
endpoints of the piece of which will be of interest, and analogously for the
other curves. The curve parameter u is indicated as going along them,
with u
0
and u
1
the initial and nal endpoints. We say that the curves form a
congruence. The variation leading from to and to can be represented
by another parameter v. It is useful to consider v as the parameter related to
another set of curves, such that each curve of the second congruence intersects
every curve of the rst, and vice versa. We have drawn only and , which
go through the endpoints of the rst family of curves. The points on all the
curves constitute a double continuum, a twoparameter domain which we
shall indicate by (u, v). It will have coordinates x

(u, v) =

(u, v). For a

xed value of v,
v
(u) = (u, v = constant) will describe a curve of the rst
family. For a xed value of u,
u
(v) = (u= constant, v) will describe a curve
of the second family. The second family is, by the way, called a variation
of any member of the rst family. There will be vectors along both families,
such as
U

=
dx

du
and V

=
dx

dv
.
64
If the connection is symmetric, the crossed absolute derivatives coincide:
DV

Du
=
DU

Dv
. (3.1)
Let us rst x v at some value, thereby choosing a xed curve of the rst
family and consider, along that curve, the function(al)
I[v] =
1
2
(u
1
u
0
)
_
u
1
u
0
g

dx

du
dx

du
du =
1
2
(u
1
u
0
)
_
u
1
u
0
g

du .
(3.2)
On the two-parameter domain, variations of that curve become simple deriva-
tives with respect to v. By conveniently adding antisymmetric and symmetric
pieces, we have the series of steps
d
dv
I[v] = (u
1
u
0
)
_
u
1
u
0
g

D
Dv
U

du = (u
1
u
0
)
_
u
1
u
0
g

DV

Du
du
= (u
1
u
0
)
_
u
1
u
0
d
du
_
g

du (u
1
u
0
)
_
u
1
u
0
g

D
Du
U

du
= (u
1
u
0
)
_
g

u
1
u
0
du (u
1
u
0
)
_
u
1
u
0
g

D
Du
U

du (3.3)
We are interested in curves with the same endpoints, and shall now collapse
the parameter v. In Figure 3.1, a
0
= b
0
= c
0
and a
1
= b
1
= c
1
. This
corresponds to putting V

= 0 at the endpoints. The rst contribution in

the righthand side above vanishes and
d
dv
I[v] = (u
1
u
0
)
_
u
1
u
0
g

D
Du
U

du . (3.4)
We now call geodesics the curves with xed endpoints for which I is station-
ary. As in the last integrand V

is arbitrary except at the endpoints, such

curves must satisfy
D
Du
U

=
dU

du
+

= 0 , (3.5)
which is the same as
geodesic
equation
d
2
x

du
2
+

dx

du
dx

du
= 0. (3.6)
65
This is the geodesic equation we have met before. It says that the vector
U

=
dx

du
has vanishing absolute derivative along the curve. This means
that it is parallel-transported along the curve. The only direction that can
be attributed to a curve is, at each point, that of its velocity U

. This
leads to a much better name for such a curve: self-parallel curve. It keeps
the same direction all along. Geodesics (solutions of the geodesic equation)
play on curved spaces the role the straight lines have on at spaces.
Comment 3.1 Why is the name self-parallel better ? There are at least two reasons:
1. the word geodesic has a very strong metric conotation; its original meaning was
that of a shortest length curve; but length, real length, only is dened by a positive-
denite metric, not our case;
2. the concept has actually no relation to metric at all; such curves can be dened for
any connection; connection is a concept quite independent of metric, and denes
parallelism; for general connections, only self-parallel makes sense.
3.3 The geodesic equation has the rst integral
g

= C. (3.7)
This is to say that, along the curve,
d
du
_
g

= 0.
Comment 3.2 Prove it for the connection given by (2.40).
3.4 Comparison with the interval expression (2.25) shows that
ds = C
1/2
du.
Eq.(3.6) is invariant under parameter changes of type
u u

= au + b . (3.8)
The constant C can, consequently, be rescaled to have only the values 1 and
0. This means that, unless C = 0, the interval parameter s can be used as
the curve parameter. That parameter is the proper time, and we recall that
66
U

=
dx

ds
is the the fourvelocity. The choice of C = 1 allows contact to be
made with Special Relativity, as (3.7) becomes
U
2
= U U = U

= g

= 1. (3.9)
The case C = 0 includes trajectories of particles with vanishing masses, in
special lightrays. As long as we keep in mind this exception, we can rewrite
the geodesic equation in the forms
D
Ds
U

=
dU

ds
+

=
d
2
x

ds
2
+

dx

ds
dx

ds
= 0. (3.10)
Once things have been interpreted in this way, we say that the velocity
is covariantly derived along . This is supported by Eq.(2.44), which now
shows the absolute derivative as the covariant derivative projected, at each
point, on the velocity.
3.5 In the case of massive particles, we can use the freedom allowed by
(3.8) to choose another parameter. Introduce s by
u u
0
= (u
1
u
0
)
1/2
s/L .
Then, comparison with(2.26) shows that I(v) =
1
2
L
2
. This leads to the
variational principle
L =
_
ds = 0 , (3.11)
which we have found before.
3.6 If we look back at what has been done in 1.15, we see that we have
there got (3.10) from (3.11). A massive particle, without additional structure
(for instance, supposing that the eect of its spin is negligible, or zero) will
follow the geodesic equation. That is actually the standard approach. It is
enough to replace, in all the discussion, the expression in the presence of a
metric by the expression in the presence of a gravitational eld. It holds,
wee see now, for massive particles, for which ds = 0.
67
3.7 Eqs.(3.6) and (3.10) are secondorder ordinary dierential equations.
Existence and unicity theorems of the theory of dierential equations state
that, given starting coordinates x
i
and velocities
dx
i
du
at a point P, there
will be a curve (u), with < u < for some > 0, which goes through P
((P) = 0) and which is unique.
3.8 The expression
A

=
dU

ds
+

(3.12)
is the covariant acceleration, the only to have a meaning in an arbitrary
coordinate system. The existence of the above rst integral is equivalent to
U A = g

= 0 . (3.13)
As in Special Relativity the acceleration is, at each point of , orthogonal
to the velocity. As the velocity is parallel (or tangent) to at each point,
we can say that the acceleration is orthogonal to . A curve is selfparallel
when its acceleration vanishes.
3.9 We have above supposed a LeviCivita connection. All the termi-
nology comes from the rst historical case, the LeviCivita connection of a
strictly Riemannian metric. We have said that any vector which is parallel
transported along a curve keeps the same angle with respect to the velocity
all along it. The concept of selfparallel curve keeps its meaning for a gen-
eral linear connection. A selfparallel curve is in that case dened as a curve
satisfying (3.10).
3.10 Let us go back to the rst integral given in Eq.(3.7). It has a very
deep meaning. As the momentum of a particle of mass m is P

= mc U

,
we rewrite it as
m
2
c
2
g

= g

= m
2
c
2
C . (3.14)
Write now
mc g

S. (3.15)
68
This means dS = mc g

dx

. This gradient form, applied to the tangent

vector U = U

x

, gives dS(U) = mc C. That reveals the meaning of S for

the present case of a massive particle: with the choice C = 1 announced in
3.4, it is the action written in terms of the sole coordinates along the curve,
as it appears in Hamilton-Jacobi theory:
P

S . (3.16)
The geodesic curve is, at each point, orthogonal to the surface S = con-
stant passing through that point. In eect, dS = mc g

dx

is (up to
the sign) just the action we have used before, but here along the curve. The
expression
Hamilton
Jacobi
equation
g

S = m
2
c
2
(3.17)
is the relativistic Hamilton-Jacobi equation for the free particle. It is possible
to recover the geodesic equation from it. Let us see how.
Consider both sides of the expression
d
du
[g

] =
d
du
U

= U

.
The left-hand side (LHS) is
d
du
(g

) = g

d
du
U

+ U

_
d
du
g

_
= g

d
du
U

+ U

.
The last piece is symmetric in the indices and , so that we can write
LHS =
d
du
(g

) = g

d
du
U

+
1
2
U

) .
Now to the right-hand side:
RHS = U

= U

S = U

S = U

= U

,
using

[g

] =

[constant] = 0 in the last step. Notice that here we

have the only contribution of Hamilton-Jacobi theory: we have used the fact
(3.16) that the momentum is a derivative to justify the exchange

. Now,

[g

] = 0 is also
0 = (

) U

+ 2 g

=
1
2
(

) U

,
69
so that
RHS = U

=
1
2
(

) U

.
Now, LHS = RHS gives
g

d
du
U

+
1
2
(

) U

= 0 ,
or
d
du
U

+
1
2
g

] U

= 0 ,
the announced result. All this can be repeated step by step, but starting
from
eikonal
equation
g

S = 0 . (3.18)
This equation has a dierent meaning: it is the eikonal equation, with S now
in the role of the eikonal. The geodesic, in this case, is the light-ray equation.
The Hamilton-Jacobi and the eikonal equations allow thus a unied view
of particle trajectories and light rays. The geodesic equation comes
from the Hamilton-Jacobi equation, as the equation of motion of a
massive particle;
from the eikonal equation, as the trajectory of a light ray.
3.11 As we have said, curves are of fundamental importance. They not
only allow testing many properties of a given space. In spacetime, every
(ideal) observer is ultimately a timelike curve.
observer
The nub of the equivalence principle is the concept of observer:
An observer is a timelike curve on spacetime, a worldline.
Such a curve represents a point-like object in 3-space, evolving in the time-
like 4-th direction. An object extended in 3-space would be necessarily
represented by a bunch of worldlines, one for each one of its points. This
mesh of curves will be necessary if, for example, the observer wishes to do
some experiment. For the time being, let us take the simplifying assumption
above, and consider only one worldline. This is an ideal, point-like observer.
If free from external forces, this line will be a geodesic.
70
And here comes the crucial point. Given a geodesic going through a
point P ((0) = P), there is always a very special system of coordinates
(Riemannian normal coordinates) in a neighborhood U of P in which the
components of the Levi-Civita connection vanish at P. The geodesic is, in
this system, a straight line: y
a
= c
a
s. This means that, as long as traverses
U, the observer will not feel gravitation: the geodesic equation reduces to the
forceless equation
du
a
ds
=
d
2
y
a
ds
2
= 0. This is an inertial observer in the absence
of external forces. If = 0, covariant derivatives reduce to usual derivatives.
If external forces are present, they will have the same expressions they have
in Special Relativity. Thus, the inertial observer will see the force equation
du
a
ds
= F
a
of Special Relativity (see Section 3.7).
3.12 There is actually more. Given any curve , it is possible to nd a
local frame in which the components of the Levi-Civita connection vanish
along . That observer would not feel the presence of gravitation.
3.13 How pointlike is a real observer ? We are used to say that an
observer can always know whether he/she is accelerated or not, by making
experiments with accelerometers and gyroscopes. The point is that all such
apparatuses are extended objects. We shall see later that a gravitational
eld is actually represented by curvature and that two geodesics are enough
to denounce its presence ( 3.46).
3.14 As repeatedly said, the principle of equivalence is a heuristic guiding
precept. It states that, as long as the dimensions involved in the denition of
an observer are negligible, an observer can choose his/hers coordinates so that
everything (s)he experiences is described by the laws of Special Relativity.
3.2 The Minimal Coupling Prescription
3.15 The equivalence principle has been used up to now in a one-way trip,
from General Relativity to Special Relativity. Some frame exists in which the
connection vanishes, so that covariant derivatives reduce to simple deriva-
tives. This can be used in the opposite sense. Given a special-relativistic
expression, how to get its version in the presence of a gravitational eld ?
71
The answer is now very simple: replace common derivatives on any tensorial
object by covariant derivatives. Symbolically, this is represented by a rule,

. (3.19)
This comma semi-colon rule is the minimal coupling prescription. To this rules
must be added the already discussed passage from at to curved metric,

ab
g

. (3.20)
Once acccepted, these two rules allow to translate special-relativistic laws,
or equations, into expressions which hold in the presence of a gravitational
eld.
3.16 Conservation of energy is one of the most important laws of
Physics. Its special-relativistic version states that the energy-momentum
tensor has vanishing divergence:

= T

,
= 0. (3.21)
In the presence of a gravitational eld, this becomes

= 0, (3.22)
being understood that every other derivative and any metric factor appear-
ing in T

are equally replaced according to (3.19) and (3.20). Notice that,

because of the metricity condition (2.39), the metric can be inserted in or
extracted from covariant derivatives at will. Equation (3.22) is usually rep-
resented in the comma-notation, as
T

;
= 0 . (3.23)
3.17 An exercise: A dust cloud (or incoherent uid) is a uid formed by
massive particles ignoring each other. It is a gas without pressure only
dust
cloud
energy is present. Special Relativity gives its energymomentum the form
T

= U

, (3.24)
where is the energy density. The 4vector U is a eld representing the
velocities of the uid streamlines.
72
Comment 3.3 The following results are immediate:
T

= U

; T

= ; g

= .
The covariant divergence of T

must vanish:
D

= U

( U

) + U

= U

( U

) +
D
Ds
U

= 0.
(3.25)
Contract this expression with U

: as U

= 1,
D

(U

) + U

D
Ds
U

= 0.
As the metric can be inserted into or extracted from the covatiant derivative
without any modication (because it is preserved by the Levi-Civita con-
nection), it follows that U

D
Ds
U

= 0. We get in consequence a continuity

equation, or energy ux conservation:
D

(U

) = 0.
Taken into (3.25), this leads to the geodesic equation for the streamlines:
D
Ds
U

= 0 .
3.18 The covariant derivative D

of a scalar eld is just the usual

derivative, D

. But the derivative

is by itself a vector. The

Laplacian operator, or (its usual name in 4 dimensions) DAlembertian, is
D

. The contracted form

of the Levi-Civita
connection has a special expression in terms of the metric. From Eq. (2.40),

=
1
2
g

} =
1
2
g

=
1
2
tr[g
1

g],
where the matrix g = (g

) and its inverse have been introduced. This is

=
1
2
tr[

ln g] =
1
2

tr[ln g].
For any matrix M, tr[ln M] = ln det M, so that

=
1
2

ln[det g] =

ln[g]
1/2
=
1

g, (3.26)
73
where g = | det g|. The DAlembertian becomes
=

+
1

g
[

g]

,
or
Laplace
Beltrami
operator
=
1

]. (3.27)
Laplaceans on curved spaces are known since long, and are called Laplace-
Beltrami operators.
3.19 Another example: the electromagnetic eld. To begin with, we
notice that, due to (3.18) and using (3.26),

;
=
1

g A

] . (3.28)
This is, by the way, the general form of the covariant divergence of a four-
covariant
divergence
vector. The eld strength F

is antisymmetric and, due to the symmetry of

the Christofell symbols, remain the same:
F

+ [

]A

.
In consequence, only the metric changes in the Lagrangean L
em
=
1
4
F
ab
F
ab
=
1
4

ac

bd
F
ab
F
cd
of Special Relativity: it becomes
L
em
=
1
4
g

. (3.29)
Maxwells equations, which in Special Relativity are written

= F
,
+ F
,
+ F
,
= 0 (3.30)

= F

,
= j

, (3.31)
take the form
F
;
+ F
;
+ F
;
= 0 (3.32)
F

;
= j

. (3.33)
74
The latter deserves to be seen in some detail. Let us rst examine the
derivative
F

;
=

.
The last term vanishes, again because

is symmetric. If another con-

nection were at work, a coupling with torsion would turn up: the last term
would be (
1
2
T

). The rst two terms give Maxwells equations in

the form
F

;
=
1

g F

] = j

. (3.34)
The energy-momentum tensor of the electromagnetic eld, which is
T
ab
=
cd
F
ac
F
bd
+
1
4

ab

ec

fd
F
ef
F
cd
in Special Relativity, becomes
T

= g

+
1
4
g

. (3.35)
Comment 3.4 Notice
g

= 0.
3.20 A uid of pressure p and energy density has an energy-momentum
tensor generalizing (3.24):
T
ab
= (p + ) U
a
U
b
p
ab
.
This changes in a subtle way. Its form is quite analogous,
energy
momentum
T

= (p + ) U

pg

, (3.36)
but the 4-velocities are U

=
dx

ds
, with ds the modied interval.
Comment 3.5 With respect to the dust cloud, only the trace changes:
T

= U

; T

= ; T = g

= 3 p .
If the gas is ultrarelativistic, the equation of state is = 3 p and, consequently, T = 0, as
for the electromagnetic eld.
75
3.21 We have seen in 3.17 that the streamlines in a dust cloud are
geodesics. The energy-momentum density (3.36) diers from that of dust by
the presence of pressure. We can repeat the procedure of that paragraph in
order to examine its eect. Here, T

;
= 0 leads to
T

;
= U

D
Ds
(p + ) + (p + )
DU

Ds
+ (p + )U

p = 0.
Contraction with U

leads now to
(U

)
;
= p U

;
,
so that the energy ux is no more conserved. Taking this expression back
into that of T

;
= 0 implies
(p + )
DU

Ds
=

p U

Dp
Ds
= (g

p . (3.37)
More will be said on this force equation in 3.50.
3.22 Suppose that T

is symmetric, as in the examples above. Dene the

quantity
L

= x

.
It is immediate that
L

;
= 0
and
U

= (x

).
3.3 Einsteins Field Equations
3.23 The Einstein tensor (2.59) is a purely geometrical second-order tensor
which has vanishing covariant derivative. It is actually possible to prove that
it is the only one. The energy-momentum tensor is a physical object with the
same property. The next stroke of genius comes here. Einstein was convinced
that some physical chacteristic of the sources of a gravitational eld should
engender the deformation in spacetime, that is, in its geometry. He looked
for a dynamical equation which gave, in the non-relativistic, classical limit,
76
the newtonian theory. This means that he had to generalize the Poisson
equation
V = 4G (3.38)
within riemannian geometry. The G

has second derivatives of the metric,

and the energy-momentum tensor contains, as one of its components, the
energy density.
He took then the bold step of equating them to each other, obtaining
what we know nowadays to be the simplest possible generalization of the
Poisson equation in a riemannian context:
eld
equation
R

1
2
g

R =
8G
c
4
T

. (3.39)
This is the Einstein equation, which xes the dynamics of a gravitational
eld. The constant in the right-hand side was at rst unknown, but he xed
it when he obtained, in the due limit, the Poisson equation of the newtonian
theory (as will be seen in 3.31 below).
Comment 3.6 A text on gravitation has always a place for the value of G. At present
time, the best experimental value for the gravitational constant is
G = 6.67390 10
11
m
3
/kg/sec
2
.
The uncertainty is 0.0014%. From that Earths and Suns masses can be obtained. The
values are M

= 5.97223(0.00008)10
24
kg and M

= 1.98843(0.00003)10
30
kg. The
apparatus used (by Jens H. Gundlach et al, University of Washington, 2000) is a moderm
version of the Cavendish torsion balance.
3.24 Contracting (3.39) with g

, we nd
R =
8G
c
4
T , (3.40)
where T = g

. This result can be inserted back into the Einstein equa-

tion, to give it the form
R

=
8G
c
4
_
T

1
2
g

. (3.41)
3.25 Consider the sourceless case, in which T

= 0. It follows from
the above equation that R

= 0 and, therefore, that R = 0. Notice that

77
this does not imply R

= 0. The Riemann tensor can be nonvanishing

even in the absence of source. Einsteins equations are non-linear and, in
consequence, the gravitational eld can engender itself. Absence of gravita-
tion is signalled by R

= 0, which means a at spacetime. This case

Minkowski spacetime is a particular solution of the sourceless equations.
Beautiful examples of solutions without any source at all are the de Sitter
spaces (see subsection 4.3.2 below).
3.26 It is usual to introduce test particles, and test elds, to probe a
gravitational eld, as we have repeatedly done when discussing geodesics.
They are meant to be just that, test objects. They are supposed to have no
inuence on the gravitational eld, which they see as a background. They
do not contribute to the source.
3.27 In reality, the Einstein tensor (2.59) is not the most general parallel-
transported purely geometrical second-order tensor which has vanishing co-
variant derivative. The metric has the same property. Consequently, it is in
principle possible to add a term g

to G

, with a constant. Equation

(3.39) becomes
R

(
1
2
R + )g

=
8G
c
4
T

. (3.42)
From the point of view of covariantly preserved objects, this equation is as
valid as (3.39). In his rst trial to apply his theory to cosmology, Einstein
looked for a static solution. He found it, but it was unstable. He then
added the term g

to make it stable, and gave to the name cosmological

constant. Later, when evidence for an expanding universe became clear, he
called this the biggest blunder in his life, and dropped the term. This
is the same as putting = 0. It was not a blunder: recent cosmological
evidence claims for = 0. Equation (3.42) is the Einsteins equation with a
cosmological term. With this extra term, Eq,(3.41) becomes
R

=
8G
c
4
_
T

1
2
g

. (3.43)
78
3.4 Action of the Gravitational Field
3.28 Einsteins equations can be derived from an action functional, the
Hilbert-Einstein action
Hilbert
action
S[g] =
_

g R d
4
x . (3.44)
It is convenient to separate the metric as soon as possible, as in R = g

.
Variations in the integration measure, which is metric-dependent, are con-
centrated in the Jacobian

g. Taking variations with respect to the metric,

S[g]
g

=
_

g
g

d
4
x+
_

g
g

d
4
x+
_

g g

d
4
x.
The rst term is

g
g

=
1
2

g
(g)
g

=
1
2

g
(exp[tr ln g])
g

=
1
2

g
exp[tr ln g]
tr ln g
g

=
1
2

g
(g) tr
ln g
g

=
1
2

g
(g) tr
_
g
1
g
g

_
=
1
2

g
(g)
_
g

_
=
1
2

g
(g)
_
g

_
=
1
2

g g

.
The last term can be shown to produce a total divergence and can be dropped.
The rst two terms give
S[g]
g

g
_
R

1
2
g

.
This is the left-hand side of (3.39). In the absence of sources, Einsteins
equation reduces to

g
_
R

1
2
g

= 0. (3.45)
3.29 A few comments:
1. variations have been taken with respect to the metric, which is the
fundamental eld
79
2. it is perhaps strange that the variation of R

, which encapsulates the

gravitational eld-strength, gives no contribution
3. If a source is present, with Lagragian density L and action
S
source
=
_

g L d
4
x, (3.46)
its contribution to the eld equation, that is, its modied energy-
momentum tensor density, will be given by

g T

=

g

S
source
=
_

g

_
g L

d
4
x
=
_

g
L
g

d
4
x +
_
L

g
g

d
4
x =

g
_
L
g

1
2
g

L
_
.
Consequently,
T

=
L
g

+
1
2
g

L . (3.47)
4. Einsteins equation (3.42) with a cosmological term comes, in and anal-
ogous way, from the action
S[g] =
_

g (R + 2) d
4
x. (3.48)
3.30 We have used in 1.15 a dimensional argument to write the action of dimensions
a test particle and of its coupling with an electromagnetic potential. Nave
dimensional analysis can be very useful in eld theories. Not only do they
help keeping trace of factors in calculations, but also lead to deep questions in
the problem of quantization. Actually, nave dimensional considerations are
enough to exhibit a fundamental dierence between gravitation and all the
other basic interactions. Indeed, all the coupling constants are dimensionless
quantities (think of the electric charge e), except that of gravitation. This
is related to the Lagrangean density

gR, whose dimension diers from
those of the other theories. Let us be nave for a while, and consider usual
mechanical dimensions in terms of mass (M), length (L) and time (T). The
dimension of a velocity will be represented as [v] = [c] = [LT
1
]; that of a
force as [F] = [MLT
2
]; an energy will have [E] = [FL] = [ML
2
T
2
]; a
80
pressure, [p] = [FL
2
] = [EL
3
]; an action, [S] = [] = [ET] = ML
2
T
1
.
Metric is dimensionless, so that [ds] = [dx] = [L]. The fact that a quantity
is dimensionless will be represented as in [e] = [g] = [0].
Field theory makes use of natural units (seemingly forbidden by interna-
tional law), in which = 1 and c = 1. In that scheme, which greatly simplify
discussions, actions and velocities are dimensionless. In consequence, [L] =
[T] = [M
1
] and only one mechanical dimension remains. Usually, the mass
M is taken as the standard reference. A new set turns out. As examples,
[F] = [2], [E] = [M], [ds] = [M
1
]. We say then that force has numeric
dimension zero, energy has dimension one, length has dimension minus one.
quantity usual dimension natural dimension numeric
mass M M + 1
length L M
1
- 1
time T M
1
- 1
velocity LT
1
M
0
0
acceleration LT
2
M + 1
force MLT
2
M
2
+ 2
Newtons G M
1
L
3
T
2
M
2
- 2
energy ML
2
T
2
M
1
+ 1
action ML
2
T
1
M
0
0
pressure ML
1
T
2
M
2
- 2
g

M
0
L
0
T
0
M
0
0
ds L M
1
- 1
charge e M
0
L
0
T
0
M
0
0
A

ML
2
T
2
M + 1
F

, E, B MLT
2
M
2
+ 2
T

ML
1
T
2
M
2
- 2
_

gF

d
4
x M
0
L
0
T
0
M
0
0

L
1
M + 1
R

, R L
2
M
2
+ 2
_

g R d
4
x L
2
M
2
- 2
81
Fields representing elementary particles have, in general, natural dimen-
sion + 1. Notice that
_

g R d
4
x has not the dimension of an action.
Actually, in order to give coherent results, it must be multiplied by a con-
stant of dimension MT
1
, actually
c
3
16G
.
3.5 Non-Relativistic Limit
3.31 A massive particle follows the geodesic equation (3.10),
dU

ds
+

=
d
2
x

ds
2
+

dx

ds
dx

ds
= 0, (3.49)
which comes from the rst term in action (1.7),
S = mc
_
ds . (3.50)
To compare with the non-relativistic case, we recall that the motion of a
particle in a gravitational eld is in that case described by the Lagrangian
L = mc
2
+
mv
2
2
mV .
This means the action
S =
_
mc
_
c
v
2
2c
+
V
c
_
dt . (3.51)
Comparison with (3.50) shows that necessarily
ds =
_
c
v
2
2c
+
V
c
_
dt
which, neglecting all the smaller terms, gives
ds
2
=
_
1 +
2V
c
2
_
c
2
dt
2
dx
2
. (3.52)
We see that, in the non-relativistic limit,
non
retativistic
metric
g
11
= g
22
= g
33
= 1 ; g
00
= 1 +
2V
c
2
. (3.53)
82
According to Special Relativity, the energy-momentum density tensor (3.36)
reduces to the sole component
T
00
= c
2
, (3.54)
where is the mass density. The trace has the same value, T = c
2
. If we
use Eq.(3.41), we have
R

=
8G
c
2

_

1
2
g

.
In particular,
R
00
=
4G
c
2
. (3.55)
The other cases are R
i=j
= 0 and R
ii
= 4GV/c
4
= in our approximation.
We can then proceed to a careful calculation of R
00
as given by (2.48), using
(3.53) and (2.40). It turns out that the terms in are of at least second
order in v/c. The derivatives with respect to x
0
= ct are of lower order if
compared with the derivatives with respect to the space coordinates x
k
, due
to the presence of the factors 1/c. What remains is R
00
=

k
00
x
k
. It is also
found that
k
00

1
2
g
kj g
00
x
j
. But this is =
1
c
2
V
x
k
. Thus,
R
00
=
1
c
2

x
k
V
x
k
=
1
c
2
V . (3.56)
Comparison with (3.55) leads then to
Poisson
recovered
V = 4G . (3.57)
This shows how Einsteins equation (3.39) reduces to the Poisson equation
(3.38) in the non-relativistic limit. By the way, the above result conrms the
value of the constant introduced in (3.39).
3.32 Notice that the non-relativistic limit corresponds to a weak gravita-
tional eld. A strong eld would accelerate the particles so that soon the
small-velocity approximation would fail.
3.33 If the cosmological constant is nonvanishing, then we should use
(3.43) to obtain R

, instead of Eq.(3.41) as above. An extra term g

00
83
appears in R
00
. The Poisson equation becomes V = 4G c
2
. To
examine the eect of this term in the non-relativistic limit, we can separate
the potential in two pieces, V = V
1
+ V
2
, with V
1
= 4G . Then, V
2
=
c
2
has for solution V
2
= c
2
r
2
/6, leading to a force F

(r) = c
2
r/3.
As > 0 by present-day evidence, this harmonic oscillator-type force is
repulsive. This is a universal eect: the cosmological constant term produces
repulsion between any two bodies.
3.34 Let us now examine what happens to the geodesic equation. First of
all, the four-velocity U has the components (, v/c) in Special Relativity,
with = 1/
_
1 v
2
/c
2
. In the non-relativistic limit,
1 +
1
2
v
2
c
2
;
U (1 +
1
2
v
2
c
2
, v/c) .
Concerning the interval ds, all 3-space distances |dx| are negligible if com-
pared with cdt, so that ds
2
c
2
dt
2
. Consequently,
dU
ds

_
d
d(ct)
_
1
2
v
2
c
2
_
,
1
c
2
d
dt
v
_
.
Now, in order to use the geodesic equation, we have to calculate the compo-
nents of the connection given in Eq.(2.40),

=
1
2
g

} .
To begin with, we notice that g
11
= g
22
= g
33
= 1 and, to the order we
are considering, g
00
= 12V/c. Given the metric (3.53), only the derivatives
with respect to the space variables of g
00
are = 0. We arrive at the general
expression

=
1
c
2
_
(
k

+
k

)
0

k
V .
The only non-vanishing components are

k
00
=

0
0k
=

0
k0
=
1
c
2

k
V .
The geodesic equations are:
84
1. for the time-like component,
0 =
dU
0
ds
+

d
d(ct)
_
1
2
v
2
c
2
_
+ 2

0
k0
U
k
U
0

d
dt
_
1
2
v
2
c
3
_
+ 2
1
c
3

k
V v
k
.
Both terms are of order 1/c
3
, equally negligible in our approximation.
2. the space-like components are more informative:
0 =
dU
k
ds
+

1
c
2
d
dt
v
k
+

k
00
U
0
U
0
=
1
c
2
d
dt
v
k
+
1
c
2

k
V ,
which is the force equation
force
equation
d
dt
v
k
=
k
V . (3.58)
3.35 In the non-relativistic limit, only g
00
remains, as well as R
00
. Ein-
steins equation greatly enlarge the the scope of the problem. What they
achieve is better stateid in the words of Mashhoon (reference of the footnote
in page 3):
the newtonian potential V is generalized to the ten compo-
nents of the metric; the acceleration of gravity is replaced by
the Christoel connection; and the tidal matrix

2
V
x
i
x
j
is replaced
by the curvature tensor; the trace of the tidal matrix is related
to the local density of matter ; the Riemann tensor is related to
the energymomentum tensor
3.6 About Time, and Space
3.6.1 Time Recovered
3.36 A gravitational eld is said to be constant when a reference frame
exists in which all the components g

are independent of the time coordi-

nate x
0
. This coordinate, by the way, is usually referred to as coordinate
time, or world time. The non-relativistic limit given above is an example
of constant gravitational eld, as the potential V in (3.53) is supposed to be
time-independent.
85
3.37 Let us now examine the relation between proper time and the coor-
dinate time x
0
= ct. Consider two close events at the same point in space.
As dx = 0, the interval between them will be ds = cd, where is time as
seen at the point. This means that ds
2
= c
2
d
2
= g
00
(dx
0
)
2
, or that
d =
1
c

g
00
dx
0
=

g
00
dt . (3.59)
As long as the coordinates are dened, the time lapse between two events at
the same point in space will be given by
=
1
c
_

g
00
dx
0
. (3.60)
This is the proper time at the point. In a constant gravitational eld,
=
1
c

g
00
x
0
. (3.61)
Once a coordinate system is established, it is the coordinate time which is
seen from abroad. Nevertheless, a test particle plunged in the eld will see
its proper time.
3.38 Consider some periodic phenomenon taking place in a constant grav-
itational eld. Its period will be given by the formula above and, as such,
will be dierent when measured in coordinate time or with a proper clock.
Its proper frequency will have, for the same reason, the values
=

0

g
00
. (3.62)
In the non-relativistic limit, (3.53) will lead to
=
0
_
1
V
c
2
_
. (3.63)
3.39 Consider a light ray goint from a point 1 to a point 2 in a weak
constant gravitational eld. If the values of the potential are V
1
and V
2
at
points 1 and 2, its proper frequency will be
1
=
0
_
1
V
1
c
2
_
at point
1 and
2
=
0
_
1
V
2
c
2
_
at point 2. It will consequently change while
86
moving from one point to the other. As the potential is a negative function,
V = |V |, this change will be given by
=
2

1
=
0
_
|V
2
| |V
1
|
c
2
_
. (3.64)
If the eld is stronger in 1 than in 2, |V
1
| > |V
2
| and
2
<
1
. The
red
shift
frequency, if in the visible region, will become redder: this is the phenomenon
of gravitational red-shift.
3.40 This eect provides one of the three classical tests of General Rel-
ativity. The others will be seen later: they are the precession of the planets
perihelia ( 4.12) and the light ray deviation ( 4.13).
3.6.2 Space
3.41 In Special Relativity, it is enough to put dx
0
= 0 in the interval ds to
get the innitesimal space distance dl. The relationship between the proper
time and the coordinate time is the same at every point. This is no more
the case in General Relativity. The standard procedure to obtain the space
interval runs as follows (see Figure 3.2). Consider two close points in space,
P = (x

) and Q = (x

+dx

). Suppose Q sends a light signal to P, at which

there is a mirror which sends it back through the same space path. With
time measured at Q, the distance between Q and P will be dl = c/2. The
interval, which vanishes for a light signal, is given by
ds
2
= g
ij
dx
i
dx
j
+ 2 g
0j
dx
0
dx
j
+ g
00
dx
0
dx
0
= 0 . (3.65)
Solving this second-order polynomial for dx
0
, we nd the two values
dx
0
(1)
=
1
g
00
_
g
0j
dx
j

_
(g
0i
g
0j
g
ij
g
00
)dx
i
dx
j
_
;
dx
0
(2)
=
1
g
00
_
g
0j
dx
j
+
_
(g
0i
g
0j
g
ij
g
00
)dx
i
dx
j
_
.
The interval in coordinate time from the emission to the reception of the
signal at Q will be
dx
0
(2)
dx
0
(1)
=
2
g
00
_
(g
0i
g
0j
g
ij
g
00
)dx
i
dx
j
.
87
The corresponding proper time is obtained by (3.59):
d =
2
c

g
00
_
(g
0i
g
0j
g
ij
g
00
)dx
i
dx
j
.
The space interval dl = c/2 is then
dl =

_
g
0i
g
0j
g
00
g
ij
_
dx
i
dx
j
,
or
dl
2
=
ij
dx
i
dx
j
, with
ij
= g
ij
+
g
0i
g
0j
g
00
. (3.66)
This is the space metric. Curiously enough, the inverse is simpler: it so
happens that

ij
= g
ij
. (3.67)
This can be seen by contracting g
ij
with
jk
.
Q P
x
0
x
0
dx
2
0
x
0
dx
1
0
x
0
x
0
+
Figure 3.2: Q sends a light sign towards a mirror at P, which sends it back.
3.42 Finite distances in space have no meaning in the general case, in
which the metric is time-dependent. If we integrate
_
dl and take the inmum
(as explained in 2.45), the result will depend on the world-lines. Only
constant gravitational elds allow nite space distances to be dened.
88
3.43 When are two events simultaneous ? In the case of Figure 3.2, the
instant x
0
of P should be simultaneous to the instant of Q which is just in
the middle between the emission and the reception of the signal, which is
x
0
+ x
0
= x
0
+
1
2
[dx
0
(2)
+ dx
0
(1)
] .
Thus, the dierence between the coordinate times of two simultaneous but
spatially distinct events is
x
0
=
g
0i
g
00
dx
i
.
Time ows dierently in dierent points of space. This relation allows one to
synchonize the clocks in a small region. Actually, that can be progressively
done along any open line. But not along a closed line: if we go along a closed
curve, we arrive back at the starting point with a non-vanishing x
0
. This
is not a property of the eld, but of the general frame. For any given eld, it
is possible to choose a system of coordinates in which the three components
g
0i
vanish. In that system, it is possible to synchonize the clocks.
3.44 Suppose that a reference system exists in which
synchronous
system
1. the synchronization condition holds:
g
0i
= 0 (3.68)
2. and further
g
00
= 1. (3.69)
Then, the coordinate time x
0
= ct in that system will represent proper time
in all points covered by the system. That system is called a synchronous
system. The interval will, in such a system, have the form
ds
2
= c
2
dt
2

ij
dx
i
dx
j
(3.70)
and

ij
= g
ij
. (3.71)
89
Direct calculation shows that

00
= 0 in such a coordinate system. This
has an important consequence: the four-velocity u

=
dx

ds
, tangent to the
world-line x
i
= 0 (i = 1, 2, 3)) has components (u
0
= 1, u
i
= 0) and satises
the geodesic equation. Such geodesics are orthogonal to the surfaces ct =
constant. In this way it is possible to conceive spacetime as a direct product
of space (the above surfaces) and time. This system is not unique: any
transformation preserving the time coordinate, or simple changing its origin,
will give another synchronous system.
3.7 Equivalence, Once Again
3.45 Consider a Levi-Civita connection

dened on a manifold M and

a point P M. Without loss of generality, we can choose around P a
coordinate system {x

} such that x

(P) = 0. Such a system will cover a

neighbourhood N of P (its coordinate neighbourhood), and will provide dual
holonomic bases {

x

, dx

}, with dx

(

x

) =

, for vector elds and 1-forms

(covector elds) on N. Any other base {e
a
, e
a
such that e
a
(e
b
) =
a
b
}, will
be given by the components of its members in terms of such initial holonomic
bases, as {e
a
= e
a

(x)

x

, e
a
= e
a

(x)dx

}. Take another coordinate system

{y
a
} around P, with a coordinate neighborhood N

such that P N
N

= . This second coordinate system will dene another holonomic base

{e
a
=

y
a
=
x

y
a

, e
a
= dy
a
=
y
a
x

dx

}. It will be enough for our purposes

to consider inside the intersection NN

a non-empty subdomain N

, small
enough to ensure that only terms up to rst order in the x

s can be retained
in the calculations.
Let us indicate by

(x) the components of in the rst holonomic

base, and by

a
b
(x) the components of referred to base {

y
a
, dy
a
}. These
components will be related by

(x) =
x

y
a

a
b
(x)
y
b
x

+
x

y
c

y
c
x

, (3.72)
or

a
b
(x) =
y
a
x

(x)
x

y
b
+
y
a
x

y
b
. (3.73)
90
Let us indicate by

(P) the value of the connection components

at the point P in the rst holonomic system. On a small enough domain N

the connection components will be approximated by

(x) =

+ x

]
P
to rst order in the coordinates x

. Choose now the second system of coor-

dinates
y
a
=
a

+
1
2

a

. (3.74)
Then,
y
a
x

=
a

+
a

;
x

y
a
=

c
y
c
.
Taken into (3.73), these expressions lead to

a
b
(x) =
a

b
_

(P)

_
x

. (3.75)
We see that at the point P, which means {x

= 0}, the connection compo-

nents in base {
a
} vanish:

a
b
(P) = 0. The curvature tensor at P, however,
is not zero:

(P) =

(P)

(P) +

.
It is thus possible to make the connection to vanish at each point by a suit-
able choice of coordinate system. The equation of force, or the geodesic
equation, acquire the expressions they have in Special Relativity. The cur-
vature, nevertheless, as a real tensor, cannot be made to vanish by a choice
of coordinates.
Furthermore, we have from Eq.(3.74)
d
ds
y
a
=
a

_
U

()
U

;
d
2
ds
2
y
a
=
a

_
d
ds
U

()
U

_
+
a

()
x

_
d
ds
U

_
.
Then, at P,
dy
a
ds
=
a

(P) ;
91
d
2
y
a
ds
2
=
a

.
Suppose a selfparallel curve goes through P in N

with velocity U

. Then,
at P,
dy
a
ds
=
a

(P) = constant = U
a
(P) and
d
2
ds
2
y
a
= 0. As by Eq.(3.74)
geodesic
system
y
a
(P) = 0, this gives for the geodesic
y
a
= U
a
(P) s
in some neighborhood around P. The geodesic is, in this system, a straight
line.
3.8 More About Curves
3.8.1 Geodesic Deviation
3.46 Curvature can be revealed by the study of two nearby geodesics. Let
us take again Eq. (3.1), rewritten in the form
U

;
= V

;
(3.76)
The deviation between two neighboring geodesics in Figure 3.1 is mea-
sured by the vector parameter

= V

V , with V a constant. It gives

the dierence between two points (such as a
1
and b
1
in that Figure) with
the same value of the parameter u. Or, if we prefer, relates two geodesics
corresponding to V and V + V .
Let us examine now the second-order derivative
D
2
V

Du
2
= (U

;
)
;
U

.
Use of Eq.(3.76) allows to write
D
2
V

Du
2
= (V

;
)
;
U

= V

;
U

;
+ V

;;
U

= U

;
V

;
+ V

;;
U

,
where once again use has been made of Eq.(3.76). Now, Eq.(2.55) leads to
D
2
V

Du
2
= (U

;
V

;
+ V

;;
) V

92
= V

(U

;
U

)
;
+

.
The rst term on the right-hand side vanishes by the geodesic equation. We
have thus, for the parameter ,
D
2

Du
2
=

. (3.77)
This the geodesic deviation equation.
As a heuristic, qualitative guide: test particles tend to close to each other
in regions of positive curvature and to part from each other in regions of
negative curvature.
3.8.2 General Observers
3.47 Let us go back to the notion of observer introduced in 3.11. An
ideal observer a timelike curve will only feel the connection, and that
can be made to vanish along a piece of that curve. Nevertheless, a real
observer will have at least two points, each one following its own timelike
curve. It will consequently feel curvature that is, the gravitational eld.
3.48 Let us go back to curves, with a mind to observers. Given a curve ,
it is convenient to attach a vector basis at each one of its points. The best
bases are those which are (pseudo-)orthogonal. This is to say that, if the
members have components h
a

, then
g

h
a

h
b

=
ab
. (3.78)
Consider a set of 4 vectors e
0
, e
1
, e
2
, e
3
at a point P on , satisfying the
following conditions:
De
0
Ds
= be
3
;
De
3
Ds
= ce
1
+ be
0
;
De
1
Ds
= de
2
ce
3
;
De
2
Ds
= d e
1
. (3.79)
These are called the FrenetSerret conditions. The choice of non-vanishing
parameters is such as to allow the 4vectors to be orthogonal. Actually, these
four vectors are furthermore required to be orthonormal at each point of :
e
0
e
0
= 1; e
1
e
1
= e
2
e
2
= e
3
e
3
= 1 . (3.80)
93
We shall always consider timelike curves, with velocity U = e
0
. In the Frenet
Serret language, e
1
, e
2
, e
3
will be the rst, second and third normals to at
P. The parameters b, c, d are real numbers, called the rst, second and third
curvatures of at P.
By what we have said above, be
3
= A, the acceleration, which is orthog-
onal to the velocity. The absolute value of the rst curvature of at P is,
thus, the acceleration modulus: |b| = |A|.
For a geodesic, b = c = d = 0. The case b = constant, c = d = 0
corresponds to an hyperbola of constant curvature and the case b = constant,
c =constant, d = 0, to a helix.
3.49 Of course, most curves are not geodesics. A geodesic represents an
observer in the absence of any external force. An observer may be settled on
a linearly accelerated rocket, or turning around Earth, or still going through
a mad spiral trajectory. An orthogonal tetrad dened as above remains an
orthogonal tetrad under parallel transport, which also preserves the compo-
nents of a vector in that tetrad. There is, however, a problem with parallel
transport: if we take, as above, e
0
as the velocity U at a point P, e
0
will not
be the velocity at other points of , unless is a geodesic.
There are other kinds of transport which preserve orthogonal tetrads
and components. There is one, in special, which corrects the mentioned
problem with the velocity:
The Fermi-Walker derivative of a vector V is dened by
Fermi
Walker
transport
D
FW
V

Ds
=
DV

Ds
b V

(e

0
e

3
e

0
e

3
) =
DV

Ds
V

(U

). (3.81)
A vector V is said to be FermiWalkertransported along a curve of
velocity eld U if its the Fermi-Walker derivative vanishes,
D
FW
V

Ds
=
DV

Ds
V

(U

) = 0 . (3.82)
We see that, applied to U,
D
FW
U

Ds
= 0, and also that
D
FW
Ds
=
D
Ds
if the
curve is a geodesic. Take two vector elds X and Y such that
D
FW
X

Ds
= 0
and
D
FW
Y

Ds
= 0. Then, it follows that the component of X along Y is
preserved along :
D(X

)
Ds
=
d(X

)
ds
= 0.
94
In particular, if e
0
= U at P, then e
0
will remain = U along the curve if it is
Fermi-Walker transported.
3.8.3 Transversality
3.50 Given a curve whose tangent velocity is U, it is interesting to
introduce a transversal metric by
h

= g

. (3.83)
Transversality is evident: h

= 0. In particular, h

= A

. A projector
is an operator P satisfying P
2
= P. Matrices (h

) = (g

) = (

) are projectors: h

= h

. They satisfy h

= h

=
h

(g

) = h

. Notice that g

= h

= h

= 3. The
energymomentum tensor (3.36) of a general uid can be rewritten as
T

= (p + )U

p g

= U

p h

. (3.84)
The transversal metric extracts the pressure:
T

= 3 p.
The Einstein equations with this energymomentum tensor give, by contrac-
tion,
R

=
4G
c
4
( + 3p) . (3.85)
We have seen in 3.21 that the streamlines of a general uid unlike those of
a dust cloud are not geodesics. The equation of force (3.37) has, actually,
the form
uid
stream
lines
(p + )
DU

Ds
= h

p . (3.86)
The pressure gradient, as usual, engenders a force. In the relativistic case the
force is always transversal to the curve. Here, it is the transversal gradient
that turns up. The equation above governs the streamlines of a general uid.
95
3.8.4 Fundamental Observers
3.51 On a pseudoRiemannian spacetime, there exists always a family of
worldlines which is preferred. They represent the motion of certain preferred
observers, the fundamental observers and the curves themselves are called the
fundamental worldlines. Proper time coincides with the line parameter, so
that he 4-velocities along these lines are U

=
dx

ds
and, consequently U
2
= U

= 1. The time derivative of a tensor T

...
...
is
d
ds
T
...
...
=
d
dx

T
...
...
dx

ds
= (T
...
...
)
,
U

. This is not covariant. The covariant time

derivative, absolute derivative, is
D
Ds
T
...
...
= (T
...
...
)
;
u

.
For example, the acceleration is
A

=
D
Ds
U

= U

;
U

.
Using the Christoel connection, it is easily seen that U

= 0. This prop-
erty is analogous to that found in Minkowski space, but here only has an
invariant sense if acceleration is covariantly dened, as above.
At each point P, under a condition given below, a fundamental observer
has a 3-dimensional space which it can consider to be hers/his own: its
restspace. Such a space is tangent to the pseudoRiemannian spacetime
and, as time runs along the fundamental worldline, orthogonal to that line
at P (orthogonal to a line means orthogonal to its tangent vector, here U

).
At each point of a worldline, that 3-space is determined by the projectors
h

.
It is convenient to introduce the notations
U
(;)
=
1
2
(U
;
+ U
;
) ; U
[;]
=
1
2
(U
;
U
;
)
for the symmetric and antisymmetric parts of U
;
. There are a few important
notions to be introduced:
the vorticity tensor

= h

U
[;]
= U
[;]
+ U
[
A
]
; (3.87)
it satises

=
[]
= -

and

= 0; it is frequently indicated
by its magnitude
2
=
1
2

0.
96
the expansion tensor

= h

U
(;)
= U
(;)
U
(
A
)
; (3.88)
its transversal trace is called the volume expansion; it is = h

= U

;
, the covariant divergence of the velocity eld; it measures
the spread of nearby lines, thereby recovering the original meaning of
the word divergence; in the Friedmann model, turns up as related to
the Hubble expansion function by = 3H(t).

1
3
h

=
()
is the symmetric tracefree shear tensor;
it satises

= 0 and

= 0 and its magnitude is dened as

2
=
1
2

0. Notice

= 2
2
+
1
3

2
.
Decomposing the covariant derivative of the 4-velocity into its symmetric
and antisymmetric parts, U
;
= U
(;)
+ U
[;]
, and using the denitions
(3.87) and (3.88), we nd
U
;
=

+ A

, (3.89)
or
U
;
=

+
1
3
h

+ A

. (3.90)
3.52 With the above characterizations of the energy density and the pres-
sure, the Einstein equations reduce to the LandauRaychaudhury equation.
Let us go back to Eq.(2.55) and take its contracted version
Landau
Raychaudhury
equation
U

;;
U

;;
=

R
U

. (3.91)
Contracting now with U

,
U

;;
U

;;
=

R
U

D
Ds
U

;
(U

;
)
;
+ U

;
U

;
+

R
U

= 0 ,

d
ds
U

;
A

;
+ U

;
U

;
+

R
U

= 0 .
To obtain U

;
U

;
, we notice that

=
1
2
_
U
;
U
;
+ U
;
U
;
A

;
97

=
1
2
_
U
;
U
;
U
;
U
;
A

.
It follows that
U

;
U

;
=

= 2
2
2
2
+
1
3

2
.
The equation acquires the aspect
d
ds
+ 2 (
2

2
) +
1
3

2
A

;
+

R
U

= 0 . (3.92)
Up to this point, only denitions have been used. Einsteins equations lead,
however, to Eq.(3.85), which allows us to put the above expression into the
form
d
ds
=
1
3

2
2
_

2
_
+ A

4G
c
4
( + 3p) + . (3.93)
The promised condition comes from a detailed examination which shows
that, actually, only when

= 0 there exists a family of 3-spaces everywhere

orthogonal to U

. In that case, there is a welldened time which is the same

over each 3-space.
From the above equation, we can see the eect of each quantity on ex-
pansion: as we proceed along a fundamental worldline, expansion
decreases (an indication of attraction) with higher values of
expansion itself
shear
energy content
increases (an indication of repulsion) with higher values of
vorticity
second-acceleration
cosmological constant .
98
3.9 An Aside: Hamilton-Jacobi
3.53 The action principle we have been using [say, as in Eq.(1.5), or (3.51)]
has a teleological character which brings forth a causal problem. When we
look for the curve which minimizes the functional
S[] =
_
t
1
t
0
Ldt =
_

PQ
Ldt ,
we seem to suppose that the behavior of a particle, starting from a point P
at instant t
0
, is somehow determined by its future, which forcibly consists
in being at a xed point Q at instant t
1
. Another notion of action exists
which avoids this diculty. Instead of as a functional S is conceived, in that
version, as a function
S(q
1
(t), q
2
(t), ..., q
n
(t), t) =
_
t
t
0
dt

L[q
1
(t

), q
2
(t

), ..., q
n
(t

), q
1
(t

), q
2
(t

), ..., q
n
(t

)] (3.94)
of the nal time t and the values of the generalized coordinates at that instant
for the real trajectory. The particle, by satisfying the Lagrange equation at
each point of its path, automatically minimizes S. In eect, taking the
variation
S =
_
t
t
0
dt

_
L
q
i
q
i
+
L
q
i
q
i
_
=
_
t
t
0
dt

_
L
q
i
q
i
+

t

_
L
q
i
q
i
_

_

t

L
q
i
_
q
i
_
=
_
L
q
i
q
i
_
t
t
0
+
_
t
t
0
dt

_
L
q
i

_

t

L
q
i
__
q
i
.
The second term vanishes by Lagranges equation. In the rst term, q
i
(t
0
) =
0, so that
S =
L
q
i
q
i
= p
i
q
i
, (3.95)
which entails
p
i
=
S
q
i
. (3.96)
99
We have been forgetting the time dependence of S. The integral (3.94) says
that L =
dS
dt
. On the other hand,
dS
dt
=
S
t
+
S
q
i
q
i
=
S
t
+ p
i
q
i
.
Consequently,
S
t
= L p
i
q
i
= H. (3.97)
Thus, the total dierential of S will be
dS = p
i
dq
i
Hdt . (3.98)
Variation of the action integral
S =
_
_
p
i
dq
i
Hdt

(3.99)
leads indeed to Hamiltons equations:
S =
_ _
p
i
dq
i
+ p
i
dq
i

H
q
i
q
i
dt
H
p
i
p
i
dt
_
=
_ _
p
i
dq
i
dt
+ p
i
dq
i
dt

H
q
i
q
i

H
p
i
p
i
_
dt
=
_ _
p
i
_
dq
i
dt

H
p
i
_

_
dp
i
dt
+
H
q
i
_
q
i
_
dt.
An integration by parts was performed to arrive at the last espression which,
to produce S = 0 for arbitrary q
i
and p
i
, enforces
dq
i
dt
=
H
p
i
;
dp
i
dt
=
H
q
i
. (3.100)
3.54 Hamiltons equations are invariant under canonical transformations
leading to new variables Q
i
= Q
i
(q
i
, p
j
, t), P
i
= P
i
(q
i
, p
j
, t), H

(Q
i
, P
j
, t).
This means that if

_
_
p
i
dq
i
Hdt

= 0 ,
then also

_
_
P
i
dQ
i
H

dt

= 0
100
must hold. In consequence, the two integrals must dier by the total dier-
ential of an arbitrary function:
p
i
dq
i
Hdt = P
i
dQ
i
H

dt + dF .
F is the generating function of the canonical transformation. It is such that
dF = p
i
dq
i
P
i
dQ
i
+ (H

H)dt .
Therefore,
p
i
=
F
q
i
; P
i
=
F
Q
i
; H

= H +
F
t
. (3.101)
In the formulas above, the generating function appears as a function of the
old and the new generalized coordinates, F = F(q, Q, t). We obtain another
generating function f = f(q, P, t), with the new momenta instead of the new
generalized coordinates, by a Legendre transformation: f = F + Q
i
P
i
, for
which
df = dF + dQ
i
P
i
+ Q
i
dP
i
= p
i
dq
i
+ Q
i
dP
i
+ (H

H)dt .
In this case,
p
i
=
f
q
i
; Q
i
=
f
P
i
; H

= H +
f
t
. (3.102)
Other generating functions, related to other choices of arguments, can of
course be chosen. For instance, the function g =

i
q
i
Q
i
generates a simple
interchange of the initial coordinates and momenta.
3.55 Let us go back to Eq. (3.97),
S
t
+ H(q
1
, q
2
, ..., q
n
, p
1
, p
2
, ..., p
n
, t) = 0 . (3.103)
By Eq. (3.96), the momenta are the gradients of the action function S.
Substituting their expressions in the above formula, we nd a rst order
partial dierential equation for S,
S
t
+ H
_
q
1
, q
2
, ..., q
n
,
S
q
1
,
S
q
2
, . . . ,
S
q
n
, t
_
= 0 . (3.104)
101
This is the Hamilton-Jacobi equation. The general solution of such an equa-
tion depends on an arbitrary function. The solution which is important for
Mechanics is not the general solution, but the so-called complete solution
(from which, by the way, the general solution can be recovered). That solu-
tion contains one arbitrary constant for each independent variable, (n+1) in
the case above. Notice that only derivatives of S appear in the equation. One
of the constants (C below) turns up, consequently, isolated. The complete
solution has the form
S = f
_
q
1
, q
2
, ..., q
n
, a
1
, a
2
, . . . , a
n
, t
_
+ C. (3.105)
We have indicated the arbitrary constants by a
1
, a
2
, . . . , a
n
and C.
The connection to the mechanical problem is made as follows. Consider
a canonical transformation with generating function f, taking the original
variables (q
1
, q
2
, ..., q
n
, p
1
, p
2
, ..., p
n
) into (Q
1
, Q
2
, ..., Q
n
, a
1
, a
2
, ..., a
n
). This is
a transformation of the type summarized in Eq. (3.102), with a
1
, a
2
, ..., a
n
as
the new momenta. From those equations and (3.105),
p
i
=
S
q
i
; Q
i
=
S
P
i
; H

= H +
S
t
= 0 . (3.106)
To get the vanishing of the last expression use has been made of Eq. (3.104).
Hamilton equations havee then the solutions Q
n
= constant, a
k
= constant.
From the equations Q
i
=
S
a
i
it is possible to obtain back
q
k
= q
k
_
Q
1
, Q
2
, ..., Q
n
, a
1
, a
2
, . . . , a
n
, t
_
,
that is, the old coordinates written in terms of 2n constants and the time.
This is the solution of the equation of motion. Summing up, the procedure
runs as follows:
given the Hamiltonian, one looks for the complete solution (3.105) of
the Hamilton-Jacobi equation (3.104);
once the solution S is obtained, one derives with respect to the con-
stants a
k
and equate the results to the new constants Q
k
; the equations
Q
i
=
S
a
i
are algebraic;
102
that set of algebraic equations are then solved to give the coordinates
q
k
(t);
the momenta are then found by using p
i
=
S
q
i
.
3.56 For conservative systems, H is time-independent. The action de-
pends on time in the form
S(q, t) = S(q, 0) E t.
It follows that
H
_
q
1
, q
2
, ..., q
n
,
S(q, 0)
q
1
,
S(q, 0)
q
2
, . . . ,
S(q, 0)
q
n
_
= E , (3.107)
which is the time-independent Hamilton-Jacobi equation.
The same happens whenever some integral of motion is known from the
start. Each constant of motion is introduced as one of the constants. For
instance, central potentials, for which the angular momentum J is a constant,
will lead to a form
S(r, t) = S(r, 0) Et + J.
These are particular cases, in which time or an angle are cyclic variables.
In eect, suppose some coordinate q
(c)
is cyclic. This means that q
(c)
does
not appear explicitly in the Hamiltonian nor, consequently in the Hamilton-
Jacobi equation. The corresponding momentum is therefore constant, p
(c)
=
S
q
(c)
= a
(c)
. It is an integral of motion, and S = S(remaining variables)
+a
(c)
q
(c)
.
3.57 Of course, a suitable choice of coordinate system is essential to isolate
a cyclic variable. For a particle in a central potential, spherical coordinates
are the obvious choice. In the Hamiltonian
H =
1
2m
_
p
2
r
+
p
2

r
2
+
p
2

r
2
sin
2

_
+ U(r) ,
the variables and , besides t, are absent. We shall use the knowledge that
the angular momentum J = mr
2

is a constant, and start from the simpler
planar Hamiltonian
H =
m
2
_
r
2
+ r
2

2
_
+ U(r) =
p
2
r
2m
+
J
2
2mr
2
+ U(r) .
103
The time-independent Hamilton-Jacobi equation is then
_
S
r
r
_
2
+
J
2
r
2
= 2m(E U(r)).
Thus,
S = E t + J +
_
dr
_
2m(E U(r))
J
2
r
2
.
Now,
S
E
= C gives
t =
_
mdr
_
2m[E U(r)]
J
2
r
2
C.
And
S
J
= C

gives
= C

+
_
J dr
r
2
_
2m[E U(r)]
J
2
r
2
.
The constants can be chosen = 0, xing simply the origins of time and angle.
The rst equation,
t =
_
m dr
_
2m[E U(r)]
J
2
r
2
, (3.108)
gives implicitly r(t). The second,
=
_
J dr
r
2
_
2m[E U(r)]
J
2
r
2
, (3.109)
gives the trajectory.
We see that what is actually at work is an eective potential U
eff
(r) =
U(r) +
J
2
2mr
2
, including the angular momentum term. The values of r for
which U
eff
(r) = E represent turning points. If the function r(t) is at rst
decreasing, it becomes increasing at that value, and vice versa.
3.58 Classical planetary motion There are two general kinds of motion:
limited (bound motion) and unlimited (scattering). We shall be concerned
here only with the rst case, in which r(t) has a nite range r
min
r(t)
104
r
max
. In one turn, that is, in the time the variable takes to vary from r
min
to r
max
and back to r
min
, the angle undergoes a change
= 2
_
r
max
r
min
J
r
2
dr
_
2m[E U(r)]
J
2
r
2
. (3.110)
The trajectory will be closed if = 2m/n. A theorem (Bertrands) says
that this can happen only for two potentials, the Kepler potential U(r) =
K/r and the harmonic oscillator potential U(r) = Kr
2
. We shall here limit
ourselves to the rst case, which describes the keplerian motion of planets
around the Sun. For U(r) = K/r, Eq.(3.109) can be integrated to give
= arccos
_
_
_
_
J
r

mK
J
_
1
_
2mE +
m
2
K
2
J
2
_
_
_
. (3.111)
This trajectory corresponds to a closed ellipse. In eect, introduce the el-
lipse parameter p =
J
2
mK
and the eccentricity e =
_
1 +
2EJ
2
mK
2
. Equation
(3.111) can then be put into the form
r =
p
1 + e cos
, (3.112)
which is the equation for the ellipse. The above choice of the integration
constants corresponds to = 0 at r = r
min
, which is the orbit perihelium.
Suppose we add another potential (for example, U

= K

/r
3
, with K

small.
The orbit, as said above, will be no more closed. Staring from = 0 at
r = r
min
, the orbit will reach the value r = r
min
at = 0 in the rst
turn, and so on at each turn. The perihelium will change at each turn.
This efect (called the perihelium precession) could come, for example, from
a non-spherical form of the Sun. A turning gas sphere can be expected to be
oblate. Observations of the Sun tend to imply that its oblateness is negligible,
in any case insucient to answer for the observed precession of Mercurys
perihelium. We shall see later (4.12) that General Relativity predicts a value
in good agreement with observations.
Comment 3.7 We recall that the equation of the ellipse in cartesian coordinates is
x
2
a
2
+
y
2
b
2
= 1, e =

a
2
b
2
a
and p = b
2
/a. For a circle, a = b and e = 0.
105
3.59 We have already seen the main interest of the Hamilton-Jacobi for-
malism: in the relativistic case, the Hamilton-Jacobi equation (3.17) for a
free particle coincides, for vanishing mass, with the eikonal equation (3.18).
The formalism allows a unied treatment of test particles and light rays.
106
Chapter 4
Solutions
Einsteins equations are a nightmare for the searcher of solutions: a system
of ten coupled non-linear partial dierential equations. Its a tribute to hu-
man ingenuity that many (almost thirty to present time) solutions have been
found. The non-linear character can be interpreted as saying that the grav-
itational eld is able to engender itself. In consequence, the equations have
non-trivial solutions even in the absence of sources. Actually, most known
solutions are of this kind. We shall only examine a few examples, divided
into two categories: small scale solutions, which are of local interest,
idealized models for stars (which gives an idea of what we mean by small)
and objects alike; and large scale solutions, of cosmological interest.
4.1 Transformations
In tackling the big task of solving so dicult a problem, it is not surprising
that the hunters have always supposed a high degree of symmetry. Let us
begin with a short comment on symmetries of spacetimes.

4.1 Let us look for the condition for a vector eld to generate a symmetry
of the metric. Consider an innitesimal transformation
x

= x

(x) ,

(x) << 1 , (4.1)

arbitrary as longs as is an arbitrary, though small, function of x. To the

A treatment of the subject is given in the chapter 13 of S. Weinberg, Gravitation and

Cosmology, J. Wiley, New York, 1972.
107
rst order in ,
x

.
Under such a tranformation, the metric components will change according
to
g

(x

) = g

(x)
x

= g

(x)
_

__

_
g

(x) + g

(x)

+ g

(x)

.
On the other hand, always keeping only the rst-order terms,
g

(x

) = g

(x + ) g

(x) +

(x)

(x) g

(x) +

(x)

(x) .
Equating both expressions,
g

(x) +

(x)

(x) = g

(x) + g

(x)

+ g

(x)

,
from which we obtain the variation of the metric components at a xed point,

(x) = g

(x) g

(x) = g

(x)

+ g

(x)

(x)
g

.
(4.2)
We now calculate the covariant derivative of

, conveniently separating
the pieces comming from the Christoel symbol:

;
=

1
2

+
1
2

) .
We see then that

;
+
;
=

. (4.3)
Therefore,

(x) =
;
+
;
. (4.4)
This gives the change in the functional form of g

under the transformation

Killing
equation
generated by the eld =

. The condition for a symmetry is

g

(x) = 0,
that is,

;
+
;
= 0 . (4.5)
108
This is the Killing equation. Fields satisfying it are called Killing elds. They
generate transformations preserving the metric, which are called isometries,
or motions. Applied to the Lorentz metric, the ten generators of the Poincare
group are found.

4.2 We shall here quote three theorems:

the rst says that the maximal number of isometries in a d-dimensional
space is d(d + 1)/2.

Consequently, a given spacetime has at most 10

isometries.
the second theorem says that this maximal number can only be attained
if the scalar curvature R is a constant. There are only three kinds of
spacetimes with 10 isometries: Minkowski spacetime, for which R = 0,
and the two families of de Sitter spacetimes, one with R > 0 and the
other with R < 0.
a third theorem states that the isometries of a given metric constitute a
group (group of isometries, or group of motions). The Poincare group
is the group of motions of Minkowski space.
4.3 The converse procedure may be useful in the search of solutions: im-
pose a certain symmetry from the start, and nd the metrics satisfying
Eq.(4.3) for the case,

. (4.6)
It should be said, however, that the Killing equation is still more useful in
the study of the simmetries of a metric given a priori.
4.4 The above procedure is a very particular case of a general and pow-
erful method. How does a transformation acts on a manifold ? We are used
to translations and rotations in Euclidean space. The same transformations,
plus boosts and time translations, are at work on Minkowski space: they
preserve the Lorentz metric. For these we use generators like, for example,

See for example W.R. Davis & G.H. Katzin, Am. J. Phys. 30 (1962) 750.

L. P. Eisenhart, Riemannian Geometry, Princeton University Press, 1949.

109
L

= x

, which is a vector eld. This can be extended to general

manifolds: given a group of transformations, each generator is represented
on the manifold by a vector eld X. This vector eld presides on the in-
nitesimal transformations undergone by every tensor eld by the so called
Lie
derivative
Lie derivative, an operation denoted L
X
. The calculation above must be
repeated for each type of tensor T: take a transformation like (4.1), compare
T

(x

) obtained from the tensor behavior with T

(x

) obtained as a Taylor
series, etc. The general result is
(L
X
T)
ab...r
ef...s
= X(T
ab...r
ef...s
) (
i
X
a
)T
ib...r
ef...s
(
i
X
b
)T
ai...r
ef...s
... (
i
X
r
)T
ab...i
ef...s
+(
e
X
i
)T
ab...r
if...s
+ (
f
X
i
)T
ab...r
ei...s
+ ... + (
s
X
i
)T
ab...r
ef...i
. (4.7)
The requirement of invariance is L
X
T = 0. For T a vector eld, L
X
T is
just the commutator:
L
X
V = [X, V ] .
The vector V is invariant with respect to the transformations engendered by
X if it commutes with X.
Equation (4.2) is just the Lie derivative of g

with respect to the eld

=

:
L

(x) = g

(x) g

(x) = g

(x)

+ g

(x)

(x)
g

.
(4.8)
We shall not go into the subject in general.

Let us only state a property

which holds when the tensor T is a dierential form. For that we need a
preliminary notion.
4.5 Given a vector eld X, the interior product of a p-form by X is that
(p-1)-form i
X
, which, for any set of elds {X
1
, X
2
, . . . , X
p1
}, satises
interior
product
i
X
(X
1
, X
2
, . . . , X
p1
) = (X, X
1
, X
2
, . . . , X
p1
). (4.9)
If is a 1-form, it gives simply its action on X:
i
X
= < , X > = (X).

A very detailed account can be found in B. Schutz,Geometrical Methods of Mathe-

matical Physics, Cambridge University Press, Cambridge, 1985.
110
The interior product of X by a 2-form is that 1-form satisfying i
X
=
(X, Y ) for any eld Y . For a form of general degree, it is enough to know
that, for a basis element,
i
X
_

1

2
. . .
p

=
p

j=1
()
j1

1

2
. . . [i
X

j
] . . .
p
.
4.6 The promised result is as follows: if is a dierential form, then its
Lie derivative has a simple expression in terms of the exterior derivative and
the interior product:
L
X
= d[i
X
] + i
X
[d] .
Notice that L
X
preserves the tensor character: it takes an r-covariant, s-
contravariant tensor into another tensor of the same type.
4.2 Small Scale Solutions
Life is much simpler when a system of coordinates can be chosen so that
invariance means just independence of some of the coordinates. In that case,
Eq.(4.8) reduces to the last term and the intuitive property holds: the metric
components are independent of those variables.
4.2.1 The Schwarzschild Solution
4.7 Suppose we look for a solution of the Einstein equations which has
spherical symmetry in the space section. This would correspond to central
potentials in Classical Mechanics. It is better, in that case, to use spherical
coordinates (x
0
, x
1
, x
2
, x
3
) = (ct, r, , ). This is one of the most studied of
all solutions, and there is a standard notation for it. The interval is written
in the form
ds
2
= e

c
2
dt
2
r
2
(d
2
+ sin
2
d
2
) e

dr
2
. (4.10)
The contravariant metric is consequently
g = (g

) =
_
_
_
_
_
e

0 0 0
0 e

0 0
0 0 r
2
0
0 0 0 r
2
sin
2

_
_
_
_
_
(4.11)
111
and its covariant counterpart,
g
1
= (g

) =
_
_
_
_
_
e

0 0 0
0 e

0 0
0 0 r
2
0
0 0 0 r
2
sin
2

_
_
_
_
_
. (4.12)
We have now to build Einsteins equations. The rst step is to calculate
the components of the Levi-Civita connection, given by Eq.(2.40). Those
which are non-vanishing are:

0
00
=
1
2
d
cdt
;
0
10
=
0
01
=
1
2
d
dr
;
0
11
=
1
2
e
d
cdt
;

1
00
=
1
2
e
d
dr
;
1
01
=
1
10
=
1
2
d
cdt
;
1
11
=
1
2
d
dr
;

1
22
= r e

;
1
33
= re

sin
2
;

2
12
=
2
21
=
1
r
;
2
33
= sin cos ;

3
13
=
3
31
=
1
r
;
3
23
=
3
32
= cot . (4.13)
As the second step, we must calculate the Ricci tensor of Eq.(2.48) and
the Einstein tensor (2.59). We list those which are non-vanishing:
G
0
0
=
1
r
2
e

_
1
r
2

1
r
d
dr

; G
1
0
=
1
r
e
d
cdt
; G
1
1
=
1
r
2
e

_
1
r
d
dr
+
1
r
2
_
G
2
2
= G
3
3
=
1
2
e

_
d
2

dr
2
+
1
2
(
d
dr
)
2
+
1
r
(
d
dr

d
dr
)
1
2
(
d
dr
d
dr
)
_
+
1
2
e

_
d
2

c
2
dt
2
+
1
2
(
d
cdt
)
2

1
2
d
cdt
d
cdt
_
. (4.14)
Each G

should now be imposed to be equal to

8G
c
4
T

. The source could

be, for example, an electromagnetic eld, in which case Eq.(3.35) would be
used. Or the uid inside a star, with the source given by Eq.(3.36) sup-
plemented by an equation of state. Notice that the symmetry requirements
made above would be satised not only by a static star, but also by a radially
112
pulsating one. We shall here consider the so-called external, or vacuum
solution for this case. We shall put T

= 0, which is the case outside the

star. The four dierential equations following from G

= 0 in (4.14) reduce
in that case to only three (see Comment 4.1 below):
G
0
0
= 0 e

_
1
r
2

1
r
d
dr

=
1
r
2
;
G
1
1
= 0 e

_
1
r
2
+
1
r
d
dr

=
1
r
2
; (4.15)
G
1
0
= 0
d
dt
= 0 .
A rst result from Eqs.(4.15) is that is time-independent. Taking the
dierence between the rst two equations shows that
d
dr
=
d
dr
. (4.16)
Substituting this back in those equations lead to
e

= 1 + r
d
dr
. (4.17)
Equation (4.16) says that + is independent of r, and is consequently
a function of time alone: + = f(t). In the interval (4.10), it is always
possible to redene the time by an arbitrary transformation t = (t

), which
corresponds to adding an arbitrary function of t to . The choice of a new
time coordinate t

=
_
t
0
e
f(t)/2
dt corresponds to changing

= +f(t).
This means that it is always possible to choose the time coordinate so as to
have + = 0.
Comment 4.1 Time-independence of entails the vanishing of the last line in (4.14).
Using (4.16) in (4.17) and taking the derivative implies that also the one-but-last line
vanishes. This shows that the equation G
2
2
= G
3
3
= 0 is indeed redundant.
Integration of the only remaining equation, which is (4.17) rewritten with
= ,
e

= 1 + r
d
dr
,
leads then to
e

= e

= 1
R
S
r
,
113
where R
S
is a constant. Far from the source, when r , we have e

=
e

1, so that the metric reduces to that of Minkowski space. Large values

of r means a weak gravitational eld. To x the constant R
S
, it is enough
to impose that, at those values of r the solution reduce to the newtonian
approximation, g
00
= 1 + 2V/c
2
, with V = GM/r and M the mass of
source body. It follows that
R
S
=
2GM
c
2
. (4.18)
The interval (4.10) is therefore
ds
2
=
_
1
2GM
c
2
r
_
c
2
dt
2
r
2
(d
2
+ sin
2
d
2
)
dr
2
1
2GM
c
2
r
(4.19)
=
_
1
R
S
r
_
c
2
dt
2
r
2
(d
2
+ sin
2
d
2
)
dr
2
1
R
S
r
. (4.20)
This is the solution found by K. Schwarzschild in 1916, soon after Eintein
had presented his nal version of General Relativity. It describes the eld
caused, outside it, by a symmetrically spherical source. We see that there is
a singularity in the metric components at the value r = R
S
. The parameter
R
S
, given in Eq.(4.18), is called the Schwarzschild radius. Its value for a
body with the mass of the Sun would be R
S
3 km. For a body with
Earths mass, R
S
0.9 cm. For such objects, of course, there exists to real
Schwarzschild radius. It would be well inside their matter distribution, where
T

= 0 and the solution is not valid.

4.8 The above solution has been obtained in the absence of the cosmolog-
ical constant. Its presence would change it to

ds
2
=
_
1
2GM
c
2
r


3
r
2
_
c
2
dt
2
r
2
(d
2
+ sin
2
d
2
)
dr
2
1
2GM
c
2
r


3
r
2
.
(4.21)

See the 96 of R.C. Tolman, Relativity, Thermodynamics and Cosmology, Dover, New
York, 1987.
114
If we compare with Eq.(3.53), we nd the potential
V =
MG
r

c
2
r
2
6
. (4.22)
Eq.(3.58) would then lead to
d
dt
v =
MG
r
2
+
1
3
c
2
r . (4.23)
We recognize Newtons law in the rst term of the righ-hand side. The extra,
cosmological term has the aspect of a harmonic oscillator but, for > 0,
produces a repulsive force.
4.9 The eld, just as in the newtonian case, depends only on the mass
M. At a large distance of any limited source, the eld will forget details
on its form and tend to have a spherical symmetry. The interval above is
approximately given, at larges distances, by
ds
2
c
2
dt
2
dr
2
r
2
(d
2
+ sin
2
d
2
)
R
S
r
_
dr
2
+ c
2
dt
2
_
. (4.24)
The last term is a correction to the Lorentz metric and the above interval
should be the asymptotic limit, for large values of r, of any eld created by
any source of limited size. We see that the Schwarzschild coordinate system
used in (4.20) is asymptotically Galilean: Schwarzschilds spacetime tends
to Minkowski spacetime when r .
As g
0j
= 0 in Eq.(3.66), the 3-dimensional space sector induced by (4.20)
will have the interval
d
2
=
dr
2
1
R
S
r
+ r
2
(d
2
+ sin
2
d
2
) , (4.25)
to be compared with the Euclidean interval
d
2
= dr
2
+ r
2
(d
2
+ sin
2
d
2
) . (4.26)
At xed and , that is, radially, the distance between two points P and
Q standing outside the Schwarzschild radius will be
_
Q
P
dr
_
1
R
S
r
> r
Q
r
P
. (4.27)
115
On the other hand, the proper time will be
d =

g
00
dt =
_
1
R
S
r
dt < dt . (4.28)
We see that d = dt when r . And we see also that, at nite distances
from the source, time marches slower than time at innity. This dierence
between proper time and coordinate time arrives at an extreme case near the
Schwarzschild radius.
4.10 We can make some checking on the results found. Given the metric
_
_
_
_
_
_
1
R
S
r
0 0 0
0
1
1
R
S
r
0 0
0 0 r
2
0
0 0 0 r
2
sin
2

_
_
_
_
_
_
,
we can proceed to the laborious computation of the Christoeln and Riemann
components. We nd, for example,
R
1
212
=
R
S
r
2
(r R
S
)
; R
1
313
=
R
S
2 r
; R
1
414
=
R
S
sin
2

2 r
R
2
323
=
R
S
2 r
; R
2
424
=
(R
S
2 r) sin
2

2 r
; R
3
434
= cos
2
+
R
S
sin
2

r
.
All components of the Ricci tensor vanish, as they should for an exterior,
T

= 0 solution. The scalar curvature, of course, vanishes also. This is a

good illustration of the statement made in 3.25, by which T

= 0 implies
R

= 0 but not necessarily R

= 0. It is also a good example of a

fundamental point of General Relativity: it is the non-vanishing of R

that indicates the presence of a gravitational eld.

4.11 We can, furthermore, examine the space section. With coordinates
(r, , ), the metric is
_
_
_
_
1
1
R
S
r
0 0
0 r
2
0
0 0 r
2
sin
2

_
_
_
_
.
116
The Christoeln form the matrices
(
1
ij
) =
_
_
_
R
S
2r(R
S
r)
0 0
0 R
S
r 0
0 0 (R
S
r) sin
2

_
_
_
(
2
ij
) =
_
_
_
0
1
r
0
1
r
0 0
0 0 sin cos
_
_
_
(
3
ij
) =
_
_
_
0 0
1
r
0 0 cot
1
r
cot 0
_
_
_
.
The Ricci tensor is given by
(R
ij
) =
_
_
_
R
S
r
2
(R
S
r)
0 0
0
R
S
2r
0
0 0
R
S
sin
2

2r
_
_
_
.
Thus, the Ricci tensor of the space sector is non-trivial. The scalar curvature,
however, is zero.
4.12 Perihelium precession The Hamilton-Jacobi equation (3.17) and
the eikonal equation (3.18) provide, as we have seen, a unied approach to
the trajectories of massive particles and light rays. Let us rst examine the
motion of a particle of mass m in the above gravitational eld. As angular
momentum is conserved, it will be a plane motion, with constant . For
reasons of simplicity, we shall choose the value = /2. With the metric
(4.19), the Hamilton-Jacobi equation g

S)(

S) = m
2
c
2
acquires the
form
_
1
R
S
r
_
1
_
S
ct
_
2

_
1
R
S
r
__
S
r
_
2

1
r
2
_
S

_
2
m
2
c
2
= 0 .
(4.29)
The solution is looked for by the Hamilton-Jacobi method described in Sec-
tion 3.9. With some constant energy E and constant angular momentum J,
we write
S = Et + J + S
r
(r). (4.30)
117
This, once inserted in (4.29), gives
S
r
=
_
dr
_
E
2
c
2
_
1
R
S
r
_
2

_
m
2
c
2
+
J
2
c
2
__
1
R
S
r
_
1
_
1/2
. (4.31)
By the method, r = r(t) is obtained from the equation
S
E
= constant,
from which comes
ct =
E
mc
2
_
dr
_
1
R
S
r
_
_
_
E
mc
2
_
2

_
1 +
J
2
m
2
c
2
r
2
_ _
1
R
S
r
_
. (4.32)
The trajectory is found from
S
J
= constant, which gives
=
_
J
r
2
dr
_
E
2
c
2

_
m
2
c
2
+
J
2
r
2
_ _
1
R
S
r
_
. (4.33)
This leads to an elliptic integral. We are putting the additive integration
constants, which merely x the origins of the coordinates ct and , equal to
zero.
We should compare the above results with their non-relativistic counter-
parts given in Eqs.(3.108), (3.109). However, in order to calculate the small
corrections given by the theory to the trajectories of the planets turning
around the Sun, it is wiser to make approximations in (4.31) before taking
the derivative
S
J
. We shall suppose radial distances very large with respect
to the Schwarzschild radius: R
S
<< r. We also change the integration vari-
able to r

r
2
rR
S
(and drop the pirmes afterwards). Writing E

for
the non-relativistic energy, we nd
S
r
=
_
dr
_
E

2
c
2
+ 2mE

+
1
r
_
2m
2
MG + 4E

MR
S
_

1
r
2
_
J
2

3
2
m
2
c
2
R
2
S
_
_
1/2
(4.34)
The term in 1/r
2
will produce a secular displacement of the orbit perihelium.
The remaining terms cause changes in the relationships between the four-
momentum of the particle and the newtonian ellipse. We shall be interested
only in the perihelium precession. The trajectory is determined by +
S
r
J
118
= constant. The variation of S
r
in one revolution is, in the approximation
considered,
S
r
= S
(0)
r

3m
2
c
2
R
2
S
4J
S
r
J
.
S
(0)
r
is the closed ellipse case. The variation of the angle in one revolution
will be
=
S
r
J
.
Taking into account that

S
(0)
r
J
=
(0)
= 2 ,
we nd
= 2 +
3m
2
c
2
R
2
S
2J
2
= 2 +
6G
2
m
2
M
2
c
2
J
2
.
The last piece gives the precession . It is usual to express it in terms of the
ellipse parameters. If the great axis is a and the eccentricity is e, we have
J
2
GMm
2
= a(1 e
2
) .
The perihelium precession is then
=
6GM
a(1 e
2
)c
2
.
For the Earth, this is a very small variation: in seconds of arc, 3.8

per
century. For Mercury, it is 43.0

per century. This is in good agreement with

measurements.
4.13 Light-ray deviation The eikonal equation (3.18) is just the Hamil-
ton-Jacobi equation (3.17) with m = 0. The trajectory will still be given
by Eq.(4.33). The interpretation is, of course, quite another. S is now the
eikonal, the energy E must be replaced by the light frequency
0
, and it is
convenient to introduce a new constant, the impact parameter = Jc/
0
.
Thus,
=
_
dr
r
2
_
1

2

1
r
2
_
1
R
S
r
_
. (4.35)
119
This gives r = / cos a straight line passing at a distance of the
coordinate origin in the non-relativistic case R
S
= 0. The procedure to
analyse the small corrections due to R
S
= 0 is analogous to that used for
m = 0. We go back to Eq. (4.31),
S
r
(r) =

0
c
_
dr
_
r
2
(r R
S
)
2

2
_
r
2
rR
S
_
1
_
1/2
. (4.36)
With the same transformations used previously, this becomes
S
r
(r) =

0
c
_
dr
_
1
2
r
2
+ 2R
S
r
1
. (4.37)
Expanding in powers of R
S
r
1
,
S
r
S
R
S
=0
r
+
R
S

0
c
_
dr
_
r
2

2
= S
R
S
=0
r
+
R
S

0
c
arccosh
r

. (4.38)
The deviation undergone by a ray coming from a large distance R down to
a distance and then again to the same distance R will be
S
r
= S
R
S
=0
r
+ 2
R
S

0
c
arccosh
r

. (4.39)
To get the variation in the angle , it is enough to take the derivative with
respect to J =
0
/c:
=
S
r
J
=
S
R
S
=0
r
J
+ 2
R
S
R

_
R
2

2
. (4.40)
The term corresponding to the straight line has = . Taking the asymp-
totic limit R ,
= + 2
R
S

. (4.41)
This gives a deviation towards the centre of an angle
= 2
R
S

=
4GM
c
2
. (4.42)
For a light ray grazing the Sun, this gives = 1.75

. This eect has been

observed by a team under the leadership of Eddington during the 1919 solar
eclipse, at Sobral. It has been considered the rst positive experimental test
of General Relativity.
120
4.14 The event horizon We can compare the energy of the particle as
seen by an observer using the Schwarzschild coordinate system with the its
energy as seen in the proper frame. The rst is, as previously seen, E =
S
t
,
and the second is E
0
=
S

= mc
2
. But
S
t
=

t
S

, so that E =

g
00
mc
2
,
or
E = mc
2
_
1
R
S
r
. (4.43)
Let us examine the case of a particle falling in purely radial motion ( = 0,
= 0, J = 0) towards the center r = 0. If it starts from a point r
0
at
the instant t
0
, its energy will be E = mc
2
_
1
R
S
r
0
. At a moment t and a
distance r, Eq.(4.32) gives
c(t t
0
) =
_
1
R
S
r
0
_
r
0
r
dr

_
1
R
S
r

_
_
R
S
r

R
S
r
0
. (4.44)
This coordinate time diverges when r R
S
. Seen from an external observer,
the particle will take an innite time to arrive at the Schwarzschild radius.
On the other hand, we can calculate the proper interval of time for the same
thing to happen: it will be
c(
0
) =
_
r
0
r
ds =
_
r
0
r
dr

g
00
c
2
_
dt
dr

_
2
+ g
11
.
From (4.44),
dt
dr
=
_
1
R
S
r
0
dr
_
1
R
S
r
_
_
R
S
r

R
S
r
0
g
00
c
2
_
dt
dr
_
2
+ g
11
=
1
R
S
r

R
S
r
0
.
Thus,
c(
0
) =
_
r
0
r
dr

_
R
S
r

R
S
r
0
. (4.45)
This is a convergent integral, giving a nite value when r R
S
. The
particle, looking at things from its own proper frame and measuring time in
121
its own clock, arrives at the Schwarzschild radius in a nite interval of time.
It will even arrive at the center r = 0 in nite time.
Suppose now a gas of particles, each one in the conditions above. They
will all fall towards the center. Each will do it in a nite interval of time in
its own frame. The gas will eventually colapse. A quite distinct picture will
be seen from the coordinate, asymptotically at frame. Seen from a distant
observer, the particles will never actually traverse the Schwarzschild radius.
All that happens inside the Schwarzschild sphere lies beyond the innity of
time. The sphere is a horizon, technically called an event horizon.
We have seen (in 3.37 and ensuing paragraphs) how coordinate time
and proper time are related. Here we have an extreme dierence, given by
Eq.(4.28). Consider a light source near the Schwarzschild radius. Seen from
a distant observer, it will be strongly red-shifted. It will actually be more
and more red-shifted as the source is closer and closer to the radius. The
red-shift will tend to innity at the radius itself: the external observer cannot
receive any signal from the sphere surface.
4.15 In consequence of all that has been said, a massive object can even-
tually fall inside its own Schwarzschild sphere. In a real star, such a gravi-
tational colapse is kept at bay by the centrifugal eects of pressure caused
by the energy production through nuclear fusion. A massive enough star
black
hole
can, however, collapse when its energy sources are exhausted. Once inside
the radius, no emission will be able to scape and reach an external observer.
Particles and radiation will go on falling down to the sphere, but nothing will
get out. Such a collapsed object, such a black hole, will only be observable by
indirect means, as the emission produced by those external particles which
are falling down the gravitational eld.
4.16 It should be said, however, that the singularity in the components
of the metric does not imply that the metric itself, an invariant tensor, be
singular. A real singularity must be independent of coordinates and should
manifest itself in the invariants obtained from the metric. For example, the
determinant is an invariant. It is g = r
4
sin
2
. This shows that the point
r = 0 is a real singularity, in which the metric is no more invertible. The
Schwarzschild sphere, however, is not. It is, as seen, an event horizon, but
122
not a singularity. This seems to have been rst noticed by Lematre in 1938,
and can be veried by transforming to other coordinate systems. Take for
instance the family of transformations
ct

= ct
_
f(r)dr
1
R
S
r
; r

= ct +
_
dr
(1
R
S
r
)f(r)
, (4.46)
involving an arbitrary function f(r) and which lead to
ds
2
=
1
R
S
r
1 f(r)
2
(c
2
dt

2
f
2
dr

2
) r
2
(d
2
+ sin
2
d
2
) .
The Schwarzschild singularity will disappear for an f such that f(
R
S
r
) = 1.
The better choice is f(r) =
_
R
S
r
, which gives a synchomous system. Notice
that there are two possible choices of the signs in (4.46). One leads to an
expanding reference system, the other to a contracting frame. In eect, the
upper signs in (4.46) give
r

ct

=
_
dr
(1 f(r)
2
)
(1
R
S
r
)f(r)
=
_
dr
f(r)
=
_
dr
r
R
S
=
2r
3/2
R
1/2
S
.
This shows that the system contracts in the old system:
r = R
1/3
S
_
3
2
(r

ct

)
_
2/3
. (4.47)
As r 0, r

ct

, the equality corresponding to the center real singularity.

The singularity would correspond to the value
3
2
(r

ct

) = R
S
.
The interval takes the form given by Lematre
Lematre
line
element
I
ds
2
= c
2
dt

dr

2
_
3
2R
S
(r

ct

)
_
2/3

_
3
2
(r

ct

)
_
4/3
R
2/3
S
(d
2
+ sin
2
d
2
) .
(4.48)
The Schwarzschild singularity does not turn up. To get a feeling, we can use
Eq.(4.47) to rewrite the interval in mixed coordinates:
ds
2
= c
2
dt

R
S
r
dr

2
r
2
(d
2
+ sin
2
d
2
) . (4.49)
123
r

=

0
r

=

R
S
t'
r'
r

=

c
o
n
s
t
a
n
t
Figure 4.1: In Lematre coordinates, a particle can go through the
Schwarzschild radius in a nite amount of time. The cones become narrower
as r

decreases.
By Eq.(4.47), to each value r = constant in the old coordinates corre-
sponds a straight line r

= a +ct

in the new system (see Figure 4.1). These

lines indicate, consequently, immobile particles in the old reference frame.
Lematres system is synchonous ( 3.44). Geodesics are represented by
vertical lines, as the dashed line in the diagram. A free particle can fall
through the Schwarzschild radius and attain the center at r = 0.
Light cones have an interesting behavior. Consider a radial (d = 0, d =
0) light ray. Equation (4.49) gives, for ds = 0, the two (future and past)
cones, one for each sign in
c
dt

dr

=
_
R
S
r
. (4.50)
We see that the cone solid angle becomes smaller for smaller values of r.
Consider again the lines r = constant, indicating immobility in the original
system of coordinates. Their inclination is c
dt

dr

= 1, and there are two

possibilities:
for r > R
S
, then

c
dt

dr

< 1: a line r = constant through the vertex lies

inside the cone; immobility is consequently causally possible;
124
for r < R
S
, then

c
dt

dr

> 1: a line r = constant through the vertex lies

outside the cone; immobility is causally impossible; any particle will
fall to the center.
The cones dened in (4.50) have one more curious behavior: they will deform
themselves progressively as they approach the center. The derivative becomes
innite at r = 0. The lines representing the cones in the Figure will meet
the line r = 0 as vertical lines.
If we choose the lower signs in (4.46), the line element will be (4.48), but
with t

:
Lematre
line
element
II
ds
2
= c
2
dt

dr

2
_
3
2R
S
(r

+ ct

)
_
2/3

_
3
2
(r

+ ct

)
_
4/3
R
2/3
S
(d
2
+ sin
2
d
2
) .
(4.51)
An analogous discussion can be done concerning the behavior of cones and
the issue of immobility. Immobility is still forbidden inside the radius. The
dierence is that, instead of falling fatally towards the center, a particle will
inexorably draw away from it. Contrary to the previous case, the system
expands in the old system:
r = R
1/3
S
_
3
2
(r

+ ct

)
_
2/3
. (4.52)
4.17 A reference frame is complete when the world line of every particle
either go to innity or stop at a true singularity. In this sense, neither of the
above coordinate systems is complete. The Schwarzschild coordinates do not
apply to the interior of the sphere. An outside particle in the contracting
Lematre system can only fall down towards the centre: initial conditions
in the opposite sense are not allowed. Just the contrary happens in the
expanding Lematre systems. Both leave some piece of space unattainable.
That a complete system of coordinates does exist was rst shown by Kruskal
and Fronsdal. We shall here only mention a few aspects of this question.

In such kind of system no singularity at all appears in the Schwarzschild

radius. The coordinates are given as implicit functions of those used above.
125
In a form given by Novikov, the metric is looked for in a form generalizing
Novikov
line
element
the mixed-coordinate interval (4.49). New time and radial variables , R are
dened so as to put the line element in the form
ds
2
= c
2
d
2

_
1 +
R
2
S
R
2
_
dr
2
r
2
(, R)(d
2
+ sin
2
d
2
) (4.53)
and supposing a dust gas as source. The coordinates , R are given implicitly
and in parametric form by
r =
R
S
2
_
1 +
R
2
S
R
2
_
(1 cos ); (4.54)
=
R
S
2
_
1 +
R
2
S
R
2
_
3/2
( + sin ). (4.55)
The parameter take values in the interval 2, 0. When it runs from 2 to
0 the time variable increases monotonically, while r increases from zero up
to a maximum value
r = R
S
_
1 +
R
2
S
R
2
_
(4.56)
and then decreases back to zero. The Kruskal diagram of Figure 4.2 suma-
Kruskal
diagram
rizes the whole thing. Coordinates , R are complete, so that all situations
are described. The Schwarzschild coordinates describe only situations exter-
nal to the line r = R
S
. Contracting Lematre coordinates cover the shaded
area, expanding coordinates cover the domain which appears shaded after
specular vertical reection. The small arrows indicate the forcible sense par-
ticles follow inside the Schwarzschild sphere: contracting in the upper side,
expanding in the lower one.
We have seen in 3.17 that dust particles follow geodesics. The system is
synchronous ( 3.44), so that such geodesics are the vertical lines R =constant
in the diagram. An example is shown as a dashed line. Starting at = 0

Details are given in an elegant form in L.D. Landau & E.M. Lifshitz, Theorie des
Champs, 4th french edition, 103.
126
r

=

0
r

=

R
S
r

=

R
S
r

=

0
1
2
3
R

Figure 4.2: Kruskal diagram.

(point 1 in the Figure), a particle attains the center r = 0 in a nite time
=
R
S
2
_
1 +
R
2
S
R
2
_
3/2
. (4.57)
Starting with outward initial conditions at the center, a situation correspond-
ing to the lower part of the diagram, a particle will follow a trajectory like
the dashed line. It will cross the Schwarzschild radius in the outward direc-
tion, attain a maximal coordinate distance given by (4.56) at = 0 (point
1 again), and then fall back, crossing the Schwarzschild radius in an inward
progression at point 2 and reaching the center at point 3.
127
4.3 Large Scale Solutions
4.18 Two of the four fundamental interactions of Nature the weak and
the strong are of very short range. Electromagnetism has a long range
but, as opposite charges attract each other and tend thereby to neutral-
ize, the eld is, so to speak, compensated at large distances. Gravitation
remains as the only uncompensated eld and, at large scales, dominates.
This is why Cosmology is deeply involved with it. We shall now examine
two of the main solutions of Einsteins equations which are of cosmologi-
cal import. The Friedmann solution provides the background for the so-
called Standard Cosmological Model, which despite some diculties
Standard
model
gives a good description of the large-scale Universe during most of its evo-
lution. It represents a spacetime where time, besides being separated from
space, is positionindependent, and space is homogeneous and isotropic at
each point. The Universe begins with a high-density singularity (the big-
bang) and evolves through two main periods, one radiation-dominated and
another matter-dominated, which lasts up to present time. There are dif-
culties both at the beginning and at present time. The rst problem is
mainly related to causality, and could be solved by a dominant cosmological
term. The second comes from the recent observation that a cosmological
term is, even at present time, dominant. It is consequently of interest to
examine models in which the cosmological term gives the only contribution.
These are the de Sitter universes, which have an additional theoretical ad-
vantage: all calculation are easily done. A more detailed account of both the
Friedmann and the de Sitter solutions, mainly concerned with cosmological
aspects, can be found in the companion notes Physical Cosmology.
4.3.1 The Friedmann Solutions
4.19 We have seen under which conditions the second term in the interval
ds
2
= g
00
dx
0
dx
0
+ g
ij
dx
i
dx
j
represents the 3-dimensional space. We shall
suppose g
00
to be spaceindependent. In that case, the coordinate x
0
can be
chosen so that the time piece is simply c
2
dt
2
, where t will be the coordinate
time. It will be a universal time, the same at every point of space.
The Universe, that is, the space part, is supposed to respect the Cos-
128
mological Principle, or Copernican Principle: it is homogeneous as a whole.
Homogeneous means looking the same at each point. Once this is accepted,
imposing isotropy around one point (for instance, that point where we are) is
enough to imply isotropy around every point. This means in particular that
space has the same curvature around each point. There are only 3 kinds of
3dimensional spaces with constant curvature:
the sphere S
3
, a closed space with constant positive curvature;
the open hyperbolic space S
2,1
, or (a pseudosphere, or sphere with
imaginary radius), whose curvature is negative; and
the open euclidean space E
3
of zero curvature (that is, at).
These three types of space are put together with the help of a parameter
k: k = +1 for S
3
, k = 1 for S
2,1
and k = 0 for E
3
. The 3dimensional line
element is then, in convenient coordinates,
dl
2
=
dr
2
1 kr
2
+ r
2
d
2
+ r
2
sin
2
d
2
. (4.58)
The last two terms are simply the line element on a 2-dimensional sphere S
2
of radius r a clear manifestation of isotropy. Notice that these symmetries
refer to space alone. The radius can be time-dependent.
4.20 The energy content is given by the energy-momentum of an ideal
uid, which in the Standard Model is supposed to be homogeneous and
isotropic:
T

= (p + c
2
) U

p g

. (4.59)
Here U

is the four-velocity and = /c

2
is the mass equivalent of the energy
density. The pressure p and the energy density are those of the matter (visible
plus dark) and radiation. When p = 0, the uid reduces to dust.
4.21 We can now put together all we have said. Instead of a time
dependent radius, it is more convenient to use xed coordinates as above,
129
and introduce an overall scale parameter a(t) for 3space, so as to have the
spacetime line element in the form
ds
2
= c
2
dt
2
a
2
(t)dl
2
. (4.60)
Thus, with the high degree of symmetry imposed, the metric is entirely xed
Friedmann
Robertson
Walker
interval
by the sole function a(t). The spacetime line element will be
ds
2
= c
2
dt
2
a
2
(t)
_
dr
2
1 kr
2
+ r
2
d
2
+ r
2
sin
2
d
2
_
. (4.61)
This is the FriedmannRobertsonWalker line element.
4.22 The above line element is a pure consequence of symmetry consid-
erations. We have now to impose the dynamical equations. Recent data,
as said, point to a non-vanishing cosmological constant. It is, consequently,
wiser to use the Einstein equations in the form (3.42). The extreme sim-
plicity of the model is reected in the fact that those 10 partial dierential
equations reduce to 2 ordinary dierential equations (in the variable t) for
the scale parameter.
In eect, once (4.59) is used, the eld equations (3.42) reduce to the two
Friedmann equations for a(t):
Friedmann
equations
a
2
=
_
2
_
4G
3
_
+
c
2
3
_
a
2
kc
2
; (4.62)
a =
_
c
2
3

4G
3
_
+
3p
c
2
__
a(t) . (4.63)
It will be convenient to absorb the length dimensionality in a(t), so that the
variable r can be seen as dimensionless. The cosmological constant has the
dimension (length)
2
. The second equation determines the concavity of the
function a(t). This has a very important qualitative consequence when
= 0. In that case, for normal sources with > 0 and p 0, a is forcibly
Big
Bang
negative for all t and the general aspect of a(t) is that of Figure 4.3. It will
consequently vanish for some time t
initial
. Distances and volumes vanish at
that time and densities become innite. This moment t
initial
is taken as the
130
beginning, the Big Bang itself. It is usual to take t
initial
as the origin of
the time coordinate: t
initial
= 0. If > 0, there is a competition between
the two terms. It may even happen that the scale parameter be = 0 for all
nite values of t.
1.2 1.4 1.6 1.8 2 2.2 2.4
t
4.5
5
5.5
6
6.5
7
a
Figure 4.3: Concavity of a(t) for = 0.
Combining both Friedmann equations, we nd
d
dt
= 3
a
a
_
+
p
c
2
_
, (4.64)
which is equivalent to
d
da
(a
3
) + 3 p a
2
= 0. (4.65)
This equation reects the energymomentum conservation: it can be alter-
natively obtained from T

;
= 0.
4.23 It is convenient to introduce the Hubble function
Hubble
function
&
constant
H(t) =
a(t)
a(t)
=
d
dt
ln a(t) , (4.66)
whose present-day value is the Hubble constant
H
0
= 100 h km s
1
Mpc
1
= 3.24 10
18
h s
1
.
131
The parameter h, of the order of unity, encapsulates the uncertainty in
present-day measurements, which is large (0.45 h 1). Another func-
tion of interest is the deceleration
q(t) =
aa
a
2
=
a
aH(t)
=
1
H
2
(t)
a
a
. (4.67)
Equivalent expressions are

H(t) = H
2
(t) (1 + q(t)) ;
d
dt
1
H(t)
= 1 + q(t) . (4.68)
Uncertainty is very large for the deceleration parameter, which is the present
day value q
0
= q(t
0
). Data seem consistent with q
0
0. Notice that we are
using what has become a standard notation, the index 0 for present-day
values: H
0
for the Hubble constant, t
0
for present time, etc.
The Hubble constant and the deceleration parameter are basically inte-
gration constants, and should be xed by initial conditions. As previously
said, the presentday values are used.
4.24 The at Universe Let us, as an exercise, examine in some detail
the particular case k = 0. The FriedmannRobertsonWalker line element is
simply
ds
2
= c
2
dt
2
a
2
(t)dl
2
, (4.69)
where dl
2
is the Euclidean 3-space interval. In this case, calculations are
much simpler in cartesian coordinates. The metric and its inverse are
(g

) =
_
_
_
_
_
1 0 0 0
0 a
2
(t) 0 0
0 0 a
2
(t) 0
0 0 0 a
2
(t)
_
_
_
_
_
;
(g

) =
_
_
_
_
_
1 0 0 0
0 a
2
(t) 0 0
0 0 a
2
(t) 0
0 0 0 a
2
(t)
_
_
_
_
_
.
132
In the Christoel symbols

, the only nonvanishing derivatives are those

with respect to x
0
. Consequently, only Christoel symbols with at least one
index equal to 0 will be nonvanishing. For example,

k
ij
= 0. Actually, the
only Christoels = 0 are:

k
0j
=
k
j
1
c
a
a
;

0
ij
=
ij
1
c
a a.
The nonvanishing components of the Ricci tensor are
R
00
=
3
c
2
a
a
= 3
H
2
(t)
c
2
q(t) ;
R
ij
=

ij
c
2
[a a + 2 a
2
] =

ij
c
2
a
2
H
2
(t)[2 q(t)] .
In consequence, the scalar curvature is
R = g
00
R
00
+ g
ij
R
ij
= 6
_
a
c
2
a
+
_
a
ca
_
2
_
= 6
H
2
(t)
c
2
[q(t) 1] .
The nonvanishing components of the Einstein tensor G

=R

1
2
Rg

are
G
00
= 3
_
a
ca
_
2
; G
ij
=

ij
c
2
[ a
2
+ 2a a] .
Let us consider the sourceless case with cosmological constant. The Einstein
equations are then
G
00
g
00
= 3
_
a
ca
_
2
= 0 ; G
ij
g
ij
=

ij
c
2
[ a
2
+ 2a a a
2
] = 0 .
Subtracting 3 times one equation from the other, we arrive at the equivalent
set
a
2

c
2
3
a
2
= 0 ; a
c
2
3
a = 0 . (4.70)
These equations are (4.62) and (4.63) for the case under consideration. Of
the two solutions, a(t) = a
0
e
H
0
(tt
0
)
, only
a(t) = a
0
e
H
0
(tt
0
)
= a
0
e
_
c
2
3
(tt
0
)
(4.71)
133
would be consistent with expansion. Expansion is a fact well established
by observation. This is enough to x the sign, and the model implies an
everlasting exponential expansion. Notice that the scalar curvature is R =
4, as is always the case in the absence of sources. Equation (4.71) is
actually a de Sitter solution. The quick growth has been called ination
and is supposed to have taken place in the very early history of the Universe.
4.25 Thermal History The presentday content of the Universe consists
of matter (visible or not) and radiation, the last constituting the cosmic mi-
crowave background. The energy density of the latter is very small, much
smaller than that of visible matter alone. Nevertheless, it comes from the
equations of state that radiation energy increases faster than matter energy
with the temperature. Thus, though matter dominates the energy content
of the Universe at present time, this dominance ceases at a turning point
time in the past. At that point radiation takes over. At about the same time,
hydrogen the most common form of matter ionizes. The photons of
the background radiation establish contact with the electrons and the whole
system is thermalized. Above that point, there exists a single temperature.
And, above the turning point, the dominating photons increase progressively
in number while their concentration grows by contraction. The opportunity
for interactions between them becomes larger and larger. When they ap-
proach the mass of an electrons, pair creation sets up as a stable process.
Radiation is now more than a gas of photons: it contains more and more
electrons and positrons. Concomitantly, nucleosynthesis stops. As we insist
in going up the temperature ladder, the photons, which are more and more
energetic, break the composite nuclei. The nucleosynthesis period is the most
remote time from which we have reasonably sure information nowadays.
The Standard model starts from present-day data and moves to the past,
taking into account the changes in the equations of state. It goes conse-
quently from a matter-dominated era through the time of hydrogen recom-
bination, then to the changeover period in the which radiation establishes
its dominance. These successive eras, the so-called thermal history of
the Universe, is analysed in the text Physical Cosmology. We shall here only
examine the de Sitter solutions also of fundamental cosmological interest
134
because, by their simplicity, they give a beautiful example in which all
calculation can be done without much ado.
4.3.2 de Sitter Solutions
4.26 de Sitter spacetimes are hyperbolic spaces of constant curvature.
They are solutions of vacuum Einsteins equation with a cosmological term.
There are two dierent kinds of them: one with positive scalar Ricci curva-
ture, and another one with negative scalar Ricci curvature. As the calcula-
tions are remarkably simple, we shall give a fairly detailed account.
We shall denote by R the de Sitter pseudo-radius, by

(, , =
0, 1, 2, 3) the Lorentz metric of the Minkowski spacetime, and
A
(A, B, . . . =
0, . . . , 4) will be the Cartesian coordinates of the pseudo-Euclidean 5spaces.
There are two types of spacetime named after de Sitter:
1. de Sitter spacetime dS(4, 1): hyperbolic 4-surface whose inclusion
in the pseudoEuclidean space E
4,1
satisfy

AB

A

B
=

4
_
2
= R
2
. (4.72)
It is a one-sheeted hyperboloid (a 4-dimensional version of the surface
seen in 2.65) with topology R
1
S
3
, and within our conven-
tions negative scalar curvature. Its group of motions is the pseudo
orthogonal group SO(4, 1)
2. antide Sitter spacetime dS(3, 2): hyperbolic 4-surface whose in-
clusion in the pseudoEuclidean space E
3,2
satisfy

AB

A

B
=

+
_

4
_
2
= R
2
. (4.73)
It is a two-sheeted hyperboloid (a 4-dimensional version of the space
seen in 2.64) with topology S
1
R
3
, and positive scalar curvature.
Its group of motions is SO(3, 2)
With the notation
44
= s, both de Sitter spacetimes can be put to-
gether in

AB

A

B
=

+ s
_

4
_
2
= sR
2
, (4.74)
where we have the following relation between s and the de Sitter spaces:
135
s = 1 for dS(4, 1)
s = +1 for dS(3, 2).
4.27 The metric Let us nd now the line element of the de Sitter spaces.
The most convenient coordinates are the stereographic conformal. The pas-
sage from the Euclidean
A
to the stereographic conformal coordinates x

(, , = 0, 1, 2, 3) is done by the transformation:

= x

;
4
= R(1 2), (4.75)
(a sign in the last expression would have no consequence for what follows)
with (x) a function of x

which we shall determine. Two expressions

preparatory to the calculation of the line element can be immediately ob-
tained by taking dierentials:
d

= x

d + dx

, (4.76)
and
_
d
4
_
2
= 4R
2
d
2
. (4.77)
Let us introduce
2
=

and rewrite the dening relation (4.74) as

2
+ s
_

4
_
2
= sR
2
. (4.78)
Equating (
4
)
2
got from (4.75) and (4.78), we nd
=
1
1 + s

2
4R
2
. (4.79)
Notice that from this expression it follows that
d = s

2
2R
2
d = s

2
4R
2
2

dx

, (4.80)
from which another preparatory result is obtained:
2

dx

= s
4R
2

2
d . (4.81)
Now, the de Sitter line element is
d
2
=
AB
d
A
d
B
=

+ s
_
d
4
_
2
,
136
or, by using (4.76),
d
2
=

(x

d + dx

)
_
x

d + dx

_
+ s
_
d
4
_
2
.
Expanding and using (4.77),
d
2
=

d
2
+ 2

dx

d +
2

dx

dx

+ s 4R
2
d
2
.
Now, using (4.81),
d
2
=
2

dx

dx

+
_

2
s
4R
2

_
d
2
+ s 4R
2
d
2
,
and then (4.79),
d
2
=
2

dx

dx

+
_
s 4R
2

d
2
+ s 4R
2
d
2
,
so that nally
d
2
= g

dx

dx

, (4.82)
where the metric g

is
g

=
2

=
1
_
1 + s

2
4R
2
_
2

. (4.83)
The de Sitter spaces are, therefore, conformally at (see 2.46), with the
conformal factor given by
2
(x).
4.28 The Christoel symbol corresponding to a conformally at metric
g

with conformal factor

2
(x) has the form

=
_

ln (x) . (4.84)
Taking derivatives in (4.79), we nd for the de Sitter spaces

= s

2R
2
_

. (4.85)
The Riemann tensor components can be found by taking the following
steps. First, take the derivative of the de Sitter connection

= s

2R
2
_

ln
137
= s

2R
2
_

sx

2R
2

= s

2R
2
_

+
x

4R
4
_

= s

2R
2
_

+
x

4R
4

2
_

= s

2R
2
_

+
1
4R
4

2
_

.
Indicating by [] the antisymmetrization (without any factor) of the in-
cluded indices, we get

= s

2R
2
_

]

[

]
_
+
1
4R
4

2
_
x

[
x
]
x

g
[
x
]
_
= s

R
2

]
+
1
4R
4

2
_
x

[
x
]
x

g
[
x
]
_
.
This is the contribution of the derivative terms.
The product terms are

=

2
4R
4
_

s
+

x
s
;

=
1
4R
4

2
_
x

[
x
]
+ x

x
[
g
]
+
2

2
g
[

.
A provisional expression for the curvature is, therefore,

= s

R
2

]
+
1
4R
4

2
_
x

[
x
]
x

g
[
x
]
_
+
1
4R
4

2
_
x

[
x
]
+ x

x
[
g
]
+
2

2
g
[

.
The rst two terms in the last line just cancel the last two in the line above
them. Therefore,

= s

R
2

]
+
1
4R
4

2
g
[

]
= s

R
2

]
+
1
4R
4

]
=
_
s

R
2
+
1
4R
4

]
= s

R
2
_
s
4R
2

2
1

]
138
Using (4.79), we nd that the bracketed term is = . Therefore, we get

= s

2
R
2

]
=
s
R
2
_

. (4.86)
The Ricci tensor will be

R
=
3 s
R
2
g

(4.87)
and the scalar curvature,

R
=
12 s
R
2
. (4.88)
The de Sitter spacetimes are spaces of constant curvature. We can now
make contact with the cosmological term. From the expressions above, we
nd that

R

1
2
g

R
+
3 s
R
2
g

= 0. (4.89)
Thus, the de Sitter spaces are solutions of the sourceless Einsteins equations
with a cosmological constant
=
3 s
R
2
=

R
/4. (4.90)
Notice the relationships to the de Sitter and the anti-de Sitter spaces:
s = 1 for the de Sitter space dS(4, 1) > 0
s = +1 for the anti-de Sitter space dS(3, 2) < 0 .
4.29 We have said ( 3.46) that positive scalar curvature tends to make
curves to close to each other, and just the contrary for negative curvature.
The relative sign in (4.90) shows that the cosmological constant has the
opposite eect: > 0 leads to diverging curves, < 0 to converging ones.
This actually depends on the initial conditions. Let us look at the geodetic
deviation equation,
D
2
X

Du
2
=

. (4.91)
Using Eq.(4.86),
139
D
2
X

Ds
2
=
s
R
2
_

=
s
R
2
[

] X

=
s
R
2
h

. (4.92)
We have recognized the transversal projector h

(see 3.50). By the

geodesic equation, the component of X along U will have vanishing con-
tributions tothe left-hand side. Consequently, only the transversal part X

will appear in the equation, which is now

D
2
X

Ds
2
+
s
R
2
X

=
D
2
X

Ds
2
+

R
12
X

=
D
2
X

Ds
2


3
X

= 0 . (4.93)
Negative leads to oscillatory solutions. Positive can lead both to con-
tracting and expanding congruences. If two lines are initially separating,
they will separate indenitely more and more.
4.30 We have been using carefully two coordinate systems. The most
convenient system for cosmological considerations is the socalled comoving
system, in which the Friedmann equations, in particular, have been writ-
ten. In that system the scale parameter appears in its utmost simplicity.
The stereographic coordinates are of interest for de Sitter spaces. We could
perform a transformation between the two systems, but that is not really
necessary: we have only taken scalar parameters from one system into the
other. The only exception, Eq. (4.89), is a tensor which vanishes in a system
and, consequently, vanishes also in the other.
The expression (4.83) for the metric is very dierent from the original
one. De Sitter has found it in another coordinate system, in the form
ds
2
=
_
1

3
r
2
_
c
2
dt
2
r
2
(d
2
+ sin
2
d
2
)
dr
2
1

3
r
2
. (4.94)
This expression is just the Schwarzschild solution in the presence of a
term, Eq.(4.21), when the source mass tends to zero. There are many other
coordinates and metric expressions of interest for dierent aims.

One of
them exhibits clearly the inationary property above discussed:
ds
2
= c
2
dt
2
e
ct/R
_
dx
2
+ dy
2
+ dz
2
_
. (4.95)

A few, included that given below, are given by R.C. Tolman, Relativity, Thermody-
namics and Cosmology, Dover, New York, 1987, 142.
140
Chapter 5
Tetrad Fields
5.1 Tetrads
For each source, set of symmetries and boundary conditions, there will be a
dierent solution of Einsteins equations. Each solution will be a spacetime.
It will be interesting to go back and review the general characteristics of
spacetimes. While in the process, we shall revisit some previously given
notions and reintroduce them in a more formal language.
5.1 A spacetime S is a fourdimensional dierentiable manifold whose
tangent space ( 2.24) at each point is a Minkowski space. Loosely speaking,
we implant a Lorentz metric
ab
on each tangent space. Bundle language
is more specic: it considers the tangent bundle TS on spacetime, an 8-
dimensional space which is locally the direct product of S and a typical
bre representing the tangent space. For a spacetime, the typical ber is
the Minkowski space M. The ber is typical because it is an ideal (in
the platonic sense) Minkowski space. The relationship between the typical
Minkowski ber and the spaces tangent to spacetime is established by tetrad
elds (see Figure 5.1). A tetrad eld will determine a copy of M on each
tangent space. M is considered not only as a at pseudo-Riemannian space,
but as a vector space as well. We are going to use letters of the latin alphabet,
a, b, c, . . . = 0, 1, 2, 3, to label components on M, and those of the greek
alphabeth, , , , . . . = 0, 1, 2, 3 for spacetime components. The rst will be
called Minkowski indices, the latter Riemann indices.
141
5.2 We shall need an initial vector basis in M. We take the simplest one,
the standard canonical basis
K
0
= (1, 0, 0, 0) K
1
= (0, 1, 0, 0)
K
2
= (0, 0, 1, 0) K
3
= (0, 0, 0, 1) .
Each K
a
is given by the entries (K
a
)
b
=
ab
. Each tetrad eld h will be a
mapping
h : M TS, h(K
a
) = h
a
.
The four vectors h
a
will constitute a vector basis on S. Actually, this is a
T S
p
p
x
x
< >
- 1
< >
- 1
h
x
a
K
a
U
(}
x
a
a
h X
= h
k
j
j
K
k
K
Minkowski
space
ideal
h (tetrad)
h
S = spacetime
ta
n
g
e
n
t M
in
k
o
w
s
k
i s
p
a
c
e
Figure 5.1: The role of a tetrad eld.
local aair: given a point p S, and around it an euclidean open set U,
the h
a
will constitute a vector basis not only for the tangent space to S at
p, denoted T
p
S, but also for all the T
q
S with q U. The extension from p
142
to U is warranted by the dierentiable structure. The dual forms h
b
, such
that h
b
(h
a
) =
b
a
, will constitute a vector basis on the cotangent space at
p, denoted T
p
S

. The dual base {h

a
} can be equally extended to all the
points of U. For example, each coordinate system {x

} on U will dene a
natural vector basis {

= /x

}, with its concomitant covector dual

basis, {dx

}. This is a very particular and convenient tetrad eld, frequently

called a trivial tetrad, given by e
a
= e(K
a
) =
a

. It is usual to x a
coordinate system around each point p from the start, and in this sense this
basis is indeed natural.
Another tetrad eld, as the above generic {h
a
} and its dual {h
b
}, can be
written in terms of {

} as
h
a
= h
a

and h
b
= h
b

dx

, (5.1)
with
h
b

h
a

=
b
a
and h
a

h
a

. (5.2)
The components h
a

(x) have one label in Minkowski space and one in space-

time, constituting a matrix with the inverse given by h
a

(x). It is usual to
designate the tetrad by the sets {h
a

(x)} or {h
a

}, but their meaning should

be clear: they represent the components of a general tetrad eld in a natural
basis previously chosen.
5.3 The tetrad Minkowski labels are vector indices, changing under Lor-
entz transformations according to
h
a

=
a

b
h
b
, (5.3)
or, in terms of components,
h
a

(x) =
a

b
(x) h
b

(x) . (5.4)
For each tangent Minkowski space,
a

b
is constant. It will, however,
depend on the point of spacetime, which we indicate by its coordinates x
= {x

}. This is better seen if we contract the last expression above with

143
h
c

(x), to obtain the Lorentz transformation in terms of the initial and nal
tetrad basis:

b
(x) = h
a

(x) h
b

(x) . (5.5)
Equation (5.4) says that each tetrad component behaves, on each Minkowski
ber, as a Lorentz vector. For each xed Riemann index , h
a

transforms
according to the vector representation of the Lorentz group. A Lorentz (ac-
tually, co-)vector on a Minkowski space has components transforming, under
a Lorentz transformation with parameters
cd
, as

=
a

b
(x)
b
=
_
exp
_
i
2

cd
J
cd
_
a

b
. (5.6)
Here, each J
cd
is a 4 4 matrix representing one of the Lorentz group gen-
erators:
[J
cd
]
a

b
= i
_

cb

db

c
_
. (5.7)
This means also that, for each xed Riemann index , the h
a

s constitute a
Lorentz basis (or frame) for M.
5.4 A tetrad eld converts tensors on M into tensors on spacetime, trans-
literating one index at a time. A general Lorentz tensor T, transforming
according to
T
a

...
=
a

c
. . . T
abc...
,
will satisfy automatically T
a

...
= h
a

h
b

h
c

. . . T
...
, which shows how
tetrad elds can mediate Lorentz transformations. As an example of that
transmutation, a tetrad will produce a eld on spacetime out of a vector in
Minkowski space by

(x) = h
a

(x)
a
(x) . (5.8)
As the tetrad, in its Minkowski label, transforms under Lorentz transfor-
mation as a vector should do,

(x) is Lorentzinvariant. Another case of

tensor transliteration is

(x) = h
a

(x)
a

b
(x) h
b

(x) =

(5.9)
144
[using (5.5)]. Thus, there is no Lorentz transformation on spacetime itself.
The Minkowski indices also called tetrad indices are lowered by
the Lorentz metric
ab
:
h
a
=
ab
h
b

.
An important consequence is that the Lorentz metric
ab
is transmuted into
the Riemannian metric
g

=
ab
h
a

h
b

. (5.10)
Of course, also g

(as any component of a spacetime tensor) is Lorentz

invariant. Here, given
ab
, the tetrad eld denes the metric g

. Dier-
ent tetrad elds transmute the same
ab
into dierent spacetime pseudo-
Riemannian metrics.
5.5 The members of a general tetrad eld {h
a
}, as vector elds, will satisfy
a Lie algebra (see 2.34) with a commutation table
[h
a
, h
b
] = c
c
ab
h
c
. (5.11)
The structure coecients c
c
ab
measure its anholonomicity they are some-
times called anoholonomicity coecients and are given by
c
c
ab
= [h
a
(h
b

) h
b
(h
a

)] h
c

. (5.12)
If {h
a
} is holonomous, c
c
ab
= 0, then h
a
= dx
a
for some coordinate system
{x
a
}, and
dx
a

=
a

b
dx
b
.
Expression (5.10) would, in that holonomic case, give just the Lorentz metric
written in another system of coordinates. This is the usual choice when we
are interested only in Minkowski space transformations, because then
a

b
= x
a

/x
b
. In this case the tetrad components can be identied with the
Lame coecients of coordinate transformations, and the metric g

will be
the Lorentz metric written in a general coordinate system.
We have up to now left quite indenite the choice of the tetrad eld
itself. In fact its choice depends on the physics under consideration. Trivial
tetrads are relevant when only coordinate transformations are considered.
A nontrivial tetrad reveals the presence of a gravitational eld, and is the
fundamental tool in the description of such a eld.
145
5.2 Linear Connections
Let us examine, in a purely descroptive way, the transformation properties
of a linear connection. A linear connection is a 1-form with values in the
linear algebra, that is, the Lie algebra of the linear group GL(4, R) of all
invertible real 4 4 matrices. This means a matrix of 1-forms. A Lorentz
connection is a 1-form with values in the Lie algebra of the Lorentz group,
which is a subgroup of GL(4, R). All connections of interest to gravitation
(the Levi-Civita in particular) are ultimately linear connections.
5.2.1 Linear Transformations
5.6 A linear transformation of N variables x
r
is an invertible transforma-
tion of the type
x
r

= M
r

s
x
s
. (5.13)
In this case (M
r

s
) is an invertible matrix. Linear transformations form
groups, which include all the groups of matrices. The set of invertible NN
matrices with real entries constitutes a group, called the real linear group
GL(N, R).
5.7 Consider the set of NN matrices. This set is, among other things, a
vector space. The simplest of such matrices will be those

whose entries
are all zero except for that of the -th row and -th column, which is 1:
(

. (5.14)
An arbitrary N N matrix K can be written K = K

.
The

s have one great quality: they are linearly independent (none

can be written as linear combinations of the other). Thus, the set {

}
constitutes a basis (the canonical basis) for the vector space of the N N
matrices. An advantage of this basis is that the components of a matrix K
as a vector written in basis {

} are the very matrix elements: (K)

=
K

.
5.8 Consider now the product of matrices: it takes each pair (A, B) of ma-
trices into another matrix AB. In our notation, matrix product is performed
146
coupling lowerright indices to higherleft indices, as in
_

=
_

, (5.15)
where in (

, is the column index.

The structure of vector space ( 2.20) includes an addition operation and
its inverse, subtraction. Once the product is given, another operation can be
dened by the commutator [A, B] = AB BA, the subtraction of two prod-
ucts. The Lie algebra of the N N real matrices with the operation dened
by the commutator is called the real N-linear algebra, denoted gl(N, R). The
set {

} is called the canonical base for gl(N, R). A Lie algebra is sum-
marized by its commutation table. For gl(N, R), the commutation relations
are
_

= f
(

)(

)
(

. (5.16)
The constants appearing in the right-hand side are the structure coecients,
whose values in the present case are
f
(

)(

)
(

)
=

. (5.17)
The choice of index positions may seem a bit awkward, but will be convenient
for use in General Relativity. There, linear connections

and Riemann
curvatures R

will play fundamental roles. These notations are quite well

established. Now, a linear connection is actually a matrix of 1-forms =

, with the components

being usual 1-forms, just

dx

. The
rst two indices refer to the linear algebra, the last to the covector character
of . A Riemann curvature is a matrix of 2-forms, R =

, where each
R

is an usual 2-form R

=
1
2
R

dx

dx

. In both cases, we talk of

algebravalued forms. In consequence, the notation we use here for the
matrices seem the best possible choice: in order to have a good notation for
the components one must sacrice somewhat the notation for the base.
The invertible N N matrices constitute, as said above, the real linear
group GL(N, R). Each member of this group can be obtained as the ex-
ponential of some K gl(N, R). GL(N, R) is thus a Lie group, of which
gl(N, R) is the Lie algebra. The generators of the Lie algebra are also called,
by extension, generators of the Lie group.
147
5.2.2 Orthogonal Transformations
5.9 A group of continuous transformations preserving an invertible real
symmetric bilinear form (see 2.44) on a vector space is an orthogonal
group (or, if the form is not positivedenite, a pseudoorthogonal group).
A symmetric bilinear form is a mapping taking two vectors into a real
number: (u, v) =

, with

. It is represented by a symmetric
matrix, which can always be diagonalized. Consequently, it is usually pre-
sented in its simplest, diagonal form in terms of some coordinates: (x, x) =

. Thus, the usual orthogonal group in E

3
is the set SO(3) of rotations
preserving (x, x) = x
2
+y
2
+z
2
; the Lorentz group is the pseudoorthogonal
group preserving the Lorentz metric of Minkowski spacetime, (x, x) = c
2
t
2
- x
2
- y
2
- z
2
. These groups are usually indicated by SO() = SO(r, s), with
(r, s) xed by the signs in the diagonalized form of . The group of rotations
in n-dimensional Euclidean space will be SO(n), the Lorentz group will be
SO(3, 1), etc.
5.10 Given the transformation x

, to say that is preserved

is to say that the distance calculated in the primed frame and the distance
calculated in the unprimed frame are the same. Take the squared distance in
the primed frame,

, and replace x

and x

by their transformation
expressions. We must have

, x . (5.18)
This is the groupdening property, a condition on the

s. When is
an Euclidean metric, the matrices are orthogonal that is, their columns
are vectors orthogonal to each other. If is the Lorentz metric, the

s
belong to the Lorentz group. We see that it is necessary that

.
The matrix form of this condition is, for each group element ,

T
= , (5.19)
where
T
is the transpose of .
148
5.11 There is a corresponding condition on the members of the group Lie
algebra. For each member A of the algebra, there will exist a group member
such that = e
A
. Taking = I+A+
1
2
A
2
+. . . and
T
= I+A
T
+
1
2
(A
T
)
2
+. . .
in the above condition and comparing order by order, we nd that A must
satisfy
A
T
=
1
A (5.20)
and will consequently have vanishing trace: tr A = tr A
T
= - tr (
1
A )
= - tr (
1
A) = - tr A tr A = 0.
5.12 If is dened on an N-dimensional space, the Lie algebras so() of
the orthogonal or pseudoorthogonal groups will be subalgebras of gl(N, R).
Given an algebra so(), both basis and entry indices can be lowered and
raised with the help of . We dene new matrices

by lowering labels
with : (

. Their commutation relations become

[

] =

. (5.21)
The generators of so() will then be J

, with commutation
relations
[J

, J

] =

. (5.22)
These are the general commutation relations for the generators of the orthog-
onal or pseudoorthogonal group related to . We shall meet many cases in
what follows. Given , the algebra is xed up to conventions. The usual
group of rotations in the 3-dimensional Euclidean space is the special or-
thogonal group, denoted by SO(3). Being special means connected to the
identity, that is, represented by 3 3 matrices of determinant = +1.
The group O(N) is formed by the orthogonal NN real matrices. SO(N)
is formed by all the matrices of O(N) which have determinant = +1. In
particular, the group O(3) is formed by the orthogonal 3 3 real matrices.
SO(3) is formed by all the matrices of O(3) which have determinant = +1.
The Lorentz group, as already said, is SO(3, 1). Its generators have just the
algebra (5.22), with the Lorentz metric.
149
5.2.3 Connections, Revisited
5.13 Suppose we are given the connection by components
a
b
, the rst
two indices being algebraic and the last a Riemann index. This supposes a
basis in the linear algebra and a basis of vector elds on the manifold. Taking
the canonical basis {
a
b
} for the algebra, and a holonomic vector basis {dx

}
on the manifold, for example, the connection is given in invariant form by
=
1
2

a
b

a
b
dx

. (5.23)
For reasons which will become clear later, the set of components {
a
b
} will
be called spin connection.
Connections have been introduced in in Section 2.3 through their behavior
under coordinate tranformations, that is, under change of holonomic tetrads.
We proceed now to a series of steps extending that presentation to general,
holonomic or not, tetrads. First, we (i) change from Minkowski indices to
Riemann indices by

a
b

= h

a

a
b
h
b

+ h

h
a

. (5.24)
This generalizes Eq.(2.33). Then, we (ii) change again through a Lorentz
transformed tetrad,

= h
a

b
+ h
a

b
, which means
that

=
a

a

a
b
_

1
_
b
b
+
a

1
_
c
b
. (5.25)
This gives the eect on
a
b
of a Lorentz transformation
a

a
= h
a

a
. In
the notation adopted,
a

a
changes V
a
into V
a

; we can write simply

c
b
for
the inverse (
1
)
c
b
, understanding that
V
c
=
c
b
V
b

=
_

1
_
c
b
V
b

.
Now, start instead with

, and (iii) change from Riemann indices to

Minkowski indices by

a
b
= h
a

b
+ h
a

b
, (5.26)
and (iv) go back to modied Riemann indices by

= h

a

a
b
h
b

+ h

h
c

.
150
Consequently,

= h

a
h
a

b
+ h

a
h
a

b
_
h
b

+ h

h
c

,
or

= B

+ B

. (5.27)
This is the eect of a change of basis given by B

= h

a
h
a

.
5.14 Vector elds transform according to (5.6):

(x) = h
e

(x)

(x) =
e

b
h
b

(x)

(x) =
e

b
(x) . (5.28)
What happens to their derivatives? Clearly, they transform in another way:

b
_

b
+
e

b
.
As the name indicates, the covariant derivative of a given object is a
derivative modied in such a way as to keep, under transformations, just the
same behavior of the object. Here, it will have to obey
D

=
e

b
D

b
. (5.29)
5.15 The way physicists introduce a connection is as a compensating
eld, an object
a
b
with a very special behavior whose action on the eld,
once added to the usual derivative, produces a covariant result. In the present
case we look for a connection such that

+
e

=
e

b
_

b
+
b
d

.
A direct calculation shows then that the required behavior is just (5.25),
which can be written also as

=
a

d
_

d
c

+
d
c

c
b
. (5.30)
5.16 As a rule, all indexed objects are tensor components and transform
accordingly. A connection, written as
a
b
, is an exception: it is tensorial
in the last index, but not in the rst two, which change in the peculiar
151
way shown above. Any connection transforming in this way will lead to a
covariant derivative of the form

a
=

a
+
a
b

b
= h
e

_
e
e

a
+
a
be

. (5.31)
Applied in particular to h
a

, it gives

h
a

h
a

+
a
b
h
b

;
applied to its inverse h
a

h
a

h
a

b
a
h
b

.
It is easily checked that
a
b
=
i
2

cd

(J
cd
)
a
b
, with J
cd
given in (5.7). It has
its values in the Lie algebra of the Lorentz group. The spin connetion
a
b
is, for this reason, said to be a Lorentz connection. We can dene the matrix

=
i
2

a
b
J
a
b
whose entries are
a
b
. Then,

a
=

a
+
i
2

cd

(J
cd
)
a
b

b
= [

]
a
.
We nd also that

= h
a

a
b
h
b

+ h
b

h
b

(5.32)
is Lorentz invariant, that is,
h
a

+
a

_
h
b

= h
a

[
a
b

+
a
b
] h
b

.
Notice that the components
a
b
and

refer to dierent spaces. Equation

(5.32) describes how changes when the algebra indices are changed into
Riemann indices. It should not be confused with (5.25), which relates the
connection components in two Lorentzrelated frames.
Comment 5.1 Equation (5.32) is frequently written as the vanishing of a total covariant
derivative of the tetrad:

h
a

h
a

+
a
c
h
c

= 0 .
152
5.17 As the two rst indices in
a
b
are not tensorial, the behavior of
is very special. On the other hand, the indices in J
ab
are tensorial. The
consequence is that the contraction

=
i
2

ab

J
ab
is not Lorentz invariant.
Actuallly,

J
a

b
=
a

db

J
a

b
+
ab

J
ab
.
Again, decomposing in terms of the tetrad, we nd
_

+ h
b

h
a

_
J
a

b
=
_

ab

+ h
b

h
a

J
ab
.
Thus, what is really invariant is

=
1
2
_

ab

+ h
b

h
a

J
ab
. (5.33)
5.18 The archaic approach to connections is more intuitive and very
suggestive. It is worth recalling, as it complements the above one. It starts
with the assumption that, under an innitesimal displacement dx

, a eld
suers a change which is proportional to its own value and to dx

. The
proportionality coecient is an ane coecient

(an oldish name !),

so that

(x) =

(x)

(x) dx

.
If we introduce the entries of the matrix J
ab
, we verify that this is the same
as

(x) =
1
2

ab

(x) (J
ab
)

(x) dx

.
Consequently, the variation in the functional form of will be

(x) =

(x)

(x) dx

(x) dx

,
which denes the covariant derivative

(x) =

(x) +

(x) .
5.19 When

= 0, we say that the eld is parallel transported. In

that case, we have

(x) = 0, or

(x) =

(x). In parallel transport,

153
the functional form of the eld does not change. To learn more on the
meaning of parallel-transport and its expression as

= 0, let us look at
the functional variation of the vector eld along a curve of tangent vector
(velocity) U. It will be the 1-form

(x) applied to the eld U = d/ds:

(x)[U] = [

] dx

[U] = U

=
D

ds
.
The purely-functional variation along the curve will vanish if

= 0. This
is what is meant when we say that

is parallel-transported along a curve :

the eld is transported in such a way that it suers no change in its functional
form. The change coming from the argument,
d

ds
=
dx

ds

= U

,
is exactly compensated by the term U

.
When

= U

itself, we have the acceleration DU

/ds. The condition

of no variation of the velocity-eld along the curve will lead to the geodesic
equation
DU

ds
= U

] =
dU

ds
+

= 0 ,
which is an equation for .
5.2.4 Back to Equivalence
5.20 We have seen in Section 3.7 how a convenient choice of coordinates
leads to the vanishing of the Levi-Civita connection at a point. Nevertheless,
we had in 3.11 introduced an observer as a timelike curve, and qualied
that notion in the ensuing paragraphs. A timelike curve is, actually, an
ideal observer, which is point-like in any local space section. Real observers
are extended in space and can always detect gravitation by comparing the
neighboring curves followed by hers/his parts ( 3.46). Ideal observers can
be arbitrary curves (subsection 3.8.2) and, eventually, can have well-dened
space-sections all along (subsection 3.8.4).
We shall now see how a tetrad {H
a
} can be chosen so that, seen from the
frame it represents, the connection can be made to vanish all along a curve.
154
Looking from that frame, an ideal observer will not feel the gravitational
eld.
5.21 Take a dierentiable curve which is an integral curve of a eld U,
with U

=
dx

ds
=
d

(s)
ds
. The condition for the connection to vanish along ,

a
b
((s)) = 0, will be
U

H
a

((s)) +

((s)) U

H
a

((s)) = 0,
that is,
d
ds
H
a

((s)) +

((s)) U

H
a

((s)) = 0. (5.34)
This is simply the requirement that the tetrad (each member H
a
of it) be
parallel-transported along . Given a curve and a linear connecion, any
vector eld can be parallel-transported along . The procedure is then very
simple: take, in the way previously discussed, a point P on the curve and
nd the corresponding trivial tetrad {H
a
(P)}; then, parallel-transport it
along the curve.
For the dual base {H
a
}, the above formula reads
d
ds
H
a

((s))

((s)) U

H
a

((s)) = 0. (5.35)
5.22 Take (5.34) in the form
H
b

(x)
d
ds
H
b

(x) =

(x)U

and contract it with U

:
U

H
b

(x)
d
ds
H
b

(x) =

(x)U

. (5.36)
Seen from the tetrad, the tangent eld will be U
b
= H
b

, and the above

formula is
U
b
d
ds
H
b

(x) +

(x)U

= 0, (5.37)
or
H
b

(x)
d
ds
U
b
=
d
ds
U

(x) +

(x)U

=
D
Ds
U

(x). (5.38)
155
This shows that, if is a self-parallel curve, then
d
ds
U
a
= 0. (5.39)
This is the equation for a geodesic, as seen from the frame {H
a
= H
a

x

}.
If an external force is present, then m
D
Ds
U

(x) = F

and
m
d
ds
U
a
= F
a
. (5.40)
This means that, looking from that tetrad, the observer will see the laws of
Physics as given by Special Relativity.
5.23 Consider now a Levi-Civita geodesic. In that case there exists a
preferred tetrad h
a
, which is not parallel-transported along the curve. Its
deviation from parallelism is measured by the spin connection. Indeed, from
(5.32),
d
ds
h
b

h
b

= h
a

a
b
U

. (5.41)
This is the same as
d
ds
h
b

= h
b

h
d

b
dc

U
c
. (5.42)
We call

U
the velocity as seen from the frame {h
a
}. Equation (5.32) is
actually a representation of the Equivalence Principle. Let us write it into
still another form,
D
Ds
h
b
=
d
ds
h
b
+
_
d
ds
_
h
b
=
a
b
_
d
ds
_
h
a
. (5.43)
This holds for any curve with tangent vector (
d
ds
). It means that the frame
{h
a
} can be parallel-transported along no curve. The spin connection forbids
it, and gives the rate of change with respect to parallel transport.
5.24 One of the versions of the Principle the etymological one says
that a gravitational eld is equivalent to an accellerated frame. Which frame
? We see now the answer: the frame equivalent to the eld represented by
the metric g

=
ab
h
a

h
b

is just the anholonomous frame {h

a
}. Another
156
piece of the Principle says that it is possible to choose a frame in which the
connection vanishes. Let us see now how to change from the equivalent
frame {h
a
} to the free-falling frame {H
a
}. Contracting (5.41) with H
a

,
H
a

d
ds
h
b

H
a

h
b

= H
a

h
c

c
b
U

,
which is the same as
d
ds
_
H
a

h
b

_
h
b

_
d
ds
H
a

H
a

_
= H
a

h
c

c
b
U

.
The second term in the left-hand side vanishes by Eq.(5.35), so that we
remain with
d
ds
_
H
a

h
b

_
H
a

h
c

c
b
U

= 0 . (5.44)
What appears here is a point-dependent Lorentz transformation relating the
metric tetrad h
a
to the frame H
a
in which vanishes:
H
a

=
a
b
h
b

, H
a
=
a
b
h
b
. (5.45)
This means

a
b
= H
a

h
b

. (5.46)
This relation holds on the common domain of denition of both tetrad elds.
Taking this into (5.44), we arrive at a relationship which holds on the inter-
section of that domain with a geodesics of the Levi-Civita connection

:
d
ds

a
b

a
c

c
bd

U
d
= 0 . (5.47)
This equation gives the change, along the metric geodesic, of the Lorentz
transformation taking the metric tetrad h
a
into the frame H
a
. The vector
formed by each row of the Lorentz matrix is parallel-transported along the
line. Contracting with the inverse Lorentz transformation
(
1
)
a
b
= h
a

H
b

, (5.48)
the expression above gives

a
bd

U
d
= (
1
)
a
c
d
ds

c
b
= (
1
d
ds
)
a
b
. (5.49)
157
This is, in the language of dierential forms,

a
bd
h
d
_
d
ds
_
= (
1
d)
a
b
_
d
ds
_
. (5.50)
Thus on the points of the curve, the connection has the form of a gauge
Lorentz vacuum:

a
b
= (
1
d)
a
b
(5.51)
It is important to stress that this is only true along a curve a one-
dimensional domain so that curvature is not aected. Curvature, the real
gravitational eld, only manisfests itself on two-dimensional domains. Seen
from the frame h
a
, the geodesic equation has the form
d
ds

U
a
+

a
bc

U
b

U
c
= 0, (5.52)
which is the same as
d
ds

U
a
+ (
1
d
ds
)
a
b

U
b
= 0.
This expression, once multiplied on the left by , gives

c
a
d
ds

U
a
+
d
ds
(
c
a
)

U
a
=
d
ds
(
c
a

U
a
) = 0. (5.53)
This is Eq. (5.39) for the present case. Summing up: at each point of the
curve/observer there is a Lorentz transformation taking the accelerated frame
{h
a
}, equivalent to the gravitational eld, into the inertial frame {H
a
}, in
which the force equation acquires the form it would have in Special Relativity.
These considerations can be enlarged to general Lorentz tensors. Take,
for instance, a second order tensor:
T
ab
=
a
c

b
d

T
cd
.
Taking
d
ds
of this expression leads to
d
ds

T
ab
+

a
c

T
cb
+

b
c

T
ac
= (
1
)
a
c
(
1
)
b
d
d
ds
T
cd
.
The covariant derivative according to the connection

, which is seen from

the frame h
a
, is the Lorentz transform of the simple derivative as seen from
the inertial frame H
a
, in which vanishes.
158
Summing up again: at each point of the curve/observer there is a Lorentz
transformation taking the accelerated frame {h
a
}, equivalent to the gravita-
tional eld, into the free falling frame {H
a
}, in which all tensorial (that is,
covariant) equations acquire the form they would have in Special Relativity.
This is the content of the Equivalence Principle.
5.25 In the gauge theories describing the other fundamental interactions,
an analogous property turns up: at each point of the curve/observer a gauge
can be found in which the potential vanishes. However, it is not the potential
but the eld strength which appears in the force equation the Lorentz
force equation (1.18) is typical. The eld strength is the curvature in these
theories, and cannot be made to vanish. This is a crucial dierence between
gravitation and the other interactions.

5.2.5 Two Gates into Gravitation

5.26 Starting with a given nontrivial tetrad eld {h
a
}, two dierent but
equivalent ways of describing gravitation are possible. In the rst, according
to (5.10), the nontrivial tetrad eld is used to dene a Riemannian metric g

,
from which we can contruct the LeviCivita connection and the correspond-
ing curvature tensor. As the starting point was a nontrivial tetrad eld, we
can say that such a tetrad is able to induce a metric structure in spacetime,
which is the structure underlying the General Relativity description of the
gravitational eld. We have seen that {h
a
} is not parallel-transported by the
Levi-Civita connection.
On the other hand, a nontrivial tetrad eld can be used to dene a very
special linear connection, called Weitzenb ock connection, with respect to
which the tetrad {h
a
} is parallel. For this reason, this kind of structure has
received the name of teleparallelism, or absolute parallelism, and is the stage-
set of the so called teleparallel description of gravitation. The important
point to be kept from these considerations is that a nontrivial tetrad eld is
able to induce in spacetime both a teleparallel and a Riemannian structure.
In what follows we will explore these structures in more detail.

Details can be found in R. Aldrovandi, P. B. Barros e J. G. Pereira, The equivalence

principle revisited, Foundations of Physics 33 (2003) 545-575 - ArXiv: gr-qc/0212034.
159
5.27 Let us consider now the covariant derivative of the metric tensor g

.
It is

,
or, by using (5.24)

= h
a

h
b

(
ab
+
ba
) . (5.54)
Therefore, the metricity condition

= 0 (5.55)
will only hold when the connection is either purely antisymmetric [(pseudo)
orthogonal],

ab
=
ba
,
or when it vanishes identically:

ab
= 0 .
The LeviCivita connection falls into the rst case and, as we are going
to see, the Weitzenbock connection into the second case. This means that
both connections preserve the metric.
160
Chapter 6
Gravitational Interaction of the
Fundamental Fields
6.1 Minimal Coupling Prescription
The interaction of a general eld with gravitation can be obtained through
the application of the so called minimal coupling prescription, according to
which the Minkowski metric must be replaced by the riemannian metric

ab
g

=
ab
h
a

h
b

, (6.1)
and all ordinary derivatives must be replaced by Fock-Ivanenko covariant
derivatives [1],

i
2

ab

S
ab
(6.2)
where
ab

is a connection assuming values in the Lie algebra of the Lorentz

group, usually called spin connection, and S
ab
is a Lorentz generator writ-
ten in a representation appropriate to the eld . This double coupling
prescription is a characteristic property of the gravitational interaction as
for all other interactions of Nature only the derivative replacement (6.2) is
necessary.
Let us explore deeper this point. The metric replacement (6.1) is a con-
sequence of the local invariance of the lagrangian under translations of the
tangentspace coordinates [2]. This part of the prescription, therefore, is
related to the coupling of the eld energymomentum to gravitation, and is
universal in the sense that it is the same for all elds. On the other hand,
161
the derivative replacement (6.2) is a consequence of the invariance of the
lagrangian under local Lorentz transformations of the tangentspace coordi-
nates. This part of the prescription, therefore, is related to the coupling of
the eld spin to gravitation, and is not universal because it depends on the
spin contents of the eld.
Another important point of the above coupling prescription is that the
metric change (6.1) is appropriate only for integerspin elds, whose la-
grangian is quadratic in the eld derivative, and consequently a metric tensor
is always present to contract these two derivatives. For half-integer spin elds,
however, the lagrangian is linear in the eld derivative, and consequently no
metric will be present to be changed. In this case, the metric change must
be replaced by the equivalent rule in terms of the tetrad eld,
e
a

h
a

, (6.3)
where e
a

is the trivial tetrad

e
a

=
x
a
x

, (6.4)
and h
a

is a nontrivial tetrad representing a true gravitational eld.

Besides being more fundamental than the metric change (6.1), the tetrad
change (6.3) allows the introduction of a full coupling prescription which en-
compasses both the metric or equivalently the tetrad and the derivative
changes. It is given by

a
D
a
h
a

= h
a

i
2

ab

S
ab
_
, (6.5)
This coupling prescription is general in the sense that it holds for both integer
and halfinteger spin elds. For the case of integer spin elds, it yields
automatically the metric replacement (6.1). For the case of halfinteger spin
elds, it yields automatically the tetrad replacement (6.3).
6.2 General Relativity Spin Connection
As is well known, a tetrad eld can be used to transform Lorentz into space-
time indices, and viceversa. For example, a Lorentz vector eld V
a
is related
to the corresponding spacetime vector V

through
V
a
= h
a

. (6.6)
162
It is important to notice that this applies to tensors only. Connections, for
example, acquire an extra vacuum term under such change [3]:

a
b
= h
a

h
b

+ h
a

h
b

. (6.7)
On the other hand, because they are used in the construction of covariant
derivatives, connections (or potentials, in physical parlance) are the most
important personages in the description of an interaction. Concerning the
specic case of the general relativity description of gravitation, the spin con-
nection, denoted here by
a
b
=

A
a
b
, is given by [4]

A
a
b
= h
a

h
b

+ h
a

h
b

h
a

h
b

. (6.8)
We see in this way that the spin connection

A
a
b
is nothing but the Levi
Civita connection

=
1
2
g

] (6.9)
rewritten in the tetrad basis. Therefore, the full coupling prescription of
general relativity is

Da
h
a

D
(6.10)
with

D
=

i
2

A
ab

S
ab
(6.11)
the general relativity FockIvanenko [1] covariant derivative operator.
Now, comes an important point. The covariant derivative (6.11) applied
to a general Lorentz tensor eld reduces to the usual Levi-Civita covariant
derivative of the corresponding spacetime tensor. For example, take again a
vector eld V
a
for which the appropriate Lorentz generator is [5]
(S
ab
)
c
d
= i (
c
a

bd

c
b

ad
) . (6.12)
It is then an easy task to verify that [6]

D
V
a
= h
a

. (6.13)
On the other hand, no LeviCivita covariant derivative can be dened for
half-integer spin elds [7]. For these elds, the only possible form of the
163
covariant derivative is that given in terms of the spin connection. For a
Dirac spinor , for example, the covariant derivative is

D
=

i
2

A
ab

S
ab
, (6.14)
where
S
ab
=
1
2

ab
=
i
4
[
a
,
b
] (6.15)
is the Lorentz spin-1/2 generator, with
a
the Dirac matrices. Therefore, we
may say that the covariant derivative (6.11), which take into account the spin
contents of the elds as dened in the tangent space, is more fundamental
than the LeviCivita covariant derivative in the sense that it is able to de-
scribe the gravitational coupling of both tensor and spinor elds. For tensor
elds it reduces to the LeviCivita covariant derivative, but for spinor elds
it remains as a FockIvanenko derivative.
6.3 Application to the Fundamental Fields
6.3.1 Scalar Field
Let us consider rst a scalar eld in a Minkowski spacetime, whose la-
grangian is
L

=
1
2
_

ab

a

b

2

, (6.16)
with
=
mc

. (6.17)
The corresponding eld equation is the so called KleinGordon equation:

a
+
2
= 0 . (6.18)
In order to get the coupling of the scalar eld with gravitation, we use
the full coupling prescription

Da
h
a

D
= h
a

i
2

A
ab

S
ab
_
. (6.19)
164
For a scalar eld, however,
S
ab
= 0 , (6.20)
and the coupling prescription in this case becomes

a
h
a

. (6.21)
Applying this prescription to the lagrangian (6.16), we get
L

g
2
_
g

. (6.22)
Then, by using the identity

g =

g
2
g

, (6.23)
it is easy to see that the corresponding eld equation is

+
2
= 0 , (6.24)
where

_
g g

_
(6.25)
is the LaplaceBeltrami operator, with

the LeviCivita covariant deriva-

tive. We notice in passing that it is completely equivalent to apply the
minimal coupling prescription to the lagrangian or to the eld equations.
Furthermore, we notice that, in a locally inertial coordinate system, the rst
derivative of the metric tensor vanishes, the LeviCivita connection vanishes
as well, and the LaplaceBeltrami becomes the freeeld dAlambertian op-
erator. This is the usual version of the (weak) equivalence principle.
6.3.2 Dirac Spinor Field
In Minkowski spacetime, the spinor eld lagrangian is
L

=
ic
2
_

a

a

_
mc
2

. (6.26)
The corresponding eld equation is the Dirac equation
i
a

a
mc = 0 . (6.27)
165
In the context of general relativity, the coupling of a Dirac spinor with
gravitation is obtained through the application of the full coupling prescrip-
tion

Da
h
a

D
= h
a

i
2

A
ab

S
ab
_
, (6.28)
where now S
ab
stands for the spin-1/2 generators of the Lorentz group, given
by
S
ab
=

ab
2
=
i
4
[
a
,
b
] . (6.29)
The spin connection, according to Eq.(6.8), is written in terms of the tetrad
eld as

A
a
b
= h
a

h
b

+ h
a

h
b

. (6.30)
Applying the above coupling prescription to the free lagrangian (6.26), we
get
L

g c
_
ic
2
_

_
mc
2

_
, (6.31)
where

= e
a

a
is the local Dirac matrix, which satisfy
{

} = 2
ab
h
a

h
b

= 2g

. (6.32)
The corresponding Dirac equation can be obtained through the use of the
Euler-Lagrange equation
L

D
L

D

)
= 0 . (6.33)
The result is the Dirac equation in a Riemann spacetime
i

D
mc = 0 . (6.34)
6.3.3 Electromagnetic Field
In Minkowski spacetime, the electromagnetic eld is described by the la-
grangian density
L
em
=
1
4
F
ab
F
ab
, (6.35)
166
where
F
ab
=
a
A
b

b
A
a
(6.36)
is the Maxwell eld strength. The corresponding eld equation is

a
F
ab
= 0 , (6.37)
which along with the Bianchi identity

a
F
bc
+
c
F
ab
+
b
F
ca
= 0 , (6.38)
constitute Maxwells equations. In the Lorentz gauge
a
A
a
= 0, the eld
equation (6.37) acquires the form

c
A
a
= 0 . (6.39)
In the framework of general relativity, the form of Maxwells equations can
be obtained through the application of the full minimal coupling prescription
(6.5), which amounts to replace

a
h
a

D
= h
a

i
2

A
ab

S
ab
_
, (6.40)
with
(S
ab
)
c
d
= i (
a
c

bd

b
c

ad
) (6.41)
the vector representation of the Lorentz generators. For the specic case of
the electromagnetic vector eld A
a
, the Fock-Ivanenko derivative acquires
the form

D
A
a
=

A
a
+

A
a
b
A
b
. (6.42)
It is important to remark once more that the FockIvanenko derivative
is concerned only to the local Lorentz indices. In other words, it ignores the
spacetime tensor character of the elds. For example, the FockIvanenko
derivative of the tetrad eld is

D
h
a

h
a

A
a
b
h
b

. (6.43)
Substituting

A
a
b
= h
a

h
b

, (6.44)
167
we get

D
h
a

h
a

. (6.45)
As a consequence, the total covariant derivative of the tetrad h
a

, that is, a
covariant derivative which takes into account both indices of h
a

, vanishes
identically:

h
a

A
a
b
h
b

h
a

= 0 . (6.46)
Now, any Lorentz vector eld A
a
can be transformed into a spacetime
vector eld A

through
A

= h
a

A
a
, (6.47)
where A

transforms as a vector under a general spacetime coordinate trans-

formation. Substituting into equation (6.42), and making use of (6.45), we
get

D
A
a
= h
a

. (6.48)
We see in this way that the FockIvanenko derivative of a Lorentz vector eld
A
c
reduces to the usual LeviCivita covariant derivative of general relativity.
This means that, for a vector eld, the minimal coupling prescription (6.40)
can be restated as

a
A
c
h
a

h
c

. (6.49)
Therefore, in the presence of gravitation, the electromagnetic eld lagrangian
acquires the form
L
em
=
1
4

g F

, (6.50)
where
F

, (6.51)
the connection terms canceling due to the symmetry of the LeviCivita con-
nection in the last two indices. The corresponding eld equation is

= 0 , (6.52)
168
or equivalently, assuming the covariant Lorentz gauge

= 0,

= 0 . (6.53)
Analogously, the Bianchi identity (6.38) can be shown to assume the form

= 0 . (6.54)
We notice in passing that the presence of gravitation does not spoil the U(1)
gauge invariance of Maxwell theory. Furthermore, like in the case of the
scalar eld, it results the same to apply the coupling prescription in the
lagrangian or in the eld equations.
169
Chapter 7
General Relativity with Matter
Fields
7.1 Global Noether Theorem
Let us start by briey reviewing the results of the global or rst
Noethers theorem [8]. As is well known, the global Noether theorem is
concerned with the invariance of the action functional under global trans-
formations. For each of such invariances, Noethers theorem determines a
conservation law. In the specic case of the invariance under a global trans-
lation of the spacetime coordinates, the corresponding Noether conserved
current is the canonical energymomentum tensor

a
b
=
L

a
b
L

, (7.1)
with L

the lagrangian of the eld .

On the other hand, in the case of the invariance under a global rota-
tion of the spacetime coordinates that is, under a Lorentz transformation
the corresponding Noether conserved current is the canonical angular
momentum tensor
J
a
bc
= M
a
bc
+S
a
bc
, (7.2)
where
M
a
bc
= x
b

a
c
x
c

a
b
, (7.3)
170
is the orbital angularmomentum, and
S
a
bc
= i
L

S
bc
(7.4)
is the spin angularmomentum, with S
bc
the generators of Lorentz transfor-
mations written in the representation appropriate for the eld . Notice that
M
a
bc
is the same for all elds, whereas S
a
bc
depends on the spin contents of
the eld .
Notice that the canonical energymomentum tensor
a
b
is not symmetric
in general. However, using the Belinfante procedure [9] it is possible to dene
a symmetric energymomentum tensor for the spinor eld,

ab
=
ab

1
2

cab
, (7.5)
where

cab
=
acb
= S
cab
+S
abc
S
bca
. (7.6)
It can be easily veried that

cab
=
ab

ba
, (7.7)
which together with (7.5) show that
ab
is in fact symmetric.
7.2 EnergyMomentum as Source of Curva-
ture
An old and controversial problem of gravitation is the conservation of energy
momentum density for both gravitational and matter elds. Concerning the
energymomentum tensor of matter elds, it becomes problematic mainly
when spinor elds are present [10]. In order to explore deeper these prob-
lems, we are going to study the denition as well as the conservation law
of the gravitational energymomentum density of a general matter eld. By
gravitational energymomentum tensor we mean the source of gravitation,
that is, the tensor appearing in the right handside of the gravitational eld
equations. For the specic case of a spinor eld, this energymomentum
tensor is sometimes believed to acquire a genuine nonsymmetric part. As
the left handside of the gravitational eld equations are always symmet-
ric, this would call for a generalization of general relativity. However, the
171
gravitational energymomentum tensor is actually always symmetric, even
for a spinor eld, which shows the consistency and completeness of general
relativity.
Let us consider the lagrangian
L = L
G
+L

, (7.8)
where
L
G
=
c
4
16G

R
(7.9)
is the EinsteinHilbert lagrangian of general relativity, and L

is the la-
grangian of a general matter eld . The functional variation of L in relation
to the metric tensor g

yields the eld equation

R

1
2
g

R
=
4G
c
4
T

, (7.10)
where
T

=
2

g
L

(7.11)
is the gravitational energymomentum tensor of the eld . The contravari-
ant components of the gravitational energymomentum tensor is
T

=
2

g
L

. (7.12)
In these expressions,
L

=
L

(7.13)
is the Lagrange functional derivative. In general relativity, therefore, energy
and momentum are the source of gravitation, or equivalently, are the source
of curvature. As the metric tensor is symmetric, the gravitational energy
momentum tensor obtained from either expression (7.11) or (7.12) is always
symmetric. These expressions yield the energymomentum tensor not only
in the case of the presence of a gravitational eld, but also in the absence.
In the absence of a gravitational eld, a transition to curvilinear coordinates
must be done before the calculation of T

. Of course, the metric tensor

in this case will not represent a true gravitational eld, but only eects of
coordinates.
172
7.3 EnergyMomentum Conservation
Let us obtain now the conservation law of the gravitational energymomen-
tum tensor of a general source eld . Denoting by L

the lagrangian of the

eld , the corresponding action integral is written in the form
S =
1
c
_
L

d
4
x . (7.14)
As a spacetime scalar, it does not change under a general transformation of
coordinates. Of course, under a coordinate transformation, the eld will
change by an amount . Due to the equation of motion satised by this
eld, the coecient of vanishes, and for this reason we are not going to
take these variations into account. For our purposes, it will be enough to
consider only the variations in the metric tensor g

. Accordingly, by using
Gauss theorem, and by considering that g

= 0 at the integration limits,

the variation of the action integral (7.14) can be written in the form [11]
S =
1
c
_
L

d
4
x =
1
c
_
L

d
4
x , (7.15)
with L

/g

the Lagrange functional derivative (7.13). But, we have al-

ready seen in section 7.2 that
L

g
2
T

, (7.16)
where T

is the gravitational energymomentum tensor of the eld . There-

fore, we have
S =
1
2c
_
T

g d
4
x =
1
2c
_
T

g d
4
x . (7.17)
On the other hand, under a spacetime general coordinate transformation
x

= x

, (7.18)
with

small quantities, the components of the metric tensor change accord-

ing to
g

(x

) g

(x

) = g

, (7.19)
where only terms linear in the transformation parameter

have been kept.

Substituting into (7.17), we get
S =
1
2c
_
T

_
g

g d
4
x . (7.20)
173
Integrating by parts the second and third terms, neglecting integrals over
hypersurfaces, and making use of the symmetry of T

, we get
S =
1
c
_ _

gT

)
1
2

g T

d
4
x , (7.21)
or equivalently
S =
1
c
_
_

gT

g
_

d
4
x . (7.22)
Then, by using the identity

g =

, (7.23)
we get
S =
1
c
_

g d
4
x . (7.24)
Therefore, from both the invariance condition S = 0 and the arbitrariness
of

, it follows that

= 0 . (7.25)
It is important to remark that this is not a true conservation law in the
sense that it does not lead to a charge conserved in time. Instead, it is
an identity satised by the gravitational energymomentum tensor, usually
called Noether identity [8]. The sum of the energymomentum of the gravi-
tational eld t

and the energymomentum of the matter eld T

is a truly
conserved quantity. In fact, this quatity satises

_
g (t

+T

= 0 , (7.26)
which, by using Gauss theorem, yields the true conservation law
dq

dt
= 0 , (7.27)
with
q

=
_
_
t
0

+T
0

_
g d
3
x (7.28)
the conserved charge.
174
7.4 Examples
7.4.1 Scalar Field
Let us take the lagrangian of a scalar eld in a Minkowski spacetime,
L

=
1
2
_

ab

a

b

2

, (7.29)
with given by (6.17). From Noethers theorem we nd that the correspond-
ing canonical energymomentum and spin tensors are given respectively by

a
b
=
a

b

a
b
L

, (7.30)
and
S
a
bc
= 0 . (7.31)
As a consequence of the vanishing of the spin tensor, the canonical energy
momentum tensor of the scalar eld is symmetric, and is conserved in the
ordinary sense:

a
b
= 0 . (7.32)
In the presence of gravitation, the scalar eld lagrangian is given by
L

g
2
_
g

. (7.33)
By using the identity

g =

g
2
g

,
the dynamical energymomentum tensor is found to be

g T

. (7.34)
Like in the free case, it is symmetric, and conserved in the covariant sense:

= 0 . (7.35)
175
7.4.2 Dirac Spinor Field
The Dirac spinor lagrangian in Minkowski spacetime is
L

=
ic
2
_

a

a

_
mc
2

. (7.36)
From the rst Noethers theorem one nds that the corresponding canonical
energymomentum and spin tensors are given respectively by

a
b
=
ic
2
_

b

b

_
, (7.37)
and
S
a
bc
=
c
2
_

a
S
bc
+

S
bc

_
, (7.38)
with
S
bc
=

bc
2
=
i
4
[
b
,
c
] . (7.39)
It should be noticed that, in contrast to the scalar eld case, the canonical
energymomentum tensor for the Dirac spinor is not symmetric. As already
discussed, however, we can use the Belinfante procedure [9] to construct a
symmetric energymomentum tensor for the spinor eld, which is given by

ab
=
ab

1
2

cab
, (7.40)
with

cab
=
acb
= S
cab
+S
abc
S
bca
. (7.41)
In the presence of gravitation, the spinor eld lagrangian is
L

g c
_
i
2
_

_
mc

_
. (7.42)
The dynamical energymomentum tensor of the spinor eld, according to the
denition (7.11), is found to be
T

1
2

, (7.43)
where

=
ic
2
_

_
(7.44)
176
is the canonical energymomentum tensor modied by the presence of grav-
itation, and

is still given by (7.6), but now written in terms of the spin

tensor modied by the presence of gravitation:
S

bc
=
c
4
_

a
h
a

bc
+

bc

a
h
a

_
. (7.45)
Equation (7.43) is a generalization of the Belinfante procedure for the
presence of gravitation. In fact, through a tedious but straightforward cal-
culation we can show that
D

= g

, (7.46)
from which we see that the dynamical energymomentum tensor T

of the
Dirac eld, that is, the EulerLagrange functional derivative of the spinor
lagrangian (7.42) with respect to the metric, is always symmetric.
7.4.3 Electromagnetic Field
In Minkowski spacetime, the electromagnetic eld is described by the la-
grangian density
L
em
=
1
4
F
ab
F
ab
. (7.47)
The corresponding canonical energymomentum and spin tensors are given
respectively by

a
b
= 4
b
A
c
F
ac
+
a
b
F
cd
F
cd
, (7.48)
and
S
a
bc
= F
a
b
A
c
F
a
c
A
b
. (7.49)
As in the spinor case, the canonical energymomentum tensor
a
b
of the
electromagnetic eld is not symmetric. By using the Belinfante procedure,
however, it is possible to dene the symmetric energymomentum tensor,

ab
=
ab

1
2

cab
, (7.50)
where

cab
=
acb
= S
cab
+S
abc
S
bca
. (7.51)
177
As can be easily veried,

cab
= 2F
ca
A
b
. (7.52)
Consequently,

ab
= 4
_
F
ac
F
b
c
+
1
4

ab
F
cd
F
cd
_
(7.53)
is in fact symmetric, and conserved in the ordinary sense:

a
b
= 0 . (7.54)
In the presence of gravitation, the electromagnetic eld lagrangian is
L
em
=
1
4

g F

, (7.55)
and the dynamical energymomentum tensor
T

=
2

g
L
em
g

(7.56)
is found to be
T

=
_
F

1
4
g

_
. (7.57)
It is symmetric and covariantly conserved:

= 0 . (7.58)
178
Chapter 8
Closing Remarks
Gravitation diers from the other three known fundamental interactions of
Nature by its more intimate relationship to spacetime. The other inter-
actions (electromagnetic, weak and strong) are also described by (gauge)
theories with a large geometrical content. However, while gravitation deals
with changes of frames, the other interactions are concerned with changes of
gauges in internal spaces. Gravitation relates to energy, while the other
interactions cope with conserved quantities (charges) which are indepen-
dent of the events on spacetime electric charge, weak isotopic spin and
hypercharge, color.
In consequence, gravitation engender forces of inertial type, quite distinct
from charge-produced forces. Hence its unique, universal character. Its pres-
ence is felt by all particles and elds in the same way as if changing the
very scene in which phenomena take place.
We hope to have given in these notes a rst glimpse into the way these
strange things happen.
179
Bibliography
[1] V. A. Fock, Z. Phys. 57, 261 (1929).
[2] V. C. de Andrade and J. G. Pereira, Phys. Rev D 56, 4689 (1997).
[3] R. Aldrovandi and J. G. Pereira, An Introduction to Geometrical
Physics (World Scientic, Singapore, 1995).
[4] P. A. M. Dirac, in: Planck Festscrift, ed. W. Frank (Deutscher Verlag
der Wissenschaften, Berlin, 1958).
[5] P. Ramond, Field Theory: A Modern Primer, 2nd edition (Addison-
Wesley, Redwood, 1989).
[6] V. C. de Andrade and J. G. Pereira, Int. J. Mod. Phys. D 8, 141 (1999).
[7] M. J. G. Veltman, Quantum Theory of Gravitation, in Methods in Field
Theory, Les Houches 1975, Ed. by R. Balian and J. Zinn-Justin (North-
Holland, Amsterdam, 1976).
[8] See, for example: N. P. Konopleva and V. N. Popov, Gauge Fields
(Harwood, New York, 1980).
[9] F. J. Belinfante, Physica 6, 687 (1939).
[10] K. Hayashi, Lett. Nuovo Cimento 5, 529 (1972).
[11] L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Perg-
amon, Oxford, 1975).
180