Visualizing curvature in General Relativity

General Relativity is a theory of gravity that is formulated in a 4-dimensional spacetime. So how do we visualize spacetimes, if they're 4D? This article will hopefully explain the mathematics of how one type of visualization, embedding diagrams, work.

Take the Schwarzschild metric:

ds2=(1Rsr)c2dt2+1(1Rsr)dr2+r2(dθ2+sin2dϕ2) ds^2 = -\left(1 - \frac{R_s}{r}\right) c^2 dt^2 + \frac{1}{\left(1 - \frac{R_s}{r}\right)} dr^2 + r^2 (d\theta^2 + \sin^2 d\phi^2)

Where:

Rs=2GMc2 R_s = \frac{2GM}{c^2}

Note that we have a problem: spacetime is 4D, but we can only visualize curvature of a function of two variables embedded in 3D space. So we need to reduce 4D spacetime down 2 dimensions to visualize it effectively. Typically, to do so, we take a 2D slice of 4D spacetime by setting t=0t = 0 and ϕ=0\phi = 0 (the ϕ\phi coordinate because the metric is spherically symmetric, so it'll look the same whichever ϕ\phi we choose, which is why we can take any arbitrary slice of it). So the metric becomes:

ds2=1(1Rsr)dr2+r2dθ2 ds^2 = \frac{1}{\left(1 - \frac{R_s}{r}\right)} dr^2 + r^2 d\theta^2

Now, we want to compare this to the metric of our embedding space, which in our case is 3D Euclidean metric in cylindrical coordinates (we chose cylindrical because it was the only coordinate system that contained a zz coordinate that we'll need in embedding but also rr and θ\theta that the Schwarzschild metric uses):

ds2=dr2+r2dθ2+dz2 ds^2 = dr^2 + r^2 d\theta^2 + dz^2

If we rearrange the terms of the Euclidean metric, we get:

ds2=dz2+dr2+r2dθ2 ds^2 = dz^2 + dr^2 + r^2 d\theta^2

Comparing the similar terms in these two metrics, we can find that:

dz2+dr2=1(1Rsr)dr2 dz^2 + dr^2 = \frac{1}{\left(1 - \frac{R_s}{r}\right)} dr^2

This is a differential equation we can solve for zz. First, note that:

dz2=1(1Rsr)dr2dr2 dz^2 = \frac{1}{\left(1 - \frac{R_s}{r}\right)} dr^2 - dr^2

If we factor out dr2dr^2, we get:

dz2=dr2[11Rsr1] dz^2 = dr^2 \left[\frac{1}{1 - \frac{R_s}{r}} -1\right]

Which means that:

(dzdr)2=[11Rsr1] \left(\frac{dz}{dr}\right)^2 = \left[\frac{1}{1 - \frac{R_s}{r}} -1\right]

This simplifies to:

(dzdr)2=RsrRs \left(\frac{dz}{dr}\right)^2 = \frac{R_s}{r - R_s}

So:

z(r,θ)=Rs1rRsdr=2RsrRs z(r, \theta) = \sqrt{R_s} \int \frac{1}{\sqrt{r - R_s}} dr = 2\sqrt{R_s} \sqrt{r - R_s}

An interactive 3D plot of this embedding diagram can be found at https://demonstrations.wolfram.com/SchwarzschildSpaceTimeEmbeddingDiagram/.

Note that the zz coordinate in the final embedding diagram has no physical meaning! It is only used as a convenient way for us to grasp the curvature of the space by embedding the 2D curvature into a higher dimension (in our case, 3D). Only the rr coordinate (which corresponds to x2+y2\sqrt{x^2 + y^2}) and θ\theta coordinate have a physical meaning.

Back to home