To the Schwarzschild Solution in General Relativity

The paper is concerned with the history of the spherically symmetric static problem solution of General Relativity found in 1916 by K. Schwarzschild   which is interpreted in modern physics as the background of the objects referred to as Black Holes. First, the modern interpretation this solution which does not exactly coincide with original solution obtained by K. Schwarzschild is discussed. Second, the basic equations of the original Schwarzschild solution are presented in modern notations allowing us to compare existing and original solutions. Finally, a modification of the Schwarzschild approach is proposed allowing us to arrive at the exact solution of the Schwarzschild problem.

Share and Cite:

Vasiliev, V. and Fedorov, L. (2018) To the Schwarzschild Solution in General Relativity. Journal of Modern Physics, 9, 2482-2494. doi: 10.4236/jmp.2018.914160.

1. Spherically Symmetric Static Problem of General Relativity

Spherically symmetric problem is one of the most discussed problems of General Relativity Theory (GRT) widely described in the literature      . This paper is concerned with the analysis of the original Schwarzschild solution of this problem in association with its modern interpretation and possible generalization.

The line element for the spherically symmetric problem is traditionally taken in the following form:

$\text{d}{s}^{2}={g}^{2}\text{d}{r}^{2}+{\rho }^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-{h}^{2}{c}^{2}\text{d}{t}^{2}$ (1)

Here $r,\theta ,\phi$ and t are space spherical and time coordinates, $g,\rho ,h$ are the metric coefficients that depend on the radial coordinate r only. The basic equations of the General Relativity Theory (GRT) link the Einstein tensor ${E}_{i}^{i}$ with the metric coefficients as 

${E}_{1}^{1}=\frac{1}{{\rho }^{2}}-\frac{1}{{g}^{2}}\frac{{\rho }^{\prime }}{\rho }\left(\frac{{\rho }^{\prime }}{\rho }+\frac{2{h}^{\prime }}{h}\right)$ (2)

${E}_{2}^{2}={E}_{3}^{3}=-\frac{1}{{g}^{2}}\left[\frac{{h}^{″}}{h}+\frac{{\rho }^{″}}{\rho }+\frac{{\rho }^{\prime }}{\rho }\left(\frac{{h}^{\prime }}{h}-\frac{{g}^{\prime }}{g}\right)-\frac{{g}^{\prime }{h}^{\prime }}{gh}\right]$ (3)

${E}_{4}^{4}=\frac{1}{{\rho }^{2}}-\frac{1}{{g}^{2}}\left[{\left(\frac{{\rho }^{\prime }}{\rho }\right)}^{2}+\frac{2{\rho }^{″}}{\rho }-\frac{2{\rho }^{\prime }{g}^{\prime }}{\rho g}\right]$ (4)

where ${\left(...\right)}^{\prime }=\text{d}\left(...\right)/\text{d}r$. The components of the Einstein tensor are proportional to the energy tensor ${T}_{i}^{i}$, i.e.

${E}_{i}^{i}=\kappa {T}_{i}^{i}$ (5)

where

$\kappa =8\text{π}\gamma /{c}^{4}$ (6)

is the GRT gravitational constant depending on the Newton constant $\gamma$ and the velocity of light c. Finally, the energy tensor (and the Einstein tensor which is proportional to it) must satisfy the following conservation equation:

${\left({T}_{1}^{1}\right)}^{\prime }-\frac{2{\rho }^{\prime }}{\rho }\left({T}_{2}^{2}-{T}_{1}^{1}\right)+\frac{{h}^{\prime }}{h}\left({T}_{1}^{1}-{T}_{4}^{4}\right)=0$ (7)

We use mixed components of the tensors E and T because for the problem under study they coincide with the corresponding physical components. The energy tensor depends on the space structure. Particularly, for the empty space

${T}_{i}^{i}=0,\left(i=1,2,3,4\right)$ (8)

and Equations (2)-(4) are homogeneous. Inside a solid sphere with radius $a$,

${T}_{1}^{1}={\sigma }_{r},\text{\hspace{0.17em}}{T}_{2}^{2}={T}_{3}^{3}={\sigma }_{\theta },\text{\hspace{0.17em}}{T}_{4}^{4}=\mu {c}^{2}$ (9)

where ${\sigma }_{r}$ and ${\sigma }_{\theta }$ are the radial and the circumferential stresses and $\mu$ is the material density.

Consider the external space $\left(r\ge a\right)$. Taking ${E}_{4}^{4}=0$, we can reduce Equation (4) to the following form:

$\frac{\text{d}}{\text{d}r}\left[\frac{{\rho }_{e}{\left({{\rho }^{\prime }}_{e}\right)}^{2}}{{g}_{e}^{2}}\right]={{\rho }^{\prime }}_{e}$ (10)

Equation (10) can be readily integrated to give

${g}_{e}^{2}=\frac{{\rho }_{e}{\left({{\rho }^{\prime }}_{e}\right)}^{2}}{{\rho }_{e}+{C}_{1}}$ (11)

where ${C}_{1}$ is the integration constant and the functions with subscript “e” correspond to the external space. Substituting this result in Equation (2), taking ${E}_{1}^{1}=0$ and integrating, we get

${h}_{e}^{2}={C}_{2}\left(1+\frac{{C}_{1}}{{\rho }_{e}}\right)$ (12)

It looks like substituting Equations (11) and (12) into Equation (3) in which ${E}_{2}^{2}=0$, we can determine the function ${\rho }_{e}\left(r\right)$. But this is not the case―under such substitution, Equation (3) is satisfied identically for any function ${\rho }_{e}\left(r\right)$. This result can be predicted―since the components of the Einstein tensor ${E}_{i}^{i}$ satisfy Equation (7), Equations (2)-(4) are not mutually independent and any solution of two of them identically satisfy the third equation.

