Schwarzschild metric

In Einstein's theory of general relativity, the Schwarzschild metric (also known as the Schwarzschild solution) is an exact solution to the Einstein field equations that describes the gravitational field outside a spherical mass, on the assumption that the electric charge of the mass, angular momentum of the mass, and universal cosmological constant are all zero. The solution is a useful approximation for describing slowly rotating astronomical objects such as many stars and planets, including Earth and the Sun. It was found by Karl Schwarzschild in 1916.

According to Birkhoff's theorem, the Schwarzschild metric is the most general spherically symmetric vacuum solution of the Einstein field equations. A Schwarzschild black hole or static black hole is a black hole that has neither electric charge nor angular momentum (non-rotating). A Schwarzschild black hole is described by the Schwarzschild metric, and cannot be distinguished from any other Schwarzschild black hole except by its mass.

The Schwarzschild black hole is characterized by a surrounding spherical boundary, called the event horizon, which is situated at the Schwarzschild radius ($r_{\text{s}}$), often called the radius of a black hole. The boundary is not a physical surface, and a person who fell through the event horizon (before being torn apart by tidal forces) would not notice any physical surface at that position; it is a mathematical surface which is significant in determining the black hole's properties. Any non-rotating and non-charged mass that is smaller than its Schwarzschild radius forms a black hole. The solution of the Einstein field equations is valid for any mass M, so in principle (within the theory of general relativity) a Schwarzschild black hole of any mass could exist if conditions became sufficiently favorable to allow for its formation.

In the vicinity of a Schwarzschild black hole, space curves so much that even light rays are deflected, and very nearby light can be deflected so much that it travels several times around the black hole.

Formulation

Further information: Derivation of the Schwarzschild solution

The Schwarzschild metric is a spherically symmetric Lorentzian metric (here, with signature convention (+, -, -, -)), defined on (a subset of) $$\mathbb {R} \times \left(E^{3}-O\right)\cong \mathbb {R} \times (0,\infty )\times S^{2}$$ where $E^{3}$ is 3 dimensional Euclidean space, and $S^{2}\subset E^{3}$ is the two sphere. The rotation group $\mathrm {SO} (3)=\mathrm {SO} (E^{3})$ acts on the $E^{3}-O$ or $S^{2}$ factor as rotations around the center $O$, while leaving the first $\mathbb {R}$ factor unchanged. The Schwarzschild metric is a solution of Einstein's field equations in empty space, meaning that it is valid only outside the gravitating body. That is, for a spherical body of radius $R$ the solution is valid for $r>R$. To describe the gravitational field both inside and outside the gravitating body the Schwarzschild solution must be matched with some suitable interior solution at ⁠$r=R$⁠, such as the interior Schwarzschild metric.

In Schwarzschild coordinates $(t,r,\theta ,\phi )$ the Schwarzschild metric (or equivalently, the line element for proper time) has the form $${ds}^{2}=c^{2}\,{d\tau }^{2}=\left(1-{\frac {r_{\mathrm {s} }}{r}}\right)c^{2}\,dt^{2}-\left(1-{\frac {r_{\mathrm {s} }}{r}}\right)^{-1}\,dr^{2}-r^{2}{d\Omega }^{2},$$ where ${d\Omega }^{2}$ is the metric on the two sphere, i.e. ⁠${d\Omega }^{2}=\left(d\theta ^{2}+\sin ^{2}\theta \,d\phi ^{2}\right)$⁠. Furthermore,

