Long years ago, Kaluza and Klein proposed the initial unification theory of the electromagnetic force and gravitational force with their revolutionary 5-dimentional metric tensor. The first metric tensor proposed by Kaluza has a 5 × 5 structure and has a scalar component, k, in the 5^th row and 5^th column as can be seen in (1-1) [1] [2].

g μ ν = ( g 00 g 01 g 02 g 03 k A 0 g 10 g 11 g 12 g 13 k A 1 g 20 g 21 g 22 g 23 k A 2 g 30 g 31 g 32 g 33 k A 3 k A 0 k A 1 k A 2 k A 3 k ) (1-1)

But unfortunately, above tensor does not satisfy the tensor condition, so Klein supplemented and modified it as can be seen in (1-2) [2] [3]. There still exists a scalar component, k, in the position of the 5^th row and 5^th column. The scalar term k is called radion or dilaton or graviscalar by the following theorists [1] [3].

g μ ν = ( g 00 + k A 0 A 0 g 01 + k A 0 A 1 g 02 + k A 0 A 2 g 03 + k A 0 A 3 k A 0 g 01 + k A 0 A 1 g 11 + k A 1 A 1 g 12 + k A 1 A 2 g 13 + k A 1 A 3 k A 1 g 02 + k A 0 A 2 g 12 + k A 1 A 2 g 22 + k A 2 A 2 g 23 + k A 2 A 3 k A 2 g 03 + k A 0 A 3 g 13 + k A 1 A 3 g 23 + k A 2 A 3 g 33 + k A 3 A 3 k A 3 k A 0 k A 1 k A 2 k A 3 k ) (1-2)

Differently from the Kaluza-Klein theory, we have shown recently that gravity and electromagnetic force can be unified by using a gauge even with a 4 × 4 metric tensor [4] [5]. However, there is no gauge term with scalar component like “k” in the Kaluza-Klein metric tensor. Although the scalar term does not directly affect the electromagnetic field, we have wondered if a scalar gauge can exist even in a 4 × 4 metric tensor and have tried to find it.

As can be shown later in this paper, although the gauge in the linear approximated space-time has degree of freedom, there is some restriction in curved space-time. For that reason, we first develop our derivation in the linear approximated gravitational equation rather than the general gravitational equation below.

R μ ν − 1 2 g μ ν R = 8 π G c 4 T μ ν (1-3)

In case of g μ ν = η μ ν + k h μ ν and k ≪ 1 . The linear approximated equation of (1-3) is acquired like followings [6],

∂ λ ∂ λ h μ ν + ∂ μ ∂ ν h − ( ∂ λ ∂ ν h μ λ + ∂ λ ∂ μ h ν λ ) − η μ ν ∂ λ ∂ λ h + η μ ν ∂ λ ∂ σ h λ σ = 8 π G c 4 T μ ν (1-4)

We go through with a mathematical deployment to see if this equation has a scalar gauge by using Nöether’s theorem [7] [8]. Since T μ ν on the right side of the Equation (1-3) is an energy-momentum tensor, the symmetry of the 4 vectors is the symmetry of space-time, and also is the conservation of an energy-momentum tensor T μ ν .

In the next chapters, we will first find out a field equation by using Nöether’s theorem [7] [8] (gauge symmetry and conservation) and conservation law in general electromagnetic fields, and then extend to the gravitational field with a similar and consistent methodology.

The only four-vectors linear in A ν that can be constructed out of A ν and no more than two of the differential operators ∂ μ are listed in (2-1) [6].

∂ μ ∂ μ A ν , ∂ ν ∂ μ A μ , A ν (2-1)

The general equation using combination of these can be constructed as follows,

∂ μ ∂ μ A ν + a ∂ ν ∂ μ A μ + b A ν = 4 π   j ν (2-2)

where A ν is rank-1 field tensor, j ν is the 4-vector source term (such as charge-current term), and a , b are unknown coefficients.

Its gauge transformation is described as following,

A ν → A ν + ∂ ν Λ (2-3)

where Λ is the gauge term of rank-1 tensor A ν .

If Equation (2-2) is symmetric with respect to the gauge transformation Equation (2-3), the result of the gauge equation must be 0 like following,

∂ μ ∂ μ ∂ ν Λ + a ∂ ν ∂ μ ∂ μ Λ + b ∂ ν Λ = 0 (2-4)

