arxiv:gr-qc/ v1 26 Jul 1999

Transcription

1 Projections and covariant divergency of the energy-momentum tensors arxiv:gr-qc/ v1 26 Jul 1999 S. Manoff Bulgarian Academy of Sciences Institute for Nuclear Research and Nuclear Energy Department of Theoretical Physics Blvd. Tzarigradsko Chaussee Sofia - Bulgaria R. Lazov Bulgarian Academy of Sciences Institute for Mathematics and Informatics G. Bonchev Str Sofia - Bulgaria address: smanov@inrne.bas.bg Abstract The invariant projections of the energy-momentum tensors of Lagrangian densities for tensor fields over differentiable manifolds with contravariant and covariant affine connections and metrics [(L n, g)-spaces] are found by the use of an non-null (non-isotropic) contravariant vector field and its corresponding projective metrics. The notions of rest mass density, momentum density, energy current density and stress tensor are introduced as generalizations of these notions from the relativistic continuum media mechanics. The energy-momentum tensors are represented by means of the introduced notions and the corresponding identities are found. The notion of covariant differential operator along a contravariant tensor field is introduced. On its basis, as a special case, the notion of contravariant metric differential operator is proposed. The properties of the operators are considered. By the use of these operators the notion of covariant divergency of a mixed tensor field is determined. The covariant divergency of tensor fields of second rank of the types 1 and 2 is found. Invariant representations of the covariant divergency of the energymomentum tensor are obtained by means of the projective metrics of a contravariant non-isotropic (non-null) vector field and the corresponding rest mass density, momentum density, and energy flux density. An invariant representation of the first Noether identity is found as well as relations between the covariant divergencies of the different energy-momentum tensors and their structures determining covariant local conserved quantities. 1

2 1 Invariant projections of a mixed tensor field of second rank In the relativistic continuum media mechanics notions are introduced as generalizations of the same notions of the classical continuum media mechanics. This has been done by means of the projections of the (canonical, symmetric of Belinfante or symmetric of Hilbert) energy-momentum tensors along or orthogonal to a non-isotropic (non-null) contravariant vector field. There are possibilities for using the projections for finding out the physical interpretations of the determined energy-momentum tensors. In an analogous way as in (pseudo) Riemannian spaces without torsion (V n -spaces), the different relations between the quantities with well known physical interpretations can be considered as well as their application for physical systems described by means of mathematical models over differentiable manifold with affine connections and metrics [(L n, g)-spaces]. The energy-momentum tensors are obtained as mixed tensor fields of second rank of type 1 [by the use of the procedure on the grounds of the method of Lagrangians with covariant derivatives (MLCD) [1]] G = G α β e β e α = G i j j dx i (1) in contrast to the mixed tensor fields of second rank of type 2 G = G β α e α e β = G j i dx i j. (2) {e α } and {e α } are non-co-ordinate (non-holonomic) covariant and contravariant basic vector fields respectively, (α, β = 1,..., n), {dx i } and { i } are coordinate (holonomic) covariant and contravariant basic vector fields respectively (i, j = 1,..., n), dim(l n, g) = n. To every covariant basic vector field a contravariant basic vector field can be juxtaposed and vice versa by the use of the contravariant and covariant metric tensor fields g(e γ ) = g αγ e α, g( j ) = g ij dx i, g(e γ ) = g αγ e α, g(dx j ) = g ij i. (3) On this basis a tensor field of the type 2 can be related to a tensor field of the type 1 by the use of the covariant and contravariant tensor fields g = g ij dx i.dx j and g = g ij i. j [dx i.dx j = (1/2)(dx i dx j + dx j dx i ), i. j = (1/2)( i j + j i )] G = g(g)g = G β α e α e β = g αγ G δ γ g δβ e α e β, (4) G β α = g αγ G δ γ g δβ, (5) G = g(g)g = G α β e β e α = g βδ G γ δ g γα e β e α, (6) 2

3 G α β = g βδ G γ δ g γα. (7) The Kronecker tensor field appears as a mixed tensor field of second rank of the type 1 Kr = g α β e α e β = g i j i dx j and can be projected by means of the non-isotropic (non-null) contravariant vector field u and its projection metrics h u and h u [h u = g 1 e g(u) g(u), h u = g 1 e u u, e = g(u, u) 0] Kr = ε Kr u g(u) + u g( Kr π) + Kr s g(u) + ( Kr S)g, ε Kr = 1 e 2 [g(u)](kr)u = 1 e 2.u α u α = 1 e 2.g αβ u α u β = = 1 e 2.f ρ β f σ f β δ g ρσ u γ u δ = 1 e 2.u i u i = 1 e 2 g ik u i u k = 1 e k, (8) k = 1 e [g(u)](kr)u = 1 e.ui u i, u i = g ij u j, (9) Kr π = 1 e [g(u)](kr)hu = 1 e g βγ uγ g β δ hδα e α = 1 e g βγ hβα u γ e α = = 1 e u β hβα e α = Kr π α e α, (10) Kr π α = 1 e g βγ uγ h βα = 1 e u β hβα, Kr π i = 1 e g kl ul h ki = 1 e u k hki, (11) Kr s = 1 e hu (g)(kr)(u) = Kr s α e α = Kr s i i, (12) Kr s α = 1 e hαβ g βδ u δ, Kr s i = 1 e hik g kl u l, (13) Kr S = h u (g)(kr)h u = Kr S αβ e α e β = Kr S ij i j, (14) Kr S αβ = h αγ g γδ h δβ, Kr S ij = h ik g kl h lj, (15) h u [g(u)] = [g(u)](h u ) = 0, h u (g)h u = h u, (Kr)g = k e u u + u Kr π + Kr s u + Kr S, (Kr)g = g αβ e α e β = g ij i j. (16) The corresponding to the Kronecker tensor field mixed tensor field of the type 2 Kr = g(kr)g = g αγ g γβ e α e β = g ik g kj dx i j (17) does not appear in general as a Kronecker tensor field. 3

4 Special case: S = C : f i j = g i j, (fα β = g α β ): Kr = g α β eβ e α = g i j dxi j, k = 1, ε Kr = 1 e, Kr π = 0, Kr s = 0, Kr S = h u, Kr = 1 e u g(u) + (hu )g = g(g), (Kr)g = g. (18) The representation of the tensor fields of the type 1 by the use of the nonisotropic (non-null) contravariant vector field u and its projective metrics h u and h u corresponds in its form to the representation of the viscosity tensor and the energy-momentum tensors in the continuum media mechanics in V 3 - or V 4 - spaces, ε G is the inner energy density, G π is the conductive momentum, e. G s is the conductive energy flux density and G S is the stress tensor density. An analogous interpretation can also be accepted for the projections of the energy-momentum tensors found by means of the method of Lagrangians with covariant derivatives (MLCD). 2 Energy-momentum tensors and the rest mass density The covariant Noether identities (generalized covariant Bianchi identities) can be considered as identities for the components of mixed tensor fields of second rank of the first type. The second covariant Noether identity can be written in the form θ α β s T α β Q α β θ s T Q, (19) θ = θ α β e β e α = θ i j j dx i, st = s T α β e β e α = s T i j j dx i, Q = Q α β e β e α = Q i j j dx i, (20) θ is the generalized canonical energy-momentum tensor (GC-EMT) of the type 1; s T is the symmetric energy-momentum tensor of Belinfante (S-EMT-B) of the type 1; Q is the variational energy-momentum tensor of Euler-Lagrange (V-EMT-EL) of the type 1. The second covariant Noether identity for the energy-momentum tensors of the type 1 is called second covariant Noether identity of type 1. By means of the non-isotropic contravariant vector field u and its corresponding projective metric the energy-momentum tensors can be represented in an analogous way as the mixed tensor fields of the type 1. The structure of the generalized canonical energy-momentum tensor and the symmetric energy-momentum tensor of Belinfante for the metric and non-metric tensor fields has similar elements and they can be written in the form G = k G L Kr, θ = k θ L Kr, st = T L Kr, (21) 4

