View source for Gravitational stress-energy tensor

'''Gravitational stress-energy tensor''' is a symmetric tensor of second valence (rank), which describes the energy density and energy flux density of gravitational field in [[Physics/Essays/Fedosin/Lorentz-invariant theory of gravitation |Lorentz-invariant theory of gravitation]]. This tensor in the [[Physics/Essays/Fedosin/Covariant theory of gravitation |covariant theory of gravitation]] is included in equation for determining metric along with [[acceleration stress-energy tensor]], [[pressure stress-energy tensor]], [[dissipation stress-energy tensor]] and stress-energy tensor of electromagnetic field. The covariant derivative of gravitational stress-energy tensor determines density of gravitational force acting on matter.

== Lorentz-invariant theory of gravitation (LITG) ==
In LITG the gravitational stress-energy tensor is determined through [[Physics/Essays/Fedosin/Gravitational tensor |gravitational tensor]] <math> ~\Phi_{ik}</math> and metric tensor <math>~ \eta^{ik}</math> in Lorentzian metrics: <ref name="fed"> [[User:Fedosin | Fedosin S.G.]] [http://lccn.loc.gov/2009457349 Fizika i filosofiia podobiia ot preonov do metagalaktik], Perm, pages 544, 1999. {{ISBN|5-8131-0012-1}}. </ref>
:<math>~ U^{ik} = - \frac{c^2} {4 \pi G }\left( -\eta^{im}\Phi_{n m}\Phi^{n k}+ \frac{1} {4} \eta^{ik}\Phi_{mr}\Phi^{mr}\right) ,</math>

where <math>~ G </math> is [[w:gravitational constant |gravitational constant]], and all gravitational quantities are defined in electromagnetic wave representation, which uses the [[w:speed of light |speed of light]] <math>~ c</math>.

After replacing <math> ~ G </math> by the [[Physics/Essays/Fedosin/Strong gravitational constant |strong gravitational constant]] <math> ~ \Gamma </math>, the gravitational stress-energy tensor can be used to describe [[Physics/Essays/Fedosin/Strong gravitation |strong gravitation]] at the level of atoms and elementary particles in [[Physics/Essays/Fedosin/Gravitational model of strong interaction |gravitational model of strong interaction]].

===Components of gravitational stress-energy tensor ===
Since gravitational tensor in LITG consists of components of vectors of [[Physics/Essays/Fedosin/Gravitational field strength |gravitational field strength]] <math> ~ \mathbf{\Gamma }</math> and [[Physics/Essays/Fedosin/Gravitational torsion field |gravitational torsion field]] <math> ~ \mathbf{\Omega}</math>, and the tensor <math>~ \eta^{ik}</math> in 4-coordinates (''ct, x,y,z'') consists of the numbers 0, 1, -1 and does not depend on coordinates and time, so the components of gravitational stress-energy tensor can be written explicitly in terms of components of the mentioned vectors:

:<math>~ U^{ik} = \begin{vmatrix} u & \frac {H_x}{c}  & \frac {H_y}{c} & \frac {H_z}{c} \\ \frac {H_x}{c} & u+ \frac{\Gamma^2_x+c^2 \Omega^2_x}{4\pi G} & \frac{\Gamma_x \Gamma_y+c^2 \Omega_x\Omega_y }{4\pi G} & \frac{\Gamma_x \Gamma_z+c^2 \Omega_x\Omega_z }{4\pi G} \\ \frac {H_y}{c} & \frac{\Gamma_x \Gamma_y+c^2 \Omega_x\Omega_y }{4\pi G} & u+\frac{\Gamma^2_y+c^2 \Omega^2_y }{4\pi G}  & \frac{\Gamma_y \Gamma_z+c^2 \Omega_y\Omega_z }{4\pi G} \\ \frac {H_z}{c} & \frac{\Gamma_x \Gamma_z+c^2 \Omega_x\Omega_z }{4\pi G}  & \frac{\Gamma_y \Gamma_z+c^2 \Omega_y\Omega_z }{4\pi G} & u+\frac{\Gamma^2_z+c^2 \Omega^2_z }{4\pi G}  \end{vmatrix}. </math>

The time-like components of the tensor denote:

1) The volumetric energy density of gravitational field, negative in value
:<math>~ U^{00} = u= -\frac{1}{8 \pi G }\left(\Gamma^2+ c^2 \Omega^2 \right).</math>