or

( 1 + a ) ∂ ν ∂ μ ∂ μ Λ + b ∂ ν Λ = 0 (2-5)

In order for Equation (2-5) to always hold, the following conditions are necessary.

1 + a = 0 (i.e. a = − 1 ) (2-6)

b = 0 (2-7)

Then Equation (2-2) becomes

∂ μ ∂ μ A ν − ∂ ν ∂ μ A μ = 4 π   j ν (2-8)

In another method, if we apply the conservation of energy-momentum ( ∂ ν j ν = 0 ) to Equation (2-2), we have

∂ ν ∂ μ ∂ μ A ν + a ∂ ν ∂ ν ∂ μ A μ + b ∂ ν A ν = 4 π ∂ ν j ν = 0 (2-9)

or

( 1 + a ) ∂ ν ∂ ν ∂ μ A μ + b ∂ ν A ν = 0 (2-10)

By the properties of dummy index, we used followings

∂ μ ∂ μ = ∂ ν ∂ ν (2-11)

∂ μ A μ = ∂ ν A ν (2-12)

From Equation (2-10), we obtain a = − 1 and b = 0 . This is the same as the result of the gauge symmetry condition. From above results, we can conclude that the gauge symmetry and energy-momentum conservation are the necessary and sufficient conditions in the rank-1 tensor equation, which is consistent with Nöether’s theorem. This rank-1 tensor field equation turns out to be the Maxwell’s electro-magnetic equation. In next chapter, we can extend to the 2^nd rank tensor linear approximated gravitational equation.

As we have done with the electromagnetic field, we can proceed with a similar procedure in the linear approximated gravitational equation.

In (3-1), we can write down all possible symmetric tensors that can be constructed out of h μ ν and no more than two gradients ∂ σ which are linear in h μ ν [6].

∂ λ ∂ λ h μ ν , ∂ μ ∂ ν h , ∂ λ ∂ ν h μ λ + ∂ λ ∂ μ h ν λ , η μ ν ∂ λ ∂ λ h , η μ ν ∂ λ ∂ σ h λ σ , h μ ν , η μ ν h (3-1)

By combining these, we can create a generalized equation as shown below.

∂ λ ∂ λ h μ ν + a ∂ μ ∂ ν h + b ( ∂ λ ∂ ν h μ λ + ∂ λ ∂ μ h ν λ ) + c η μ ν ∂ λ ∂ λ h   + d η μ ν ∂ λ ∂ σ h λ σ + e h μ ν + f η μ ν h = 8 π G c 4 T μ ν (3-2)

where h μ ν is rank-2 field tensor, T μ ν is the rank-2 source tensor (energy-momentum tensor), and a , b , c , d , e , f are unsolved coefficients.

Its gauge transformation can be described as following [6] [9].

h μ ν → h μ ν + ∂ ν Λ μ + ∂ μ Λ ν (3-3)

where ∂ ν Λ μ + ∂ μ Λ ν is the gauge term of rank-2 tensor h μ ν . From now on, we will call this gauge a vector gauge. We name a vector gauge because the gauge is made from vector.

We have referenced the previous trial to find the unknown coefficients, and the following 3 sequential conditions (titled Old-1,2,3) were used in that reference [6].

Old-1. The term e h μ ν + f η μ ν h must vanish empirically because it reduces the potential exponentially as the distance increases, so that e = f = 0 .

Old-2. The results that a = − c , b = − 1 , d = 1 , e = 0 , and f = 0 are acquired by using the conservation condition ( ∂ μ T μ ν = 0 ).

Old-3. To find the unsolved a or c , they tried to replace h μ ν to h μ ν + C η μ ν h (where C is constant) in Equation (3-1) and then acquired a = 1 , c = − 1 even if C = 0 .

But Old-3 step is somewhat weird and need more explanation [6].

In the other way, we try to obtain exact solutions using 2 different conditions (titled New-1,2), gauge symmetry and energy-momentum conservation, which is similar to the work in chapter 2.

New-1. Gauge symmetry condition:

The linear equation of rank-2 tensor of h μ ν and T μ ν , as the energy momentum source, is given by Equation (2-1). To see the gauge symmetry, we try to apply ∂ ν Λ μ + ∂ μ Λ ν . To reduce the calculation, we only put the gauge term into Equation (3-2). Then the result becomes zero like (3-4).