5 kθ = k θ α β e β e α = k θ i j j dx i. (22) On the analogy of the notions of the continuum media mechanics k G is called viscosity tensor field. G and k G can be written by means of u, h u and h u in the form G = ε G u g(u) + u g( G π) + G s g(u) + ( G S)g, kg = ε k u g(u) + u g( k π) + k s g(u) + ( k S)g. (23) From the relation (21) the relations between the different projections of G and k G follow. If we introduce the abbreviation ε G = ρ G, then ε k = ρ G + 1 e L k, k π = G π + L Kr π, k s = G s + L Kr s, k S = G S + L Kr S, (24) and G can be written by means of (21) in the form G = (ρ G + 1 e L k) u g(u) L Kr + u g(k π) + k s g(u) + ( k S)g, (25) ρ G = 1 e 2 [g(u)](g)(u) is the rest mass density of the energy-momentum tensor G of the type 1. This type of representation of a given energy-momentum tensor G by means of the projective metrics of u and ρ G is called representation of G by means of the projective metrics of the contravariant non-isotropic (non-null) vector field u and the rest mass density ρ G. There are other possibilities for representation of G by means of u and its corresponding projective metrics. If we introduce the abbreviation p G = ρ G u + G π, (26) p G is the momentum density of the energy-momentum tensor G of the type 1, then G can be written in the form G = u g(ρ G u + G π) + G s g(u) + ( G S)g, G = u g(p G ) + G s g(u) + ( G S)g. (27) The representation of G by means of the last relation is called representation of G by means of the projective metrics of the contravariant non-isotropic contravariant vector field u and the momentum density p G. By the use of the relations g( G π, u) = 0, ( G S)[g(u)] = 0, (28) valid (because of their constructions) for every energy-momentum tensor G and the definitions e G = G(u) = (G)(u) = e (ρ G u + G s), g(u, u) = e 0, (29) 5

6 e G is the energy flux density of the energy-momentum tensor G of the type 1, the tensor field G can be written in the form G = (ρ G u + G s) g(u) + u g( G π) + ( G S)g, G = 1 e e G g(u) + u g( G π) + ( G S)g. (30) The representation of G by means of the last expression is called representation of G by means of the projective metrics of the contravariant non-isotropic vector field u and the energy flux density e G. The generalized canonical energy-momentum tensor θ can be represented, in accordance with the above described procedure, by the use of the projective metrics of u and the rest mass density ρ θ θ = k θ L Kr, kθ = θ + L Kr, (31) θ = (ρ θ + 1 e L k) u g(u) L Kr + u g(θ π) + θ s g(u) + ( θ S)g, (32) kθ = k θ α β e β e α, kθ α β = t α β K α β W α βγ γ + L g β α, (33) ρ θ = 1 e 2 [g(u)](θ)(u), k = 1 [g(u)](kr)(u), (34) e ρ θ = 1 e g 2 αβ u β θ α γ u γ = 1 e θ α 2 γ u α u γ = = 1 e g 2 ij u j θ i k u k = 1 e θ i 2 k u i u k, (35) ε k θ = ρ θ + 1 e L k, (36) θ π = 1 e [g(u)]( kθ)h u = θ π α e α = θ π i i, (37) θ π α = 1 e kθ γ β u β h γα, θ π i = 1 e kθ l k u k h li, (38) θ s = 1 e hu (g)( k θ)(u) = θ s α e α = θ s i i, (39) θ s α = 1 e hαβ g βγ kθ δ γ u δ, θ s i = 1 e hij g jk kθ l k u l, (40) θ S = h u (g)( k θ)h u = θ S αβ e α e β = θ S ij i j, (41) θ S αβ = h αγ g γδ kθ κ δ h κβ, θ S ij = h ik g kl kθ m l h mj, (42) (θ)g = (ρ θ + 1 e L k) u u L Kr(g) + u θ π + θ s u + θ S = = θ αβ e α e β = θ ij i j, (43) θ αβ = g βγ θ γ α = θ γ α g βγ, θ ij = θ k i g kj, (44) g(θ) = (ρ θ + 1 e L k) g(u) g(u) L g(kr) + g(u) g(θ π)+ +g( θ s) g(u) + g( θ S)g, (45) 6

7 g(θ) = θ αβ e α e β = θ ij dx i dx j, θ αβ = g αγ θ β γ, θ ij = g ik θ j k, (46) g(kr) = g αβ e α e β = g ij dx i dx j, g αβ = g αγ f γ β = g γα f γ β = g βα, g ij = g ji, (Kr)g = Kr(g) = g αβ e α e β = g ij i j. (47) In a co-ordinate basis θ, (θ)g and g(θ) can be represented in the forms θ j i = (ρ θ + 1 e L k) u i u j L g j i + θ π i u j + u i θs j + g ik θs jk, (48) θ ij = θ k i g kj = (ρ θ + 1 e L k) ui u j L g ij +u i θπ j + θ s i u j + θ S ij, (49) θ ij = g ik θ j k = (ρ θ + 1 e L k) u i u j L g ij +u i θπ j + θ s i u j + g ik θs kl g lj, (50) u i = g ij u j, θ π i = g ik θπ k, θ s i = g il θs l, θ S ij = g ik θs kl g lj. (51) The symmetric energy-momentum tensor of Belinfante s T can be represented in an analogous way by the use of the projective metrics h u and h u and the rest mass density ρ T in the form st = (ρ T + 1 e L k).u g(u) L Kr + u g(t π) + T s g(u) + ( T S)g, (52) st = s T α β e β e α = s T i j j dx i, st = sk T L Kr, skt = s T + L Kr = T, (53) ε T = ρ T + 1 e L k, (54) ρ T = 1 e 2 [g(u)]( s T)(u) = 1 e 2 g αβ u β s T γ α u γ = 1 e 2 s T γ α u α u γ = = 1 e 2 g ij u j s T k i u k = 1 e 2 st k i u i u k, (55) T π = 1 e [g(u)](t)hu = 1 e [g(u)]( skt)h u = T π α e α = T π i i, (56) T π α = 1 e T γ β u β h γα, T π i = 1 e T l k u k h li, (57) T s = 1 e hu (g)(t)(u) = T s α e α = T s i i, (58) T s α = 1 e hαβ g βγ T δ γ u δ, T s i = 1 e hij g jk T l k u l, (59) T S = h u (g)(t)h u = T S αβ e α e β = T S ij i j, (60) T S αβ = h αγ g γδ T κ δ h κβ, T S ij = h ik g kl T m l h mj, (61) 7