2) The vector of energy flux density of gravitational field or the [[Physics/Essays/Fedosin/Heaviside vector |Heaviside vector]] 
:<math>~\mathbf{H} =-\frac{ c^2 }{4 \pi G }[\mathbf{\Gamma }\times \mathbf{\Omega }].</math>

The components of the vector <math>~\mathbf{H} </math> are part of the corresponding tensor components <math>~ U^{01}, U^{02}, U^{03}</math>, and due to the symmetry of the tensor indices <math>~ U^{01}= U^{10}, U^{02}= U^{20}, U^{03}= U^{30}</math>.

According to the Heaviside theorem, the relation holds:
:<math>~ \nabla \cdot \mathbf{H} = - \frac{ \partial U^{00}}{\partial t} - \mathbf{\Gamma }\cdot \mathbf{J},</math>

where <math>~\mathbf{J} </math> is the 3-vector of mass current density.

3) The space-like components of the tensor form a submatrix 3 x 3, which is the 3-dimensional gravitational stress tensor, taken with a minus sign. Gravitational stress tensor can be written as <ref name="fed"/>

:<math>~ \sigma^{p q} = \frac {1}{4 \pi G } \left( - \Gamma^p \Gamma^q - c^2 \Omega^p \Omega^q + \frac {1}{2} \delta^{pq} (\Gamma^2 + c^2 \Omega^2 ) \right) ,</math>

where <math>~p,q =1,2,3, </math> <math>~\Gamma^1=\Gamma_x, </math> <math>~\Gamma^2=\Gamma_y, </math> <math>~\Gamma^3=\Gamma_z, </math>  <math>~\Omega^1=\Omega_x, </math> <math>~\Omega^2=\Omega_y, </math> <math>~\Omega^3=\Omega_z, </math> <math>~\delta^{pq}</math> is the [[w:Kronecker delta |Kronecker delta]], <math>~\delta^{pq}=1 </math> if <math>~p=q, </math> and <math>~\delta^{pq}=0 </math> if <math>~p \not=q. </math>

Calculation of three-dimensional divergence of gravitational stress tensor gives:
:<math>~ \partial_q \sigma^{p q} = f^p +\frac {1}{c^2 } \frac{ \partial H^p}{\partial t}, </math>
where <math>~ f^p </math> denote the components of three-dimensional density of gravitational force, <math>~ H^p </math> – components of Heaviside vector.

=== Gravitational force ===
The gravitational stress-energy tensor has such form that it allows us to find the 4-vector of gravitational force density <math>~ f^\alpha </math> by differentiation in four-dimensional space:
:<math>~ f^\alpha = -\partial_\beta U^{\alpha \beta} = {\Phi^\alpha}_{i} J^i . \qquad (1) </math>

As we can see from formula (1), the 4-vector of gravitational force density can be calculated in a different way, through the gravitational tensor with mixed indices <math> {\Phi^\alpha}_{i} </math> and the 4-vector of mass current density <math> ~J^i</math>. This is due to the fact that in LITG the gravitational field equations have the form:
:<math>~ \partial_n \Phi_{ik} + \partial_i \Phi_{kn} + \partial_k \Phi_{ni}=0, </math>

:<math>~\partial_k \Phi^{ik} = \frac {4 \pi G }{c^2} J^i .</math>

Expressing from the latter equation <math> ~J^i</math> in terms of <math> \Phi^{ik} </math> and substituting in (1) and also using definition of gravitational stress-energy tensor, we can prove the validity of equation (1). The components of 4-vector of gravitational force density are as follows:
:<math>~ f^\alpha = (\frac {\mathbf{\Gamma } \cdot \mathbf{J} }{c }, \mathbf{f} ),</math>