∂ λ ∂ λ ( ∂ ν Λ μ + ∂ μ Λ ν ) + a ∂ μ ∂ ν ( ∂ σ Λ σ + ∂ λ Λ λ ) + b ( ∂ λ ∂ ν ( ∂ μ Λ λ + ∂ λ Λ μ )   + ∂ λ ∂ μ ( ∂ ν Λ λ + ∂ λ Λ ν ) ) + c η μ ν ∂ λ ∂ λ ( ∂ σ Λ σ + ∂ λ Λ λ ) + d η μ ν ∂ λ ∂ σ ( ∂ λ Λ σ   + ∂ σ Λ λ ) + e ( ∂ ν Λ μ + ∂ μ Λ ν ) + f η μ ν ( ∂ σ Λ σ + ∂ λ Λ λ ) = 0 (3-4)

Multiplying both sides by η μ ν

η μ ν ∂ λ ∂ λ ( ∂ ν Λ μ + ∂ μ Λ ν ) + a η μ ν ∂ μ ∂ ν ( ∂ σ Λ σ + ∂ λ Λ λ )   + b η μ ν ( ∂ λ ∂ ν ( ∂ μ Λ λ + ∂ λ Λ μ ) + ∂ λ ∂ μ ( ∂ ν Λ λ + ∂ λ Λ ν ) )   + 4 c ∂ λ ∂ λ ( ∂ σ Λ σ + ∂ λ Λ λ ) + 4 d ∂ λ ∂ σ ( ∂ λ Λ σ + ∂ σ Λ λ )   + 4 e ( ∂ μ Λ μ + ∂ ν Λ ν ) + 4 f ( ∂ σ Λ σ + ∂ λ Λ λ ) = 0 (3-5)

The last two terms on the left side are similar terms according to the following properties.

∂ μ Λ μ + ∂ ν Λ ν = ∂ σ Λ σ + ∂ λ Λ λ (3-6)

In this case, the following conditions are required for Equation (3-5) to hold always.

1 + b = 0 (3-7)

a + b = 0 (3-8)

c + d = 0 (3-9)

e + f = 0 (3-10)

Unlike electromagnetic equation, in gravity, e and f (the coefficients of non-derivative terms) are alive even if the gauge symmetry is satisfied. From Equation (3-7) to Equation (3-10), we have that b = − 1 and a = 1 , but c , d , e and f remain unknown. This means that all the coefficients are not determined by only ∂ ν Λ μ + ∂ μ Λ ν . For the gravitational equation, unlike the electromagnetic equation, the gauge symmetry does not guarantee the conservation of energy-momentum.

New-2. Energy-momentum conservation:

When we apply the condition for conservation of energy-momentum ( ∂ μ T μ ν = 0 ) to Equation (3-2), then

∂ μ ∂ λ ∂ λ h μ ν + a ∂ μ ∂ μ ∂ ν h + b ∂ μ ( ∂ λ ∂ ν h μ λ + ∂ λ ∂ μ h ν λ ) + c ∂ ν ∂ λ ∂ λ h   + d ∂ ν ∂ λ ∂ σ h λ σ + e ∂ μ h μ ν + f η μ ν ∂ μ h = 8 π G c 4 ∂ μ T μ ν = 0 (3-11)

Notice that similar terms have the same colors. In order for Equation (3-11) to always hold, the following conditions are necessary.

1 + b = 0 (3-12)

a + c = 0 (3-13)

b + d = 0 (3-14)

e ∂ μ h μ ν + f η μ ν ∂ μ h = 0 (3-15)

Therefore

b = − 1 (3-16)

a = − c (3-17)

d = 1 (3-18)

Since ∂ μ h μ ν and η μ ν ∂ μ h are different terms, Equation (3-11) can be established by setting both e = 0 and f = 0 . But in the case of Equation (3-15), we treat the 2 terms like a similar term, because when ∂ μ h μ ν and η μ ν ∂ μ h are multiplied by η μ ν , both become terms of h. For this reason there is a possibility other than e = f = 0 .

According to Equations (3-16)-(3-18), a , c , e and f are still unknown.

