# Landau-Lifschitz gravitational energy-momentum pseudotensor

## Energy conservation principle in general relativity

Two of the most popular definitions of the combined Energy and Momentum Conservation Principle (EMCP) are:

- the energy and momentum of a closed isolated system are conserved
- the stress-energy-momentum (SEM) tensor of a closed isolated system is a constant quantity

The second definition is preferred in theoretical physics. The *SEM tensor*, *T _{i}^{k}*, which combines energy, mass, momentum, their fluxes and stresses (also called "super energy-momentum tensor",

*etc.*), is the

**content**of the Universe. It is both the matter and the electromagnetic fields. EMCP requires that the SEM tensor is constant. Imagine the closed isolated system as a volume (Universe, black hole), enclosed by an impenetrable envelope such that no matter or radiation can pass through. The SEM tensor within the enclosed volume is constant: there is no tensor flow through the envelope. Tensor flow, or tensor divergence, is the sum of partial derivatives of tensor components with respect to coordinates (similar to vector divergence). Divergence should be zero according to EMCP (LL32,4):

∂T_{i}^{k}/ ∂x^{ k} ≡ T_{i}^{k}_{,k} = 0 (EMCP).

Presence of gravitation spoils this clear picture in several ways. Gravitation curves space. In curved space, tensor (and vector) flows are ill-defined because they change direction in each point. Parallel transfer, which is essential for definition of derivatives, is very complicated on curved paths. To cut it short, derivatives are obtained only after making corrections for the curvature – one correction for each tensor index. These corrected derivatives are known as *covariant derivatives*. [1] In curved space, ordinary derivatives are replaced by covariant ones and (EMCP) becomes

∂T_{i}^{k}/ ∂x^{ k} − Γ^{ l}_{ik}T_{l}^{k} + Γ^{ k}_{kl}T_{i}^{l} ≡
T_{i}^{k}_{;k} = 0. (85,13) [2].

The *Γ ^{l}_{ik}* terms above are the corrections for spacetime curvature: a negative correction for the lower index in

*T*and a positive correction for the upper index.

_{i}^{k}*T*is a symmetric tensor and according to 86,11

_{i}^{k}T_{i}^{k}_{;k} = (√−g)^{ −1} (T_{i}^{k}√−g)_{ , k} − ½g_{kl,i}T^{kl} (96,1).

The first term in the right hand side can be made an ordinary derivative of the SEM tensor by pulling out the constant *√−g* from inside the differential operator. However, the real difficulty lies in the second term *−½g _{kl,i}T^{kl}*. The metric tensor

*g*and its derivatives can be anything depending on the spacetime topology. Taking up the above analogy, the volume that contains the SEM tensor itself stretches, shrinks, and moves in every imaginable way. It cannot be bound in an envelope with zero volume. The spacetime itself generates tensor flows. Do what you want, you can't just ignore or remove the second term. Failure to do so means nothing less than the startling conclusion:

_{kl}*in the curved spacetime of general relativity, the EMCP is no longer valid.*Whether this is real or apparent remains subject to interpretation. It is easy to see, however, that this result can be expected from what we know from elementary physics. In the heart of the concept of the energy conservation principle is the assumption that space is isotropic while the momentum conservation principle requires symmetric space. Both isotropy and symmetry are no longer assumed in curved spacetime. The energy and momentum conservation principles fail in principle. The very basis on which they are built sinks into a quicksand. Is there a way to salvage something?

## In search of symmetry

Many efforts have been and are being made to sidetrack this important unsolvable problem. Landau-Lifshitz turned their attention to a couple of possible loopholes.

- the immediate neighborhood of every point is locally flat; this means that the tiny volume around this point can be enclosed in an envelope into which EMCP is valid
- the identity EMCP does not fix the SEM tensor exactly but leaves a quantity undefined. Namely, according to LL(32,7), one can always find an antisymmetric tensor
*η*such that if^{ikl}= −η^{ilk}*T*then^{ ik}= &eta^{ikl}_{,l}*T*≡ 0. [3]^{ik}_{,k}= η^{ikl}_{,k,l}

The possibility to locally fulfil EMCP gives the idea to find a SEM tensor that locally has a zero divergence and then correct it to make it valid globally, in curved spacetime. An inkling to the nature of this correction may be taken from the analogous case with the covariant derivatives. There all corrections contain Christoffel symbols which are composed from first derivatives of the metric tensor. The way to make the required SEM tensor is already obvious: it must be such antisymmetric expression that upon differention becomes a second derivative which is identically zero. When constructing it, we must throw out all first metric derivatives that would prevent the expression from turning to zero.