where <math>~ \mathbf{f}= \rho \mathbf{\Gamma } + [\mathbf{J} \times \mathbf{\Omega} ]</math>  is 3-vector of gravitational force density, <math>~\rho</math> is density of moving matter, <math>~\mathbf{J} =\rho \mathbf{v} </math> is 3-vector of mass current density, <math>~\mathbf{v} </math> is 3-vector of particle velocity.

The integral of (1) over the three-dimensional volume of a small particle or a matter unit, calculated in the reference frame co-moving with the particle, gives the gravitational [[four-force]]:
:<math>~ F^\alpha = {\Phi^\alpha}_{i} M u^i= {\Phi^\alpha}_{i} p^i = (\frac {\mathbf{\Gamma } \cdot \mathbf{p} }{c}, \frac{E}{ c^2 } \mathbf{\Gamma }+ [\mathbf{p} \times \mathbf{\Omega}]) . </math>

In the integration it was taken into account that <math> ~J^i = \rho_0 u^i</math>, where <math> ~ \rho_0</math> is mass density in the co-moving reference frame, <math> ~ M</math> is the invariant mass, <math> ~ u^i</math> is the 4-velocity of the particle, <math> ~ p^i</math> is the 4-momentum of the particle, <math> ~\mathbf{p} </math> is the relativistic momentum, <math> ~E </math> is the relativistic energy of the particle. It is also assumed that the mass densities <math> ~ \rho_0</math> and <math> ~ \rho</math> include contributions from the mass-energy of the proper gravitational field and the electromagnetic field of the particle. The obtained 4-force is acting on the particle in the presence of the gravitational field with the tensor <math> ~ {\Phi^\alpha}_{i} </math>, and in some cases we can neglect the proper gravitational field of the particle and consider its motion only in the external field.

=== Integral vector ===
The gravitational stress-energy tensor contains time-like components <math>~ U^{0k}</math>, integrating which over the moving volume we can calculate integral vector of energy flux of gravitational field:
:<math>~ Q^k = \int {\frac { U^{0k}}{ c } dV} = (\frac {U}{c}, \mathbf {Q}) , </math>

where <math>~ U = \int { U^{00} dV}</math> is the energy of gravitational field, <math>~ \mathbf {Q} = \frac {1}{c^2 } \int { \mathbf {H}dV}</math> is a three-dimensional vector having the dimension of momentum.