8 ( s T)g = (ρ T + 1 e L k) u u L Kr(g) + u T π + T s u + T S = = s T αβ e α e β = s T ij i j, (62) st αβ = g βγ st γ α = s T γ α g βγ, st ij = s T k i g kj, (63) g( s T) = (ρ T + 1 e L k) g(u) g(u) L g(kr) + g(u) g(t π)+ +g( T s) g(u) + g( T S)g, g( s T) = s T αβ e α e β = s T ij dx i dx j, st αβ = g αγ st β γ, st ij = g ik st j k. In a co-ordinate basis s T, ( s T)g and g( s T) can be represented in the forms (64) (65) st i j = (ρ T + 1 e L k) u i u j L g j i + T π i u j + u i T s j + g ik TS jk, (66) st ij = s T i g kj k = (ρ T + 1 e L k) ui u j L g ij +u i Tπ j + T s i u j + T S ij, (67) st ij = g ik st k j = (ρ T + 1 e L k) u i u j L g ij +u i Tπ j + T s i u j + g ik TS kl g lj, (68) T π i = g ik Tπ k, T s i = g il Ts l, T S ij = g ik TS kl g lj. The variational energy-momentum tensor of Euler-Lagrange Q can be represented in the standard manner by the use of the projective metrics h u, h u and the rest mass density ρ Q in the form Q = ρ Q u g(u) u g( Q π) Q s g(u) ( Q S)g, (69) ρ Q = 1 [g(u)](q)(u), (70) e2 ρ Q = 1 e 2 g αβ u β Q γ α u γ = 1 e 2 Q γ α u α u γ = = 1 e 2 g ij u j Q k i u k = 1 e 2 Q k i u i u k, (71) Q π = 1 e [g(u)](q)hu = Q π α e α = Q π i i, (72) Q π α = 1 e Q γ β u β h γα, Q π i = 1 e Q l k u k h li, (73) Q s = 1 e hu (g)(q)(u) = Q s α e α = Q s i i, (74) Q s α = 1 e hαβ g βγ Q δ γ u δ, Q s i = 1 e hij g jk Q l k u l, (75) Q S = h u (g)(q)h u = Q S αβ e α e β = Q S ij i j, (76) Q S αβ = h αγ g γδ Q κ δ h κβ, Q S ij = h ik g kl Q m l h mj, (77) 8

9 (Q)g = ρ Q u u u Q π Q s u Q S = = Q αβ e α e β = Q ij i j, (78) Q αβ = g βγ Q γ α = Q γ α g βγ, Q ij = Q k i g kj, g(q) = ρ Q g(u) g(u) g(u) g( Q π) g( Q s) g(u) g( Q S)g, (79) g(q) = Q αβ e α e β = Q ij dx i dx j, Q αβ = g αγ Q β γ, Q ij = g ik Q j k. (80) In a co-ordinate basis Q, (Q)g and g(q) can be represented in the forms Q i j = ρ Q u i u j Q π i u j u i Qs j g ik QS jk, (81) Q ij = Q k i g kj = ρ Q u i u j u i Qπ j Q s i u j Q S ij, (82) Q ij = g ik Q j k = ρ Q u i u j u i Qπ j Q s i u j g ik QS kl g lj, (83) Q π i = g ik Qπ k, Q s i = g il Qs l, Q S ij = g ik QS kl g lj. The introduced abbreviations for the different projections of the energymomentum tensors have their analogous forms in V 3 - and V 4 -spaces, their physical interpretations have been proposed [2], [3] (S ). The stress tensor in V 3 -spaces has been generalized to the energy-momentum tensor s T in V 4 -spaces. The viscosity stress tensor sk T appears as the tensor T in the structure of the symmetric energy-momentum tensor of Belinfante s T. On the analogy of the physical interpretation of the different projections, the following definitions can be proposed for the quantities in the representations of the different energy-momentum tensors: A. Generalized canonical energy-momentum tensor of the type θ (a) Generalized viscous energy-momentum tensor of the type k θ (b) Rest mass density of the generalized canonical energy-momentum tensor θ....ρ θ (c) Conductive momentum density of the generalized canonical energy-momentum tensor θ.... θ π (d) Conductive energy flux density of the generalized canonical energy-momentum tensor θ....e θs (e) Stress tensor of the generalized canonical energy-momentum tensor θ.... θ S B. Symmetric energy-momentum tensor of Belinfante of the type s T (a) Symmetric viscous energy-momentum tensor of the type 1....T (b) Rest mass density of the symmetric energy-momentum tensor of Belinfante s T....ρ T (c) Conductive momentum density of the symmetric energy-momentum tensor of Belinfante s T.... T π 9

10 (d) Conductive energy flux density of the symmetric energy-momentum tensor of Belinfante s T....e Ts (e) Stress tensor of the symmetric energy-momentum tensor of Belinfante s T.... T S C. Variational (active) energy-momentum tensor of Euler-Lagrange... Q (a) Rest mass density of the variational energy-momentum tensor of Euler-Lagrange Q... ρ Q (b) Conductive momentum density of the variational energy-momentum tensor of Euler-Lagrange Q... Q π (c) Conductive energy flux density of the variational energy-momentum tensor of Euler-Lagrange Q... e Q s (d) Stress tensor of the variational energy-momentum tensor of Euler-Lagrange Q.... Q S The projections of the energy-momentum tensors have properties which are due to their construction, the orthogonality of the projective metrics h u and h u correspondingly to the vector fields u and g(u) [h u (u) = 0, h u [g(u)] = 0] g(u, θ π) = g( θ π, u) = 0, g(u, T π) = 0, g(u, Q π) = 0, (84) g(u, θ s) = g( θ s, u) = 0, g(u, T s) = 0, g(u, Q s) = 0, (85) g(u)( θ S) = 0, ( θ S)g(u) = 0, g(u)( T S) = 0, ( T S)g(u) = 0, g(u)( Q S) = 0, ( Q S)g(u) = 0. (86) From the properties of the different projections it follows that the conductive momentum density π (or π) is a contravariant vector field orthogonal to the vector field u. The conductive energy flux density e s (or e s) is also a contravariant vector field orthogonal to u. The stress tensor S (or S) is orthogonal to u independently of the side of the projection by means of the vector field u. The second covariant Noether identity θ s T Q can be written by the use of the projections of the energy-momentum tensors in the form (ρ θ ρ T + ρ Q ) u g(u) + u g( θ π T π + Q π)+ + ( θ s T s + Q s) g(u) + ( θ S T S + Q S)g 0. (87) After contraction of the last expression consistently with u and g and taking into account the properties (84) (86) the second covariant Noether identity disintegrates in identities for the different projections of the energy-momentum tensors ρ θ ρ T ρ Q, θ π T π Q π, θ s T s Q s, θ S T S Q S. (88) If the covariant Euler-Lagrange equations of the type δ v L/δV A B = 0 are fulfilled for the non-metric tensor fields of a Lagrangian system and g Q = 0, then the variational energy-momentum of Euler-Lagrange Q = v Q+ g Q is equal to zero. This fact leads to vanishing the invariant projections of Q (ρ Q = 0, 10

11 Q π = 0, Q s = 0, Q S = 0). The equality which follows between θ and s T has as corollaries the identities ρ θ ρ T, θ π T π, θ s T s, θ S T S. (89) From the first identity (ρ θ ρ T ) and the identity (88) for ρ, it follows that the covariant Euler-Lagrange equations of the type δ v L/δV A B = 0 for non-metric fields V and g Q = 0 appear as sufficient conditions for the unique determination of the notion of rest mass density ρ for a given Lagrangian system. Proposition 1 The necessary and sufficient condition for the equality ρ θ = ρ T (90) is the condition ρ Q = 0. Proof. It follows immediately from the first identity in (88). The condition ρ Q 0 leads to the violation of the unique determination of the notion of rest mass density and to the appearance of three different notions of rest mass density corresponding to the three different energy-momentum tensors for a Lagrangian system. Therefore, the violation of the covariant Euler- Lagrange equations δ v L/δV A B = 0 for the non-metric tensor fields or the existence of metric tensor fields in a Lagrangian density with g Q 0 induce a new rest mass density (a new rest mass respectively) for which the identity (88) is fulfilled. The identity ρ θ ρ T ρ Q can be related to the physical hypotheses about the inertial, passive and active gravitational rest mass densities in models for describing the gravitational interaction. To every energy-momentum tensor a non-null rest mass density corresponds. The existence of the variational energymomentum tensor of Euler-Lagrange is connected with the existence of the gravitational interaction in a Lagrangian system in Einstein s theory of gravitation [4] and therefore, with the existence of a non-null active gravitational rest mass density. When a Lagrangian system does not interact gravitationaly, the active gravitational rest mass density is equal to zero and the principle of equivalence between the inertial and the passive rest mass density is fulfilled [5] - [8]. From the second covariant Noether identity of the type 1. by means of the relations G = g(g)g, G = g(g)g, one can find the second covariant Noether identity of the type 2. for the energymomentum tensors of the type 2. in the form θ s T Q, (91) θ = g(θ)g, st = g( s T)g, Q = g(q)g. (92) 11