Thus, we have only two Einstein equations for three unknown functions $g\left(r\right)$, $\rho \left(r\right)$ and $h\left(r\right)$. The fact that the set of GTR equations is not complete and must be supplemented with some coordinate conditions was first mentioned by D. Hilbert  . By now, the general form of these conditions has not been developed, though some particular forms (e.g., the so-called De-Donder-Fock harmonic coordinate conditions) have been used in spherically symmetric problem  .

Not knowing the function ${\rho }_{e}\left(r\right)$, we can make some qualitative conclusions about its behavior. It is natural to suppose the for $r\to \infty$ we should have ${\rho }_{e}\to r$ and Equations (11) and (12) should reduce to the corresponding results of the Newton gravitation theory, according to which 

${g}_{\infty }^{2}=1-\frac{2\phi }{{c}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{h}_{\infty }^{2}=1+\frac{2\phi }{{c}^{2}}$ (13)

where $\phi =-\gamma m/r$ is the Newton gravitational potential and m is the mass inducing the gravitation field. Equations (13) allows us to determine the constants in Equations (11) and (12) and to present the result as

${g}_{e}^{2}=\frac{{\left({{\rho }^{\prime }}_{e}\right)}^{2}}{1-{r}_{g}/{\rho }_{e}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{h}_{e}^{2}=1-\frac{{r}_{g}}{{\rho }_{e}}$ (14)

Here,

${r}_{g}=2m\gamma /{c}^{2}$ (15)

is the so-called gravitational radius sometimes referred to as the Schwarzschild radius (though K. Schwarzschild did not use this term).

Consider the internal space $\left(0\le r\le a\right)$. For an elastic sphere with known density, we have totally four equations, i.e., Equations (2)-(4) in which the left-hand parts are specified by Equations (5) and (9) and Equation (7) in which the energy tensor should be expressed with the aid of Equation (9). These equations include five unknown functions―the metric coefficients $g,\rho ,h$ and the stresses ${\sigma }_{r}$ and ${\sigma }_{\theta }$. To solve the problem, we must supplement the GRT equations with the equations for stresses similar to the compatibility equations of the Theory of Elasticity. Such equations can be derived  , but we restrict ourselves to a particular solution obtained by K. Schwarzschild for a sphere consisting of a perfect incompressible fluid. In this case, ${\sigma }_{r}={\sigma }_{\theta }=-p\left(r\right)$ in which p is the pressure in the fluid and the fluid density ${\mu }_{0}$ does not depend on p and r. Taking ${E}_{4}^{4}=\kappa {\mu }_{0}{c}^{2}$ in Equation (4), we can reduce it to the form 

$\frac{\text{d}}{\text{d}r}\left[\frac{{\rho }_{i}}{{g}_{i}^{2}}\left({{\rho }^{\prime }}_{i}\right)\right]=\left(1-\kappa {\mu }_{0}{c}^{2}{\rho }_{i}^{2}\right){{\rho }^{\prime }}_{i}$ (16)

Equation (16) can be readily integrated to give

${g}_{i}^{2}=\frac{{\left({\rho }^{\prime }\right)}^{2}}{1-{\lambda }^{2}{\rho }_{i}^{2}-{C}_{3}/{\rho }_{i}}$ (17)

where

${\lambda }^{2}=\frac{1}{3}\kappa {\mu }_{0}{c}^{2}$ (18)

and the functions with subscript “i” correspond to the internal space. We do not know the function ${\rho }_{i}\left(r\right)$, but can make some reasonable prediction concerning its behavior. Particularly, it is natural to suppose that at the sphere center ${\rho }_{i}=r=0$. Then, in accordance with the symmetry condition at the sphere center, we should take ${С}_{3}=0$ and Equation (17) becomes

${g}_{i}^{2}=\frac{{\left({{\rho }^{\prime }}_{i}\right)}^{2}}{1-{\lambda }^{2}{\rho }_{i}^{2}}$ (19)

Continue the derivation and consider Equations (2) and (7). Taking ${E}_{1}^{1}=-\kappa p$ and ${T}_{1}^{1}={T}_{2}^{2}=-p$, ${T}_{4}^{4}={\mu }_{0}{c}^{2}$, we arrive at

$\frac{1}{{\rho }_{i}^{2}}-\frac{{{\rho }^{\prime }}_{i}}{{g}_{i}^{2}{\rho }_{i}}\left(\frac{{{\rho }^{\prime }}_{i}}{{\rho }_{i}}-\frac{2{{h}^{\prime }}_{i}}{{h}_{i}}\right)=-\kappa p,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{p}^{\prime }+\frac{{{h}^{\prime }}_{i}}{h}\left(p+{\mu }_{0}{c}^{2}\right)=0$ (20)

Consider the first Equation. Substituting Equation (19), express the time metric coefficient ${h}_{i}$, i.e.,

$\frac{\text{d}{h}_{i}}{{h}_{i}}=\frac{\text{d}{\rho }_{i}}{2}\left[\frac{{\rho }_{i}}{1-{\lambda }^{2}{\rho }_{i}^{2}}\left(\kappa p+\frac{1}{{\rho }_{i}^{2}}\right)-\frac{1}{{\rho }_{i}}\right]$ (21)

and rewrite the second equation of Equations (20) as

$\text{d}p+\frac{\text{d}{h}_{i}}{{h}_{i}}\left(p+{\mu }_{0}{c}^{2}\right)=0$ (22)

Substituting Equation (21) and using Equation (18), we arrive at the following equation for the pressure:

$\text{d}p+\frac{{\lambda }^{2}{\rho }_{i}\text{d}{\rho }_{i}}{2\left(1-{\lambda }^{2}{\rho }_{i}^{2}\right)}\left(p+{\mu }_{0}{c}^{2}\right)\left(1+\frac{3p}{{\mu }_{0}{c}^{2}}\right)=0$ (23)

The solution of Equation (23) that satisfies the boundary condition on the free sphere surface $r=a$ according to which $p\left({\rho }_{i}={\rho }_{a}\right)=0$ is

$p=-{\mu }_{0}{c}^{2}\frac{\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}-\sqrt{1-{\lambda }^{2}{\rho }_{a}^{2}}}{\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}-3\sqrt{1-{\lambda }^{2}{\rho }_{a}^{2}}}$ (24)

