To the Complete Set of Equations for a Static Problem of General Relativity

The paper is concerned with the formulation of the static problem of general relativity. As known, this problem is reduced to ten equations for the compo-nents of the Einstein tensor and the solution of these equations is associated with two principal problems. First, since the components of the Einstein tensor identically satisfy four conservation equations, only six of these equations are mutually independent. So, the set of the Einstein equations actually contains six independent equations for ten components of the metric tensor and should be supplemented with four additional equations which are missing in the original theory. Second, for a deformable solid the Einstein tensor is associated with the energy tensor which is expressed in terms of six stresses induced by gravitation. These stresses are not known and the relativity theory does not propose any equations for them. Thus, the static problem of general relativity cannot be properly formulated because the set of governing equations is not complete. In the paper, the problem of completeness of the general relativity governing set of equations is analyzed in application to the spherically symmetric static problem and the proposed approach is further described for the general case. As an example, linearized axisymmetric problem is considered.

Share and Cite:

Vasiliev, V. and Fedorov, L. (2019) To the Complete Set of Equations for a Static Problem of General Relativity. Journal of Modern Physics, 10, 1401-1415. doi: 10.4236/jmp.2019.1012093.

1. Introduction. General Relativity Equations

The Einstein equation which specifies the Einstein tensor has the following form:

${E}_{i}^{j}={R}_{i}^{j}-\frac{1}{2}{\delta }_{i}^{j}R$, $\left(i,j=1,2,3,4\right)$ (1)

in which ${R}_{i}^{j}$ ( $R={R}_{i}^{i}$ ) are the components of the Ricci curvature tensor (we use mixed components because for the spherically symmetric problem considered further they coincide with the physical components). The Einstein tensor is associated with the energy tensor as

${E}_{i}^{j}=\chi {T}_{i}^{j}$ (2)

where

$\chi =8\pi G/{c}^{4}$ (3)

is the relativity gravitational constant expressed in terms of the classical constant G and the velocity of light c. The energy tensor expressed with the aid of Equation (1) and Equation (2) identically satisfies the conservation equations

${T}_{i,k}^{k}=0$ (4)

For the static problem,

${T}_{i}^{j}={\sigma }_{i}^{j}$, ${T}_{i}^{4}=0$ $\left(i,j=1,2,3\right)$, ${T}_{4}^{4}=\mu {c}^{2}$ (5)

where ${\sigma }_{i}^{j}$ is the stress tensor and $\mu$ is the density.

Consider two problems associated with the formulation of the general relativity static problem. First, substituting Equations (2) in Equation (1), we arrive at ten equations for ten components of the metric tensor in four-dimensional Riemannian space. However, the right parts of these equations identically satisfy Equations (4) which means that only six of ten Equation (1) are mutually independent. Thus, we have six equations for ten unknown functions. The additional equations which are usually referred to as coordinate conditions should be imposed on the metric tensor. As known, the metric tensor of the Euclidean space must satisfy the Lame equations. Such equations do not exist for the Riemannian space. However, we can suppose that the Riemannian space induced by gravitation cannot be arbitrary and must be somehow restricted, e.g., by coordinate conditions. The necessity of such conditions was first mentioned by D. Hilbert . By now the widely recognized general form of these conditions has not been proposed. Existing particular conditions are discussed further in application to the spherically symmetric problem.

Second, changing E to T and then to $\sigma$ in Equations (1), we arrive at the set of equations containing stresses on the left sides. These stresses are not known. Traditionally, Equation (1) is used to study the gravitation in the empty space for which ${T}_{i}^{j}=0$ and Equation (1) are homogeneous. For solids the stresses are not zero and the equations cannot be solved. The solution can be obtained if the solid is simulated with perfect fluid. In this case the nonzero stresses ${\sigma }_{1}^{1}={\sigma }_{2}^{2}={\sigma }_{3}^{3}=-p$, where p is the pressure that can be found from the corresponding equation of Equation (4). However, in equilibrium state, the perfect fluid can have only spherical shape and the solution can be obtained for this particular case only. For the general case, consider the analogy with the theory of elasticity. Conservation equations, Equation (4), correspond to equilibrium equations of this theory, whereas Equation (1) is similar to the equations which allow us to express the stresses in terms of stress functions and to satisfy the equilibrium equations. Thus, the metric tensor of general relativity is analogous to the tensor of stress functions in theory of elasticity. In this theory, the stress functions are found from the compatibility equations which postulate that the geometry of the stressed solid is Euclidean. In general relativity, the geometry is Riemannian and the compatibility equations of the theory of elasticity cannot be directly applied.

Thus, the traditional set of the general relativity equations is not complete and should be supplemented with some additional equations that are discussed further. The proposed approach is demonstrated in application to the spherically symmetric problem which has the exact solution.

2. Spherically Symmetric Static Problem

2.1. Classical Linear Solution

For comparison with the general relativity solutions that are discussed further, consider the problem of the theory of elasticity for a linear elastic isotropic solid sphere loaded with gravitation forces following from the Newton theory. For a sphere with constant density $\mu$, the gravitational potential $\phi$ is the solution of the Poisson equation

$\Delta \phi ={\phi }^{″}+\frac{2}{r}{\phi }^{\prime }=4\pi G\mu$ (6)

Here, ${\left(\cdot \cdot \cdot \right)}^{\prime }=\text{d}\left(\cdot \cdot \cdot \right)/\text{d}r$ and r is the radial coordinate ( $0\le r\le R$ ). For the external space ( $r\ge R$, index “e”), $\mu =0$ and the solution of Equation (6) is ${\phi }_{e}=-Gm/r$ in which

$m=\frac{4}{3}\pi {R}^{3}\mu$ (7)

is the sphere mass. Introduce the so-called gravitational radius

${r}_{g}=\frac{2Gm}{{c}^{2}}$ (8)

Then, ${\phi }_{e}=-{r}_{g}{c}^{2}/2r$. For the internal space ( $0\le r\le R$, index “i”), the regular solution of Equation (6) is

${\phi }_{i}=\frac{2}{3}\pi G\mu {r}^{2}+C$

Determining constant C from the boundary condition ${\phi }_{i}\left(R\right)={\phi }_{e}\left(R\right)$ and using Equation (8), we get

