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:
Here and t are space spherical and time coordinates, 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 with the metric coefficients as 
where . The components of the Einstein tensor are proportional to the energy tensor , i.e.
is the GRT gravitational constant depending on the Newton constant 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:
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
and Equations (2)-(4) are homogeneous. Inside a solid sphere with radius ,
where and are the radial and the circumferential stresses and is the material density.
Consider the external space . Taking , we can reduce Equation (4) to the following form:
Equation (10) can be readily integrated to give
where is the integration constant and the functions with subscript “e” correspond to the external space. Substituting this result in Equation (2), taking and integrating, we get
It looks like substituting Equations (11) and (12) into Equation (3) in which , we can determine the function . But this is not the case―under such substitution, Equation (3) is satisfied identically for any function . This result can be predicted―since the components of the Einstein tensor 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 , and . 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 , we can make some qualitative conclusions about its behavior. It is natural to suppose the for we should have and Equations (11) and (12) should reduce to the corresponding results of the Newton gravitation theory, according to which 
where 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
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 . 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 and the stresses and . 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, in which p is the pressure in the fluid and the fluid density does not depend on p and r. Taking in Equation (4), we can reduce it to the form 
Equation (16) can be readily integrated to give
and the functions with subscript “i” correspond to the internal space. We do not know the function , but can make some reasonable prediction concerning its behavior. Particularly, it is natural to suppose that at the sphere center . Then, in accordance with the symmetry condition at the sphere center, we should take and Equation (17) becomes
Continue the derivation and consider Equations (2) and (7). Taking and , , we arrive at
Consider the first Equation. Substituting Equation (19), express the time metric coefficient , i.e.,
and rewrite the second equation of Equations (20) as
Substituting Equation (21) and using Equation (18), we arrive at the following equation for the pressure:
The solution of Equation (23) that satisfies the boundary condition on the free sphere surface according to which is
where . To determine , integrate Equation (22) to get in which is the integration constant that can be found from the boundary condition on the sphere surface according to which . Using the second equation in Equations (14), we finally have
It should be noted that substitution of the obtained solutions and in the Einstein Equation (3) where , does not allow us to find the function because Equation (3) is identically satisfied for any function . 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
To fulfill the solution, we need to find two functions and such that allow us to satisfy the boundary conditions at the sphere surface , i.e.,
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
for ten components of the metric tensor . 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 . 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
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
To fulfill the solution, we need to satisfy the boundary conditions (28) on the sphere surface . The first condition is satisfied automatically, whereas the second one yields
However, the parameters and 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:
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
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:
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 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 . To overcome this problem, K. Schwarzschild introduced new variables such that
In new coordinates, the volume element becomes and the line element takes the form
Three functions can be found from Equations (2) and (4) supplemented with equation which follows from Equation (36). The final solution is
, , (39)
in which and 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 which can be a source of singularity. Equating the denominator to zero and using Equation (37) for we get
Let singularity take place at the origin . Then, Equation (40) yields . Introducing the new variable
and using Equations (38) and (39), K. Schwarzschild arrived at the following final result of his paper:
This form formally coincides with Equation (31), but it should be taken into account that R is not the radial coordinate r. The constant 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 or . Thus, the original Schwarzschild solution has only one singular point― .
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
As can be proved, the first term of this equation becomes zero at which cannot be true. The origin of the mistake is in Equation (40) from which it follows that .
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
As can be seen, the first term has the proper behavior at if we take (not as in the Schwarzschild solution). The resulting expression
completely corresponds to the modern interpretation of the Schwarzschild solution in Equation (31), if we apply the asymptotic Equations (13) which give .
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 is introduced in accordance with the following equation:
where is specified by Equation (18) (in the original equation ) and is the integration constant. The final original result for the space part of the metric form is
To discuss the result obtained by K. Schwarzschild, change to r using Equation (46). Then, Equation (47) becomes
The first coefficient becomes zero at , which cannot be true. To obtain the realistic metric, we must take and arrive at the expression
In the closure of his paper, K. Schwarzschild analyzed the obtained solution. Particularly, the sphere mass was found in the form
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 . This radius coincides with the gravitational radius 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 cannot be called the Schwarzschild radius.
Now return to Equation (44) which specifies the Schwarzschild solution for the external space. Taking , and performing asymptotic analysis, we can prove that the metric coefficients in Equation (44) reduce to Equations (13) if we take . Then Equation (44) becomes
In contrast to the traditional Equation (31), this equation contains one more integration constant― . 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 because of the behavior of the first coefficient in Equation (50) at and reduced Equation (44) to Equation (45). But now we study the external space of a fluid sphere for which . 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 . Thus, the final form of Equation (50) for the external space of a fluid sphere
coincides with the traditional Equation (31). The first and the third coefficients of Equations (48) and (51) are continuous at if 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 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 and , so that . Using Equations (14) and (19), we get for the external and internal spaces
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
Consider the variational equations providing the minimum values of the functional in Equations (53) and (54), i.e.,
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 which means that the space tends to become homogeneous with respect to the space density d. The condition 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 and in Equations (52), we arrive at two differential equations for functions and , i.e.,
Consider the second equation. Since , the sphere mass becomes
and coincides with the Euclidean mass in Equation (34) which means that the condition in Equation (33) is valid and . Then, the second equation in Equations (56) reduces to
The solution of Equation (57) which satisfies the boundary condition is 
Recall that at the sphere surface . Taking and in Equation (58), we get
The general solution of the first equation in Equations (56) is 
The integration constant can be found from the boundary condition on the so here surface according to which . Then,
Thus, the functions and are specified by Equations (58) and (61). The first boundary condition in Equations (28) according to which is satisfied which follows from Equations (60) and (62). The second boundary condition in Equations (28) according to which is satisfied because the obtained solution is based on the condition from which it follows that if the function is continuous at , the function is also continuous. For low levels of gravitation, i.e. for , Equations (58) and (61) yield .
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 . Otherwise, the solution becomes imaginary. The minimum possible value of is . Assume that this minimum value corresponds to the sphere radius . Then, substituting in Equation (60), we get
The solution of Equation (63) is . Thus, the critical radius is larger than the gravitational radius. For , the solution is not singular and gives finite values for the metric coefficients. Particularly, for we get and . For , the solution becomes imaginary which means that GRT is not valid for such high level of gravitation. For the sphere with radius , 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  .
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.