where ${\rho }_{a}={\rho }_{i}\left(r=a\right)={\rho }_{e}\left(r=a\right)$. To determine ${h}_{i}$, integrate Equation (22) to get ${h}_{i}\left(p+{\mu }_{0}{c}^{2}\right)={C}_{4}$ in which ${C}_{4}$ is the integration constant that can be found from the boundary condition on the sphere surface $r=a$ according to which ${h}_{e}\left({\rho }_{e}={\rho }_{a}\right)={h}_{i}\left({\rho }_{i}={\rho }_{a}\right)$. Using the second equation in Equations (14), we finally have

${h}_{i}=\frac{\sqrt{1-{r}_{g}/{\rho }_{a}}}{2\sqrt{1-{\lambda }^{2}{\rho }_{a}^{2}}}\left(3\sqrt{1-{\lambda }^{2}{\rho }_{0}^{2}}-\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}\right)$ (25)

It should be noted that substitution of the obtained solutions ${g}_{i}\left({\rho }_{i}\right),{h}_{i}\left({\rho }_{i}\right)$ and $p\left({\rho }_{i}\right)$ in the Einstein Equation (3) where ${E}_{2}^{2}=-\kappa p$, does not allow us to find the function ${\rho }_{i}\left(r\right)$ because Equation (3) is identically satisfied for any function ${\rho }_{i}\left(r\right)$. The situation is similar to the external space and for the same reason―since the right-hand parts of Equations (2) and (3) satisfy Equation (7), only three of Equations (2)-(4) and (7) are mutually independent.

Thus, Equation (1) which specifies the metric forms of the external and internal spaces of the fluid sphere can be presented as

$\text{d}{s}_{e}^{2}=\frac{{\left({{\rho }^{\prime }}_{e}\right)}^{2}\text{d}{r}^{2}}{1-{r}_{g}/{\rho }_{e}}+{\rho }_{e}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-\left(1-\frac{{r}_{g}}{{\rho }_{e}}\right){c}^{2}\text{d}{t}^{2}$ (26)

$\begin{array}{c}\text{d}{s}_{i}^{2}=\frac{\left({\rho }_{i}^{2}\right)\text{d}{r}^{2}}{1-{\lambda }^{2}{\rho }_{i}^{2}}+{\rho }_{i}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1-{r}_{g}/{\rho }_{a}}{4\left(1-{\lambda }^{2}{\rho }_{a}^{2}\right)}{\left(3\sqrt{1-{\lambda }^{2}{\rho }_{a}^{2}}-\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}\right)}^{2}{c}^{2}\text{d}{t}^{2}\end{array}$ (27)

To fulfill the solution, we need to find two functions ${\rho }_{e}\left(r\right)$ and ${\rho }_{i}\left(r\right)$ such that allow us to satisfy the boundary conditions at the sphere surface $r=a$, i.e.,

${\rho }_{e}\left(a\right)={\rho }_{i}\left(a\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{g}_{e}\left(a\right)={g}_{i}\left(a\right)$ (28)

However, the equations allowing us to determine these functions are missing in GRT. The same problem exists in the general case―as known, the set of Einstein equations is not complete. In the four-dimensional Riemannian space, this set consists of 10 equations

${E}^{ij}=-\kappa {T}^{ij}$ (29)

for ten components of the metric tensor ${g}^{ij}$. However the Einstein tensor satisfies equations which are analogous to Equations (7). As a result, only six of Equations (29) are mutually independent and to determine the metric tensor we should supplement Equations (29) with four coordinate conditions for ${g}^{ij}$. Some authors declare that these conditions cannot be covariant because there forms depend on the particular coordinate frame   . Consider some particular cases.

2. Modern Interpretation of the Schwarzschild Solution

Traditional description of the Schwarzschild solution can be found elsewhere  . The coordinate condition mentioned in the closure of the previous section is taken in the form

${\rho }_{e}\left(r\right)={\rho }_{i}\left(r\right)=r$ (30)

though K. Schwarzschild did not use directly this relationship. Applying Equations (30) to Equations (14) and (19), we can specify the metric coefficients for this case and present the metric form in Equation (1) as

$\text{d}{s}_{e}^{2}=\frac{\text{d}{r}^{2}}{1-{r}_{g}/r}+{r}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-\left(1-\frac{{r}_{g}}{r}\right){c}^{2}\text{d}{t}^{2}$ (31)

$\begin{array}{c}\text{d}{s}_{i}^{2}=\frac{\text{d}{r}^{2}}{1-{\lambda }^{2}{r}^{2}}+{r}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1-{r}_{g}/a}{4\left(1-{\lambda }^{2}{a}^{2}\right)}{\left(3\sqrt{1-{\lambda }^{2}{a}^{2}}-\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}\right)}^{2}{c}^{2}\text{d}{t}^{2}\end{array}$ (32)