${\phi }_{i}=-\frac{{r}_{g}{c}^{2}}{4R}\left(3-\frac{{r}^{2}}{{R}^{2}}\right)$ (9)

The equilibrium equation for the sphere under the action of gravitational body forces ${f}_{g}=-\mu {{\phi }^{\prime }}_{i}$ is

${{\sigma }^{\prime }}_{r}+\frac{2}{r}\left({\sigma }_{r}-{\sigma }_{\theta }\right)-\frac{\mu {r}_{g}{c}^{2}}{2{R}^{3}}r=0$ (10)

where ${\sigma }_{r}$ and ${\sigma }_{\theta }$ are the radial and the circumferential stresses. The second equation for the stresses follows from the compatibility equation. There are two ways to derive this equation. First, introduce the corresponding strains expressed in terms of the radial displacement u as

${\epsilon }_{r}={u}^{\prime }$, ${\epsilon }_{\theta }=u/r$ (11)

The compatibility equation follows from these equations and has the form

${\left(r{\epsilon }_{\theta }\right)}^{\prime }={\epsilon }_{r}$ (12)

Express the strains in terms of stresses with the aid of Hooke’s law

${\epsilon }_{r}=\frac{1}{E}\left({\sigma }_{r}-2\nu {\sigma }_{\theta }\right)$, ${\epsilon }_{\theta }=\frac{1}{E}\left[\left(1-\nu \right){\sigma }_{\theta }-\nu {\sigma }_{r}\right]$ (13)

in which E is the elastic modulus and $\nu$ is the Poisson’s ratio. Substituting the strains in Equation (12), we finally get

$r\left[\left(1-\nu \right){{\sigma }^{\prime }}_{\theta }-\nu {\sigma }_{r}\right]+\left(1+\nu \right)\left({\sigma }_{\theta }-{\sigma }_{r}\right)=0$ (14)

Equation (12) means that the geometry of the deformed sphere is Euclidean. In general relativity, the geometry is Riemannian, the displacement u and Equations (11) do not exist. However, there is the second way to obtain Equation (14) not attracting Equations (11). This approach is based on the principle of minimum of the complementary energy

$U=\frac{2\pi }{E}\underset{0}{\overset{R}{\int }}\left[{\sigma }_{r}^{2}+2\left(1-\nu \right){\sigma }_{\theta }^{2}-4\nu {\sigma }_{r}{\sigma }_{\theta }\right]{r}^{2}\text{d}r$

under the condition that the stresses satisfy the equilibrium equation, Equation (10). Introducing this equation with the aid of the Lagrange multiplier $\lambda$, construct the augmented functional

$U=\frac{2\pi }{E}\underset{0}{\overset{R}{\int }}F\text{d}r$,

$F=\left[{\sigma }_{r}^{2}+2\left(1-\nu \right){\sigma }_{\theta }^{2}-4\nu {\sigma }_{r}{\sigma }_{\theta }\right]{r}^{2}+\lambda \left(r\right)\left[{{\sigma }^{\prime }}_{r}+\frac{2}{r}\left({\sigma }_{r}-{\sigma }_{\theta }\right)-\frac{\mu {r}_{g}{c}^{2}}{2{R}^{3}}r\right]$ (15)

The Euler equations providing $\delta U=0$

$\frac{\partial F}{\partial {\sigma }_{r}}-\frac{\text{d}}{\text{d}r}\left(\frac{\partial F}{\partial {{\sigma }^{\prime }}_{r}}\right)=0$, $\frac{\partial F}{\partial {\sigma }_{\theta }}=0$ (16)

take the form

$2\left({\sigma }_{r}-2\nu {\sigma }_{\theta }\right){r}^{2}+\frac{2}{r}\lambda -{\lambda }^{\prime }=0$, $2\left(1-\nu \right){\sigma }_{\theta }-2\nu {\sigma }_{r}-\frac{\lambda }{{r}^{3}}=0$ (17)

Expressing $\lambda$ from the second equation and substituting in the first equation, we arrive at the compatibility Equation (14).

Thus, we get two equations, Equation (10) and Equation (14) for two stresses. The final solution which satisfies the boundary condition ${\sigma }_{r}\left(R\right)=0$ is

${\stackrel{¯}{\sigma }}_{r}=-k{\stackrel{¯}{r}}_{g}\left(1-{\stackrel{¯}{r}}^{2}\right)$, ${\stackrel{¯}{\sigma }}_{\theta }=-k{\stackrel{¯}{r}}_{g}\left(1+\frac{1+\nu }{3-\nu }{\stackrel{¯}{r}}^{2}\right)$, $k=\frac{3-\nu }{20\left(1-\nu \right)}$ (18)

Here,

${\stackrel{¯}{\sigma }}_{r}=\frac{{\sigma }_{r}}{\mu {c}^{2}}$, ${\stackrel{¯}{\sigma }}_{\theta }=\frac{{\sigma }_{\theta }}{\mu {c}^{2}}$, ${\stackrel{¯}{r}}_{g}=\frac{{r}_{g}}{R}$, $\stackrel{¯}{r}=\frac{r}{R}$ (19)

For a sphere of perfect fluid, ${\sigma }_{r}={\sigma }_{\theta }=-p$ and the pressure p can be found from the equilibrium Equation (10) not attracting the compatibility Equation (14). The result is

$\stackrel{¯}{p}=\frac{{\stackrel{¯}{r}}_{g}}{4}\left(1-{\stackrel{¯}{r}}^{2}\right)$ (20)

In general relativity, the space geometry is Riemannian and the line element in spherical coordinates $r,\theta ,\phi$ is taken in the form

$\text{d}{s}^{2}={g}_{11}\text{d}{r}^{2}+{g}_{22}\left(\text{d}{\theta }^{2}+{\mathrm{sin}}^{2}\theta \text{d}{\phi }^{2}\right)-{g}_{44}{c}^{2}\text{d}{t}^{2}$ (21)

The components of the metric tensor depend on the radial coordinate only. For the foregoing linear solution, these components are 

${g}_{11}=1+\frac{{\stackrel{¯}{r}}_{g}}{\stackrel{¯}{r}}$, ${g}_{22}={r}^{2}$ (22)