• $d\tau ^{2}$ is positive for timelike curves, in which case $\tau$ is the proper time (time measured by a clock moving along the same world line with a test particle),
• $c$ is the speed of light,
• $t$ is, for ⁠$r>r_{\text{s}}$⁠, the time coordinate (measured by a clock located infinitely far from the massive body and stationary with respect to it),
• $r$ is, for ⁠$r>r_{\text{s}}$⁠, the radial coordinate (measured as the circumference, divided by 2π, of a sphere centered around the massive body),
• $\Omega$ is a point on the two sphere ⁠$S^{2}$⁠,
• $\theta$ is the colatitude of $\Omega$ (angle from north, in units of radians) defined after arbitrarily choosing a z-axis,
• $\phi$ is the longitude of $\Omega$ (also in radians) around the chosen z-axis, and
• $r_{\text{s}}$ is the Schwarzschild radius of the massive body, a scale factor which is related to its mass $M$ by ⁠$r_{\text{s}}={2GM}/{c^{2}}$⁠, where $G$ is the gravitational constant.

The Schwarzschild metric has a singularity for r = 0, which is an intrinsic curvature singularity. It also seems to have a singularity on the event horizon r = r_s. Depending on the point of view, the metric is therefore defined only on the exterior region $r>r_{\text{s}}$, only on the interior region $r<r_{\text{s}}$ or their disjoint union. However, the metric is actually non-singular across the event horizon, as one sees in suitable coordinates (see below). For ⁠$r\gg r_{\text{s}}$⁠, the Schwarzschild metric is asymptotic to the standard Lorentz metric on Minkowski space. For almost all astrophysical objects, the ratio ${\frac {r_{\text{s}}}{R}}$ is extremely small. For example, the Schwarzschild radius $r_{\text{s}}^{({\text{Earth}})}$ of the Earth is roughly 8.9 mm, while the Sun, which is 3.3×10⁵ times as massive has a Schwarzschild radius $r_{\text{s}}^{({\text{Sun}})}$ of approximately 3.0 km. The ratio becomes large only in close proximity to black holes and other ultra-dense objects such as neutron stars.

The radial coordinate turns out to have physical significance as the "proper distance between two events that occur simultaneously relative to the radially moving geodesic clocks, the two events lying on the same radial coordinate line".

The Schwarzschild solution is analogous to a classical Newtonian theory of gravity that corresponds to the gravitational field around a point particle. Even at the surface of the Earth, the corrections to Newtonian gravity are only one part in a billion.

History

The Schwarzschild solution is named in honour of Karl Schwarzschild, who found the exact solution in 1915 and published it in January 1916, a little more than a month after the publication of Einstein's theory of general relativity. It was the first exact solution of the Einstein field equations other than the trivial flat space solution. Schwarzschild died shortly after his paper was published, as a result of a disease (thought to be pemphigus) he developed while serving in the German army during World War I.

Johannes Droste in 1916 independently produced the same solution as Schwarzschild, using a simpler, more direct derivation.

In the early years of general relativity there was a lot of confusion about the nature of the singularities found in the Schwarzschild and other solutions of the Einstein field equations. In Schwarzschild's original paper, he put what we now call the event horizon at the origin of his coordinate system. In this paper he also introduced what is now known as the Schwarzschild radial coordinate (r in the equations above), as an auxiliary variable. In his equations, Schwarzschild was using a different radial coordinate that was zero at the Schwarzschild radius.

A more complete analysis of the singularity structure was given by David Hilbert in the following year, identifying the singularities both at r = 0 and r = r_s. Although there was general consensus that the singularity at r = 0 was a 'genuine' physical singularity, the nature of the singularity at r = r_s remained unclear.

In 1921, Paul Painlevé and in 1922 Allvar Gullstrand independently produced a metric, a spherically symmetric solution of Einstein's equations, which we now know is coordinate transformation of the Schwarzschild metric, Gullstrand–Painlevé coordinates, in which there was no singularity at r = r_s. They, however, did not recognize that their solutions were just coordinate transforms, and in fact used their solution to argue that Einstein's theory was wrong. In 1924 Arthur Eddington produced the first coordinate transformation (Eddington–Finkelstein coordinates) that showed that the singularity at r = r_s was a coordinate artifact, although he also seems to have been unaware of the significance of this discovery. Later, in 1932, Georges Lemaître gave a different coordinate transformation (Lemaître coordinates) to the same effect and was the first to recognize that this implied that the singularity at r = r_s was not physical. In 1939 Howard Robertson showed that a free falling observer descending in the Schwarzschild metric would cross the r = r_s singularity in a finite amount of proper time even though this would take an infinite amount of time in terms of coordinate time t.