To fulfill the solution, we need to satisfy the boundary conditions (28) on the sphere surface $r=a$. The first condition is satisfied automatically, whereas the second one yields

${r}_{g}={\lambda }^{2}{a}^{3}$ (33)

However, the parameters ${r}_{g}$ and ${\lambda }^{2}$ are specified by Equations (15) and (18) and are known. So, Equation (33) cannot be satisfied in the general case and the second boundary condition in Equations (28) is violated. Substituting formally Equations (15) and (18) in Equation (33), we arrive at the following expression:

$m={m}_{e}=\frac{4}{3}\text{π}{\mu }_{0}{a}^{3}$ (34)

which specifies the mass of a homogeneous solid sphere in the Euclidean space. However, the space in GRT is not Euclidean and the mass of the sphere with the metric coefficients corresponding to Equations (19) is

$\begin{array}{c}m=4\text{π}{\mu }_{0}\underset{0}{\overset{a}{\int }}{g}_{i}{r}^{2}\text{d}r=\frac{2}{{\lambda }^{2}}\text{π}a{\mu }_{0}\left(\frac{1}{\lambda a}{\mathrm{sin}}^{-1}\lambda a-\sqrt{1-{\lambda }^{2}{a}^{2}}\right)\\ \approx {m}_{e}\left(1+\frac{15}{128}{\lambda }^{2}{a}^{2}+\cdot \cdot \cdot \right)\end{array}$ (35)

As can be seen, the obtained result does not coincide with Equation (34) and the second boundary condition in Equations (28) is not satisfied. The reason for this discrepancy is associated with Equations (30). Equation (4), being originally of the second order, under transformation in accordance with Equation (30) reduces to the equation of the first order. As a result, the solution does not contain the proper number of integration constants allowing us to satisfy the complete set of the boundary conditions.

Thus, the coordinate conditions in Equations (30) do not look suitable for the problem under study.

3. Original Schwarzschild’s Solution

In 1916 K. Schwarzschild presented the solution of the external spherically symmetric problem  . He did not use the final version of the Einstein equations, however the field equations that he applied can be now associated with Equations (2)-(4). We can suppose that he understood that only two of these equations were mutually independent because he attracted for the analysis only two equations, i.e. Equations (2) and (4), and ignored Equation (3). The third equation which is necessary to solve the problem, was obtained under the following condition imposed on the determinant of the metric tensor:

$|{g}_{mn}|=1$ (36)

Introducing this equation, K. Schwarzschild followed A. Einstein who used it in general theory to specify the coordinate frame  . Governing equations of GRT contain symbols ${\Gamma }_{ij}^{k}$ which include derivatives of the determinant and become zero under condition (33) thus simplifying the equations. However, Equation (33) cannot be directly applied in spherical coordinates in which the volume element in the three-dimensional Euclidean space is $\text{d}v={r}^{2}\mathrm{sin}\theta \text{d}r\text{d}\theta \text{d}\phi$. To overcome this problem, K. Schwarzschild introduced new variables ${x}_{i}$ such that

${x}_{1}={r}^{3}/3,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{2}=-\mathrm{cos}\theta ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{3}=\phi ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{4}=t$ (37)

In new coordinates, the volume element becomes $\text{d}v=\text{d}{x}_{1}\text{d}{x}_{2}\text{d}{x}_{3}$ and the line element takes the form

$\text{d}{s}^{2}={f}_{1}\text{d}{x}_{1}^{2}+{f}_{2}\left[\frac{\text{d}{x}_{2}^{2}}{1-{x}_{2}^{2}}+\left(1-{x}_{2}^{2}\right)\text{d}{x}_{3}^{2}\right]-{f}_{4}\text{d}{x}_{4}^{2}$ (38)

Three functions ${f}_{1},{f}_{2},{f}_{4}$ can be found from Equations (2) and (4) supplemented with equation ${f}_{1}{f}_{2}^{2}{f}_{4}=1$ which follows from Equation (36). The final solution is

${f}_{1}=\frac{{\left(3{x}_{1}+{\beta }^{3}\right)}^{-4/3}}{1-\alpha {\left(3{x}_{1}+{\beta }^{3}\right)}^{-1/3}}$, ${f}_{2}={\left(3{x}_{1}+{\beta }^{3}\right)}^{2/3}$, ${f}_{4}=1-\alpha \left(3{x}_{1}+{\beta }^{3}\right){}^{-1/3}$ (39)

in which $\alpha$ and ${\beta }^{3}$ are the integration constants. As stated in the Schwarzschild paper, this solution identically satisfies Equation (3) (which should be the case).

The final part of the paper can hardly be understood. Directly following K. Schwarzschild, consider the function ${f}_{1}$ which can be a source of singularity. Equating the denominator to zero and using Equation (37) for ${x}_{1}$ we get

$\alpha {\left({r}^{3}+{\beta }^{3}\right)}^{-1/3}=1$ (40)

Let singularity take place at the origin $r=0$. Then, Equation (40) yields $\beta =\alpha$. Introducing the new variable

$R={\left({r}^{3}+{\alpha }^{3}\right)}^{1/3}$ (41)

and using Equations (38) and (39), K. Schwarzschild arrived at the following final result of his paper:

$\text{d}{s}^{2}=\frac{\text{d}{R}^{2}}{1-\alpha /R}+{R}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-\left(1-\alpha /R\right){с}^{2}\text{d}{t}^{2}$ (42)

This form formally coincides with Equation (31), but it should be taken into account that R is not the radial coordinate r. The constant $\alpha$ is declared to depend on the mass located at the origin, but is not found.