12 The invariant representation of the energy-momentum tensors by means of the projective metrics h u, h u and the rest mass density allows a comparison of these tensors with the well known energy-momentum tensors from the continuum media mechanics (for instance, with the energy-momentum tensor of an ideal liquid in V 4 -spaces: st i j = (ρ + 1 e p) u i u j p g j i, e = const. 0, k = 1 ). (93) It follows from the comparison that the Lagrangian invariant L can be interpreted as the pressure p = L characterizing the Lagrangian system. This possibility for an other physical interpretation than the usual one (in the mechanics L is interpreted as the difference between the kinetic and the potential energy) allows a description of Lagrangian systems on the basis of phenomenological investigations determining the dependence of the pressure on other dynamical characteristics of the system. If these relations are given, then by the use of the method of Lagrangians with covariant derivatives (MLCD) the corresponding covariant Euler-Lagrange equations can be found as well as the energy-momentum tensors. The use of a contravariant non-isotropic (non-null) vector field u and its corresponding projective metrics h u and h u is analogous to the application of a non-isotropic (time-like) vector field in the s. c. monad formalism [(3 + 1)- decomposition] in V 4 -spaces for description of dynamical systems in Einstein s theory of gravitation (ETG) [9] - [17], [3], [18]. The contravariant time-like vector field has been interpreted as a tangential vector field at the world line of an observer determining the frame of reference (the reference frame) in the space-time. By the use of this reference frame a given physical system is observed and described. The characteristics of the vector field determine the properties of the reference frame. Moreover, the vector field is assumed to be an absolute element in the scheme for describing the physical processes, i. e. the vector field is not an element of the model of the physical system. It is introduced as a priory given vector field which does not depend on the Lagrangian system. In fact, the physical interpretation of the contravariant non-isotropic vector field u can be related to two different approaches analogous to the method of Lagrange and the method of Euler in describing the motion of liquids in the hydrodynamics [19]. In the method of Lagrange, the object with the considered motion appears as a point (particle) of the liquid. The motion of this point is given by means of equations for the vector field u interpreted as the velocity of the particle. The solutions of these equations give the trajectories of the points in the liquid as basic characteristics of the physical system. In this case, the vector field u appears as an element of the model of the system. It is connected with the motions of the system s elements. Therefore, a Lagrangian system (and respectively its Lagrangian density) could contain as an internal characteristic a contravariant vector field u obeying equations of the type of the Euler-Lagrange equations and describing the evolution of the system. 12

13 In the method of Euler, the object with the considered motion appears as the model of the continuum media. Instead of the investigation of the motion of every (fixed by its velocity and position) point (particle), the kinematic characteristics in every immovable point in the space are considered as well as the change of these characteristics after moving on to an other space point. The motion is assumed to be described if the vector field u is considered as a given (or known) velocity vector field. The (n 1) + 1 decomposition (monad formalism) can be related to the method of Euler or to the method of Lagrange: (a) Method of Euler. The vector field u is interpreted as the velocity vector field of an observer who describes a physical system with respect to his vector field (his velocity). This physical system is characterized by means of a Lagrangian system (Lagrangian density). The motion of the observer (his velocity vector field) is given independently of the motion of the considered Lagrangian system. (b) Method of Lagrange. The vector field u is interpreted as the velocity vector field of a continuum media with a co-moving with it observer. The last assumption means that the velocity vector field of the observer is identical with the velocity vector field of the media he is situated. The motion of the observer (his velocity vector field) is determined by the characteristics of the (Lagrangian) system. Its velocity vector field is, on the other side, determined by the dynamical characteristics of the system by means of equations of the type of the Euler-Lagrange equations. 3 Energy-momentum tensors and the momentum density In the first section of this paper the notion of momentum density of an energymomentum tensor G has been introduced as p G = ρ G u + G π, ρ G.u is the convective momentum density of the energy-momentum tensor G; G π is the conductive momentum density of the energy-momentum tensor G. The tensor field G can be represented by means of the projective metrics h u and h u in the form G = u g(p G ) + G s g(u) + ( G S)g. At the same time the relations are fulfilled (G)g = u p G + G s u + G S = G αβ e α e β = G ij i j, (94) g(g) = g(u) g(p G ) + g( G s) g(u) + g( G S)g, (95) 13

14 p G = 1 [g(u)](g)(g). (96) e The Kronecker tensor can be represented in the form Kr = u g(p Kr ) + Kr s g(u) + ( Kr S)g, (97) p Kr = 1 e k u+ Kr π, L Kr = u g(l p Kr )+L Kr s g(u)+(l Kr S)g. (98) The generalized canonical energy-momentum tensor θ will have then the form θ = u g(p θ ) + θ s g(u) + ( θ S)g, (99) p θ = ρ θ u + θ π = ρ θ u + θ π L Kr π. (100) In a co-ordinate basis θ, (θ)g and g(θ) can be written in the forms θ i j = p θi u j + u i θs j + g ik θs jk, (101) p θ = p i θ i, p θj = g jk p k θ, θ ij = θ k i g kj = u i p j θ + θ s i u j + θ S ij, (102) θ ij = g ik θ j k = u i p θj + θ s i u j + θ S ij. (103) The symmetric energy-momentum tensor of Belinfante s T can be represented by means of the projective metrics h u, h u and the momentum density p T in the form st = u g(p T ) + T s g(u) + ( T S)g = = s T β α e β e α = s T j i j dx i (104), p T = ρ T u + T π = ρ T u + T π L Kr π, (105) ( s T)g = u p T + T s u + T S = s T αβ e α e β = s T ij i j, (106) g( s T) = g(u) g(p T ) + g( T s) g(u) + g( T S)g = = s T αβ e α e β = s T ij dx i dx j, (107) g(u, p T ) = ρ T e, g(u, T π) = 0, h u (p T ) = h u ( T π), (108) 14