In 1950, John Synge produced a paper that showed the maximal analytic extension of the Schwarzschild metric, again showing that the singularity at r = r_s was a coordinate artifact and that it represented two horizons. A similar result was later rediscovered by George Szekeres, and independently Martin Kruskal. The new coordinates nowadays known as Kruskal–Szekeres coordinates were much simpler than Synge's but both provided a single set of coordinates that covered the entire spacetime. However, perhaps due to the obscurity of the journals in which the papers of Lemaître and Synge were published their conclusions went unnoticed, with many of the major players in the field including Einstein believing that the singularity at the Schwarzschild radius was physical. Synge's later derivation of the Kruskal–Szekeres metric solution, which was motivated by a desire to avoid "using 'bad' [Schwarzschild] coordinates to obtain 'good' [Kruskal–Szekeres] coordinates", has been generally under-appreciated in the literature, but was adopted by Chandrasekhar in his black hole monograph.

Real progress was made in the 1960s when the mathematically rigorous formulation cast in terms of differential geometry entered the field of general relativity, allowing more exact definitions of what it means for a Lorentzian manifold to be singular. This led to definitive identification of the r = r_s singularity in the Schwarzschild metric as an event horizon, i.e., a hypersurface in spacetime that can be crossed in only one direction.

Singularities and black holes

The Schwarzschild solution appears to have singularities at r = 0 and r = r_s; some of the metric components "blow up" (entail division by zero or multiplication by infinity) at these radii. Since the Schwarzschild metric is expected to be valid only for those radii larger than the radius R of the gravitating body, there is no problem as long as R > r_s. For ordinary stars and planets this is always the case. For example, the radius of the Sun is approximately 700000 km, while its Schwarzschild radius is only 3 km.

The singularity at r = r_s divides the Schwarzschild coordinates in two disconnected patches. The exterior Schwarzschild solution with r > r_s is the one that is related to the gravitational fields of stars and planets. The interior Schwarzschild solution with 0 ≤ r < r_s, which contains the singularity at r = 0, is completely separated from the outer patch by the singularity at r = r_s. The Schwarzschild coordinates therefore give no physical connection between the two patches, which may be viewed as separate solutions. The singularity at r = r_s is an illusion however; it is an instance of what is called a coordinate singularity. As the name implies, the singularity arises from a bad choice of coordinates or coordinate conditions. When changing to a different coordinate system (for example Lemaître coordinates, Eddington–Finkelstein coordinates, Kruskal–Szekeres coordinates, Novikov coordinates, or Gullstrand–Painlevé coordinates) the metric becomes regular at r = r_s and can extend the external patch to values of r smaller than r_s. Using a different coordinate transformation one can then relate the extended external patch to the inner patch.

The case r = 0 is different, however. If one asks that the solution be valid for all r one runs into a true physical singularity, or gravitational singularity, at the origin. To see that this is a true singularity one must look at quantities that are independent of the choice of coordinates. One such important quantity is the Kretschmann invariant, which is given by

$R^{\alpha \beta \gamma \delta }R_{\alpha \beta \gamma \delta }={\frac {12r_{\mathrm {s} }^{2}}{r^{6}}}={\frac {48G^{2}M^{2}}{c^{4}r^{6}}}\,.$

At r = 0 the curvature becomes infinite, indicating the presence of a singularity. At this point the metric cannot be extended in a smooth manner (the Kretschmann invariant involves second derivatives of the metric), spacetime itself is then no longer well-defined. Furthermore, Sbierski showed the metric cannot be extended even in a continuous manner. For a long time it was thought that such a solution was non-physical. However, a greater understanding of general relativity led to the realization that such singularities were a generic feature of the theory and not just an exotic special case.