As can be seen, the first term in Equation (42) becomes singular if $R=\alpha$ or $r=0$. Thus, the original Schwarzschild solution has only one singular point― $r=0$.

However, it looks like Equation (42) is not correct. To show this, change R to r in Equation (42) with the aid of Equation (41) to get

$\begin{array}{c}\text{d}{s}^{2}=\frac{{\left({r}^{3}+{\alpha }^{3}\right)}^{-4/3}{r}^{4}\text{d}{r}^{2}}{1-\alpha {\left({r}^{3}+{\alpha }^{3}\right)}^{-1/3}}+{\left({r}^{3}+{\alpha }^{3}\right)}^{2/3}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\left[1-\alpha {\left({r}^{3}+{\alpha }^{3}\right)}^{-1/3}\right]{с}^{2}\text{d}{t}^{2}\end{array}$ (43)

As can be proved, the first term of this equation becomes zero at $r=0$ which cannot be true. The origin of the mistake is in Equation (40) from which it follows that $\beta =\alpha$.

To demonstrate the alternative approach, substitute Equations (39) in Equation (38) and return to spherical coordinates with the aid of Equations (37). The resulting equation is

$\begin{array}{c}\text{d}{s}^{2}=\frac{{\left({r}^{3}+{\beta }^{3}\right)}^{-4/3}{r}^{4}d{r}^{2}}{1-\alpha {\left({r}^{3}+{\beta }^{3}\right)}^{-1/3}}+{\left({r}^{3}+{\beta }^{3}\right)}^{2/3}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\left[1-\alpha {\left({r}^{3}+{\beta }^{3}\right)}^{-1/3}\right]{с}^{2}\text{d}{t}^{2}\end{array}$ (44)

As can be seen, the first term has the proper behavior at $r=0$ if we take $\beta =0$ (not $\beta =\alpha$ as in the Schwarzschild solution). The resulting expression

$\text{d}{s}^{2}=\frac{\text{d}{r}^{2}}{1-\alpha /r}+{r}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-\left(1-\frac{\alpha }{r}\right){с}^{2}\text{d}{t}^{2}$ (45)

completely corresponds to the modern interpretation of the Schwarzschild solution in Equation (31), if we apply the asymptotic Equations (13) which give $\alpha ={r}_{g}$.

Consider the solution of the internal problem that was published by K. Schwarzschild in 1916  . This solution was not supported by A. Einstein  because the concept of an incompressible fluid involves infinitely high velocity of the wave in the fluid which does not correspond to the basic GRT concept. However, the solution for compressible fluids does not demonstrate qualitative deviation from the Schwarzschild solution  which is discussed below

The method of the solution is the same that for the external problem, i.e., the new variables in Equations (37) are introduced and the field equations are supplemented with Equation (36). Further, another variable $\chi$ is introduced in accordance with the following equation:

$\mathrm{sin}\chi =\lambda {\left({r}^{3}+{\eta }^{3}\right)}^{1/3}$ (46)

where $\lambda$ is specified by Equation (18) (in the original equation $с=1$ ) and ${\eta }^{3}$ is the integration constant. The final original result for the space part of the metric form is

$\text{d}{s}_{i}^{2}=\frac{1}{\lambda }\left[\text{d}{\chi }^{2}+{\mathrm{sin}}^{2}\chi \left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\right]-{\left(\frac{3}{2}\mathrm{cos}{\chi }_{a}-\frac{1}{2}\mathrm{cos}\chi \right)}^{2}{с}^{2}\text{d}{t}^{2}$ (47)

where ${\chi }_{a}=\chi \left(r=a\right)$.

To discuss the result obtained by K. Schwarzschild, change $\chi$ to r using Equation (46). Then, Equation (47) becomes

$\begin{array}{c}\text{d}{s}_{i}^{2}=\frac{{\left({r}^{3}+{\eta }^{3}\right)}^{-4/3}{r}^{4}\text{d}{r}^{2}}{1-{\lambda }^{2}{\left({r}^{3}+{\eta }^{3}\right)}^{2/3}}+{\left({r}^{3}+{\eta }^{3}\right)}^{2/3}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{4}{\left[3\sqrt{1-{\lambda }^{2}{\left({a}^{3}+{\eta }^{3}\right)}^{2/3}}-\sqrt{1-{\lambda }^{2}{\left({r}^{3}+{\eta }^{3}\right)}^{2/3}}\right]}^{2}{с}^{2}\text{d}{t}^{2}\end{array}$

The first coefficient becomes zero at $r=0$, which cannot be true. To obtain the realistic metric, we must take $\eta =0$ and arrive at the expression

$\text{d}{s}_{i}^{2}=\frac{\text{d}{r}^{2}}{1-{\lambda }^{2}{r}^{2}}+{r}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-\frac{1}{4}{\left(3\sqrt{1-{\lambda }^{2}{a}^{2}}-\sqrt{1-{\lambda }^{2}{r}^{2}}\right)}^{2}{с}^{2}\text{d}{t}^{2}$ (48)

In the closure of his paper, K. Schwarzschild analyzed the obtained solution. Particularly, the sphere mass was found in the form

$m=\frac{6\text{π}}{\kappa \lambda }\left({\chi }_{a}-\frac{1}{2}\mathrm{sin}2{\chi }_{a}\right)$ (49)

which coincides with Equation (35) after transformation with the aid of Equation (46). The first term in Equation (48) becomes singular at some critical sphere radius ${a}_{g}=1/\lambda$. This radius coincides with the gravitational radius ${r}_{g}$ only under the condition imposed by Equation (33). If this equation is valid, the metric form in Equation (48) coincides with Equation (32). But then, the sphere mass is specified by Equation (34) which corresponds to the Euclidean space. However, the mass found by K. Schwarzschild is given by Equation (49) and corresponds to the Riemannian space. Thus, the radius ${r}_{g}$ cannot be called the Schwarzschild radius.