In case ${\stackrel{¯}{r}}_{g}=0$, the space is Euclidean and gravitation vanishes. For real objects, the ratio ${\stackrel{¯}{r}}_{g}$, as a rule, is extremely small. For example, for Earth ${\stackrel{¯}{r}}_{g}=1.4×{10}^{-6}$, for Sun ${\stackrel{¯}{r}}_{g}=4.3×{10}^{-6}$, for the largest of the observed visible stars—red supergiant UI Scutti ( $R=11.9×{10}^{11}\text{m}$, $m=64×{10}^{30}\text{kg}$  ) ${\stackrel{¯}{r}}_{g}=8×{10}^{-9}$.

2.2. General Relativity Solution

For a spherically symmetric problem, the field equations following from Equation (1), Equation (2) and Equation (5) reduce to

$\chi {\sigma }_{r}=\frac{1}{{g}_{22}}-\frac{1}{{g}_{11}}\left[\frac{1}{4}{\left(\frac{{{g}^{\prime }}_{22}}{{g}_{22}}\right)}^{2}+\frac{{{g}^{\prime }}_{22}{{g}^{\prime }}_{44}}{2{g}_{22}{g}_{44}}\right]$ (23)

$\begin{array}{c}\chi {\sigma }_{\theta }=-\frac{1}{2{g}_{11}}\left[\frac{{{g}^{″}}_{44}}{{g}_{44}}-\frac{1}{2}{\left(\frac{{{g}^{\prime }}_{44}}{{g}_{44}}\right)}^{2}+\frac{{{g}^{″}}_{22}}{{g}_{22}}-\frac{1}{2}{\left(\frac{{{g}^{\prime }}_{22}}{{g}_{22}}\right)}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{{{g}^{\prime }}_{22}}{2{g}_{22}}\left(\frac{{{g}^{\prime }}_{44}}{{g}_{44}}-\frac{{{g}^{\prime }}_{11}}{{g}_{11}}\right)-\frac{{{g}^{\prime }}_{11}{{g}^{\prime }}_{44}}{2{g}_{11}{g}_{44}}\right]\end{array}$ (24)

$\chi \mu {c}^{2}=\frac{1}{{g}_{22}}-\frac{1}{{g}_{11}}\left[\frac{{{g}^{″}}_{22}}{{g}_{22}}-\frac{1}{4}{\left(\frac{{{g}^{\prime }}_{22}}{{g}_{22}}\right)}^{2}-\frac{{{g}^{\prime }}_{11}{{g}^{\prime }}_{22}}{2{g}_{11}{g}_{22}}\right]$ (25)

The only one conservation equation, Equation (4), becomes

${{\sigma }^{\prime }}_{r}+\frac{{{g}^{\prime }}_{22}}{{g}_{22}}\left({\sigma }_{r}-{\sigma }_{\theta }\right)+\frac{{{g}^{\prime }}_{44}}{2{g}_{44}}\left({\sigma }_{r}-\mu {c}^{2}\right)=0$ (26)

The solution of the external ( $r\ge R$ ) problem must satisfy the asymptotic conditions and to reduce to Equation (22) for $r\to \infty$. The solution for the internal ( $0\le r\le R$ ) problem must satisfy the symmetry condition at the sphere center according to which ${g}_{11}\left(0\right)=1$, ${g}_{22}\left(0\right)=0$. Both solutions must satisfy the boundary conditions on the sphere surface, i.e.

${g}_{11}^{e}\left(R\right)={g}_{11}^{i}\left(r\right)$, ${g}_{22}^{e}\left(R\right)={g}_{22}^{i}\left(R\right)$, ${g}_{44}^{e}\left(r\right)={g}_{44}^{i}\left(R\right)$ (27)

As in the general case (Section 1), substitution of the left parts of Equations (23)-(25) in Equation (26) identically satisfies this equation. So, only three of four Equations (23)-(26) are mutually independent. Traditionally , the simplest set of equations including Equation (23), Equation (25) and Equation (26) is used. The obtained solution identically satisfies Equation (24).

To solve the problem, we should supplement Equation (23), Equation (25) and Equation (26) which include three components of the metric tensor and two stresses with one coordinate condition for the metric tensor and one equation for the stresses. The first coordinate condition was proposed by K. Schwarzchild  who changed the spherical coordinates to ${x}_{1}={r}^{3}/3$, ${x}_{2}=-\mathrm{cos}\theta$, ${x}_{3}=\phi$, ${x}_{4}=t$ and applied the condition $g=1$, where g is the determinant of the metric tensor components in coordinates ${x}_{i}$. This condition is equivalent to ${g}_{22}={r}^{2}$  and reduces the order of Equation (25). As a result, the solution does not contain the proper number of integration constants and the first boundary condition in Equation (27) cannot be satisfied . The internal problem was solved for a sphere of perfect fluid  and did not require the additional equation. The other way involves the application of the so-called harmonic coordinate conditions which in the general case have the following form 

$\frac{\partial }{\partial {x}^{i}}\left(\sqrt{g}{g}^{ij}\right)=0$

External spherically symmetric problem was solved with the harmonic coordinate condition by V. Fock . Internal problem and boundary conditions were not considered.

To obtain the general solution of the spherically symmetric static problem, apply the set of Equation (23), Equation (25) and Equation (26). To simplify these equations, introduce new notations for the components of the metric tensor, i.e., put ${g}_{11}={q}^{2}$, ${g}_{22}={\rho }^{2}$, ${g}_{44}={h}^{2}$. Then, Equation (23), Equation (25) and Equation (26) reduce to

$\chi {\sigma }_{r}=\frac{1}{{\rho }^{2}}-\frac{{\rho }^{\prime }}{{q}^{2}\rho }\left(\frac{{\rho }^{\prime }}{\rho }+\frac{2{h}^{\prime }}{h}\right)$, $\chi \mu {c}^{2}=\frac{1}{{\rho }^{2}}-\frac{1}{{q}^{2}}\left[{\left(\frac{{\rho }^{\prime }}{\rho }\right)}^{2}+\frac{2{\rho }^{″}}{\rho }-\frac{2{q}^{\prime }{\rho }^{\prime }}{q\rho }\right]$ (28)