The Schwarzschild solution, taken to be valid for all r > 0, is called a Schwarzschild black hole. It is a perfectly valid solution of the Einstein field equations, although (like other black holes) it has rather bizarre properties. For r < r_s the Schwarzschild radial coordinate r becomes timelike and the time coordinate t becomes spacelike. A curve at constant r is no longer a possible worldline of a particle or observer, not even if a force is exerted to try to keep it there; this occurs because spacetime has been curved so much that the direction of cause and effect (the particle's future light cone) points into the singularity. The surface r = r_s demarcates what is called the event horizon of the black hole. It represents the point past which light can no longer escape the gravitational field. Any physical object whose radius R becomes less than or equal to the Schwarzschild radius has undergone gravitational collapse and become a black hole.

Alternative coordinates

The Schwarzschild solution can be expressed in a range of different choices of coordinates besides the Schwarzschild coordinates used above. Different choices tend to highlight different features of the solution. The table below shows some popular choices.

Alternative coordinates

Coordinates	Line element	Notes	Features
Eddington–Finkelstein coordinates (ingoing)	$\left(1-{\frac {r_{\mathrm {s} }}{r}}\right)\,dv^{2}-2\,dv\,dr-r^{2}\,g_{\Omega }$		regular at future horizon past horizon is at v = −∞
Eddington–Finkelstein coordinates (outgoing)	$\left(1-{\frac {r_{\mathrm {s} }}{r}}\right)\,du^{2}+2\,du\,dr-r^{2}g_{\Omega }$		regular at past horizon extends across past horizon future horizon at u = ∞
Gullstrand–Painlevé coordinates	$\left(1-{\frac {r_{\mathrm {s} }}{r}}\right)\,dT^{2}\pm 2{\sqrt {\frac {r_{\mathrm {s} }}{r}}}\,dT\,dr-dr^{2}-r^{2}\,g_{\Omega }$		regular at past and future horizons
Isotropic coordinates	${\frac {\left(1-{\frac {r_{\mathrm {s} }}{4R}}\right)^{2}}{\left(1+{\frac {r_{\mathrm {s} }}{4R}}\right)^{2}}}\,{dt}^{2}-\left(1+{\frac {r_{\mathrm {s} }}{4R}}\right)^{4}\,\left(dx^{2}+dy^{2}+dz^{2}\right)$	$R={\sqrt {x^{2}+y^{2}+z^{2}}}$ Valid only outside the event horizon: $R>r_{\text{s}}/4$	isotropic lightcones on constant time slices
Kruskal–Szekeres coordinates	${\frac {4r_{\mathrm {s} }^{3}}{r}}e^{-{\frac {r}{r_{\mathrm {s} }}}}\,\left(dT^{2}-dR^{2}\right)-r^{2}\,g_{\Omega }$	$T^{2}-R^{2}=\left(1-{\frac {r}{r_{\mathrm {s} }}}\right)e^{\frac {r}{r_{\mathrm {s} }}}$	regular at horizon; maximally extends to full spacetime
Lemaître coordinates	$dT^{2}-{\frac {r_{\mathrm {s} }}{r}}\,dR^{2}-r^{2}\,g_{\Omega }$	$r=\left({\tfrac {3}{2}}(R\pm T)\right)^{\frac {2}{3}}r_{\mathrm {s} }^{\frac {1}{3}}$	regular at either past or future horizon
Harmonic coordinates	${\frac {\rho -r_{\mathrm {s} }/2}{\rho +r_{\mathrm {s} }/2}}dt^{2}-{\frac {\rho +r_{\mathrm {s} }/2}{\rho -r_{\mathrm {s} }/2}}d\rho ^{2}-(\rho +r_{\mathrm {s} }/2)^{2}g_{\Omega }$	$\rho =r-r_{\mathrm {s} }/2$

In table above, some shorthand has been introduced for brevity. The speed of light c has been set to one. The notation

$g_{\Omega }=d\theta ^{2}+\sin ^{2}\theta \,d\varphi ^{2}$

is used for the metric of a unit radius 2-dimensional sphere. Moreover, in each entry R and T denote alternative choices of radial and time coordinate for the particular coordinates. Note, the R or T may vary from entry to entry.

The Kruskal–Szekeres coordinates have the form to which the Belinski–Zakharov transform can be applied. This implies that the Schwarzschild black hole is a form of gravitational soliton.

Flamm's paraboloid

A plot of Flamm's paraboloid. It should not be confused with the unrelated concept of a gravity well.

The spatial curvature of the Schwarzschild solution for r > r_s can be visualized as the graphic shows. Consider a constant time equatorial slice H through the Schwarzschild solution by fixing θ = ⁠π/2⁠, t = constant, and letting the remaining Schwarzschild coordinates (r, φ) vary. Imagine now that there is an additional Euclidean dimension w, which has no physical reality (it is not part of spacetime). Then replace the (r, φ) plane with a surface dimpled in the w direction according to the equation (Flamm's paraboloid)

$w=2{\sqrt {r_{\mathrm {s} }\left(r-r_{\mathrm {s} }\right)}}.$

This surface has the property that distances measured within it match distances in the Schwarzschild metric, because with the definition of w above,

$dw^{2}+dr^{2}+r^{2}\,d\varphi ^{2}={\frac {dr^{2}}{1-{\frac {r_{\mathrm {s} }}{r}}}}+r^{2}\,d\varphi ^{2}=-ds^{2}$

Thus, Flamm's paraboloid is useful for visualizing the spatial curvature of the Schwarzschild metric. It should not, however, be confused with a gravity well. No ordinary (massive or massless) particle can have a worldline lying on the paraboloid, since all distances on it are spacelike (this is a cross-section at one moment of time, so any particle moving on it would have an infinite velocity). A tachyon could have a spacelike worldline that lies entirely on a single paraboloid. However, even in that case its geodesic path is not the trajectory one gets through a "rubber sheet" analogy of gravitational well: in particular, if the dimple is drawn pointing upward rather than downward, the tachyon's geodesic path still curves toward the central mass, not away. See the gravity well article for more information.

Flamm's paraboloid may be derived as follows. The Euclidean metric in the cylindrical coordinates (r, φ, w) is written

$-ds^{2}=dw^{2}+dr^{2}+r^{2}\,d\varphi ^{2}\,.$

Letting the surface be described by the function w = w(r), the Euclidean metric can be written as

$-ds^{2}=\left(1+\left({\frac {dw}{dr}}\right)^{2}\right)\,dr^{2}+r^{2}\,d\varphi ^{2}\,,$

Comparing this with the Schwarzschild metric in the equatorial plane (θ = ⁠π/2⁠) at a fixed time (t = constant, dt = 0)

$-ds^{2}=\left(1-{\frac {r_{\mathrm {s} }}{r}}\right)^{-1}\,dr^{2}+r^{2}\,d\varphi ^{2}\,,$

yields an integral expression for w(r):

$w(r)=\int {\frac {dr}{\sqrt {{\frac {r}{r_{\mathrm {s} }}}-1}}}=2r_{\mathrm {s} }{\sqrt {{\frac {r}{r_{\mathrm {s} }}}-1}}+{\mbox{constant}}$

whose solution is Flamm's paraboloid.

Orbital motion

Comparison between the orbit of a test particle in Newtonian (left) and Schwarzschild (right) spacetime; note the apsidal precession on the right.

Further information: Schwarzschild geodesics

A particle orbiting in the Schwarzschild metric can have a stable circular orbit with r > 3r_s. Circular orbits with r between 1.5r_s and 3r_s are unstable, and no circular orbits exist for r < 1.5r_s. The circular orbit of minimum radius 1.5r_s corresponds to an orbital velocity approaching the speed of light. It is possible for a particle to have a constant value of r between r_s and 1.5r_s, but only if some force acts to keep it there.

Noncircular orbits, such as Mercury's, dwell longer at small radii than would be expected in Newtonian gravity. This can be seen as a less extreme version of the more dramatic case in which a particle passes through the event horizon and dwells inside it forever. Intermediate between the case of Mercury and the case of an object falling past the event horizon, there are exotic possibilities such as knife-edge orbits, in which the satellite can be made to execute an arbitrarily large number of nearly circular orbits, after which it flies back outward.

Symmetries

The isometry group of the Schwarzchild metric is ⁠$\mathbb {R} \times \mathrm {O} (3)\times \{\pm 1\}$⁠, where $\mathrm {O} (3)$ is the orthogonal group of rotations and reflections in three dimensions, $\mathbb {R}$ comprises the time translations, and $\{\pm 1\}$ is the group generated by time reversal.

This is thus the subgroup of the ten-dimensional Poincaré group which takes the time axis (trajectory of the star) to itself. It omits the spatial translations (three dimensions) and boosts (three dimensions). It retains the time translations (one dimension) and rotations (three dimensions). Thus it has four dimensions. Like the Poincaré group, it has four connected components: the component of the identity; the time reversed component; the spatial inversion component; and the component which is both time reversed and spatially inverted.

Curvatures

The Ricci curvature scalar and the Ricci curvature tensor are both zero. Non-zero components of the Riemann curvature tensor are given by

$-R^{t}{}_{rtr}=2R^{\theta }{}_{r\theta r}=2R^{\phi }{}_{r\phi r}={\frac {r_{\text{s}}}{r^{2}(r_{\text{s}}-r)}},$

$2R^{t}{}_{\theta t\theta }=2R^{r}{}_{\theta r\theta }=-R^{\phi }{}_{\theta \phi \theta }=-{\frac {r_{\text{s}}}{r}},$

$2R^{t}{}_{\phi t\phi }=2R^{r}{}_{\phi r\phi }=-R^{\theta }{}_{\phi \theta \phi }=-{\frac {r_{\text{s}}\sin ^{2}(\theta )}{r}},$

$R^{r}{}_{trt}=-2R^{\theta }{}_{t\theta t}=-2R^{\phi }{}_{t\phi t}=c^{2}{\frac {r_{\text{s}}(r_{\text{s}}-r)}{r^{4}}},$

from which one can see that $R^{\gamma }{}_{\alpha \gamma \beta }=0$. Six of these formulas are Eq. 5.13 in Carroll and imply the other 6 by $R^{\alpha }{}_{\beta \gamma \delta }=g^{\alpha \kappa }g_{\beta \lambda }R^{\lambda }{}_{\kappa \delta \gamma }$. Components which are obtainable by other symmetries of the Riemann tensor are not displayed.

To understand the physical meaning of these quantities, it is useful to express the curvature tensor in an orthonormal basis. In an orthonormal basis of an observer the non-zero components in geometric units are

$R^{\hat {r}}{}_{{\hat {t}}{\hat {r}}{\hat {t}}}=-R^{\hat {\theta }}{}_{{\hat {\phi }}{\hat {\theta }}{\hat {\phi }}}=-{\frac {r_{\text{s}}}{r^{3}}},$

$R^{\hat {\theta }}{}_{{\hat {t}}{\hat {\theta }}{\hat {t}}}=R^{\hat {\phi }}{}_{{\hat {t}}{\hat {\phi }}{\hat {t}}}=-R^{\hat {r}}{}_{{\hat {\theta }}{\hat {r}}{\hat {\theta }}}=-R^{\hat {r}}{}_{{\hat {\phi }}{\hat {r}}{\hat {\phi }}}={\frac {r_{\text{s}}}{2r^{3}}}.$

Again, components which are obtainable by the symmetries of the Riemann tensor are not displayed. These results are invariant to any Lorentz boost, thus the components do not change for non-static observers. The geodesic deviation equation shows that the tidal acceleration between two observers separated by $\xi ^{\hat {j}}$ is $D^{2}\xi ^{\hat {j}}/D\tau ^{2}=-R^{\hat {j}}{}_{{\hat {t}}{\hat {k}}{\hat {t}}}\xi ^{\hat {k}}$, so a body of length $L$ is stretched in the radial direction by an apparent acceleration $(r_{\text{s}}/r^{3})c^{2}L$ and squeezed in the perpendicular directions by $-(r_{\text{s}}/(2r^{3}))c^{2}L$.

See also

• Derivation of the Schwarzschild solution
• Reissner–Nordström metric (charged, non-rotating solution)
• Kerr metric (uncharged, rotating solution)
• Kerr–Newman metric (charged, rotating solution)
• Black hole, a general review
• Schwarzschild coordinates
• Kruskal–Szekeres coordinates
• Eddington–Finkelstein coordinates
• Gullstrand–Painlevé coordinates
• Lemaître coordinates (Schwarzschild solution in synchronous coordinates)
• Frame fields in general relativity (Lemaître observers in the Schwarzschild vacuum)
• Tolman–Oppenheimer–Volkoff equation (metric and pressure equations of a static and spherically symmetric body of isotropic material)
• Planck length

References

• Schwarzschild, K. (1916). "Über das Gravitationsfeld eines Massenpunktes nach der Einsteinschen Theorie". Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften. 7: 189–196. Bibcode:1916AbhKP1916..189S.

• Text of the original paper, in Wikisource
• Translation: Antoci, S.; Loinger, A. (1999). "On the gravitational field of a mass point according to Einstein's theory". arXiv:physics/9905030.
• A commentary on the paper, giving a simpler derivation: Bel, L. (2007). "Über das Gravitationsfeld eines Massenpunktesnach der Einsteinschen Theorie". arXiv:0709.2257 [gr-qc].

• Schwarzschild, K. (1916). "Über das Gravitationsfeld einer Kugel aus inkompressibler Flüssigkeit". Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften. 1: 424.

• Text of the original paper, in Wikisource
• Translation: Antoci, S. (1999). "On the gravitational field of a sphere of incompressible fluid according to Einstein's theory". arXiv:physics/9912033.

• Flamm, L. (1916). "Beiträge zur Einstein'schen Gravitationstheorie". Physikalische Zeitschrift. 17: 448.
• Adler, R.; Bazin, M.; Schiffer, M. (1975). Introduction to General Relativity (2nd ed.). McGraw-Hill. Chapter 6. ISBN 0-07-000423-4.
• Landau, L. D.; Lifshitz, E. M. (1951). The Classical Theory of Fields. Course of Theoretical Physics. Vol. 2 (4th Revised English ed.). Pergamon Press. Chapter 12. ISBN 0-08-025072-6.
• Misner, C. W.; Thorne, K. S.; Wheeler, J. A. (1970). Gravitation. W.H. Freeman. Chapters 31 and 32. ISBN 0-7167-0344-0.
• Weinberg, S. (1972). Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity. John Wiley & Sons. Chapter 8. ISBN 0-471-92567-5.
• Taylor, E. F.; Wheeler, J. A. (2000). Exploring Black Holes: Introduction to General Relativity. Addison-Wesley. ISBN 0-201-38423-X.
• Heinzle, J. M.; Steinbauer, R. (2002). "Remarks on the distributional Schwarzschild geometry". Journal of Mathematical Physics. 43 (3): 1493–1508. arXiv:gr-qc/0112047. Bibcode:2002JMP....43.1493H. doi:10.1063/1.1448684. S2CID 119677857.