Now return to Equation (44) which specifies the Schwarzschild solution for the external space. Taking $r\to \infty$, and performing asymptotic analysis, we can prove that the metric coefficients in Equation (44) reduce to Equations (13) if we take $\alpha ={r}_{g}$. Then Equation (44) becomes

$\begin{array}{c}\text{d}{s}^{2}=\frac{{\left({r}^{3}+{\beta }^{3}\right)}^{-4/3}{r}^{4}\text{d}{r}^{2}}{1-{r}_{g}{\left({r}^{3}+{\beta }^{3}\right)}^{-1/3}}+{\left({r}^{3}+{\beta }^{3}\right)}^{2/3}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\left[1-{r}_{g}{\left({r}^{3}+{\beta }^{3}\right)}^{-1/3}\right]{с}^{2}\text{d}{t}^{2}\end{array}$ (50)

In contrast to the traditional Equation (31), this equation contains one more integration constant― $\beta$. This result looks natural because K. Schwarzschild did not use Equations (30) and, hence, did not reduce the order of Equation (4). Considering the space with a point mass, we took $\beta =0$ because of the behavior of the first coefficient in Equation (50) at $r=0$ and reduced Equation (44) to Equation (45). But now we study the external space of a fluid sphere for which $r\ge a$. So, we can try to use this constant to obtain the continuous solution in Equations (48) and (50) at the sphere surface. However, matching equations (48) and (50), we can conclude that the second terms can be continuous only if $\beta =0$. Thus, the final form of Equation (50) for the external space of a fluid sphere

$\text{d}{s}^{2}=\frac{\text{d}{r}^{2}}{1-{r}_{g}/r}+{r}^{2}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-\left(1-\frac{{r}_{g}}{r}\right){с}^{2}\text{d}{t}^{2}$ (51)

coincides with the traditional Equation (31). The first and the third coefficients of Equations (48) and (51) are continuous at $r=a$ if ${r}_{g}={\lambda }^{2}{a}^{3}$ which coincides with Equation (33). Thus, Equation (36) applied by K. Schwarzschild is not the proper coordinate condition and actually gives the same results that the conditions in Equations (30). The original Schwarzschild solution, as well as it modern interpretation, does not provide the solution which satisfies GRT equations and all asymptotic and boundary conditions for a fluid sphere.

4. New Model of Space and Spherically Symmetric Problem

Traditionally GRT is associated with Riemannian geometry which describes the so-called curved space. A three-dimensional curved space can be hardly imagined. This space can be formally embedded into traditional Euclidean space. However such space has six dimensions which do not provide better understanding of the problem. The proposed interpretation of the Riemannian space is based on the following assumptions. First, we assume that the space is not an object of geometry, but is a material substance (ether, physical vacuum or whatever else). Second, we think that the curved space does not exist in reality and the Riemannian geometry is only a mathematical model of a special Euclidean space. This space is not homogeneous and is characterized with so-called space density that is a function of the coordinates to which the space is referred. The space density $d=\text{d}{v}_{R}/\text{d}{v}_{E}$ is the ratio of the three-dimensional volume elements corresponding to the Riemannian and to the Euclidean geometries in the same coordinate frame   . For spherical coordinates and the metric form in Equation (1), we have $\text{d}{v}_{R}=g{\rho }^{2}\mathrm{sin}\theta \text{d}r\text{d}\theta \text{d}\phi$ and $\text{d}{v}_{E}={r}^{2}\mathrm{sin}\theta \text{d}r\text{d}\theta \text{d}\phi$, so that $d=g{\left(\rho /r\right)}^{2}$. Using Equations (14) and (19), we get for the external and internal spaces

${d}_{e}=\frac{\rho {}_{e}{}^{2}g{}_{e}}{{r}^{2}}=\frac{{{\rho }^{\prime }}_{e}{\rho }_{e}^{2}}{{r}^{2}\sqrt{1-{r}_{g}/{\rho }_{e}}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{d}_{i}=\frac{{\rho }_{i}^{2}{g}_{i}}{{r}^{2}}=\frac{{{\rho }^{\prime }}_{i}{\rho }_{i}^{2}}{{r}^{2}\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}}$ (52)

The space densities in Equations (52) are characterized with some specific properties. Consider the total density of the external and internal spaces for the fluid sphere with radius a

${D}_{e}=4\text{π}\underset{a}{\overset{\infty }{\int }}{d}_{e}{r}^{2}\text{d}r=4\text{π}\underset{a}{\overset{\infty }{\int }}F\left({\rho }_{e},{{\rho }^{\prime }}_{e}\right)\text{d}r,\text{\hspace{0.17em}}\text{\hspace{0.17em}}F\left({\rho }_{e},{{\rho }^{\prime }}_{e}\right)=\frac{{{\rho }^{\prime }}_{e}{\rho }_{e}^{2}}{\sqrt{1-{r}_{g}/{\rho }_{e}}}$ (53)

${D}_{i}=4\text{π}\underset{0}{\overset{a}{\int }}{d}_{i}{r}^{2}\text{d}r=4\text{π}\underset{0}{\overset{a}{\int }}F\left({\rho }_{i},{{\rho }^{\prime }}_{i}\right)\text{d}r,\text{\hspace{0.17em}}\text{\hspace{0.17em}}F\left({\rho }_{i},{{\rho }^{\prime }}_{i}\right)=\frac{{{\rho }^{\prime }}_{i}{\rho }_{i}^{2}}{\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}}$ (54)