This looks that Nöether’s theorem saying gauge symmetry is conservation, does not hold. Because we do not know all the variables yet, we need more condition to find out the unknown variables. To obtain additional equations, the step Old-3 uses the fact that the field equation should be maintained even if the trace reversed form like (3-19) is used [6]. But Old-3 calculation looks a little weird, and we show that calculation is inconsistent. After that we suggest the way to resolve the inconsistency.

At first, let’s define h-bar like (3-19).

h ¯ μ ν = h μ ν + C η μ ν h (where C is arbitrary) (3-19)

If we substitute Equation (3-19) into Equation (3-2), then we have

∂ λ ∂ λ ( h μ ν + C η μ ν h ) + a ∂ μ ∂ ν ( h + 4 C h )   − ( ∂ λ ∂ ν ( h λ μ + C η λ μ h ) + ∂ λ ∂ μ ( h λ ν + C η λ ν h ) )   − a η μ ν ∂ λ ∂ λ ( h + 4 C h ) + η μ ν ∂ λ ∂ σ ( h λ σ + C η λ σ h ) = 8 π G c 4 T μ ν (3-20)

Separating the terms of h μ ν and C η μ ν h

∂ λ ∂ λ ( h μ ν ) + a ∂ μ ∂ ν ( h ) − ( ∂ λ ∂ ν ( h λ μ ) + ∂ λ ∂ μ ( h λ ν ) ) − a η μ ν ∂ λ ∂ λ ( h )   + η μ ν ∂ λ ∂ σ ( h λ σ ) + C ( − ∂ λ ∂ ν ( η λ μ h ) − a ∂ μ ∂ ν ( 4 h ) + ∂ λ ∂ λ ( η μ ν h )   + a η μ ν ∂ λ ∂ λ ( 4 h ) − η μ ν ∂ λ ∂ σ ( h λ σ ) ) = 8 π G c 4 T μ ν (3-21)

Rearranging

∂ λ ∂ λ ( h μ ν ) + a ∂ μ ∂ ν ( h ) − ( ∂ λ ∂ ν ( h λ μ ) + ∂ λ ∂ μ ( h λ ν ) ) − a η μ ν ∂ λ ∂ λ ( h )   + η μ ν ∂ λ ∂ σ ( h λ σ ) − ∂ λ ∂ ν ( C η λ μ h ) − a ∂ μ ∂ ν ( C 4 h ) + C ( ∂ μ ∂ ν h + ∂ μ ∂ ν h )   + a C η μ ν ∂ λ ∂ λ ( 4 h ) − C η μ ν ∂ λ ∂ σ h = 8 π G c 4 T μ ν (3-22)

The same color is a similar term. When the blue term is defined as a ′ , then

a ′ = a ( 1 − 4 C ) + 2 C (3-23)

or

C = a ′ − a 2 − 4 a (3-24)

On the other hand, the green term in Equation (3-22) should disappear so the coefficients should be 0 as shown below.

4 C a − 2 C = 0 (3-25)

or

2 C ( − 2 a + 1 ) = 0 (3-26)

The solution of Equation (3-26) is either C = 0 or a = 1 / 2 . But because a = 1 / 2 cannot maintain the field equation, C must be 0, and a ′ = a in Equation (3-24). Even though a cannot be determined in this approach, a and a ′ are set to 1 without justifiable explanation. The requirement is to find the value of a . But choosing C = 0 doesn’t satisfy the initial request of independent C. In the next chapter, we suggest other way to prevent this contradiction.

In Equation (3-26), in order for a to be an integer, the value of trace should be 2, instead of 4 (i.e. η μ ν η μ ν = 2 ). Then, Equation (3-23) becomes

a ′ = a ( 1 − 2 C ) + 2 C (4-1)

Comparing Equation (4-1) with Equation (3-23), 4 is changed to 2. a ′ must be equal to a to maintain the field equation form, thus we have

a ′ = a ( 1 − 2 C ) + 2 C = a (4-2)

Rearranging

0 = − 2 a C + 2 C = 2 C ( 1 − a ) (4-3)

For equation (4-3) to hold independent of C , a must be 1. This is what we hope for. In Equation (3-19), h plays as only a scalar, so we can generalize it and change h to a general scalar Λ . Therefore, we propose a new gauge as shown below and try to acquire all the unknown coefficients.

h μ ν → h μ ν + η ′ μ ν Λ (4-4)