15 h u (h u ) = h u (g) = g(h u ), h u ( G π) = ( G π)h u = g( G π) = ( G π)g, (109) h u (p G ) = g( G π), G π = g[h u (p G )], (110) g( T π) = h u (p T ), T π = g[h u (p T )]. In a co-ordinate basis s T, ( s T)g and g( s T) can be written in the forms st i j = p Ti u j + u i Ts j + g ik TS jk, (111) p T = p i T i, p Tj = g jk p k T, st ij = s T k i g kj = u i p j T + T s i u j + T S ij, (112) st ij = g ik st j k = u i p Tj + T s i u j + T S ij. (113) The variational energy-momentum tensor of Euler-Lagrange Q can be represented by means of the projective metrics h u, h u and the momentum density p Q in the forms Q = u g(p Q ) Q s g(u) ( Q S)g = = Q α β e β e α = Q i j j dx i, (114) (Q)g = u p Q Q s u Q S = Q αβ e α e β = Q ij i j, (115) g(q) = g(u) g(p Q ) g( Q s) g(u) g( Q S)g = = Q αβ e α e β = Q ij dx i dx j, h u ( Q π) = g( Q π), g(u, Q π) = 0, h u (p Q ) = g( Q π), p Q = ρ Q u + Q π, Q π = g[h u (p Q )], h u ( Q s) = g( Q s), h u (h u ) = h u (h u ) = h u (g) = g(h u ). (116) (117) In a co-ordinate basis Q, (Q)g and g(q) will have the forms Q i j = p Qi u j u i Qs j g ik QS jk, (118) Q ij = Q k i g kj = u i p j Q Q s i u j Q S ij, (119) Q ij = g ik Q j k = u i p Qj Q s i u j Q S ij, (120) 15

16 p Qi = g ik p k Q, Q s i = g il Q s l, Q S ij = g ik Q S kl g lj. On the analogy of the notions in the continuum media mechanics one can introduce the following definitions connected with the notion momentum density of a given energy-momentum tensor: A. Generalized canonical energy-momentum tensor θ (a) Momentum density of θ p θ = 1 e [g(u)](θ)(g) = ρ θ u + θ π. (b) Convective momentum density of θ... ρ θ u. B. Symmetric energy-momentum tensor of Belinfante s T (a) Momentum density of s T p T = 1 e [g(u)]( st)(g) = ρ T u + T π. u. (b) Convective momentum density of s T... ρ T C. Variational energy-momentum tensor of Euler-Lagrange Q (a) Momentum density of Q p Q = 1 e [g(u)](q)(g) = ρ Q u Q π. (b) Convective momentum density of Q... ρ Q u. From the covariant Noether identities for the rest mass and the conductive momentum densities, it follows the identity for the momentum density of the different energy-momentum tensors p θ p T p Q. (121) If the Euler-Lagrangian equations δ v L/δV A B = 0 for the non-metric tensor fields are fulfilled and g Q = 0, then p θ = p T. (122) The representations of the energy-momentum tensors of the type 2 is analogous to the representations of the energy-momentum tensors of the type 1 by the use of h u, h u and the momentum density. 4 Energy-momentum tensors and the energy flux density In the first section the notion of energy flux density has been introduced for a given energy-momentum tensor G as e G = G(u) = e (ρ G u + G s). 16

17 The energy-momentum tensors can be now represented by the use of the projective metrics h u, h u and the energy flux density in the forms G = 1 e e G g(u) + u g( G π) + ( G S)g = = G α β e β e α = G i j j dx i, (G)g = 1 e e G u + u G π + G S = G αβ e α e β = G ij i j, (123) g(g) = 1 e g(e G) g(u) + g(u) g( G π) + g( G S)g = = G αβ e α e β = G ij dx i dx j. (124) The Kronecker tensor Kr can also be represented by the use of e Kr, e Kr = Kr(u) = k u + e Kr s = e ( 1 e k u + Kr s). (125) The structure of e G allows the introduction of the abbreviations: (a) Convective energy flux density of the energy-momentum tensor G...ρ G e u. (b) Conductive energy flux density of the energy-momentum tensor G...e. s. The generalized canonical energy-momentum tensor θ can be represented in the form θ = 1 e e θ g(u) + u g( θ π) + ( θ S)g = = θ β α e β e α = θ j i j dx i (126), or by means of the forms e θ = θ(u) = e (ρ θ.u + θ s) = e α θ e α = e i θ i, (θ)g = 1 e e θ u + u θ π + θ S = θ αβ e α e β = θ ij i j, (127) g(θ) = 1 e g(e θ) g(u) + g(u) g( θ π) + g( θ S)g = = θ αβ e α e β = θ ij dx i dx j. In a co-ordinate basis θ, (θ)g and g(θ) will have the forms (128) θ i j = 1 e u i e j θ + θ π i u j + g ik θs jk, (129) θ ij = θ k i g kj = 1 e ei θ u j + u i θπ j + θ S ij, (130) θ ij = g ik θ j k = 1 e e θi u j + u i θπ j + θ S ij, (131) e θi = g ik e k θ, θ S ij = g ik θs kl.g lj. 17

18 The symmetric energy-momentum of Belinfante s T can be represented by the use of e T in the form or by means of the forms st = 1 e e T g(u) + u g( T π) + ( T S)g = = s T α β e β e α = s T i j j dx i, (132) e T = s T(u) = e (ρ T u + T s) = e α T e α = e i T i, (133) ( s T)g = 1 e e T u + u T π + T S = s T αβ e α e β = s T ij i j, (134) g( s T) = 1 e g(e T) g(u) + g(u) g( T π) + g( T S)g = = s T αβ e α e β = s T ij dx i dx j. In a co-ordinate basis s T, ( s T)g and g( s T) will have the forms (135) st i j = 1 e u i e j T + T π i u j + g ik TS jk, (136) st ij = s T k i g kj = 1 e ei T u j + u i Tπ j + T S ij, (137) st ij = g ik st j k = 1 e e Ti u j + u i Tπ j + T S ij, (138) e Ti = g ik e k T, T S ij = g ik TS kl g lj. The variational energy-momentum tensor of Euler-Lagrange Q can be represented in an analogous way by the use of the energy flux density e Q in the forms or by means of the forms Q = 1 e e Q g(u) u g( Q π) ( Q S)g = = Q α β e β e α = Q i j j dx i, (139) e Q = Q(u) = e (ρ Q u + Q s) = e α Q e α = e i Q i, (140) (Q)g = 1 e e Q u u Q π Q S = Q αβ e α e β = Q ij i j, (141) g(q) = 1 e g(e Q) g(u) g(u) g( Q π) g( Q S)g = = Q αβ e α e β = Q ij dx i dx j. In a co-ordinate basis Q, (Q)g and g(q) will have the forms (142) Q i j = 1 e u i e j Q Q π i u j g ik QS jk, 18