Consider the variational equations providing the minimum values of the functional in Equations (53) and (54), i.e.,

$\frac{\partial F}{\partial \rho }-\frac{\partial }{\partial r}\frac{\partial F}{\partial {\rho }^{\prime }}=0$ (55)

As can be readily proved, Equation (55) is satisfied identically for both functions F in Equations (53) and (54). Thus, the space densities in Equations (52) provide the minimum total density of the space. However, the space density is caused by gravitation and is minimum in the absence of gravitation, i.e., if the space geometry is Euclidean or if $d=1$ which means that the space tends to become homogeneous with respect to the space density d. The condition $d=1$ looks slightly similar to Equation (36) applied by K. Schwarzschild. It has a simple physical and geometrical meaning―gravitation, changing the space geometry, does not affect the space volume.

Taking ${d}_{e}=1$ and ${d}_{i}=1$ in Equations (52), we arrive at two differential equations for functions ${\rho }_{e}\left(r\right)$ and ${\rho }_{i}\left(r\right)$, i.e.,

${{\rho }^{\prime }}_{e}{\rho }_{e}^{2}={r}^{2}\sqrt{1-{r}_{g}/{\rho }_{e}}$, ${{\rho }^{\prime }}_{i}{\rho }_{i}^{2}={r}^{2}\sqrt{1-{\lambda }^{2}{\rho }_{i}^{2}}$ (56)

Consider the second equation. Since ${d}_{i}={g}_{i}{\rho }_{i}^{2}/{r}^{2}=1$, the sphere mass becomes

$m=4\text{π}{\mu }_{0}\underset{0}{\overset{a}{\int }}{g}_{i}{\rho }_{i}^{2}\text{d}r=4\text{π}{\mu }_{0}\underset{0}{\overset{a}{\int }}{r}^{2}\text{d}r=\frac{4}{3}\text{π}{\mu }_{0}{a}^{3}$

and coincides with the Euclidean mass in Equation (34) which means that the condition in Equation (33) is valid and ${\lambda }^{2}={r}_{g}/{a}^{3}$. Then, the second equation in Equations (56) reduces to

${{\rho }^{\prime }}_{i}{\rho }_{i}^{2}={r}^{2}\sqrt{1-{r}_{g}{\rho }_{i}^{2}/{a}^{3}}$ (57)

The solution of Equation (57) which satisfies the boundary condition ${\rho }_{i}\left(r=0\right)=0$ is 

$\frac{1}{\sqrt{{\stackrel{¯}{r}}_{g}}}{\mathrm{sin}}^{-1}\left({\stackrel{¯}{\rho }}_{i}\sqrt{{\stackrel{¯}{r}}_{g}}\right)-{\stackrel{¯}{\rho }}_{i}\sqrt{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}}=\frac{2}{3}{\stackrel{¯}{r}}_{g}{\stackrel{¯}{r}}^{3}$ (58)

where

$\stackrel{¯}{r}=r/a,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{\rho }=\rho /a$ (59)

Recall that at the sphere surface ${\rho }_{i}\left(r=a\right)={\rho }_{a}$. Taking $\stackrel{¯}{r}=1$ and ${\stackrel{¯}{\rho }}_{i}={\stackrel{¯}{\rho }}_{a}$ in Equation (58), we get

$\frac{1}{\sqrt{{\stackrel{¯}{r}}_{g}}}{\mathrm{sin}}^{-1}\left({\stackrel{¯}{\rho }}_{a}\sqrt{{\stackrel{¯}{r}}_{g}}\right)-{\stackrel{¯}{\rho }}_{a}\sqrt{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{a}^{2}}=\frac{2}{3}{\stackrel{¯}{r}}_{g}$ (60)

The general solution of the first equation in Equations (56) is 