Landau-Lifshitz start from the Einstein field equations to put *T ^{ik}* into the required form

*T ^{ik} =
η^{ikl}_{,l}* (etaform).

The field equations in contravariant form are

T^{ik}
= (c^{4} / 8πk) (R^{ik} −
½g^{ik}R). (95,5)

All tensors in 95,5 are decomposed in metric derivatives and only the second derivatives are taken. The detailed transformation of the right hand side of the last expression (explanation of each step is provided in notes) is:

(16πk / c^{4}) T^{ik}
note 1 = 2 (R^{ik} −
½g^{ik}R) = 2 (g^{im}g^{kp}g^{ln} − ½g^{ik}g^{ln}g^{mp}) R_{lmnp}
note 2 ≈ (g^{im}g^{kp}g^{ln} − ½g^{ik}g^{ln}g^{mp})(g_{lp,m,n} + g_{mn,l,p} − g_{ln,m,p} − g_{mp,l,n})
note 3 = g^{im}g^{kp}g^{ln}g_{lp,m,n} + g^{im}g^{kp}g^{ln}g_{mn,l,p
} − g^{im}g^{kp}g^{ln}g_{ln,m,p} −
g^{im}g^{kp}g^{ln}g_{mp,l,n
} −
½ g^{ik}g^{ln}g^{mp}g_{lp,m,n
} +
½ g^{ik}g^{ln}g^{mp}g_{mn,l,p} −
½ g^{ik}g^{ln}g^{mp}g_{ln,m,p
} −
½ g^{ik}g^{ln}g^{mp}g_{mp,l,n}
note 4 = g^{im}g^{kp}g^{ln}g_{np,m,l} +
g^{ip}g^{km}g^{ln}g_{np,m,l} −
g^{im}g^{kl}g^{np}g_{np,m,l} −
g^{in}g^{kp}g^{lm}g_{np,m,l} −
½ g^{ik}g^{ln}g^{mp}g_{np,m,l} +
½
g^{ik}g^{ln}g^{mp}g_{np,m,l} −
½ g^{ik}g^{np}g^{ml}g_{np,m,l} −
½ g^{ik}g^{lm}g^{np}g_{np,m,l}
note 5 = g_{np,m,l} (g^{im}g^{kp}g^{ln
} +
g^{ip}g^{km}g^{ln
} −
g^{im}g^{kl}g^{np} − g^{in}g^{kp}g^{lm} −
½g^{ik}g^{ln}g^{mp} − ½g^{ik}g^{ln}g^{mp
} +
½g^{ik}g^{np}g^{ml} + ½g^{ik}g^{lm}g^{np})
note 6 = g_{np,m,l} (g^{im}g^{kp}g^{ln} + g^{ip}g^{km}g^{ln} −
g^{im}g^{kl}g^{np} − g^{in}g^{kp}g^{lm} − g^{ik}g^{ln}g^{mp}
+
g^{ik}g^{np}g^{ml})
note 7 = [g_{np,m}(g^{im}g^{kp}g^{ln} + g^{ip}g^{km}g^{ln} −
g^{im}g^{kl}g^{np} − g^{in}g^{kp}g^{lm} −
g^{ik}g^{ln}g^{mp} + g^{ik}g^{np}g^{ml})]_{ , l} −
g_{np,m} (g^{im}g^{kp}g^{ln} +
g^{ip}g^{km}g^{ln} −
g^{im}g^{kl}g^{np} −
g^{in}g^{kp}g^{lm} −
g^{ik}g^{ln}g^{mp}
+ g^{ik}g^{np}g^{ml})_{,l}
note 8 ≈ [g_{np,m}(g^{im}g^{kp}g^{ln} +
g^{ip}g^{km}g^{ln} −
g^{im}g^{kl}g^{np} −
g^{in}g^{kp}g^{lm} −
g^{ik}g^{ln}g^{mp} +
g^{ik}g^{np}g^{ml})]_{,l
}
note 9 = (g^{im}g^{kp}g^{ln}g_{np}_{,m} +
g^{ip}g^{km}g^{ln}g_{np}_{,m} −
g^{im}g^{kl}g^{np}g_{np}_{,m} −
g^{in}g^{kp}g^{lm}g_{np}_{,m} −
g^{ik}g^{ln}g^{mp}g_{np}_{,m} +
g^{ik}g^{np}g^{ml}g_{np}_{,m})_{,l}
note 10 = [−g^{im}g^{kp}g_{np}g^{ln}_{,m} −
g^{ip}g^{km}g_{np}g^{ln}_{,m} −
g^{im}g^{kl}(g_{,m}/g) +
g^{in}g^{lm}g_{np}g^{kp}_{,m} +
g^{ik}g^{ln}g_{np}g^{mp}_{,m} +
g^{ik}g^{ml}(g_{,m}/g)]_{,l}
note 11 = [−δ^{k}_{n}g^{im}g^{ln}_{,m} − δ^{i}_{n}g^{km}g^{ln}_{,m
} −
g^{im}g^{kl}(g_{,m}/g) + δ^{i}_{p}g^{lm}g^{kp}_{,m} +
δ^{l}_{p}g^{ik}g^{mp}_{,m} +
g^{ik}g^{ml}(g_{,m}/g)]_{,l}
note 12 = [− g^{im}g^{kl}_{,m} −
g^{km}g^{il}_{,m} −
g^{im}g^{kl}(g_{,m}/g) +
g^{lm}g^{ik}_{,m} +
g^{ik}g^{lm}_{,m} +
g^{ik}g^{ml}(g_{,m}/g)]_{,l}
note 13 = [−g^{i}^{l}g^{k}^{m}_{,m} −
g^{km}g^{il}_{,m} − g^{i}^{l}g^{k}^{m}(g_{,m}/g) +
g^{lm}g^{ik}_{,m} +
g^{ik}g^{lm}_{,m} +
g^{ik}g^{lm}(g_{,m}/g)]_{,l}
note 14 = [(g^{lm}g^{ik}_{,m} +
g^{ik}g^{lm}_{,m})
− (g^{il}g^{km}_{,m} +
g^{km}g^{il}_{,m}) +
(g_{,m}/g)(g^{ik}g^{lm} − g^{il}g^{km})]_{,l}
note 15 = [(g^{ik}g^{lm})_{,m} −
(g^{il}g^{km})_{,m
} + (g_{,m}/g)(g^{ik}g^{lm} −
g^{il}g^{km})]_{,l}
note 16 = [(g^{ik}g^{lm} −
g^{il}g^{km})_{,m} +
(g_{,m}/g)(g^{ik}g^{lm} −
g^{il}g^{km})]_{, l}
note 17 = [(−g/− g)(g^{ik}g^{lm} −
g^{il}g^{km})_{,m
} +
(−1/−1)(g_{,m}/g)(g^{ik}g^{lm} −
g^{il}g^{km})]_{, l}
note 18 = {1/(−g)[(−g)(g^{ik}g^{lm} −
g^{il}g^{km})_{,m
} +
(−g)_{,m}(g^{ik}g^{lm} − g^{il}g^{km})]},_{l }
note 19 = {1/(−g)[(−g)(g^{ik}g^{lm} −
g^{il}g^{km})]_{,m}}_{,l}
note 20 = 1/(−g)[(−g)(g^{ik}g^{lm} −
g^{il}g^{km})]_{,m,l}
note 21.