If in this volume there is matter as a source of proper gravitational (electromagnetic) field, we should consider total integral vector of energy flux <math>~ p^k  </math>, which includes contributions from all the fields in the given volume, including [[acceleration field]] and [[pressure field]].
In particular, for a uniform spherical body with the radius <math>~ R </math> the integral vector in view of proper gravitational field of the body has the form: <ref> Fedosin S.G. [http://journals.yu.edu.jo/jjp/Vol8No1Contents2015.html  Relativistic Energy and Mass in the Weak Field Limit.] Jordan Journal of Physics. Vol. 8 (No. 1), pp. 1-16, (2015).</ref>
:<math>~ p^k = (m_g + \frac {U_0}{2c^2}) u^k =M u^k ,</math>

where <math>~ U_0 = - \frac {3 G m^2_g }{5R}</math> is total gravitational energy of the body in the reference frame in which the body is at rest, <math>~ M</math> is invariant inertial mass, which includes contributions of mass-energy from all the fields.
<math>~ m_g </math> is gravitational mass, which is equal to total mass <math>~ m_b </math> of all particles in the system, found by integrating over volume of the body the product of invariant mass density and Lorentz factor of the moving particles.

It is assumed that the mass <math> ~ M </math> equal to the mass <math> ~ m'</math> of matter particles without their gravitational binding energy, and when the particles combine in the whole body the bulk of their gravitational energy is compensated by internal energy of the particles motion and the energy of the body pressure. Since the energy <math>~ U_0</math> is negative, then the condition holds <math>~ m_b =m_g > M =m' </math>,  that is as long as the scattered matter is undergoing gravitational contraction into the body of finite size, the gravitational mass is increasing. We can find also the [[invariant energy]] of the physical system:
:<math>~ E_0 =Mc^2= c \sqrt {g_{ik}p^i p^k } . </math>

A more accurate analysis of masses of the system shows that the following relationship holds:<ref> Fedosin S.G. [https://www.bpasjournals.com/physics/chapter-details.php?chap_id=2367&issue_type=177&volume=117&journal=7 The Mass Hierarchy in the Relativistic Uniform System]. Bulletin of Pure and Applied Sciences, Vol. 38 D (Physics), No. 2, pp. 73-80 (2019). http://dx.doi.org/10.5958/2320-3218.2019.00012.5. </ref> 
:<math>~m'  < M < m < m_b = m_g ,</math>

where the gauge mass <math>~m' </math> is related to cosmological constant and represents mass-energy of the matter’s particles in four-potentials of the system’s fields; <math>~M </math> is the inertial mass; the auxiliary mass <math>~m </math> is equal to the product of the particles’ invariant mass density by the volume of the body.

This can be compared with the approach of the [[general relativity]], in which such total mass of baryons <math>~ m_b </math> is used that when we add to it negative mass from the gravitational field energy, we obtain the relativistic mass: <math>~ M = m_b + \frac {U_0}{c^2}</math>, and there is a relation: <math>~ m_g = M < m_b  <  m' </math>.

== Covariant theory of gravitation (CTG) ==
CTG is a generalization of LITG to any reference frames and phenomena that occur in the presence of fields and accelerations of acting forces. In CTG all deviations from LITG relations are described by the [[w:metric tensor |metric tensor]] <math>~ g^{ik}</math> which becomes a function of coordinates and time. In addition, in equations the operator of 4-gradient <math>~ \partial_k</math> is replaced by the [[w:covariant derivative |covariant derivative]] <math>~ \nabla_k</math>. After replacing <math>~ \eta^{ik}</math> by <math>~ g^{ik}</math> the gravitational stress-energy tensor becomes in the following form:
:<math>~ U^{ik} = - \frac{c^2} {4 \pi G }\left( -g^{im}\Phi_{n m}\Phi^{n k}+ \frac{1} {4} g^{ik}\Phi_{mr}\Phi^{mr}\right) .</math>

Gravitational field is considered at the same time as a component of [[Physics/Essays/Fedosin/General field |general field]]. <ref>Fedosin S.G. Two components of the macroscopic general field. Reports in Advances of Physical Sciences, Vol. 1, No. 2, 1750002, 9 pages (2017). http://dx.doi.org/10.1142/S2424942417500025. </ref>

Transforming contravariant indices in [[Physics/Essays/Fedosin/Gravitational tensor |gravitational tensor]] <math>~ \Phi^{mr}</math> into covariant indices with the help of metric tensor, and interchanging some indices, which are summed up, we obtain:
:<math>~ U^{ik} = - \frac{c^2} {4 \pi G }g^{ms} \left( - g^{ir} g^{n k} + \frac{1} {4} g^{ik} g^{n r} \right) \Phi_{m r}\Phi_{s n}.</math>

Since gravitational tensor <math>~ \Phi_{mr}</math> with covariant indices consists of components of vector of gravitational field strength <math> ~ \mathbf{\Gamma }</math>, divided by the speed of light <math> ~ c</math>, and components of vector of gravitational torsion field <math> ~ \mathbf{\Omega}</math>, the formula shows that in curved spacetime the gravitational stress-energy tensor is a sum of products of components of these vectors with the corresponding coefficients of components of metric tensor. And it turns out that the energy density of gravitational field <math>~ U^{00}</math> contains mixed products of the form <math>~ \Gamma_x \Omega_y </math>, etc. There are no such products in the flat Minkowski space, which leads to the fact that in spacetime of special relativity the energy associated with the strength of gravitational field is not mixed with the energy of torsion field. The same situation takes place in electromagnetism: in Minkowski space the energy of electric field is calculated separately from the energy of magnetic field, but in curved spacetime in the energy density of electromagnetic field there is additional energy from mixed components with the products of components of strengths of electric and magnetic fields.

Due to use of covariant differentiation in four-dimensional space in CTG the gravitational field equations are changed as well as the expression for 4-vector of gravitational force density (1), while the expression for gravitational 4-force remains the same: <ref name="fs"> Fedosin S.G. [http://sergf.ru/knen.htm Fizicheskie teorii i beskonechnaia vlozhennost’ materii.] – Perm, 2009, 844 pages, Tabl. 21, Pic. 41, Ref. 289. {{ISBN|978-5-9901951-1-0}}. (in Russian).</ref> 

:<math>~ \nabla_n \Phi_{ik} + \nabla_i \Phi_{kn} + \nabla_k \Phi_{ni}=0, </math>

:<math>~\nabla_k \Phi^{ik} = \frac {4 \pi G }{c^2} J^i ,</math>

:<math>~ f^\alpha = -\nabla_\beta U^{\alpha \beta} = {\Phi^\alpha}_{i} J^i , </math>

:<math>~ F^\alpha = {\Phi^\alpha}_{i} M u^i= {\Phi^\alpha}_{i} p^i. </math>

=== Equation for metric ===
In CTG the metric tensor is determined by solving the equation similar to Hilbert-Einstein equation. In covariant indices this equation can be written as follows: <ref name="fs1"> Fedosin S.G. [http://journals.yu.edu.jo/jjp/Vol9No1Contents2016.html About the cosmological constant, acceleration field, pressure field and energy.] Jordan Journal of Physics. Vol. 9 (No. 1), pp. 1-30 (2016). </ref>
:<math>~ R_{ik} - \frac{1} {4 }g_{ik}R = \frac{8 \pi G \beta  }{ c^4} \left( B_{ik}+ P_{ik}+ U_{ik}+ W_{ik} \right), </math>

where <math>~ R_{ik}</math> is the Ricci tensor, <math>~ R=R_{ik}g^{ik}</math> is the scalar curvature, <math>~ \beta </math> is the coefficient to be determined, <math>~ B_{ik}</math>, <math>~ P_{ik}</math>, <math>~ U_{ik}</math> and <math>~ W_{ik}</math> are the stress-energy tensors of the acceleration field, pressure field, gravitational and electromagnetic fields, respectively.

In contrast to the general theory of relativity, in this equation, there is no cosmological constant <math>~ \Lambda </math>, there is an additional constant <math>~ \beta </math> and the metric is dependent on the gravitational stress-energy tensor. The latter is the consequence of the fact that in CTG gravitation is an independent physical force as well as the electromagnetic force, and therefore participates in determining the metric according to the principles of the [[Physics/Essays/Fedosin/Metric theory of relativity |metric theory of relativity]].

=== Equation of motion ===
Using the principle of least action allows us to deduce not only the formula for the gravitational stress-energy tensor, but also gives the equation of motion written in tensor form: <ref> Fedosin S.G. Lagrangian formalism in the theory of relativistic vector fields. International Journal of Modern Physics A, Vol. 40, No. 02, 2450163 (2025). https://doi.org/10.1142/S0217751X2450163X. </ref>
:<math>~ \nabla_\beta \left( {B_\alpha}^\beta + {U_\alpha}^\beta +{W_\alpha}^\beta + {P_\alpha}^\beta \right) = 0. \qquad (2)</math>

The covariant derivative of the [[acceleration stress-energy tensor]] defines up to the sign the density of the 4-force acting on the field from the matter. At the same time the [[operator of proper-time-derivative]]  is applied to the 4-potential <math> ~ U_\alpha </math> of acceleration field in the Riemannian space: <ref> Fedosin S.G. Equations of Motion in the Theory of Relativistic Vector Fields. International Letters of Chemistry, Physics and Astronomy, Vol. 83, pp. 12-30 (2019). https://doi.org/10.18052/www.scipress.com/ILCPA.83.12. </ref>

:<math> ~ f_\alpha =  \nabla_\beta {B_\alpha}^\beta = - u_{\alpha k} J^k = \rho_0 \frac {DU_\alpha }{D \tau}- J^k \nabla_\alpha U_k = \rho_0 \frac {dU_\alpha }{d \tau} - J^k \partial_\alpha U_k ,</math>

where <math> ~ u_{\alpha k}</math> is the [[acceleration tensor]], <math> ~\tau </math> is the proper time of the particle in comoving reference frame.

The total density of the 4-force of the gravitational and electromagnetic fields and pressure field is determined by transfer of the stress-energy tensors of the fields to the right side of the equation of motion (2) and then applying the covariant derivative:
: <math> ~ f_\alpha = - \nabla_\beta \left( {U_\alpha}^\beta + {W_\alpha}^\beta + {P_\alpha}^\beta \right) =  \Phi_{\alpha k} J^k + F_{\alpha k} j^k  + f_{\alpha k} J^k  , </math>

where <math> ~ F_{\alpha k}</math> is the [[w:electromagnetic tensor |electromagnetic tensor]], <math> ~ f_{\alpha k}</math> is the [[pressure field tensor]], <math>~j^k = \rho_{0q} u^k </math> is the 4-vector of the electromagnetic current density, <math>~\rho_{0q}</math> is the density of electric charge of the matter unit in the reference frame at rest.

=== Conservation laws===
In the weak field limit, when the covariant derivative can be replaced by the partial derivative, for the time-like component in (2), which has the index <math>~ i=0</math>, the local conservation law of energy and energy flux in matter and gravitational and electromagnetic fields can be written as follows: 
<ref name="f"> Fedosin S.G. [http://vixra.org/abs/1403.0973 The Integral Energy-Momentum 4-Vector and Analysis of 4/3 Problem Based on the Pressure Field and Acceleration Field.] American Journal of Modern Physics. Vol. 3, No. 4, pp. 152-167 (2014). http://dx.doi.org/10.11648/j.ajmp.20140304.12.</ref> <ref name="gen"> Fedosin S.G. The generalized Poynting theorem for the general field and solution of the 4/3 problem. International Frontier Science Letters, Vol. 14, pp. 19-40 (2019). https://doi.org/10.18052/www.scipress.com/IFSL.14.19. </ref>
:<math>~ \nabla \cdot (\mathbf{ K }+ \mathbf{H}+\mathbf{P}+ \mathbf{F} ) = -\frac{\partial (B^{00}+U^{00}+W^{00}+P^{00} )}{\partial t},</math>

where <math>~ \mathbf{ K }</math> is the vector of the acceleration field energy flux density,  <math>~ \mathbf{H}</math> is the [[Physics/Essays/Fedosin/Heaviside vector |Heaviside vector]], <math>~ \mathbf{ P }</math> is the [[w:Poynting vector |Poynting vector]], <math>~ \mathbf{F}</math> is the vector of the pressure field energy flux density, which are determined in the [[special relativity]].

According to this law, the work of the field to accelerate the masses and charges is compensated by the work of the matter to create the field. As a result, the change in time of the total energy in a certain volume is possible only due to the inflow of energy fluxes into this volume.

The analysis of equation (2) for the space-like components with the index <math>~ i=1,2,3</math> in a weak field shows that, taking into account the equation of motion of the matter in the field, all the force densities and the change rates of the momentums of the matter and fields are cancelled.

However in general case, when the spacetime is significantly curved by the existing fields and matter, in (2) we should consider the contributions with the additional non-zero components of the metric tensor and their derivatives, which are absent in the special theory of relativity. This follows from the fact that the covariant derivative is expressed through the partial derivative and the Christoffel coefficients. Since in CTG the purpose of using the metric tensor is correction of the motion equations in order to take into account the dependence on the fields of the time intervals and spatial distances, measured by the electromagnetic (gravitational) waves, then such correction changes notation of many physical quantities in the form of 4-vectors and tensors. In particular, if we consider equation (2) as a local conservation law of energy and energy flux in matter and gravitational and electromagnetic fields, then appearing of additional contributions with the components of the metric tensor leads to specification of the theory for the case of the curved space. And the physical meaning of the results obtained for the flat spacetime and the weak field remains unchanged.

As another example we can consider the integral of the left side of (2) over the entire four-dimensional space. In flat Minkowski space the covariant divergence of the sum of the tensors becomes an ordinary 4-divergence, to which we can apply the divergence theorem. By this theorem the integral of the 4-divergence of some tensor over the 4-space can be replaced by the integral of the tensor over the hypersurface surrounding the 4-volume, over which the integration is done. If we choose a projection of this hypersurface on the hyperplane of the constant time in the form of a three-dimensional volume the integral of the left side of (2) is transformed into the integral of the sum of time components of the tensors in (2) over the volume, which must equal to conserved integral vector of the physical system under consideration:
:<math>~ \mathbb{Q}^i= \int{ \left( B^{i0}+ U^{i0} +W^{i0}+P^{i0} \right) dV }. </math>

For these tensor components the integral vector <math> ~ \mathbb {Q}^i </math> vanishes. <ref name="f"/> Vanishing of the four-vector allows us to explain the 4/3 problem, according to which the mass-energy of gravitational or electromagnetic field in the flux of field of moving system in 4/3 more than in the field energy of fixed system. On the other hand, according to, <ref name="gen"/> the generalized Poynting theorem and the integral vector should be considered differently inside the matter and beyond its limits. As a result, the occurrence of the 4/3 problem is associated with the fact that the time components of the stress-energy tensors do not form four-vectors, and therefore they cannot define the same mass in the fields’ energy and flux energy in principle. It follows that the integral vector <math>~ \mathbb{Q}^i</math> cannot be considered as a four-momentum of the system, determining only the distribution of energy and energy fluxes of the fields of the system.<ref> Fedosin S.G. [http://www.bpasjournals.com/physics/chapter-details.php?chap_id=2121&issue_type=158&volume=110&journal=7 The covariant additive integrals of motion in the theory of relativistic vector fields]. Bulletin of Pure and Applied Sciences, Vol. 37 D (Physics), No. 2, pp. 64-87 (2018). http://dx.doi.org/10.5958/2320-3218.2018.00013.1. </ref>
This is shown in more detail in the papers describing the generalized 4-momentum <ref>Fedosin S.G. Generalized Four-momentum for Continuously Distributed Materials. Gazi University Journal of Science, Vol. 37, Issue 3, pp. 1509-1538 (2024). https://doi.org/10.35378/gujs.1231793.</ref> and the total 4-momentum of a physical system.<ref>Fedosin S.G. What should we understand by the four-momentum of physical system? Physica Scripta, Vol. 99, No. 5, 055034 (2024). https://doi.org/10.1088/1402-4896/ad3b45.? </ref>

== General theory of relativity (GTR) ==
In GTR, there is a problem of determining the gravitational stress-energy tensor. The reason is that due to the applied principle of geometrization of physics, all manifestations of gravitation are completely replaced by the geometric effect – the spacetime curvature. Thus, the gravitational field is reduced to the metric field, set by the metric tensor and its derivatives with respect to coordinates and time. Since in each reference frame there is its own metric, so gravitation is not independent physical interaction but is assumed as the consequence of the metric of the reference frame. This leads to the fact that instead of the gravitational stress-energy tensor in GTR there is a pseudotensor which depends on the metric. The feature of this pseudotensor is that locally it can be made equal to zero at any point by choosing the corresponding reference frame. As a result, GTR have to refuse from the possibility to accurately determine the localization of gravitational energy and the momentum of gravitational field in a physical system, which significantly impedes understanding the physics of gravitation on a fundamental level, and describing the phenomena in a classical form.
The purpose of Hilbert-Einstein equations in GTR is to find the metric:
:<math>~ R_{ik} - \frac{1} {2 }g_{ik}R + g_{ik} \Lambda = \frac{8 \pi G }{ c^4} \left( t_{ik} + W_{ik} \right). </math>

These equations do not contain the gravitational stress-energy tensor, which is present in the right side of the equations for the metric in CTG. The metric in GTR depends only on the matter and the electromagnetic field in the considered reference frame. After finding such metric it is used to find the pseudotensor of the energy-momentum of gravitational field, and this pseudotensor must depend only on the metric tensor, be symmetrical with respect to the indices, and when added to the stress-energy tensors of matter <math>~ t_{ik} </math> and electromagnetic field <math>~ W_{ik} </math> it must give zero divergence so that the conservation law of energy-momentum of the matter and field would be satisfied. This pseudotensor in a weak field, that is, in the special theory of relativity, must vanish in order to ensure the principle of equivalence of free fall in the gravitational field and inertial motion.

The Landau-Lifshitz pseudotensor corresponds to the mentioned criteria. <ref>  Landau L.D., Lifshitz E.M. (1975). The Classical Theory of Fields. Vol. 2 (4th ed.). Butterworth-Heinemann. {{ISBN|978-0-750-62768-9}}.</ref> There is also an Einstein pseudotensor, but it is not symmetrical with respect to the indices. <ref> Albert Einstein. Das hamiltonisches Prinzip und allgemeine Relativitätstheorie. (The Hamiltonian principle and general relativity). Sitzungsber. preuss. Acad. Wiss. 1916, 2, 1111–1116. </ref> Note that in contrast to GTR, in mathematics under pseudotensor we understand something different, more precisely a tensor quantity, which changes its sign when it is transformed into the coordinate reference frame with the opposite orientation of the basis.

Since in GTR there is no gravitational stress-energy tensor, the gravitational energy of a body is determined indirectly. For example, one method is when by the given mass density and body size the mass of the body is calculated, first in the absence of the metric’s influence, and then taking into account the metric’s influence and the corresponding change in the volume differential in the integral of the mass under influence of the gravitational field. The difference of the mentioned masses is equated to the mass-energy of gravitational field as to the manifestation of the metric field. <ref>  Abhas Mitra. [http://arxiv.org/abs/gr-qc/0605066 Why Gravitational Contraction Must be Accompanied by Emission of Radiation both in Newtonian and Einstein Gravity.] Phys. Rev. D74, (2006), 024010. </ref> In this case, the components of the metric tensor are used, found in GTR from the Hilbert-Einstein equation for the metric.

== See also ==
* [[w:Electromagnetic stress-energy tensor |Electromagnetic stress-energy tensor]]
* [[Pressure stress-energy tensor]]
* [[Acceleration stress-energy tensor]]
* [[Dissipation stress-energy tensor]]
* [[Physics/Essays/Fedosin/Gravitational tensor |Gravitational tensor]]
* [[w:Einstein tensor |Einstein tensor]]
* [[w:Four-momentum |Four-momentum]]
* [[Physics/Essays/Fedosin/General field |General field]]
* [[Dissipation field]]
* [[Acceleration field]]
* [[Pressure field]]

== References ==
<references/>

==External links ==
* [http://www.wikiznanie.ru/ru-wz/index.php/%D0%A2%D0%B5%D0%BD%D0%B7%D0%BE%D1%80_%D1%8D%D0%BD%D0%B5%D1%80%D0%B3%D0%B8%D0%B8-%D0%B8%D0%BC%D0%BF%D1%83%D0%BB%D1%8C%D1%81%D0%B0_%D0%B3%D1%80%D0%B0%D0%B2%D0%B8%D1%82%D0%B0%D1%86%D0%B8%D0%BE%D0%BD%D0%BD%D0%BE%D0%B3%D0%BE_%D0%BF%D0%BE%D0%BB%D1%8F Gravitational stress-energy tensor in Russian]

[[Category:Tensors]] 
[[Category:Theories of gravitation]]
[[Category:Covariant theory of gravitation]]