${{\sigma }^{\prime }}_{r}+\frac{2{\rho }^{\prime }}{\rho }\left({\sigma }_{r}-{\sigma }_{\theta }\right)+\frac{{h}^{\prime }}{h}\left({\sigma }_{r}-\mu {c}^{2}\right)=0$ (29)

For the external space ( ${\sigma }_{r}={\sigma }_{\theta }=0,\mu =0$ ), the solution of Equation (28) which satisfies the asymptotic conditions for $r\to \infty$ is 

${q}_{e}^{2}=\frac{{\rho }_{e}{\left({{\rho }^{\prime }}_{e}\right)}^{2}}{{\rho }_{e}-{r}_{g}}$, ${h}_{e}^{2}=1-\frac{{r}_{g}}{{\rho }_{e}}$ (30)

To determine ${\rho }_{e}\left(r\right)$, we need to add a coordinate condition. Introduce a new interpretation of the Riemannian space   according to which it is a mathematical model of the actual nonhomogeneous Euclidean space characterized with the space density $d=\sqrt{{g}_{R}/{g}_{E}}$ in which ${g}_{R}$ and ${g}_{E}$ are the determinants of the metric tensors in Riemannian and Euclidean three-dimensional spaces in the same coordinates. Assume that in space coordinates ${x}^{1},{x}^{2},{x}^{3}$ the space density satisfies the following variational equation:

$\delta D=0$, $D=\iiint d\sqrt{{g}_{E}}\text{d}{x}^{1}\text{d}{x}^{2}\text{d}{x}^{3}=\iiint \sqrt{{g}_{R}}\text{d}{x}^{1}\text{d}{x}^{2}\text{d}{x}^{3}$ (31)

Equation (31), written for a four-dimensional space, is known in general relativity as a possible way to derive the field equations . However, if the variations $\delta {g}_{ij}$ are mutually independent, $\delta D\ne 0$. The situation becomes different if these variations are not independent. For spherically symmetric problem,

${d}_{e}=\frac{{q}_{e}{\rho }_{e}^{2}}{{r}^{2}}=\frac{{{\rho }^{\prime }}_{e}{\rho }_{e}^{2}}{{r}^{2}\sqrt{1-{r}_{g}/{\rho }_{e}}}$ (32)

The Euler equation similar to Equation (16) is satisfied identically which means that $\delta D=0$ for any function ${\rho }_{e}\left(r\right)$. To specify this function, we can minimize ${d}_{e}$ taking ${d}_{e}=1$. Thus, it is assumed that gravitation transforming Euclidean space into Riemannian does not change the density of the metric tensor. For ${d}_{e}=1$, Equation (32) yields the following equation for ${\rho }_{e}\left(r\right)$ :

${{\rho }^{\prime }}_{e}{\rho }_{e}^{2}={r}^{2}\sqrt{1-{r}_{g}/{\rho }_{e}}$ (33)

For the internal space ( $\mu =const$ ), the solution of the second equation in Equation (28) which satisfies the regularity condition at the sphere center is 

${q}_{i}^{2}=\frac{{\left({{\rho }^{\prime }}_{i}\right)}^{2}}{1-u{\rho }_{i}^{2}}$, $u=\frac{1}{3}\chi \mu {c}^{2}$ (34)

As for the external space, take Equation (31) as the coordinate condition in which

${d}_{i}=\frac{{q}_{i}{\rho }_{i}^{2}}{{r}^{2}}=\frac{{{\rho }^{\prime }}_{i}{\rho }_{i}^{2}}{{r}^{2}\sqrt{1-u{\rho }_{i}^{2}}}$ (35)

where ${d}_{i}$ is the space density. The Euler equation is satisfied identically and we can take ${d}_{i}=1$. This means that that the sphere mass is not affected by gravitation and is specified by Equation (7). Using this equation in conjunction with Equation (3) and Equation (8) for $\chi$ and ${r}_{g}$, we get $u={r}_{g}/{R}^{3}$. Then, Equation (34) and Equation (35) take the following final forms:

${q}_{i}^{2}=\frac{{\left({{\rho }^{\prime }}_{i}\right)}^{2}}{1-{r}_{g}{\rho }_{i}^{2}/{R}^{3}}$, ${d}_{i}=\frac{{{\rho }^{\prime }}_{i}{\rho }_{i}^{2}}{{r}^{2}\sqrt{1-{r}_{g}{\rho }_{i}^{2}/{R}^{3}}}$

Putting ${d}_{i}=1$, we arrive at the following equation for ${\rho }_{i}\left(r\right)$ :

${{\rho }^{\prime }}_{i}{\rho }_{i}^{2}={r}^{2}\sqrt{1-{r}_{g}{\rho }_{i}^{2}/{R}^{3}}$

The solution of this equation 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}$ (36)

where in addition to Equation (19) $\stackrel{¯}{\rho }=\rho /R$. For the sphere surface $\stackrel{¯}{r}=1$, ${\stackrel{¯}{\rho }}_{i}={\stackrel{¯}{\rho }}_{1}$, and Equation (36) yields

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

The general solution of Equation (33) 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$ (38)

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

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

For ${\stackrel{¯}{r}}_{g}\ll 1$, Equation (36) and Equation (38) reduce to ${\stackrel{¯}{\rho }}_{i}={\stackrel{¯}{\rho }}_{e}=\stackrel{¯}{r}$.

Thus, the functions ${\rho }_{i}\left(r\right)$ and ${\rho }_{e}\left(r\right)$ are specified by Equation (36) and Equation (38). The obtained solution satisfies the asymptotic and the boundary conditions in Equation (27) . As follows from Equation (39), the solution exists if ${\stackrel{¯}{\rho }}_{1}\ge {\stackrel{¯}{r}}_{g}$. Otherwise, the solution becomes imaginary. The minimum possible value of ${\stackrel{¯}{\rho }}_{1}$ is ${\stackrel{¯}{r}}_{g}$. Assume that this minimum value correspond to the sphere radius ${R}_{g}$. Then, substituting ${\stackrel{¯}{\rho }}_{1}={\rho }_{1}/{R}_{g}={r}_{g}/{R}_{g}$ in Equation (37), we get

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