where Λ is an arbitrary scalar and η ′ μ ν is a metric tensor with trace value 2. From now on, we call this gauge a scalar gauge compared with a vector gauge in chapter 3. We replace h μ ν with η ′ μ ν Λ in Equation (3-2) and this gauge equation should be 0.

∂ λ ∂ λ η ′ μ ν Λ + a ∂ μ ∂ ν ( 2 Λ ) + b ( ∂ λ ∂ ν η ′ μ λ Λ + ∂ λ ∂ μ η ′ ν λ Λ ) + c η μ ν ∂ λ ∂ λ ( 2 Λ )   + d η μ ν ∂ λ ∂ σ η ′ σ λ Λ + e η ′ μ ν Λ + f η μ ν ( 2 Λ ) = 0 (4-5)

Multiplying both sides by η ′ μ ν and using the properties for η ′ μ ν η ′ μ ν = 2 , η ′ μ ν η μ ν = 2 , we obtain

∂ λ ∂ λ 2 Λ + a η ′ μ ν ∂ μ ∂ ν ( 2 Λ ) + b ( ∂ λ ∂ ν 2 Λ + ∂ λ ∂ μ 2 Λ ) + 2 c ∂ λ ∂ λ ( 2 Λ )   + 2 d ∂ λ ∂ σ η ′ σ λ Λ + 2 e Λ + 2 f ( 2 Λ ) = 0 (4-6)

Similar terms have the same colors. Rearranging Equation (4-6),

2 ( 1 + 2 c + d ) ∂ λ ∂ λ Λ + ( 2 a + 2 b ) ∂ μ ∂ ν η ′ μ ν Λ + 2 ( e + 2 f ) Λ = 0 (4-7)

The following conditions are necessary for the above expression to hold always.

1 + 2 c + d = 0 (4-8)

a + b = 0 (4-9)

e + 2 f = 0 (4-10)

As we can see Equation (4-8)-(4-10), coefficients a , b , c , d , e and f are not obtained by η ′ μ ν Λ only. When we use two gauge conditions for the vector and the scalar gauge, we use Equation (3-7)-(3-10) and Equation (4-8)-(4-10), thus all coefficients can be found like (4-11).

a = 1 ,   b = − 1 ,   c = − 1 ,   d = 1 ,   e = f = 0 (4-11)

This is the purpose of our trial for searching an additional new gauge named a scalar gauge. Applying Equation (4-11) to Equation (3-2), we have

∂ λ ∂ λ h μ ν + ∂ μ ∂ ν h − ( ∂ λ ∂ ν h μ λ + ∂ λ ∂ μ h ν λ ) − η μ ν ∂ λ ∂ λ h + η μ ν ∂ λ ∂ σ h λ σ = 8 π G c 4 T μ ν (4-12)

We can write the full gauge transformation as follows.

h μ ν → h μ ν + ∂ ν Λ μ + ∂ μ Λ ν + η ′ μ ν Λ (4-13)

In Equation (4-13) ∂ ν Λ μ + ∂ μ Λ ν is fully 4 × 4 form, on the other hand η ′ μ ν Λ is only a 2 × 2 form included in the 4 × 4 equation. The equations of coefficients that can be obtained through gauge symmetry and conservation can be summarized in Table 1.

It was shown that all coefficients were determined using case (1) and case (2). This shows that symmetry is conservation, which is equal to case (3). The combination for case (1), (2), (3) is left for the our future study.

In this chapter, we want to see whether the scalar gauge obtained in the previous chapter is also symmetric in the curved space-time like the case of approximated gravitational equation by putting it into the Bianchi identities which represent the curved space-time equations.

First, let’s check the Bianchi identity.

R i j k l + R i k l j + R i l j k = 0 (5-1)

Here R i j k l is written as follows [10].

R i j k l = ∂ k Γ i j l − ∂ l Γ i j k + Γ j l p Γ i p k − Γ j k p Γ i p l (5-2)

Here, Λ is very small, so if we ignore the 3^rd and 4^th on the right side of Equation (5-2),

R i j k l ≈ ∂ k Γ i j l − ∂ l Γ i j k (5-3)

= ∂ k ( ∂ j g i l + ∂ l g i j − ∂ i g j l ) − ∂ l ( ∂ j g i k + ∂ k g i j − ∂ i g j k ) (5-4)

The gauge to be substituted here is as follows.