**Note 1.** Both sides are multiplied by *16πk/c ^{4}* to get rid of the constant in the right hand side. The factor 2 on the right should subsequently neutralise the ½ in the decomposition.

**Note 2.** The Ricci tensor is

*R ^{ik} = g^{im}g^{kp}g^{ln}R_{lmnp}* (92,6).

The term with the scalar curvature, *g ^{ik}R*, is

g^{ik}R
= g^{ik}g^{ln}g^{mp}R_{lmnp} (92,9).

**Note 3.** The Riemann tensor, *R _{lmnp}*, in which all first metric derivatives in Christoffel symbols are ignored, is approximated by

R_{lmnp} ≈ ½ (g_{lp,m,n} + g_{mn,l,p} − g_{ln,m,p} − g_{mp,l,n}) (92,1).

**Note 4.** All parentheses are opened resulting in 2 x 4 = 8 terms.

**Note 5.** Pairs of dummies are exchanged so as to make all secondary derivatives, *g _{np,m,l}* with uniform indexes.

**Note 6.** The second derivatives *g _{np,m,l}* are collected from all terms. Note that the last four terms are pairwise equal when we take into account that the metric tensor is symmetric, e.g.,

*g*. The pairwise equal terms are marked in different colors: red and blue.

^{ml}= g^{lm}**Note 7.** The last four terms that are pairwise equal are added; now we have a second derivative of the metric multiplied by a combination of 6 metric terms that do not contain the factor ½.