19 Q ij = Q k i g kj = 1 e ei Q uj u i Q π j Q S ij, Q ij = g ik Q j k = 1 e e Qi u j u i Qπ j Q S ij, e Qi = g ik e k Q, Q S ij = g ik QS kl g lj. On the analogy of the determined notions, the following abbreviations can be introduced for a given energy-momentum tensor: A. Generalized canonical energy-momentum tensor θ. (a) Energy flux density of θ...e θ. (b) Convective energy flux density of θ...ρ θ e u. B. Symmetric energy-momentum tensor of Belinfante s T. (a) Energy flux density of s T...e T. (b) Convective energy flux density of s T...ρ T e u. C. Variational energy-momentum tensor of Euler-Lagrange Q. (a) Energy flux density of Q...e Q. (b) Convective energy flux density of Q...ρ Q e u. By means of the covariant Noether identities for the rest mass density and for the conductive energy flux density the identity for the energy flux density follows in the form e θ e T e Q. (143) The relation between the momentum density and the energy flux density follows from the structure of their definitions p G = ρ G u + G π, e G = (ρ G u + G s) e, G (θ, s T, Q). From the last expressions the relations p G = 1 e e G + G π G s, (144) e G = (p G + G s G π) e (145) follow. By the use of the different representations of the energy-momentum tensors the different physical processes can be investigated. The application of an representation will depend on the role of the considered quantity (rest mass density, momentum density or energy flux density) in the dynamical process. The physical interpretation of the introduced notions has been used for describing Lagrangian systems in V 4 -spaces, the contravariant vector field u has been interpreted as a time-like vector field tangential to the trajectories of the moving in the space-time particles. A V 4 -space is considered as a model of the space-time. 19

20 5 Covariant divergency of a mixed tensor field The operation of the covariant differentiation along a contravariant vector field can be extended to covariant differentiation along a contravariant tensor field. The Lie derivative ξ u of a contravariant vector field u along a contravariant vector field ξ can be expressed by the use of the covariant differential operators ξ and u in the form ξ u = ξ u u ξ T(ξ, u), ξ, u T(M), T(ξ, u) is the contravariant torsion vector field T(ξ, u) = T αβ γ ξ α u β e γ = T ij k ξ i u j k, constructed by means of the components T αβ γ (or T ij k ) of the contravariant torsion tensor field T. The Lie derivative ξ V of a contravariant tensor field V = V A e A = V A A m (M) along a contravariant vector field ξ can be written on the analogy of the relation for ξ u and by the use of the covariant differential operator ξ and an operator V in the form ξ V = ξ V V ξ T(ξ, V ), ξ T(M), V m (M), V ξ = ξ α /β S Bα Aβ V B e A = ξ i ;j S Ci Aj V C A, T(ξ, V ) = T Bγ A ξ γ V B e A = T Ck A ξ k V C A, T Bγ A = T βγ α S Bα Aβ, T Ck A = T jk i S Ci Aj. V ξ appears as a definition of the action of the operator V on the vector field ξ. Let we now consider more closely this operator and its properties. Let a mixed tensor field K k l(m) be given in a non-co-ordinate (or co-ordinate) basis K = K C D e C e D = K C1α D e C1 e α e D, e C = e C1 e α, e D = e α1... e α l. The action of the operator V on the mixed tensor field K can be defined on the analogy of the action of V on a contravariant vector field ξ V K = K C1α D/β S Bα Aβ V B e C1 e D e A, (146) V is the covariant differential operator along a contravariant tensor field V V : K V K, K k l(m), V m (M), V K k 1+m l(m). Remark. There is an other possibility for a generalization of the action of V on a mixed tensor field V K = k m=1 Kα1...αm...α k D/β.S Bαm Aβ.V B.e α1... e αm 1 e αm+1... e αk e D e A. 20

21 The result of the action of this operator on a contravariant vector field ξ is identical with the action of the above defined operator V. Remark. A covariant differential operator along a contravariant tensor field can also be defined through its action on mixed tensor fields in the form V K = K C D/β.V A1β.e C e D e A1. The operator V differs from V in its action on a contravariant vector field and does not appear as a generalization of the already defined operator by its action on a contravariant vector field. It appears as a new differential operator acting on mixed tensor fields. The covariant differential operator V has the properties: (a) Linear operator V (α K 1 + β K 2 ) = α V K 1 + β V K 2, α, β R (or C), K 1, K 2 k l(m). The proof of this property follows immediately from (146) and the linear property of the covariant differential operator along a basic contravariant vector field. (b) Differential operator (not obeying the Leibniz rule) V (K S) = eβ K S β + K V S, K = K A B e A e B, eβ K = K A B/β e A e B, S = S C D e C e D = S C1α D e C1 e α e D, S β = S C1α D S Eα Fβ V E e C1 e D e F. The proof of this property follows from the action of the defined in (146) operator V and the properties of the covariant derivative of the product of the components of the tensor fields K and S. If the tensor field V is given as a contravariant metric tensor field g, then the covariant differential operator V (V = g) will have additional properties connected with the properties of the contravariant metric tensor field. Definition 2 Contravariant metric differential operator g. Covariant differential operator V for V = g. By means of the relations S Bα Aβ g B e A = (g σ α g βκ + g κ α g βσ ).e σ e κ = = (g σ α gβκ + g κ α gβσ ).e σ.e κ, (147) e σ.e κ = 1 2 (e σ e κ + e κ e σ ), K C1α D/β g σ α = KC1σ D/β, (148) 21

22 the action of the contravariant metric differential operator on a mixed tensor field K can be represented in the form g K = (K C1σ D/β g βκ + K C1κ D/β g βσ ) e C1 e D e σ e κ = (K C1σ D/β g βκ + K C1κ D/β g βσ ).e C1 e D e σ.e κ. (149) The properties of the operator g are determined additionally by the properties of the contravariant metric tensor field of second rank: (a) g : K g K, K k l(m), g K k+1 l(m). (b) Linear operator g (α K 1 + β K 2 ) = α g K 1 + β g K 2. (c) Differential operator (not obeying the Leibniz rule) g (K S) = eβ K S β + K g S, K k l(m), S m r(m), K = K A B e A e B, eβ K = K A B/β e A e B, S = S C D e C e D = S C1α D e C1 e α e D, S β = ( S C1σ D g βκ + S C1κ D g βσ ) e C1 e D e σ.e κ. (150) Remark. The definition of g in (149) differs from the definition in [6], g g, i. e. the contravariant metric differential operator is defined in the last case as a special case of the covariant differential operator V for V = g. The notion of covariant divergency of a mixed tensor field has been used in V 4 -spaces for the determination of conditions for the existence of local conserved quantities and in identities of the type of the first covariant Noether identity. Usually, the covariant divergency of a contravariant or mixed tensor field has been given in co-ordinate or non-co-ordinate basis in the form δk = K Aβ B/β e A e B = K Ci D;i C dx D, (151) K Aβ B/β = K Aβ B/γ g γ β, KCi D;i = K Ci D;k g k i. (152) For full anti-symmetric covariant tensor fields (differential forms) the covariant divergency (called also codifferential) δ is defined by means of the Hodge operator, its reverse operator 1 and the external differential operator a d in the form [20], [21] (pp ) δ = 1 a d. (153) Remark. The Hodge operator is constructed by means of the permutation (Levi-Chivita) symbols. It maps a full covariant anti-symmetric tensor of rank (0, p) a p M) Λ p (M) in a full covariant anti-symmetric tensor of rank (0, n p) a n p (M) Λ n p (M), dim M = n, : a A a A, aa Λ p (M), a A Λ n p (M), 22