The solution of this equation is ${R}_{g}=1.115{r}_{g}$. Thus, the obtained solution gives the critical radius which is larger than the gravitational radius. In contrast to the Schwarzchild solution, for the sphere with the critical radius ${R}_{g}$ the solution is not singular and gives finite values for the metric coefficients. Particularly, for $\rho ={\rho }_{1}$ we get ${q}_{e}={q}_{i}=1.243$ and ${\rho }_{e}={\rho }_{i}=0.8968R$. For $R<{R}_{g}$, the solution becomes imaginary which means that the general relativity is not valid for such high levels of gravitation. Dependences of the space metric coefficients on the radial coordinate for the sphere with the critical radius ${R}_{g}$ is shown in Figure 1.

As can be seen, $g\left(r\to \infty \right)=1$ and $\rho \left(r\to \infty \right)=r$ (dashed line in Figure 1).

Consider the propagation of light from the sphere surface. The trajectory of light in the equatorial ( $\theta =\pi /2$ ) plane is specified by the following equations  :

$\frac{\text{d}r}{\text{d}t}=c\frac{{h}_{e}}{{q}_{e}}\sqrt{1-{\left(\frac{s{h}_{e}}{{\rho }_{e}}\right)}^{2}}$, $\frac{\text{d}\phi }{\text{d}t}=cs{\left(\frac{{h}_{e}}{{\rho }_{e}}\right)}^{2}$ (40)

in which ${q}_{e}$ and ${h}_{e}$ are given by Equation (30), ${\rho }_{e}\left(r\right)$ is the solution of Equation (38) and s is the integration constant that can be found from the initial condition. Physical components of the velocity are

Figure 1. Dependences of the space metric coefficients on the radial coordinate for the sphere with the critical radius.

${v}_{r}=\frac{{g}_{e}}{{h}_{e}}\frac{\text{d}r}{\text{d}t}=c\sqrt{1-{\left(\frac{s{h}_{e}}{{\rho }_{e}}\right)}^{2}}$, ${v}_{\phi }=c\frac{s{h}_{e}}{{\rho }_{e}}$ (41)

so that ${v}_{r}^{2}+{v}_{\phi }^{2}={c}^{2}$. Assume that light propagates from point A on the sphere surface ( $\stackrel{¯}{\rho }={\stackrel{¯}{\rho }}_{1}$, $\stackrel{¯}{r}=1$ ) at angle $\alpha$ with respect to the radius (Figure 2).

The initial conditions are ${v}_{r}=c\mathrm{cos}\alpha$, ${v}_{\phi }=c\mathrm{sin}\alpha$ and Equations (41) yield $s={\rho }_{1}\mathrm{sin}\alpha /{h}_{1}$, where ${h}_{1}={h}_{e}\left({\rho }_{e}={\rho }_{1}\right)$. Consider the case $\alpha =0$ for which $s=0$ and ${v}_{r}=c,{v}_{\phi }=0$. This result allows us to conclude that in the radial direction light propagates for any spherical object with velocity c. However, if $\alpha \ne 0$, the situation can be different. Using Equation (30) and Equation (40), we can obtain the following equation for the light trajectory:

$\frac{\text{d}{\rho }_{e}}{\text{d}\phi }=\frac{\text{d}r}{\text{d}\phi }\frac{\text{d}{\rho }_{e}}{\text{d}r}={\rho }_{e}^{2}\sqrt{\frac{1}{{\rho }_{1}^{2}{\mathrm{sin}}^{2}\alpha }\left(1-\frac{{r}_{g}}{{\rho }_{1}}\right)-\frac{1}{{\rho }_{e}^{2}}\left(1-\frac{{r}_{g}}{{\rho }_{e}}\right)}$ (42)

Numerical integration of this equation allows us to plot ${\rho }_{e}\left(\phi \right)$. Using further Equation (38), we change ${\rho }_{e}$ to r. The resulting trajectory $\stackrel{¯}{r}\left(\phi \right)$ for $\alpha =\pi /2$ and ${\stackrel{¯}{r}}_{g}=0.5$ is shown in Figure 2 with the dashed line. As follows from Equation (42), for the sphere with the critical radius, ${\rho }_{1}={r}_{g}$ the trajectory becomes imaginary and light does not propagate.

Return to the internal problem and determine the stresses. Consider the first equation in Equation (28). Substitute ${q}_{i}$ from Equation (34) and express

$\frac{{{h}^{\prime }}_{i}}{{h}_{i}}=\frac{{\rho }_{i}{{\rho }^{\prime }}_{i}\left(u-\chi {\sigma }_{r}\right)}{2\left(1-u{\rho }_{i}^{2}\right)}$ (43)

Here, $u=\chi \mu {c}^{2}/3={r}_{g}/{R}^{3}$. Substitute Equation (43) in Equation (29) to get

${{\sigma }^{\prime }}_{r}-\frac{2{{\rho }^{\prime }}_{i}}{{\rho }_{i}}\left({\sigma }_{r}-{\sigma }_{\theta }\right)+\frac{{\rho }_{i}{{\rho }^{\prime }}_{i}}{2\left(1-u{\rho }_{i}^{2}\right)}\left(u-\chi {\sigma }_{r}\right)\left({\sigma }_{r}-\mu {c}^{2}\right)=0$ (44)

The first equation for the stresses follows from this equation if we change variable r to $\rho$, use notations (19), and Equation (3), Equation (7) for $\chi$ and ${r}_{g}$, i.e.,

Figure 2. Propagation of light from the sphere surface.

$\frac{\text{d}{\stackrel{¯}{\sigma }}_{r}}{\text{d}{\stackrel{¯}{\rho }}_{i}}+\frac{2}{{\stackrel{¯}{\rho }}_{i}}\left({\stackrel{¯}{\sigma }}_{r}-{\stackrel{¯}{\sigma }}_{\theta }\right)-\frac{{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}}{2\left(1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}\right)}\left(1-3{\stackrel{¯}{\sigma }}_{r}\right)\left(1-{\stackrel{¯}{\sigma }}_{r}\right)=0$ (45)