**Note 8.** Let *g _{np} = u* and

*g*. Using the formula for a derivative of a product:

^{im}g^{kp}g^{ln}+ g^{ip}g^{km}g^{ln}− g^{im}g^{kl}g^{np}− g^{in}g^{kp}g^{lm}− g^{ik}g^{ln}g^{mp}+ g^{ik}g^{np}g^{ml}= v*u*we get 2 terms: the first (red) includes second derivatives while the second (blue) includes only first derivatives

_{,m,l}v = (u_{,m}v)_{,l}− u_{,m}v_{,l}**Note 9.** The second term is dropped because it contains only first derivatives; the resulting expression is approximation valid in the local point.

**Note 10.** Brackets are expanded by distributing *g _{np,m}* among terms.

**Note 11.** Indexes shown in red are dummies that were exchanged. We have two cases: first, where we have the same upper
and lower dummies, LL(86,7) was used in the form: *g ^{ln}g_{np}_{,m}
= −g^{np}g_{ln}_{,m}* ;
second, where both upper and lower dummies are the same, as in

*g*, the ubiquitous LL(86,4) was used in the form:

^{np}g_{np}_{,m}*g*so that

_{,m}= gg^{np}g_{np,m}*g*

^{np}g_{np,m}= (g_{,m}/g).**Note 12.** Products of metric tensors that contain the same dummy indexes are equal to Kronecker deltas according to LL(83,10): *g ^{kp}g_{np} = δ^{k}_{n}*, and so on.

**Note 13.** Kronecker deltas when multiplied with metric tensors change the indexes of the latter in the following way: *δ ^{k}_{n}g^{im}g^{ln}_{,m} = g^{im}g^{kl}_{,m}*; making these changes in the above expression results in some shortening.

**Note 14.** In some terms, *m* and *l* indexes are exchanged, using the symmetry of the differential operator.

**Note 15.** Terms are grouped conveniently ...

**Note 16.** ... so that the formula for differentiation of a product: *g ^{lm}g^{ik}_{,m} + g^{ik}g^{lm}_{,m} = (g^{ik}g^{lm})_{,m}* can be used to present the terms in the first two parentheses as derivatives of products.

**Note 17.** The first two parentheses are factored by the *,m* differentiation (associative property).

**Note 18.** The first term is multiplied by *(−g/−g)* = 1 and the second term by (−1/−1) = 1