23 with aa = A [i1...i p].dx i1... dx ip, a A = A [j1...j n p].dx j1... dx jn p, A [j1...j n p] = 1 p!.ε i 1...i pj 1...j n p.a [i1...ip], A [i1...ip] = g i1k1...g ipkp.a [k1...k p], 1 = ( 1) p.(n p).. By the use of the contravariant metric differential operator g, the covariant metric tensor field g and the contraction operator one can introduce the notion of covariant divergency of a mixed tensor field K with finite rank. Definition 3 Covariant divergency δk of a mixed tensor field K δk = 1 2.[ gk]g = K Aβ B/β.e A e B = K Ci D;i. C dx D, K = K Aβ B e A e β e B = K Ci D C i dx D, K k l(m), k 1. δ is called operator of the covariant divergency δ : K δk, K k l(m), δk k 1 l(m), k 1. Remark. The symbol δ has also been introduced for the variation operator. Both operators are different from each other and can easily be distinguished. Ambiguity would occur only if the symbol δ is used out of the context. In such a case, the definition of the symbol δ is necessary. The properties of the operator of the covariant divergency δ are determined by the properties of the contravariant metric differential operator, the contraction operator and the metric tensor fields g and g (a) The operator of the covariant divergency δ is a linear operator δ(α K 1 + β K 2 ) = α δk 1 + β δk 2, α, β R (or C), K 1, K 2 k l(m). (154) The proof of this property follows immediately from the definition of the covariant divergency. (b) Action on a tensor product of tensor fields δ(k S) = S K + K δs, (155) K = K A B e A e B, S = S Cβ D e C e β e D, S K = K A B/β S Cβ D e A e B e C e D (see above V ). (156) The proof of this property follows from the properties of g and from the relations 1 2 [ e β K S β ]g = K A B/β S Cβ D e A e B e C e D, (157) 23

24 K 1 2 [K gs]g = K 1 2 [ gs]g = K δs. (158) (c) Action on a contravariant vector field u δu = 1 2 [ gu]g = u β /β = u i ;i. (159) (d) Action on the tensor product of two contravariant vector fields u and v δ(u v) = v u + δu v, v u = v u. (160) (e) Action of the product of an invariant function L and a mixed tensor field δ(l.k) = K L + L δk, (161) K L = L /β K Aβ B e A e B, δk = K Aβ B/β e A e B, L /β = e β L, L ;i = L,i, K = K Aβ B e A e β e B k l(m). (162) Special case: Action of the product of an invariant function L and the contravariant metric tensor g: δ(l.g) = (L /β g αβ + L g αβ /β) = (L,j g ij + L g ij ;j). (163) (f) Action on an anti-symmetric tensor product of two contravariant vector fields u and v δ(u v) = 1 2 ( vu u v + δv u δu v) = = 1 2 [ uv + T(u, v) + δu v δv u]. (164) (g) Action on a full anti-symmetric contravariant tensor field A of second rank δa = 1 2 (Aαβ A βα ) /β e α = 1 2 (Aij A ji ) ;j i, A = A αβ e α e β = A ij i j, A αβ = A βα (165). (h) Action on a tensor product of a contravariant vector field u, multiplied with an invariant function, and a covariant vector field g(v) with the contravariant vector field v δ(ε u g(v)) = (uε) g(v) + ε [δu g(v) + ( u g)(v) + g( u v)] = = [uε + ε δu] g(v) + ε [( u g)(v) + g( u v)], ε C r (M), ε (x k ) = ε(x k ), u, v T(M). (166) Special case: v u: δ(ε u g(u)) = [uε + ε δu] g(u) + ε [( u g)(u) + g(a)], u u = a. (167) Special case: ε = 1: δ(u g(v)) = δu g(v) + ( u g)(v) + g( u v). (168) 24

25 5.1 Covariant divergency of a mixed tensor field of second rank From the definition of the covariant divergency δk of a mixed tensor field K, the explicit form of the covariant divergency of tensor fields of second rank of the type 1. or 2. follows as δg = 1 2 [ gg]g = G α β /β e α = G i j ;j dx i, (169) δg = 1 2 [ gg]g = G β α/β e α = G j i;j dx i. (170) By the use of the relations (159) (164), (166) (168), and the expression [see (155) (158)] v (g(u)) = v (g(u)) = ( v g)(u) + g( v u), (171) δ(g(u) v) = δv g(u) + ( v g)(u) + g( v u), (172) δ(( G S)g) = (g αγ GS βγ ) /β e α, (173) the covariant divergency of the representation of G by means the rest mass density ρ G (ε G = ρ G ) G = ρ G u g(u) + u g( G π) + G s g(u) + ( G S)g, can be found in the form ( u u = a) δg = ρ G g(a) + (uρ G + ρ G δu + δ G s) g(u) + δu g( G π) + g( u G π)+ + g( G su) + ρ G ( u g)(u) + ( u g)( G π) + ( G sg)(u) + δ(( G S)g). (174) g(δg) will have the form g(δg) = ρ G a + (uρ G + ρ G δu + δ G s) u + δu Gπ + u G π+ + G su + ρ G g( u g)(u) + g( u g)( G π) + g( G sg)(u) + g(δ(( G S)g)). (175) In a co-ordinate basis δg and g(δg) will have the forms j G i ;j = ρ G a i + (ρ G,j u j + ρ G u j ;j + G s j ;j) u i + u j ;j Gπ i + + g ij ( G π j ;k u k + u j ;k Gs k ) + g ij;k (ρ G u k u j + u k Gπ j + G s k u j )+ + (g ik G S jk ) ;j, (176) a i = g ij a j, a i = u i ;j u j, ρ G;i = ρ G,i, u i = g ik u k, G π i = g ij G π j, g ik G k j ;j = ρ G a i + (ρ G,j u j + ρ G u j ;j + G s j ;j) u i + u j ;j Gπ i + + G π i ;j u j + u i ;j Gs j + g il g lj;k (ρ G u k u j + u k Gπ j + G s k u j )+ + g il (g lk GS jk ) ;j. (177) 25