$\left(\frac{1}{3}{\stackrel{¯}{\rho }}_{e}^{2}+\frac{5}{12}{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{e}+\frac{5}{8}{\stackrel{¯}{r}}_{g}^{2}\right)\sqrt{{\stackrel{¯}{\rho }}_{e}\left({\stackrel{¯}{\rho }}_{e}-{\stackrel{¯}{r}}_{g}\right)}+\frac{5}{8}{\stackrel{¯}{r}}_{g}^{3}\mathrm{ln}\left(\sqrt{{\stackrel{¯}{\rho }}_{e}}+\sqrt{{\stackrel{¯}{\rho }}_{e}-{\stackrel{¯}{r}}_{g}}\right)=\frac{1}{3}{\stackrel{¯}{r}}^{3}+{C}_{5}$ (61)

The integration constant can be found from the boundary condition on the so here surface according to which ${\stackrel{¯}{\rho }}_{e}\left(\stackrel{¯}{r}=1\right)={\stackrel{¯}{\rho }}_{a}$. Then,

${C}_{5}=\left(\frac{1}{3}{\stackrel{¯}{\rho }}_{a}^{2}+\frac{5}{12}{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{a}+\frac{5}{8}{\stackrel{¯}{r}}_{g}^{2}\right)\sqrt{{\stackrel{¯}{\rho }}_{a}\left({\stackrel{¯}{\rho }}_{a}-{\stackrel{¯}{r}}_{g}\right)}+\frac{5}{8}{\stackrel{¯}{r}}_{g}^{3}\mathrm{ln}\left(\sqrt{{\stackrel{¯}{\rho }}_{a}}+\sqrt{{\stackrel{¯}{\rho }}_{a}-{\stackrel{¯}{r}}_{g}}\right)-\frac{1}{3}$ (62)

Thus, the functions ${\rho }_{i}\left(r\right)$ and ${\rho }_{e}\left(r\right)$ are specified by Equations (58) and (61). The first boundary condition in Equations (28) according to which ${\rho }_{i}\left(a\right)={\rho }_{e}\left(a\right)={\rho }_{a}$ is satisfied which follows from Equations (60) and (62). The second boundary condition in Equations (28) according to which ${g}_{i}\left(a\right)={g}_{e}\left(a\right)$ is satisfied because the obtained solution is based on the condition $g{}_{i}\rho {}_{i}^{2}={g}_{e}{\rho }_{e}^{2}$ from which it follows that if the function $\rho \left(r\right)$ is continuous at $r=a$, the function $g\left(r\right)$ is also continuous. For low levels of gravitation, i.e. for ${\stackrel{¯}{r}}_{g}\ll 1$, Equations (58) and (61) yield ${\rho }_{i}={\rho }_{e}=r$.

Consider Equations (60) and (62) which allow us to satisfy the boundary conditions, i.e., to solve the problem that turned out to be critical for the solutions discussed above. As follows from Equation (62), the solution exists if ${\stackrel{¯}{\rho }}_{a}\ge {\stackrel{¯}{r}}_{g}$. Otherwise, the solution becomes imaginary. The minimum possible value of ${\stackrel{¯}{\rho }}_{a}$ is ${\stackrel{¯}{r}}_{g}$. Assume that this minimum value corresponds to the sphere radius ${a}_{g}$. Then, substituting ${\stackrel{¯}{\rho }}_{a}={\rho }_{a}/{a}_{g}={r}_{g}/{a}_{g}$ in Equation (60), we get

$\sqrt{\frac{{a}_{g}}{{r}_{g}}}{\mathrm{sin}}^{-1}\left(\frac{r{}_{g}}{{a}_{g}}\sqrt{\frac{{r}_{g}}{{a}_{g}}}\right)-\frac{{r}_{g}}{{a}_{g}}\sqrt{1-{\left(\frac{{r}_{g}}{{a}_{g}}\right)}^{3}}=\frac{2{r}_{g}}{3{a}_{g}}$ (63)

The solution of Equation (63) is ${a}_{g}=1.115{r}_{g}$. Thus, the critical radius is larger than the gravitational radius. For $a={a}_{g}$, the solution is not singular and gives finite values for the metric coefficients. Particularly, for $\rho ={\rho }_{a}$ we get ${g}_{e}={g}_{i}=1.243$ and ${\rho }_{e}={\rho }_{i}=0.8968a$. For $a<{a}_{g}$, the solution becomes imaginary which means that GRT is not valid for such high level of gravitation. For the sphere with radius ${a}_{g}$, the escape velocity is equal to the velocity of light and such sphere is invisible  . More results concerning the discussed solution can be found elsewhere  .

5. Conclusion

As follows from the foregoing analysis, the Schwarzschild solution after some minor correction and reconstruction coincides with the traditional  interpretation of this solution. Both solutions do not satisfy the boundary condition on the fluid sphere surface for the radial space metric coefficient. A proposed model of the Riemannian space as the Euclidean space of variable density allows us to obtain the solution which satisfies equations GRT and all boundary conditions for the spherically symmetric problem for a fluid sphere. In future, the authors plan to generalize the approach discussed in Section 4 to the axially symmetric problem of GRT.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

  Schwarzschild, K. (1916) Sitz. Preuss. Akad. Wiss., 1916, 189-207.  Schwarzschild, K. (1916) Sitz. Preuss. Akad. Wiss., 1916, 424-432.  Synge, J.L. (1960) Relativity: the General Theory. Amsterdam, North Holland.  Weinberg, S. (1972) Gravitation and Cosmology. Amsterdam, North Holland.  Misner, S.W., Thorne, K.S. and Wheeler, G.A. (1973) Gravitation. W.H. Freeman and Co., San Francisco.  Chandrasecar, S. (1983) The Mathematical Theory of Black Holes. Oxford University Press, Oxford.  Frolov, V.P. and Zelnikov, A. (2011) Introduction to Black Hole Physics. Oxford University Press, Oxford. https://doi.org/10.1093/acprof:oso/9780199692293.001.0001  Logunov, A.A., Mestvirishvili, M.A. and Petrov, V.A. (2004) Physics-Uspekhi, 47, 607-621. (In Russian) https://doi.org/10.1070/PU2004v047n06ABEH001817  Fock, V. (1959) The Theory of Space, Time and Gravitation. Pergamon Press, London.  Landau, L.D. and Lifshits, E.M. (1988) Field Theory. Nauka, Moscow. (In Russian)  Vasiliev, V.V. and Fedorov, L.V. (2017) Applied Physics Research, 6, 8-13. https://doi.org/10.5539/apr.v9n6p8  Vasiliev, V.V. (2017) Journal of Modern Physics, 8, 1087-1100. https://doi.org/10.4236/jmp.2017.87070  Einstein, A. (1916) Annalen der Physik, 49, 769-819. https://doi.org/10.1002/andp.19163540702  Einstein, A. (1939) Annals of Mathematics, 40, 922-936. https://doi.org/10.2307/1968902  Vasiliev, V.V. and Fedorov, L.V. (2012) Applied Physics Research, 2, 166-174.  Vasiliev, V.V. (1989) Mechanics of Solids, 5, 30-34.  Vasiliev, V.V. and Fedorov, L.V. (2013) Applied Physics Research, 1, 115-125.  Vasiliev, V.V. and Fedorov, L.V. (2015) Applied Physics Research, 3, 84-92. 