g μ ν → g μ ν + η ′ μ ν Λ (5-5)

Let’s substitute the gauge term only and check if there is symmetry.

R i j k l = ∂ k ( η ′ i l ∂ j Λ + η ′ i j ∂ l Λ − η ′ j l ∂ i Λ ) − ∂ l ( η ′ i k ∂ j Λ + η ′ i j ∂ k Λ − η ′ j k ∂ i Λ ) (5-6)

R i k l j = ∂ l ( η ′ i k ∂ j Λ + η ′ i j ∂ k Λ − η ′ j k ∂ i Λ ) − ∂ j ( η ′ i l ∂ k Λ + η ′ i k ∂ l Λ − η ′ k l ∂ i Λ ) (5-7)

R i l j k = ∂ j ( η ′ i k ∂ l Λ + η ′ i l ∂ k Λ − η ′ k l ∂ i Λ ) − ∂ k ( η ′ i l ∂ j Λ + η ′ i j ∂ l Λ − η ′ j l ∂ i Λ ) (5-8)

It can be seen that they exactly cancel each other and have symmetry.

Now let’s check the 2^nd Bianchi identity [10]

∇ l R m i j k + ∇ j R m i k l + ∇ k R m i l j = 0 (5-9)

Substituting Equation (5-5) here

∇ l R m i j k = ∇ l ∂ j ( η ′ m k ∂ i Λ + η ′ m i ∂ k Λ − η ′ k i ∂ m Λ ) − ∇ l ∂ k ( η ′ m i ∂ j Λ + η ′ m j ∂ i Λ − η ′ j i ∂ m Λ ) (5-10)

∇ j R m i k l = ∇ j ∂ k ( η ′ m l ∂ i Λ + η ′ m i ∂ l Λ − η ′ l i ∂ m Λ ) − ∇ j ∂ l ( η ′ m i ∂ k Λ + η ′ m k ∂ i Λ − η ′ k i ∂ m Λ ) (5-11)

∇ k R m i l j = ∇ k ∂ l ( η ′ m j ∂ i Λ + η ′ m i ∂ j Λ − η ′ j i ∂ m Λ ) − ∇ k ∂ j ( η ′ m i ∂ l Λ + η ′ m l ∂ i Λ − η ′ l i ∂ m Λ ) (5-12)

In general, since ∇ l ∂ j ∂ i ≠ ∇ j ∂ l ∂ k , the gauge symmetry does not hold in the above case. However, the special case where Λ = constant is possible to satisfy the gauge symmetry. In other words, only the global gauge has symmetry for the scalar gauge. Since the space is curved, it is considered that the symmetry is more limited than in the case of linear approximation. It is also unusual that there is no such distinction in the Kaluza-klein theory. We summarize all gauge terms in Table 2.

The expression that describes the gauge as the electromagnetic field A μ in Table 2, refers to previous papers [4] [5]. Even there, space-time is studied dividing by both linear approximated and curved [4] [5].

Where μ = 0 , 1 , 2 , 3 , ν = 0 , 1 , 2 , 3 .

We have revealed that there is a gauge including a scalar quantity even in the four dimensions. We have found that the previously known vector gauge alone can not satisfy the law of conservation of symmetry in the linear approximated gravitational equation, and an additional gauge is needed. So, by introducing a scalar gauge, we could solve this problem completely. In Kaluza-Klein theory, the scalar component compared with the scalar gauge proposed in this paper exists only in the fifth row and column of the five-dimentional metric tensor. But in this paper, we have shown that scalar gauges (trace value of η ′ μ ν as 2) exist at two points on the diagonal of a metric tensor. It is expected that physical interpretation with respect to the scalar gauge will be studied in the future. It is also found that, although the gauge in the linear approximated space-time has a degree of freedom, the gauge is constant in curved space-time. That is to say, it does not have a local gauge, but a global gauge in curved space-time. It looks natural that space-time is more symmetrical when it is more linear and vice versa. We may conclude that the limit of symmety depends on the limit of curvatue of the space-time and the symmtry term becomes a constant from a variable as space-time becomes curved.

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

Yun, Y.H. and Jang, K. (2021) Proposal of a New Scalar Gauge in Relativity Field Equation. Journal of High Energy Physics, Gravitation and Cosmology, 7, 1524-1534. https://doi.org/10.4236/jhepgc.2021.74093