**Note 19.** 1/*(−g)* is collected in front of the brackets. Inside the brackets we can distinguish products of two terms: *(−g) = u* and *(g ^{ik}g^{lm}
− g^{il}g^{km}) = v*.

**Note 20.** The property of a derivative of product *u v _{,m} + u_{,m} v = (u v)_{,m}* is again used to factor the second differentiation operator

*,m*.

**Note 21.** 1/(*−g*) is pulled out of derivatives by *,l* which is possible because 1/(*−g*) is constant in the special point. This means that 1/(*−g*) goes in front of the braces which are then dropped as there is no need for them – the expression that remains in the brackets is doubly differentiated.

To summarize, the final expression for the SEM tensor as a result of the above transformations is:

(16πk/c^{4})T^{ik}
= 1/(−g)[(−g)(g^{ik}g^{lm}
− g^{il}g^{km})]_{,m,l}.

For the following discussion we need this expression somewhat rearranged. Namely, the constant *(c ^{4}/16πk)* is put under

*,l*(but not under

*,m*) and both sides are multiplied by

*(−g)*:

(−g)T^{
ik} = (c^{4}/16πk)[(−g)(g^{ik}g^{lm} − g^{il}g^{km})]_{,m,l} = {(c^{4}/16πk)[(−g)(g^{ik}g^{lm} − g^{il}g^{km})]_{,m}}_{,l}.

After these "simple transformations" (according to Landau-Lifshitz), we obtained a SEM tensor in a completely antisymmetric form:

(−g)T^{ik} = {(c^{4}/16πk)[(−g)(g^{ik}g^{lm} − g^{il}g^{km})]_{,m}}_{,l}.

The quantity in the braces corresponds to *η ^{ikl}* in (etaform):

η^{ikl}
= (c^{4}/16πk)[(−g)(g^{ik}g^{lm} − g^{il}g^{km})]_{,m
} (96,3).

According to [3], *(−g)T ^{ik}* vanishes upon differentiation by

*,k*(check). Furthermore,

*(−g)T*is symmetric:

^{ik}*(−g)T*(check).

^{ik}= (−g)T^{ki}## Landau-Lifshitz pseudotensor as first metric derivatives

The nice and pliable SEM tensor obtained above whose divergence disappears automatically just by dint of symmetry (or, more precisely, antisymmetry) so that it fulfils EMCP, does not work outside the very small flat locality of the point that we chose. Remember that in its derivation above we threw out first derivatives of the metric tensor at a few steps. First derivatives usually combine to form connections (Christoffel symbols) and these last are directly connected to curvature. So we may expect that outside our flat neighbourhood first derivatives will reappear to play a prominent role as corrections for curvature in the same way as they act in covariant derivatives.

We have every reason to expect, then, that *(−g)T ^{ik}* will no longer be equal to

*η*

^{ikl}_{,l}
*η ^{ikl}_{,l} − (−g)T^{ik}* ≠ 0.

There will be some non-zero difference that is the quantity that corrects for the transition from local to global geometry, from flat to curved spacetime. Let's define this difference as *(−g)t ^{ik}*, then by definition

*(−g)(T ^{ik} + t^{ik}) = η^{ikl}_{,l}* (96,5).

We can't tell much about the quantities *t ^{ik}* except that they have absorbed all the first derivatives which we threw out before, and also that they are symmetric:

*t ^{ik} = t^{ ki}* (96,6).

The latter stems from the definition of *t ^{ik}*, since they are composed of

*T*and

^{ik}*η*which are symmetric as we proved in [5].

^{ikl}_{,l}The method to find *t ^{ik}* that Landau-Lifshitz propose, is direct but excessively long and tedious: " ... after a rather lenghty calculation we find the following expression ... " One seldom finds in their book such a phrase at the place of the gory detail of the intermediate calculations though it is an understatement as we see in the derivation of LL96,8 that follows.

To be more specific, we shall articulate the plan of calculations in several steps:

*t*is expressed according to (96,5) as a difference between the Landau-Lifsitz pseudotensor in the special point and the SEM tensor^{ik}- the two differentiations in the Landau-Lifsitz pseudotensor in the special point are carried out
- derivatives and double derivatives of the metric determinant are transformed into derivatives and double derivatives of the metric tensor according to (86,4)
- all double derivatives are written as covariant metric tensors rather than as contravariant ones using when necessary (86,4) in the covariant form
- the SEM tensor is expressed by the Einstein equations and transformed sequentially into Riemann tensor (through (92,6) and (92,9)), Christoffel symbols (through (92,1)), and metric tensor derivatives (through (86,3))
- the SEM tensor in metric derivatives is subtracted from the Landau-Lifshitz pseudotensor in the special point, expressed also in metric derivatives
- the second metric tensor derivatives in the resultant expression are presented in such a manner as to cancel and derive an expression containing only first derivatives of the metric tensor (base form of the Landau-Lifshitz pseudotensor)
- the base form (first metric derivatives) is transformed into (96,8) (Christoffel symbols) using (86,8) and (86,5) or
- the base form (first metric derivatives) is transformed into (96,9) (metric tensor densities) using (86,6) As before, the intermediate calculations are provided with detailed notes.

t^{ik} = 1/(−g)η^{ikl}_{,l} − T^{ik}
note 22 = (c^{4}/16πk){1/(−g)[(−g)(g^{ik}g^{lm} −
g^{il}g^{km})]}_{,m,l} − (c^{4}/8πk)(R^{ik} −
½g^{ik}R)
note 23 = (c^{4}/16πk){(1/g)[g(g^{ik}g^{lm} −
g^{il}g^{km})]_{,m,l} − (2R^{ik} −
g^{ik}R)
note 24 = (c^{4}/16πk){(1/g)[(g^{ik}g^{lm} −
g^{il}g^{km})g_{,m,l} +
g_{,m} (g^{lm}g^{ik}_{,l} − g^{km}g^{il}_{,l} −
g^{il}g^{km}_{,l} + g^{ik}g^{lm}_{,l} ) +
g_{,l} (g^{lm}g^{ik}_{,m} − g^{km}g^{il}_{,m} −
g^{il}g^{km}_{,m} + g^{ik}g^{lm}_{,m} ) + g(−g^{il}_{,m}
g^{km}_{,l} − g^{il}_{,l} g^{km}_{,m} + g^{ik}_{,m} g^{lm}_{,l} +
g^{ik}_{,l} g^{lm}_{,m} + g^{lm}g^{ik}_{,m,l} −
g^{km}g^{il}_{,m,l} &minus g^{il}g^{km}_{,m,l} +
g^{ik}g^{lm}_{,m,l} )] −
(2g^{il}g^{km} − g^{ik}g^{lm})R_{lm}}
note 25
= (c^{4}/16πk){(1/g)[(g^{ik}g^{lm} −
g^{il}g^{km})g_{,m,l} +
g_{,l} (
g^{lm}g^{ik}
_{,m} − g^{kl}g^{im}_{,m} −
g^{im}g^{kl}_{,m} + g^{ik}g^{lm}_{,m} ) +
g_{,l} (g^{lm}g^{ik}_{,m} − g^{km}g^{il}_{,m} −
g^{il}g^{km}_{,m} + g^{ik}g^{lm}_{,m} ) + g(−g^{il}_{,m}
g^{km}_{,l} − g^{il}_{,l} g^{km}_{,m} + 2g^{ik}_{,l} g^{lm}_{,m} + g^{lm}g^{ik}_{,m,l} −
g^{km}g^{il}_{,m,l} &minus g^{il}g^{km}_{,m,l} +
g^{ik}g^{lm}_{,m,l} )] −
(2g^{il}g^{km} − g^{ik}g^{lm})
(Γ^{ n}_{lm,n} − Γ^{ n}_{ln,m} +
Γ^{ n}_{lm}Γ^{ p}_{np} − Γ^{ p}_{ln}Γ^{ n}_{mp})}
note 26
= (c^{4}/16πk){(1/g)[
2g(g^{ik}g^{lm} −
g^{il}g^{km})(2Γ^{ n}_{ln}
Γ^{ p}_{mp} +
Γ^{ n}_{ln,m} ) +
2gΓ^{ n}_{ln} (
2g^{lm}g^{ik}_{,m} − g^{kl}g^{im}_{,m} −
g^{km}g^{il}_{,m} −
g^{im}g^{kl}_{,m} −
g^{il}g^{km}_{,m} +
2g^{ik}g^{lm}_{,m} ) + g(−g^{il}_{,m}
g^{km}_{,l} − g^{il}_{,l} g^{km}_{,m} + 2g^{ik}_{,l} g^{lm}_{,m} + g^{lm}g^{ik}_{,m,l} −
g^{km}g^{il}_{,m,l} &minus g^{il}g^{km}_{,m,l} +
g^{ik}g^{lm}_{,m,l} )] −
(2g^{il}g^{km} − g^{ik}g^{lm})
(Γ^{ n}_{lm,n} − Γ^{ n}_{ln,m} +
Γ^{ n}_{lm}Γ^{ p}_{np} − Γ^{ p}_{ln}Γ^{ n}_{mp})}
note 27
= (c^{4}/16πk)
[−g^{il}_{,m}
g^{km}_{,l} − g^{il}_{,l} g^{km}_{,m} +
2g^{ik}_{,l} g^{lm}_{,m} + g^{lm}g^{ik}_{,m,l} −
g^{km}g^{il}_{,m,l} &minus g^{il}g^{km}_{,m,l} +
g^{ik}g^{lm}_{,m,l} +
4g^{lm}g^{ik}_{,m}
Γ^{ n}_{ln} − 2g^{kl}g^{im}_{,m}
Γ^{ n}_{ln} −
2g^{km}g^{il}_{,m}
Γ^{ n}_{ln} −
2g^{im}g^{kl}_{,m}
Γ^{ n}_{ln} −
2g^{il}g^{km}_{,m}
Γ^{ n}_{ln} +
4g^{ik}g^{lm}_{,m}
Γ^{ n}_{ln} +
4g^{ik}g^{lm}Γ^{ n}_{ln}
Γ^{ p}_{mp} +
2g^{ik}g^{lm}Γ^{ n}_{ln,m} −
4g^{il}g^{km}
Γ^{ n}_{ln}
Γ^{ p}_{mp} −
2g^{il}g^{km}
Γ^{ n}_{ln,m} +
g^{ik}g^{lm}Γ^{ n}_{lm,n} −
g^{ik}g^{lm}Γ^{ n}_{ln,m} +
g^{ik}g^{lm}Γ^{ n}_{lm}Γ^{ p}_{np} −
g^{ik}g^{lm}Γ^{ p}_{ln}Γ^{ n}_{mp} − 2g^{il}g^{km}Γ^{ n}_{lm,n} + 2g^{il}g^{km}Γ^{ n}_{ln,m} −
2g^{il}g^{km}Γ^{ n}_{lm}Γ^{ p}_{np} + 2g^{il}g^{km}Γ^{ p}_{ln}Γ^{ n}_{mp} ]
note 28
= (c^{4}/16πk)
[−
(−g^{ln}Γ^{
i}_{mn}
−g^{in}Γ^{ l}_{mn})
(−g^{mp}Γ^{
k}_{lp}
−g^{kp}Γ^{ m}_{lp})
−
(−g^{lp}Γ^{
i}_{lp}
−g^{ip}Γ^{ l}_{lp})
(−g^{mn}Γ^{
k}_{mn}
−g^{kn}Γ^{
m}_{mn}) +
2
(−g^{kn}Γ^{
i}_{ln}
−g^{in}Γ^{ k}_{ln})
(−g^{mp}Γ^{
l}_{mp}
−g^{lp}Γ^{
m}_{mp}) + g^{lm}g^{ik}_{,m,l} −
g^{km}g^{il}_{,m,l} &minus g^{il}g^{km}_{,m,l} +
g^{ik}g^{lm}_{,m,l} +
4g^{lm}
(−g^{kp}Γ^{
i}_{mp}
−g^{ip}Γ^{ k}_{mp})
Γ^{ n}_{ln} − 2g^{kl}
(−g^{mp}Γ^{
i}_{mp}
−g^{ip}Γ^{ m}_{mp})
Γ^{ n}_{ln} −
2g^{km}
(−g^{lp}Γ^{
i}_{mp}
−g^{ip}Γ^{ l}_{mp})
Γ^{ n}_{ln} −
2g^{im}
(−g^{lp}Γ^{
k}_{mp}
−g^{kp}Γ^{ l}_{mp})
Γ^{ n}_{ln} −
2g^{il}
(−g^{mp}Γ^{
k}_{mp}
−g^{kp}Γ^{ m}_{mp})
Γ^{ n}_{ln} +
4g^{ik}
(−g^{mp}Γ^{
l}_{mp}
−g^{lp}Γ^{ m}_{mp})
Γ^{ n}_{ln} +
4g^{ik}g^{lm}Γ^{ n}_{ln}
Γ^{ p}_{mp} +
g^{ik}g^{lm}Γ^{ n}_{ln,m} −
4g^{il}g^{km}
Γ^{ n}_{ln}
Γ^{ p}_{mp} +
g^{ik}g^{lm}Γ^{ n}_{lm,n} +
g^{ik}g^{lm}Γ^{ n}_{lm}Γ^{ p}_{np} −
g^{ik}g^{lm}Γ^{ p}_{ln}Γ^{ n}_{mp} − 2g^{il}g^{km}Γ^{ n}_{lm,n} −
2g^{il}g^{km}Γ^{ n}_{lm}Γ^{ p}_{np} + 2g^{il}g^{km}Γ^{ p}_{ln}Γ^{ n}_{mp} ]
note 29
= (c^{4}/16πk)
[−
(g^{lm}
g^{np}Γ^{
i}_{ln}Γ^{k}_{mp} +
g^{il}g^{mn}Γ^{
k}_{np}
Γ^{ p}_{lm} + g^{kl}g^{mn}Γ^{
i}_{np}
Γ^{
p}_{lm} +
g^{il}g^{km}Γ^{
n}_{lp}Γ^{ p}_{mn})
− (g^{lm}g^{np}Γ^{
i}_{ln}
Γ^{
k}_{mp}
+ g^{kl}g^{mn}Γ^{
i}_{mn}Γ^{ p}_{lp} + g^{il}g^{mn}Γ^{
k}_{mn}Γ^{
p}_{lp}
+ g^{il}g^{km}Γ^{ n}_{mn}Γ^{
p}_{lp}) +
2
(−g^{kn}Γ^{
i}_{ln}
−g^{in}Γ^{ k}_{ln})
(−g^{mp}Γ^{
l}_{mp}
−g^{lp}Γ^{
m}_{mp}) + g^{lm}g^{ik}_{,m,l} −
g^{km}g^{il}_{,m,l} &minus g^{il}g^{km}_{,m,l} +
g^{ik}g^{lm}_{,m,l} +
4g^{lm}
(−g^{kp}Γ^{
i}_{mp}
−g^{ip}Γ^{ k}_{mp})
Γ^{ n}_{ln} − 2g^{kl}
(−g^{mp}Γ^{
i}_{mp}
−g^{ip}Γ^{ m}_{mp})
Γ^{ n}_{ln} −
2g^{km}
(−g^{lp}Γ^{
i}_{mp}
−g^{ip}Γ^{ l}_{mp})
Γ^{ n}_{ln} −
2g^{im}
(−g^{lp}Γ^{
k}_{mp}
−g^{kp}Γ^{ l}_{mp})
Γ^{ n}_{ln} −
2g^{il}
(−g^{mp}Γ^{
k}_{mp}
−g^{kp}Γ^{ m}_{mp})
Γ^{ n}_{ln} +
4g^{ik}
(−g^{mp}Γ^{
l}_{mp}
−g^{lp}Γ^{ m}_{mp})
Γ^{ n}_{ln} +
4g^{ik}g^{lm}Γ^{ n}_{ln}
Γ^{ p}_{mp} +
g^{ik}g^{lm}Γ^{ n}_{ln,m} −
4g^{il}g^{km}
Γ^{ n}_{ln}
Γ^{ p}_{mp} +
g^{ik}g^{lm}Γ^{ n}_{lm,n} +
g^{ik}g^{lm}Γ^{ n}_{lm}Γ^{ p}_{np} −
g^{ik}g^{lm}Γ^{ p}_{ln}Γ^{ n}_{mp} − 2g^{il}g^{km}Γ^{ n}_{lm,n} −
2g^{il}g^{km}Γ^{ n}_{lm}Γ^{ p}_{np} + 2g^{il}g^{km}Γ^{ p}_{ln}Γ^{ n}_{mp} ]
note 30
=

**Note 22.** *t ^{ik}* is obtained by isolating it in (96,5).

**Note 23.** *η ^{ikl}_{,l}* is obtained from (96,3)

*η*and

^{ikl}_{,l}= (c^{4}/16πk) [(−g)(g^{ik}g^{lm}− g^{il}g^{km})]_{,m,l}*T*is expressed with

^{ik}*R*and

^{ik}*g*through the Einstein field equations (95,5).

^{ik}R**Note 24.** *(c ^{4}/16πk)* is pulled out in front; as a result, the

*R*and

^{ik}*g*terms are multiplied by 2. The minus sign from

^{ik}R*1/(−g)*is put under the derivative where it cancels with the minus sign of

*(−g)*.

**Note 25.** The expression in brackets is doubly differentiated by *,m* and *,l*. Terms that become identical on exchange of *l* and *m* indexes are summed (in red). Now in curved spacetime, the determinant *g* is no longer constant so that its derivatives are generally not zeros. This fact greatly contributes to the complexity of the resulting expressions. The SEM is expressed entirely through the covariant Ricci tensor *R _{lm}* according to LL(92,6)

*R*and LL(92,9)

^{ik}= g^{il}g^{km}R_{lm}*g*.

^{ik}R = g^{ik}g^{lm}R_{lm}**Note 24.** Differention of the expression in the first brackets (former *η ^{ikl}_{,l}* ) is carried out to give a mixture of first and second derivatives. Christoffel symbols are not decomposed in metric tensors here because our first aim is to obtain LL(96,8) which is the Landau-Lifshitz pseudotensor expressed in Christoffel symbols.

**Note 25.** Derivatives of metric determinant are replaced with derivatives of metric tensor with the help of the following variants of LL(86,4): *g _{,m} = gg^{np}g_{np,m} ; g_{,l} = gg^{np}g_{np,l}* and

*g*(replacements are in red).

_{,m,l}= (gg^{np}g_{np,m})_{,l}= gg^{np}g^{qr}g_{np,m}g_{qr,l}+ gg^{np}_{,l}g_{np,m}+ gg^{np}g_{np,m,l}**Note 26.** Expression is fully expanded to obtain 32 terms overall. All metric tensor determinants cancel and disappear. Some terms can be made equal by exchanging pairs of dummy indexes *l, m, n, p*. Pairs of equal terms (shown in various colors) are added together, reducing the number of terms to 27.