For a sphere of perfect fluid, ${\sigma }_{r}={\sigma }_{\theta }=-p$ and Equation (45) reduces to

$\frac{\text{d}\stackrel{¯}{p}}{\text{d}{\stackrel{¯}{\rho }}_{i}}+\frac{{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}}{2\left(1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}\right)}\left(1+3\stackrel{¯}{p}\right)\left(1+\stackrel{¯}{p}\right)=0$

The solution of this equation which satisfies the boundary condition $\stackrel{¯}{p}\left({\stackrel{¯}{\rho }}_{i}={\stackrel{¯}{\rho }}_{1}\right)=0$ is 

$\stackrel{¯}{p}=\frac{\sqrt{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}}-\sqrt{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{1}^{2}}}{\sqrt{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}}-3\sqrt{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{1}^{2}}}$ (46)

In contrast to the Schwarzchild solution, the pressure is not singular. For ${\stackrel{¯}{r}}_{g}\ll 1$, Equation (46) degenerates into Equation (20). For the fluid sphere, Equation (43), being transformed to variables $\stackrel{¯}{p},{\stackrel{¯}{\rho }}_{i},{\stackrel{¯}{r}}_{g}$, becomes

$\frac{1}{{h}_{i}}\frac{\text{d}{h}_{i}}{\text{d}{\stackrel{¯}{\rho }}_{i}}=\frac{{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}\left(1+3\stackrel{¯}{p}\right)}{2\left(1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}\right)}$

Integrating, we can find ${h}_{i}$ for the fluid sphere. The integration constant allows us to satisfy the last boundary condition in Equation (27).

To obtain the stresses, we need to supplement Equation (45) with an additional equation. To derive this equation, we minimize the functional in Equation (15), where with regard to Equation (44)

$\begin{array}{c}F=\left[{\sigma }_{r}^{2}+2\left(1-\nu \right){\sigma }_{\theta }^{2}-4\nu {\sigma }_{r}{\sigma }_{\theta }\right]{q}_{i}{\rho }_{i}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\lambda \left[{{\sigma }^{\prime }}_{r}+\frac{2{{\rho }^{\prime }}_{i}}{{\rho }_{i}}\left({\sigma }_{r}-{\sigma }_{\theta }\right)+\frac{{\rho }_{i}{{\rho }^{\prime }}_{i}\left(u-\chi {\sigma }_{r}\right)}{2\left(1-u{\rho }_{i}^{2}\right)}\left({\sigma }_{r}-\mu {c}^{2}\right)\right]\end{array}$

Using Equation (34) for ${q}_{i}$, we can present the Euler equations, Equation (16), as