26 The relation between δg and δg follows from the relation between G and G G = g(g)g : δg = δ(g(g)g) = 1 2 [ g(g(g)g)]. (178) The covariant divergency of the Kronecker tensor field can be found in an analogous way as the covariant divergency of a tensor field of second rank of the type 1, since Kr = g α β e α e β = g i j i dx j δkr = 1 2 [ gkr]g = gα β /β e α = g j i;j dxi. (179) If we use the representation of Kr Kr = 1 e k u g(u) + u g(kr π) + Kr s g(u) + ( Kr S)g, then the covariant divergency δkr can be written in the form δkr = 1 e k g(a) + [u(1 e k) + 1 e k δu + δkr s] g(u)+ +δu g( Kr π) + g( u Kr π) + g( Kr su)+ + 1 e k ( ug)(u) + ( u g)( Kr π) + ( Kr sg)(u) + δ(( Kr S)g), (180) or in the form g(δkr) = 1 e k a + [u(1 e k) + 1 e k δu + δkr s] u+ + δu Kr π + u Kr π + Kr su+ + 1 e k g( ug)(u) + g( u g)( Kr π) + g( Kr sg)(u) + g(δ(( Kr S)g)). (181) In a co-ordinate basis δkr and g(δkr) will have the forms g j i ;j = 1 e k a i + [( 1 e k),j u j + 1 e k uj ;j + Kr s j ;j] u i + + u j ;j Kr π i + g ij ( Kr π j ;k u k + u j ;k Kr s k )+ + g ij;k ( 1 e k uj u k + Kr π j u k + u j Kr s k ) + (g ik Kr S jk ) ;j, (182) g ik g j k ;j = 1 e k ai + [( 1 e k),j u j + 1 e k uj ;j + Kr s j ;j] u i + + u j ;j Kr π i + Kr π i ;j u j + u i ;j Kr s j + + g il g lj;k ( 1 e k uj u k + Kr π j u k + u j Kr s k ) + g il (g lk Kr S jk ) ;j. (183) 6 Covariant divergency of the energy-momentum tensors and the rest mass density The covariant divergency of the energy-momentum tensors can be represented by the use of the projective metrics h u, h u of the contravariant vector field u and the rest mass density for the corresponding energy-momentum tensor. In this case the representation of the energy-momentum tensor is in the form G = (ρ G + 1 e L k) u g(u) L Kr + u g(k π) + k s g(u) + ( k S)g, 26

27 k π = G π, k s = G s, k S = G S. (184) By the use of the relations u(ρ G + 1 e L k) = uρ G + L u( 1 e k) + 1 e k (ul) = = [ρ G/α + L ( 1 e k) /α + 1 e k L /α] u α = (185) = [ρ G,j + L ( 1 e k),j + 1 e k L,j] u j, KrL = L /α e α = L,i dx i = Kr L, (186) δ(l Kr) = KrL + L δkr, (187) δ(l Kr) = 1 2 [ g(l.kr)]g = (L.g β α ) /β e α = (L g j i ) ;j dx i = = (L /β g β α + L gβ α /β) e α = (L,i + L g j i ;j) dx i, (188) δ(u g( G π)) = δu g( G π) + g( u G π) + ( u g)( G π), (189) δ( G s g(u)) = δ G s g(u) + g( G su) + ( G sg)(u), (190) δ(( G S)g) = (g αγ GS βγ ) /β e α [see (173)], δg and g(δg) can be found in the forms δg = (ρ G + 1 e L k) g(a)+ +[u(ρ G + 1 e L k) + (ρ G + 1 e L k) δu + δg s] g(u) KrL L δkr + δu g( G π) + g( G u π) + g( G su)+ + (ρ G + 1 e L k) ( ug)(u) + ( u g)( G π) + ( G sg)(u)+ + δ(( G S)g), g(δg) = (ρ G + 1 e L k) a+ +[u(ρ G + 1 e L k) + (ρ G + 1 e L k) δu + δg s] u g(krl) L.g(δKr) + δu Gπ + G u π + su+ G + (ρ G + 1 e L k) g( ug)(u) + g( u g)( G π) + sg)(u)+ g( G + g(δ(( G S)g)). In a co-ordinate basis δg and g(δg) will have the forms G i j ;j = (ρ G + 1 e L k) a i+ + [(ρ G + 1 e L k),j u j + (ρ G + 1 e L k) uj ;j + G s j ;j] u i L,i L g j i;j + uj ;j Gπ i + g ij ( G π j ;k u k + u j ;k Gs k )+ + g ij;k [(ρ G + 1 e L k) uj u k + G π j u k + u j Gs k ]+ + (g ik GS jk ) ;j, g ik G k j ;j = (ρ G + 1 e L k) ai + + [(ρ G + 1 e L k),j u j + (ρ G + 1 e L k) uj ;j + G s j ;j] u i L,j g ij L g ik g j k ;j + uj ;j Gπ i + G π i ;j u j + u i ;j Gs j + + g il g lj;k [(ρ G + 1 e L k) uj u k + G π j u k + u j Gs k ]+ + g il (g lk GS jk ) ;j. (191) (192) (193) (194) 27

28 On the grounds of the representations of δg and g(δg) the representation of the different energy-momentum tensors θ, s T and Q can be found. The covariant divergency δθ of the generalized canonical energy-momentum tensor (GC-EMT) θ can be written in the form δθ = (ρ θ + 1 e L k) g(a)+ +[u(ρ θ + 1 e L k) + (ρ θ + 1 e L k) δu + δθ s] g(u) KrL L.δKr + δu g( θ π) + g( θ u π) + su)+ g( θ + (ρ θ + 1 e L k) ( ug)(u) + ( u g)( θ π) + sg)(u)+ ( θ + δ(( θ S)g), (195) or in the form g(δθ) = (ρ θ + 1 e L k) a+ +[u(ρ θ + 1 e L k) + (ρ θ + 1 e L k) δu + δθ s] u g(krl) L g(δkr) + δu θπ + θ u π + θ su+ + (ρ θ + 1 e L k).g( ug)(u) + g( u g)( θ π) + g( θ sg)(u)+ + g(δ(( θ S)g)). (196) In a co-ordinate basis δθ and g(δθ) will have the forms θ i j ;j = (ρ θ + 1 e L k) a i+ + [(ρ θ + 1 e L k),j u j + (ρ θ + 1 e L k) uj ;j + θ s j ;j] u i L,i L g j i;j + uj ;j θπ i + g ij ( θ π j ;k u k + u j ;k θs k )+ + g ij;k [(ρ θ + 1 e L k) uj u k + θ π j u k + u j θs k ]+ + (g ik θs jk ) ;j, (197) g ik θ k j ;j = (ρ θ + 1 e L k) ai + + [(ρ θ + 1 e L k),j u j + (ρ θ + 1 e L k) uj ;j + θ s j ;j] u i L,j g ij L g ik g j k ;j + uj ;j θπ i + θ π i ;j u j + u i ;j θs j + + g il g lj;k [(ρ θ + 1 e L k) uj u k + θ π j u k + u j θs k ]+ + g il (g lk θs jk ) ;j. (198) The covariant divergency δ s T of the symmetric energy-momentum tensor of Belinfante (S-EMT-B) s T, represented in the form st = (ρ T + 1 e L k).u g(u) L Kr + u g(t π) + T s g(u) + ( T S)g, can be found in the form δ s T = (ρ T + 1 e L k).g(a)+ +[u(ρ T + 1 e L k) + (ρ T + 1 e L k) δu + δt s] g(u) KrL L δkr + δu g( T π) + g( u T π) + g( T su)+ + (ρ T + 1 e L k) ( ug)(u) + ( u g)( T π) + ( T sg)(u)+ + δ(( T S)g), (199) 28