In this post I’ll be deriving one of the simplest spacetime metrics: the Schwarzschild Metric. It is named after Karl Schwarzschild who published the solution in 1916, just a little more than a month after the publication of Einstein’s Theory of General Relativity.
It is a metric for a mass that’s:
Spherically symmetric
Non-rotating
Non-charged
and it describes the metric for:
Slowly rotating planets/stars
Non-rotating Black Holes
and it predicts:
Gravitational time dilation around celestial bodies
Gravitational doppler effect
Bending of light due to gravity
Shift in perihelion of orbits
Existence of Black Holes with the event horizon the size of Schwarzschild radius rs
So let’s get into it!
Einstein Field Equations
To find the metric, we first have to solve the Einstein Field Equations.
Rμν−21Rgμν−Λgμν=c48πGTμν
We input the energy-momentum tensor and we get out the metric tensor. That is: we give the description of the mass-energy and momentum and we get the description of the geometry of spacetime.
Using the metric, we can find out the connection coefficients:
Γμνσ=21gσα(∂νgαμ+∂μgαν−∂αgμν)
whence we can derive the geodesic equation which gives us the paths of masses and light:
dλ2d2xσ+Γμνσdλdxμ⋅dλdxν=0
Tμν=0 inside the body, so let’s find the metric outside of our body where Tμν=0. We can also assume that it’s at non-cosmological scales, i.e. Λ=0. So that gives us:
for vacuum outside the body. This property is called Ricci flat: there’s no immediate volume changes in test particles, there’s only squishing due to curvature of spacetime.
We’re going to assume that as we go farther away from the mass, the effects of gravity are going to be negligible (i.e. described by the Minkowski metric).
We’re also going to assume that the spacetime is static, which means that:
The metric doesn’t depend on time: ∂tgμν=0
t→−t doesn’t change the metric (this guarantees that the black hole is non-rotating)
Since we want spherical symmetry, the θ and ϕ components should use the metric for a sphere of radius r. Hence, to not violate spherical symmetry, the metric becomes:
ds2=−A(r~)c2dt2+B(r~)dr~2+r~2(dθ2+sin2θdϕ2)
where C(r~) is a function that scales the values and is independent of θ and ϕ.
To make things look cleaner, let’s define:
C(r~)r~≡r
To find A(r) and B(r), we first have to find the metric tensor gμν, which we can use to find the connection coefficients Γμνσ, which we can use to find the Ricci tensor Rμν. Then we can force the metric to match Newtonian gravity at weak field and low velocity, hence solve for A(r) and B(r).
We can calculate the connection coefficients using:
Γμνσ=21gσα(∂νgαμ+∂μgαν−∂αgμν)
Hence we can find the non-zero connection coefficients to be:
Since Rμν=0, that implies R00=R11=R22=0. We can use this to solve for A and B.
The Riemann tensor is:
Rσμνρ=∂μΓνσρ−∂νΓμσρ+ΓνσαΓμαρ−ΓμσβΓνβρ
Hence, the Ricci tensor is obtained by contraction:
Rνδ=Rνμδμ
Calculating R00
R00=R0μ0μ==0 only for Γ001∂μΓ00μ−∂0=0Γμ0μ+=0 only for α=1Γ00αΓμαμ−(Γμ0βΓ0βμ)=∂1Γ001−0+Γ001Γμ1μ−(Γμ00Γ00μ+Γμ01Γ01μ)=∂1Γ001+Γ001Γ010+Γ001Γ111+Γ212=Γ313Γ001Γ212+Γ001Γ313−Γ100Γ001−Γ001Γ010=∂1Γ001+Γ001Γ111+2Γ001Γ212−Γ100Γ001
Solving Rμν=0
From above, we are working with the static, spherically symmetric ansatz
ds2=−A(r)c2dt2+B(r)dr2+r2(dθ2+sin2θdϕ2)
where A(r) and B(r) are functions of r only, and we want the vacuum equations
Rμν=0
It is standard to write derivatives with respect to r using a prime:
A′=drdA,B′=drdB,A′′=dr2d2A
The Ricci components
Using the connection coefficients (Christoffels) computed earlier, the nontrivial Ricci components reduce to three independent equations. A convenient set is R00=0, R11=0, and R22=0.
A standard simplification is to work with the mixed structure of these equations and write them in terms of A,B and their derivatives. The results are:
The R00=0 equation
R00=0
is equivalent to
AA′′−2AA′(AA′+BB′)+rA2A′=0
The R11=0 equation
R11=0
is equivalent to
−AA′′+2AA′(AA′+BB′)+rB2B′=0
The R22=0 equation
R22=0
is equivalent to
1−B1+2Br(−AA′+BB′)=0
(And R33=0 gives the same condition as R22=0, up to the sin2θ factor.)
First key simplification: A(r)B(r)=constant
Add the R00=0 and R11=0 equations. The A′′ terms cancel, and the messy middle term cancels too, leaving a very simple relation:
rA2A′+rB2B′=0
Multiply both sides by r/2:
AA′+BB′=0
This means
drdln(AB)=0
so
A(r)B(r)=C
for some constant C.
Now we impose asymptotic flatness: far away from the mass, we want Minkowski space, so
A(r)→1,B(r)→1as r→∞
Therefore
C=1
and we get the important relation
B(r)=A(r)1
So the metric becomes a 1-function ansatz:
ds2=−A(r)c2dt2+A(r)1dr2+r2(dθ2+sin2θdϕ2)
Solve for A(r) using R22=0
Take the equation
1−B1+2Br(−AA′+BB′)=0
and substitute B=1/A.
First note:
B1=A
and also
BB′=drdlnB=drdln(A1)=−AA′
So the bracket becomes
−AA′+BB′=−AA′−AA′=−A2A′
Plugging into R22=0 gives
1−A+2Br(−A2A′)=0
But B=1/A, so B1=A and therefore 2B1=2A. Hence
1−A+r(2A)(−A2A′)=0
The A cancels nicely:
1−A−rA′=0
Rewrite it as
rA′+A=1
This is a first-order ODE:
drd(rA)=1
Integrate:
rA=r+K
so
A(r)=1+rK
Asymptotic flatness already ensured the 1 is correct. The remaining constant K must be related to the mass.
Match to Newtonian gravity (weak-field limit)
In the weak field, the metric should reduce to Newtonian gravity. In particular, the time-time component behaves like
gtt≈−(1+2Φ/c2)c2
where Φ is the Newtonian potential
Φ(r)=−rGM
In our metric,
gtt=−A(r)c2
So for large r,
A(r)≈1+c22Φ=1−c2r2GM
Comparing with
A(r)=1+rK
we identify
K=−c22GM
Therefore
A(r)=1−c2r2GM
and since B=1/A,
B(r)=(1−c2r2GM)−1
The Schwarzschild metric
Substituting A(r) and B(r) back into the line element, we obtain the Schwarzschild metric:
The coordinate singularity at r=rs corresponds to the event horizon of a non-rotating, uncharged black hole.
Beyond Schwarzschild: Other Black Hole Solutions
The Schwarzschild solution describes a very special case: a black hole that is non-rotating and uncharged.
In reality, more general black holes can carry electric charge, angular momentum, or both.
Remarkably, General Relativity admits exact analytic solutions for all these cases.
We will now derive (or at least motivate) the remaining classical black hole metrics.
Charged Black Holes: The Reissner–Nordström Metric
Physical setup
We now allow the black hole to carry an electric charge Q, while still assuming:
spherical symmetry
no rotation
vacuum outside the body (except for the electromagnetic field)
This means the stress–energy tensor is not zero, but instead given by the electromagnetic field.
Einstein–Maxwell equations
The Einstein equations become:
Rμν−21Rgμν=c48πGTμν(EM)
where the electromagnetic stress–energy tensor is
Tμν(EM)=μ01(FμαFνα−41gμνFαβFαβ)
with Fμν the electromagnetic field tensor.
For a static, spherically symmetric charge distribution, the only nonzero component is the radial electric field:
E(r)=4πε0r2Q
Metric ansatz
Spherical symmetry again forces the metric to take the form
ds2=−A(r)c2dt2+A(r)1dr2+r2(dθ2+sin2θdϕ2)
Solving the coupled Einstein–Maxwell equations yields