$2\left({\sigma }_{r}-2\nu {\sigma }_{\theta }\right)\frac{{\rho }_{i}^{2}{{\rho }^{\prime }}_{i}}{\sqrt{1-u{\rho }_{i}^{2}}}+\lambda \frac{{{\rho }^{\prime }}_{i}}{{\rho }_{i}}\left[2+\frac{{\rho }_{i}^{2}}{2\left(1-u{\rho }_{i}^{2}\right)}\left(u-2\chi {\sigma }_{r}+\chi \mu {c}^{2}\right]-{\lambda }^{\prime }=0$

$2\left[\left(1-\nu \right){\sigma }_{\theta }-\nu {\sigma }_{r}\right]\frac{{\rho }_{i}^{2}}{\sqrt{1-u{\rho }_{i}^{2}}}-\frac{\lambda }{{\rho }_{i}}=0$

Expressing $\lambda$ from the second equation, substituting in the first equation and reducing the resulting equation to the form similar to Equation (45), we arrive at

${\stackrel{¯}{\rho }}_{i}\frac{\text{d}\stackrel{¯}{\sigma }}{\text{d}{\stackrel{¯}{\rho }}_{i}}+{\stackrel{¯}{\sigma }}_{r}-2\nu {\stackrel{¯}{\sigma }}_{\theta }+\stackrel{¯}{\sigma }\left\{2+\frac{1}{1-{\stackrel{¯}{r}}_{g}{\stackrel{¯}{\rho }}_{i}^{2}}\left[{\stackrel{¯}{r}}_{g}{\rho }_{i}^{2}\left(4-3{\stackrel{¯}{\sigma }}_{r}\right)-3\right]\right\}=0$ (47)

where $\stackrel{¯}{\sigma }=\left(1-\nu \right){\stackrel{¯}{\sigma }}_{\theta }-\nu {\stackrel{¯}{\sigma }}_{r}$. For ${r}_{g}=0$ and $\rho =r$, this equation coincides with Equation (14). Numerical integration of Equation (45) and Equation (47) under the boundary conditions ${\sigma }_{r}\left({\rho }_{i}=0\right)={\sigma }_{\theta }\left({\rho }_{i}=0\right)$ and ${\sigma }_{r}\left({\rho }_{i}={\rho }_{1}\right)=0$ allows us to obtain the dependences of stresses on ${\stackrel{¯}{\rho }}_{i}$ which can be changed to $\stackrel{¯}{r}$ with the aid of Equation (38). The dependences ${\stackrel{¯}{\sigma }}_{r}\left(\stackrel{¯}{r}\right)$ and ${\stackrel{¯}{\sigma }}_{\theta }\left(\stackrel{¯}{r}\right)$ corresponding to ${\stackrel{¯}{r}}_{g}=0.25$ and $\nu =0$ are shown in Figure 3 with solid lines. Dashed lines correspond to the linear classical solution in Equation (18).

3. The General Theory

Return to Section 1 and consider the general case. Ten Einstein’s equations

${R}^{ij}-\frac{1}{2}{g}^{ij}R=\chi {T}^{ij}$ (48)

in which the energy tensor

${T}^{ij}={\sigma }^{ij}$ $\left(i,j=1,2,3\right)$, ${T}^{i4}=0$ $\left(i=1,2,3\right)$, ${T}^{44}=\mu {c}^{2}$

satisfies the conservation equations

Figure 3. Dependences of the normalized stresses on the radial coordinate corresponding to the obtained solution (¾¾) and the classical linear solution (- - - -).

${\partial }_{k}{\sigma }^{ik}+{\Gamma }_{mk}^{i}{\sigma }^{mk}+{\Gamma }_{kn}^{n}{\sigma }^{ik}+{\Gamma }_{44}^{i}\mu {c}^{2}=0$ (49)

and includes 10 components of the metric tensor ${g}^{ij}$. Because of Equation (49), only six of Equation (48) are mutually independent and we have six equations for 16 functions, i.e., 10 coefficients ${g}^{ij}$ and six stresses ${\sigma }^{ij}$. Assume that we use six independent equations of Equation (48) to express six metric coefficients in terms of four. To derive four additional equations, we propose to use the variational equation, Equation (31), in which six metric coefficients are expressed in terms of four. Variation with respect to these coefficients allows us to write four Euler’s equations and to obtain the set of 10 equations for the metric tensor.

To derive the equations for stresses, introduce the strain energy

$U=\frac{1}{2}\iiint {\sigma }^{ij}{\epsilon }_{ij}\sqrt{{g}_{R}}\text{d}{x}^{1}\text{d}{x}^{2}\text{d}{x}^{3}$

in which ${g}_{R}$ is the determinant of the metric tensor in the Riemannian three-dimensional space. Expressing strains in terms of stresses through Hooke’s la

${\epsilon }_{mn}={c}_{mnij}{\sigma }^{ij}$

in which ${c}_{mnij}$ is the compliance tensor and introducing Equation (49) with aid of the Lagrange multipliers, construct the augmented functional

$U=\frac{1}{2}\iiint F\text{d}{x}^{1}\text{d}{x}^{2}\text{d}{x}^{3},$

$F={c}_{mnij}{\sigma }^{mn}{\sigma }^{ij}\sqrt{{g}_{R}}+{\lambda }_{i}\left({\partial }_{k}{\sigma }^{ik}+{\Gamma }_{mk}^{i}{\sigma }^{mk}+{\Gamma }_{kn}^{n}{\sigma }^{ik}+{\Gamma }_{44}^{i}\mu {c}^{2}\right)$

Minimization with respect to the stresses and $\lambda$ -multipliers yields 10 equations for six stresses and four multipliers . Thus, we have arrived at the complete set of 20 equations for 10 metric coefficients 6 stresses and 4 multipliers.

4. Linearized Axisymmetric Problem

Spherically symmetric problem discussed above requires only one coordinate condition. To demonstrate a more complicated case, consider an axisymmetric problem for which we need two conditions. Since the general problem can hardly be solved because the equations are too complicated, obtain the linearized solution for the external space. The line element in cylindrical coordinates $r,\phi ,z$ can be presented as

$\text{d}{s}^{2}=\left(1+{f}_{1}{}_{1}\right)\text{d}{r}^{2}+{r}^{2}\left(1+{f}_{22}\right)\text{d}{\phi }^{2}+\left(1+{f}_{33}\right)\text{d}{z}^{2}+{f}_{13}\text{d}r\text{d}z-\left(1+{f}_{44}\right){c}^{2}\text{d}{t}^{2}$

Assume that functions ${f}_{mn}\left(r,z\right)$ are small in comparison with unity. For the external space with zero stresses and density, Equations (1), i.e. ${E}_{11}={E}_{13}={E}_{22}={E}_{33}={E}_{44}=0$, reduce to

$\begin{array}{l}r\frac{{\partial }^{2}}{\partial {z}^{2}}\left({f}_{22}+{f}_{44}\right)+\frac{\partial }{\partial r}\left({f}_{33}+{f}_{44}\right)-2\frac{\partial {f}_{13}}{\partial z}=0\\ r\frac{{\partial }^{2}}{\partial r\partial z}\left({f}_{22}+{f}_{44}\right)-\frac{\partial }{\partial z}\left({f}_{11}-{f}_{22}\right)=0\end{array}$ (50)

$\begin{array}{l}2\frac{{\partial }^{2}{f}_{13}}{\partial r\partial z}-\frac{{\partial }^{2}}{\partial {r}^{2}}\left({f}_{33}+{f}_{44}\right)-\frac{{\partial }^{2}}{\partial {z}^{2}}\left({f}_{11}+{f}_{44}\right)=0\\ r\frac{{\partial }^{2}}{\partial {r}^{2}}\left({f}_{22}+{f}_{44}\right)-\frac{\partial }{\partial r}\left({f}_{11}-2{f}_{22}-{f}_{44}\right)=0\end{array}$ (51)

$\begin{array}{l}\frac{{\partial }^{2}}{\partial {r}^{2}}\left({f}_{22}+{f}_{33}\right)+\frac{{\partial }^{2}}{\partial {z}^{2}}\left({f}_{11}+{f}_{22}\right)-2\frac{{\partial }^{2}{f}_{13}}{\partial r\partial z}\\ -\frac{\partial }{\partial r}\left({f}_{11}-2{f}_{22}-{f}_{33}\right)-2\frac{\partial {f}_{13}}{\partial z}=0\end{array}$ (52)

For the axially symmetric problem, we have two conservation equations, so only three of five Equations (50)-(52) are mutually independent. Consider Equation (52) and subtract from it the first equation in Equation (50) and the first equation in Equation (51). The resulting equation

$\Delta {f}_{44}=\frac{{\partial }^{2}{f}_{44}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial {f}_{44}}{\partial r}+\frac{{\partial }^{2}{f}_{44}}{\partial {z}^{2}}=0$ (53)

allows us to conclude that ${f}_{44}\left(r,z\right)$ is the classical gravitational potential which can be found in terms of exponential functions with respect to z and Bessel functions with respect to r . For the external problem, the solution must satisfy the asymptotic conditions according to which ${f}_{44}\to 0$ for $r\to \infty$ and $z\to \infty$. Proceeding, express ${f}_{13}$ from the first equation in Equations (50), i.e.,

$\frac{\partial {f}_{13}}{\partial z}=\frac{1}{2}\frac{\partial }{\partial r}\left({f}_{33}+{f}_{44}\right)+\frac{r}{2}\frac{{\partial }^{2}}{\partial {z}^{2}}\left({f}_{22}+{f}_{44}\right)$ (54)

and substitute this result in the first equation in Equation (51) to get

$r\frac{{\partial }^{3}}{\partial r\partial {z}^{2}}\left({f}_{22}+{f}_{44}\right)-\frac{{\partial }^{2}}{\partial {z}^{2}}\left({f}_{11}-{f}_{22}\right)=0$

This equation can be ignored because it follows from the second equation in Equation (50). Integration of Equation (50) yields

${f}_{11}={f}_{22}+r\frac{\partial }{\partial r}\left({f}_{22}-{f}_{44}\right)+{\phi }_{1}\left(r\right)$ (55)

where ${\phi }_{1}\left(r\right)$ is the integration function. Substituting Equation (55) in the second equation of Equation (51), we get ${{\phi }^{\prime }}_{1}=0$, so ${\phi }_{1}=C$. Thus, Equation (50) and Equation (51) allow us to express ${f}_{13}$ and ${f}_{11}$ in terms of two unknown functions— ${f}_{22}$ and ${f}_{33}$. To proceed, we need to introduce two coordinate conditions. As earlier, apply Equation (31), i.e., $\delta D=0$. The linearized space density is $d=1+{f}_{11}+{f}_{22}+{f}_{33}$. Using Equation (55), we can construct the following functional:

$D=\iiint F\text{d}r\text{d}\theta \text{d}z$, $F=r\left[1+2{f}_{22}+r\frac{\partial }{\partial r}\left({f}_{22}+{f}_{44}\right)+{f}_{33}+C\right]$

The condition $\delta D=0$ is satisfied if ${f}_{33}=0$ and ${f}_{22}$ is an arbitrary function. To identify this function, we use, as earlier, the minimum condition $d=1$. Then, ${f}_{22}=-{f}_{11}$ and Equation (55) yields

$r\frac{\partial {f}_{22}}{\partial r}+2{f}_{22}=r\frac{\partial {f}_{44}}{\partial r}-C$

The solution of this equation is

${f}_{22}=-{f}_{11}=\frac{1}{{r}^{2}}\left[{f}_{44}-C\mathrm{ln}r+{\phi }_{2}\left(z\right)\right]$

Integration of Equation (54) allows us to find the last metric coefficient, i.e.,

${f}_{13}=r\frac{\partial {f}_{44}}{\partial z}+\int \frac{\partial {f}_{44}}{\partial r}\text{d}z+\frac{1}{{r}^{3}}{{\phi }^{\prime }}_{2}\left(z\right)+{\phi }_{3}\left(r\right)$

Using the asymptotic conditions, we can conclude that constant C and functions ${\phi }_{2},{\phi }_{3}$ are zero. Thus, the solution is

${f}_{11}=-{f}_{22}=-\frac{{f}_{44}}{{r}^{2}}$, ${f}_{13}=r\frac{\partial {f}_{44}}{\partial z}+\int \frac{\partial {f}_{44}}{\partial r}\text{d}z$, ${f}_{33}=0$

where ${f}_{44}$ is the solution of Equation (53). As follows from the foregoing derivation, the space density is uniform in the external space.

5. Gravitation and Space Density

The space density introduced in Section 2.2 allows us to propose the new interpretation of gravitation. As follows from the foregoing discussion, the isolated object in space can be in equilibrium under the action of gravitation and stresses induced by gravitation. It is important that the gradient of the space density outside the object is zero. Two objects in space cannot be in equilibrium and it is natural to suppose that the space density between them is not uniform. To take the equilibrium state and to reduce the gradient of the space density between the objects, they should move towards each other. Two situations are possible resulting in stable equilibrium or stable motion. First, the collision and the merge into one object can take place resulting in the equilibrium of the new object and zero gradient of the space density. Second case can take place if the trajectories of the moving objects are affected by perturbations induced by other objects in space. In this case, the collision does not occur and the objects orbit in elliptical paths.

6. Conclusion

The general relativity equations are supplemented with the coordinate conditions following from the stationarity condition of the three-dimensional metric tensor density and equations for the stresses similar to compatibility equations of the theory of elasticity. The solution of the obtained complete set of equations is demonstrated for linearized and general spherically symmetric problems and linearized axially symmetric problem. The space density which is the ratio of the three-dimensional metric tensor densities in Riemannian and Euclidean spaces in the same coordinates is introduced and used to explain the attraction of objects under the action of gravitation.

Conflicts of Interest

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

  Logunov, A.A., Mestvirishvili, M.A. and Petrov, V.A. (2004) Physics-Uspekhi, 47, 607-621. https://doi.org/10.1070/PU2004v047n06ABEH001817  Landau, L.D. and Lifshits, E.M. (1988) Field Theory. Nauka, Moscow. (In Russian)  Arroyo-Torres, D., Wittkovski, M, Marcaide, J.M. and Hauschildt, P.H. (2013) Astronomy and Astrophysics, 554, A76. https://doi.org/10.1051/0004-6361/201220920  Singe, J.L. (1960) Relativity: The General Theory. North-Holland Publishing Company, Amsterdam.  Schwarzschild, K. (1916) Sitz Preuss. Akad. Wiss., 189-207.  Vasiliev, V.V. and Fedorov, L.V. (2018) Journal of Modern Physics, 9, 2482-2494. https://doi.org/10.4236/jmp.2018.914160  Schwarzschild, K. (1916) Sitz Preuss. Akad. Wiss., 424-432.  Weinberg, S. (1972) Gravitation and Cosmology. John Wiley and Sons, Inc., New York.  Fock, V. (1959) The Theory of Space, Time and Gravitation. Pergamon Press, London.  Vasiliev, V.V. (2017) Journal of Modern Physics, 8, 1087-1100. https://doi.org/10.4236/jmp.2017.87070  Vasiliev, V.V. (1989) Mechanics of Solids, 5, 30-34.  Schrodinger, E. (1950) Space-Time Structure. The University Press, Cambridge.  Logunov, A.A. (2006) Relativistic Theory of Gravitation. Nauka, Moscow. (In Russian)  Vasiliev, V.V. and Fedorov, L.V. (2018) Mechanics of Solids, 53, 256-261. https://doi.org/10.3103/S0025654418070038  Polyanin, A.D. (2001) Handbook on Linear Equations of Mathematical Physics. Fizmatlit, Moscow. (In Russian) 