New Way to Calculate Ricci Tensor and Ricci Scalar

Abed El Karim S. Abou Layla

Independent Researcher, Gaza City, Palestine.

DOI: 10.4236/jhepgc.2019.53042                

Abstract

In the theory of general relativity, the finding of the Einstein Field Equation happens in a complex mathematical operation, a process we don’t need any more. Through a new theory in vector analysis, we’ll see that we can calculate the components of the Ricci tensor, Ricci scalar, and Einstein Field Equation directly in an easy way without the need to use general relativity theory hypotheses, principles, and symbols. Formulating the general relativity theory through another theory will make it easier to understand this relativity theory and will help combining it with electromagnetic theory and quantum mechanics easily.

Keywords

General Relativity, Ricci Tensor, Ricci Scalar, Einstein Field Equation, Stress-Energy Tensor, Robertson-Walker Metric, Schwarzschild Metric

1. Introduction

All the books that explain general relativity theory present it in a conventional and complex way, an issue made it impossible to unit this theory with the electromagnetic theory and the quantum mechanics.

I resolved this problem through redefining the relativity theory in a different way and tie it with a new theory in the vector analysis; at present time I publish my research papers in this theory leading to unify the three theories in one.

To define my new theory (called the extended field theory), I published my first research papers, Old Mechanics, Gravity, Electromagnetics and Relativity in One Theory: “Part I”, in “Journal of High Energy Physics, Gravitation and Cosmology (JHEPGC)”, I implemented some principles and derived some of physics equations that can apply to both electric and gravitational fields, without distinguishing in this theory.

My second research published in the same journal titled “General relativity without Space-Time (₲ R)” in which, I calculated the stress-energy tensor of the electric and gravitational fields in General relativity theory without the need of curved Space-Time concept.

In this research, we will calculate the components of the Ricci tensor and Ricci scalar for many metrics, such as The Robertson-Walker metric and Schwarzschild metric directly, means that we don’t need the hypotheses, principles, and symbols of general relativity theory to get stress energy tensor, and this is important to annul the difference between the electromagnetic and gravitational, to better understanding of the universe, gravity, and black holes.

We still have to point that there is a general state to calculate the components of the Ricci tensor and Ricci scalar, which I will publish in next researches with all proofs and details.

2. Basic Notions

2.1. The Equation of Interval Angular Vector

Let’s assume the orthonormal vector z as the following

$z=z_0{e}_0+z_1{e}_1+z_2{e}_2+z_3{e}_3$

where $e_0,e_1,e_2,e_3$ 4-orthogonal unit vectors

assume that:

$z=\text{d}s=h_0\text{d}{x}_0{e}_0+h_1\text{d}{x}_1{e}_1+h_2\text{d}{x}_2{e}_2+h_3\text{d}{x}_3{e}_3$

And

$\frac{1}{z}\equiv \frac{1}{z_0}{e}_0+\frac{1}{z_1}{e}_1+\frac{1}{z_2}{e}_2+\frac{1}{z_3}{e}_3$

We conclude that

$z_0=h_0\text{d}{x}_0,z_1=h_1\text{d}{x}_1,z_2=h_2\text{d}{x}_2,z_3=h_3\text{d}{x}_3$

We can write the equation of angular velocity $\stackrel{\to }{\omega }$ as the following

$\stackrel{\to }{\omega }=\frac{1}{2}\left(\nabla \times V\right)$

In our paper, we will use the last equation as the form

$\stackrel{\to }{\omega }=\left(\nabla \times V\right)$

where,

$\stackrel{\to }{\omega }$ : 4-angular velocity

$V$ : the particle 4-velocity, $V=\frac{\text{d}s}{\text{d}{x}_0}$

$\nabla$ : bold gradient operator, $\nabla =\frac{1}{h_0}\frac{\partial }{\partial x_0}{e}_0+\frac{1}{h_1}\frac{\partial }{\partial x_1}{e}_1+\frac{1}{h_2}\frac{\partial }{\partial x_2}{e}_2+\frac{1}{h_3}\frac{\partial }{\partial x_3}{e}_3$

Now we can find the equation of 4-interval angular $\text{d}\stackrel{\to }{\phi }$ as the following,

$\text{d}\stackrel{\to }{\phi }=\stackrel{\to }{\omega }\text{d}{x}_0=\left(\nabla \times V\right)\text{d}{x}_0=\left(\nabla \times \frac{\text{d}s}{\text{d}{x}_0}\right)\text{d}{x}_0=\nabla \times z$

The cross product of two vectors $\nabla$ and $z$ as shown in my first paper [1] is given by

$\nabla \times z=\left|\begin{array}{cccc}{e}_0& e_1& e_2& e_3\\ \frac{1}{h_0}\frac{\partial }{\partial x_0}& \frac{1}{h_1}\frac{\partial }{\partial x_1}& \frac{1}{h_2}\frac{\partial }{\partial x_2}& \frac{1}{h_3}\frac{\partial }{\partial x_3}\\ z_0& z_1& z_2& z_3\end{array}\right|$

then we get

$\begin{array}{c}\text{d}{\phi }_{ij}={\left(\nabla \times z\right)}_{ij}=\left(\frac{1}{h_i}\frac{\partial }{\partial x_i}{z}_j-\frac{1}{h_j}\frac{\partial }{\partial x_j}{z}_i\right)e_i\times e_j\\ =\left(\frac{\text{d}{x}_j}{h_i}\frac{\partial }{\partial x_i}{h}_j-\frac{\text{d}{x}_i}{h_j}\frac{\partial }{\partial x_j}{h}_i\right)v_{ij}\end{array}$

where

$v_{ij}$ is a complex orthogonal unit (explained in my first paper)

$\text{d}{x}_{i,j}$ as a constant value

2.2. The Inverse of the Interval Angular Vector $d\stackrel{\leftarrow }{\phi }$

The equation of inverse angular $\text{d}\stackrel{\leftarrow }{\phi }$ is defined as

$\text{d}\stackrel{\leftarrow }{\phi }=\frac{1}{z}\times \nabla$

Or as the following

$\text{d}{\stackrel{\leftarrow }{\phi }}_{ij}=\left(\frac{1}{z_j{h}_i}\frac{\partial }{\partial x_i}-\frac{1}{z_i{h}_j}\frac{\partial }{\partial x_j}\right)v_{ij}$

3. The Details of Abou Layla’s Methods

In our new way, there are two kinds of field equations K and P, defined by

$\left(k\right)=D_{\sigma \gamma }{F}_{\sigma \gamma }$ (3.1)

$\left(P\right)={\tau }_{\sigma \gamma }{F}_{\sigma \gamma }{F}_{\sigma \gamma }$ (3.2)

where:

$D=\eta \text{d}\stackrel{\leftarrow }{\phi }$ (3.3)

$F=\eta \text{d}\stackrel{\to }{\phi }$ (3.4)

${\tau }_{\sigma \gamma }=\eta \left(\frac{1}{z_{\sigma }{z}_{\gamma }}\right)$ (3.5)

Or as the form,

${\tau }_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{1}{h_0{h}_1\text{d}{x}_0\text{d}{x}_1}& -\frac{1}{h_0{h}_2\text{d}{x}_0\text{d}{x}_2}& -\frac{1}{h_0{h}_3\text{d}{x}_0\text{d}{x}_3}\\ \frac{1}{h_0{h}_1\text{d}{x}_0\text{d}{x}_1}& 0& \frac{1}{h_1{h}_2\text{d}{x}_1\text{d}{x}_2}& \frac{1}{h_1{h}_3\text{d}{x}_1\text{d}{x}_3}\\ \frac{1}{h_0{h}_2\text{d}{x}_0\text{d}{x}_2}& \frac{1}{h_1{h}_2\text{d}{x}_1\text{d}{x}_2}& 0& \frac{1}{h_2{h}_3\text{d}{x}_2\text{d}{x}_3}\\ \frac{1}{h_0{h}_3\text{d}{x}_0\text{d}{x}_3}& \frac{1}{h_1{h}_3\text{d}{x}_1\text{d}{x}_3}& \frac{1}{h_2{h}_3\text{d}{x}_2\text{d}{x}_3}& 0\end{array}\right)$

(k) kinetic field matrix in the local frame

(P) potential field matrix in the local frame

${\eta }_{ab}=\text{diag}\left(-1,1,1,1\right)$

3.1. Define Ricci Tensor and Ricci Scalar

3.1.1. Ricci Tensor

The Ricci tensor in the local frame ${\hat{R}}_{\sigma \gamma }$ is equal the total field matrix as the following,

${\hat{R}}_{\sigma \gamma }=k+P$ (3.6)

The Ricci tensor $R_{\sigma \gamma }$ in the coordinate basis can be written as

$R_{\sigma \gamma }=g_{\sigma \gamma }{\hat{R}}_{\sigma \gamma }$ (3.7)

where

$g_{\sigma \gamma }$ : The metric of the space

3.1.2. Ricci Scalar

The Ricci scalar R is equal the sum of the elements on the main diagonal of ${\hat{R}}_{\sigma \gamma }$

$R=tr\left({\hat{R}}_{\sigma \gamma }\right)$ (3.8)

$tr\left({\hat{R}}_{\sigma \gamma }\right)$ : is the trace of a matrix ${\hat{R}}_{\sigma \gamma }$

3.2. Einstein Field Equation

In general relativity, the Einstein Field Equation is defined by [2]

$G_{\sigma \gamma }=R_{\sigma \gamma }-\frac{1}{2}R{g}_{\sigma \gamma }$

And the equation of stress-energy tensor $T_{\sigma \gamma }$ is

$G_{\sigma \gamma }=\kappa T_{\sigma \gamma }$

where:

$\kappa$ : is a constant, $\kappa =\frac{8\text{π}G}{c^4}$

İn our paper we suppose that

${\hat{G}}_{\sigma \gamma }={\hat{R}}_{\sigma \gamma }-\frac{1}{2}RI$ (3.9)

The Einstein Field Equation G in the coordinate basis can be written as

$G=g_{\sigma \gamma }\hat{G}$

3.3. Important Note

3.3.1. Options of the Potential Matrix

To find the potential matrix, there are two possibility option

Option 1

$\left(P\right)={\tau }_{\sigma \gamma }{F}_{\sigma \gamma }{F}_{\sigma \gamma }=\left({\tau }_{\sigma \gamma }{F}_{\sigma \gamma }\right)F_{\sigma \gamma }=S{F}_{\sigma \gamma }$ (3.10)

$\stackrel{\to }{S}={\tau }_{\sigma \gamma }{F}_{\sigma \gamma }$ (3.11)

Option 2

$\left(P\right)={\tau }_{\sigma \gamma }{F}_{\sigma \gamma }{F}_{\sigma \gamma }={\tau }_{\sigma \gamma }\left(F_{\sigma \gamma }{F}_{\sigma \gamma }\right)={\tau }_{\sigma \gamma }{F}_{\sigma \gamma }^2$ (3.12)

3.3.2. The Easiest Way to Find the Components of Potential Matrix

We can easily find the components of potential matrix by Assuming that,

${\delta }_{ij}=\frac{\text{d}{x}_i}{\text{d}{x}_j}$

Therefore

${\delta }_{ij}=1$ if $i=j$ , ${\delta }_{ij}=0$ if $i\ne j$

${\delta }_{ij}{\delta }_{kl}=\frac{\text{d}{x}_i}{\text{d}{x}_j}\frac{\text{d}{x}_k}{\text{d}{x}_l}=\frac{1}{\text{d}{x}_j}\frac{1}{\text{d}{x}_l}\text{d}{x}_i\text{d}{x}_k$

${\delta }_{ij}{\delta }_{ik}=\frac{\text{d}{x}_i}{\text{d}{x}_j}\frac{\text{d}{x}_i}{\text{d}{x}_k}=\frac{1}{\text{d}{x}_j}\frac{1}{\text{d}{x}_k}\text{d}{x}_i\text{d}{x}_i=0$ , it looks like $\boxed{\text{d}{x}_i\text{d}{x}_i\equiv 0}$ because k, j always different $k\ne j$

For $i\ne j\ne k$ , $\text{d}{x}_l\ne 0$ we find

${\delta }_{ij}{\delta }_{lk}=\frac{\text{d}{x}_i}{\text{d}{x}_j}\frac{\text{d}{x}_l}{\text{d}{x}_k}=\left(\frac{1}{\text{d}{x}_j\text{d}{x}_k}\text{d}{x}_i\right)\text{d}{x}_l=0$ , it looks like $\boxed{\frac{1}{\text{d}{x}_j\text{d}{x}_k}\text{d}{x}_i\equiv 0}$

4. Discussion

In this discussion, we will calculate the components of the Ricci tensor, Ricci scalar and Einstein field equation for many metrics directly by using our new way, we will see that all results will be correct and consistent with the calculations founded in General Relativity Books

4.1. Calculate the Components of the Ricci Tensor

EXAMPLE 1

Consider the metric

$\text{d}{s}^{\text{2}}=-\text{d}{t}^{\text{2}}+{\text{e}}^{\text{2}b}\left(t,r\right)\text{d}{r}^{\text{2}}+R^2\left(t,r\right)\text{d}{\varphi }^{\text{2}}$

and use Abou Layla’s structure equations to find the components of Ricci tensor in the local frame

SOLUTION

Looking at the metric, we define the orthonormal vector z as the following

$z=\text{d}tL+{\text{e}}^b\left(t,r\right)\text{d}ri+0\times h_{\theta }\text{d}\theta j+R\left(t,r\right)\text{d}\varphi k$

Then we have

$\nabla =\frac{\partial }{\partial t}L+\frac{1}{{\text{e}}^b}\frac{\partial }{\partial r}i+0\frac{\partial }{\partial \theta }j+\frac{1}{R}\frac{\partial }{\partial \varphi }k$

The equation of 4-interval angular $\text{d}\stackrel{\to }{\phi }$ is given by

$\nabla \times z=\left|\begin{array}{cccc}L& i& j& k\\ \frac{\partial }{\partial t}& \frac{1}{{\text{e}}^b}\frac{\partial }{\partial r}& 0\frac{\partial }{\partial \theta }& \frac{1}{R}\frac{\partial }{\partial \varphi }\\ 1\text{d}t& {\text{e}}^b\text{d}r& 0\text{d}\theta & R\text{d}\varphi \end{array}\right|$ , ( $\text{d}{x}_i$ as a constant value)

then we get

$\begin{array}{c}\text{d}\stackrel{\to }{\phi }=\nabla \times z\\ =\left(\frac{\partial }{\partial t}{\text{e}}^b\text{d}r-\frac{1}{{\text{e}}^b}\frac{\partial }{\partial r}1\text{d}t\right)L\times i+\left(\frac{\partial }{\partial t}0\text{d}\theta -0\frac{\partial }{\partial \theta }1\text{d}t\right)L\times j\\ +\left(\frac{\partial }{\partial t}R\text{d}\varphi -\frac{1}{R}\frac{\partial }{\partial \varphi }1\text{d}t\right)L\times k+\left(\frac{1}{{\text{e}}^b}\frac{\partial }{\partial r}0\text{d}\theta -0\frac{\partial }{\partial \theta }{\text{e}}^b\text{d}r\right)i\times j\\ +\left(\frac{1}{{\text{e}}^b}\frac{\partial }{\partial r}R\text{d}\varphi -\frac{1}{R}\frac{\partial }{\partial \varphi }{\text{e}}^b\text{d}r\right)i\times k+\left(0\frac{\partial }{\partial \theta }R\text{d}\varphi -\frac{1}{R}\frac{\partial }{\partial \varphi }0\text{d}\theta \right)j\times k\end{array}$

$\begin{array}{c}\text{d}\stackrel{\to }{\phi }=\left(\frac{\partial }{\partial t}{\text{e}}^b\text{d}r\right)v^{01}+\left(0\right)v^{02}+\left(\frac{\partial }{\partial t}R\text{d}\varphi \right)v^{03}+\left(0\right)v^{12}\\ +\left({\text{e}}^{-b}\frac{\partial }{\partial r}R\text{d}\varphi \right)v^{13}+\left(0\right)v^{23}\end{array}$

Therefore, we have

$\begin{array}{c}\text{d}{\phi }_{\sigma \gamma }=\left(\begin{array}{cccc}0& \frac{\partial }{\partial t}\left({\text{e}}^b\text{d}r\right)& 0& \frac{\partial }{\partial t}\left(R\text{d}\varphi \right)\\ -\frac{\partial }{\partial t}\left({\text{e}}^b\text{d}r\right)& 0& 0& {\text{e}}^{-b}\frac{\partial }{\partial r}\left(R\text{d}\varphi \right)\\ 0& 0& 0& 0\\ -\frac{\partial }{\partial t}\left(R\text{d}\varphi \right)& -{\text{e}}^{-b}\frac{\partial }{\partial r}\left(R\text{d}\varphi \right)& 0& 0\end{array}\right)\\ =\left(\begin{array}{cccc}0& \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& \frac{\partial R}{\partial t}\text{d}\varphi \\ -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& 0& {\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi \\ 0& 0& 0& 0\\ -\frac{\partial R}{\partial t}\text{d}\varphi & -{\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi & 0& 0\end{array}\right)\end{array}$

Thus,

$F_{\sigma \gamma }=\eta \text{d}{\phi }_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& -\frac{\partial R}{\partial t}\text{d}\varphi \\ -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& 0& {\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi \\ 0& 0& 0& 0\\ -\frac{\partial R}{\partial t}\text{d}\varphi & -{\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi & 0& 0\end{array}\right)$

We can find the components of matrix $\text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }$ directly from matrix $\text{d}{\stackrel{\to }{\phi }}_{\sigma \gamma }$ as the following

$\text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }=\left(\begin{array}{cccc}0& \frac{1}{\left({\text{e}}^b\text{d}r\right)}\frac{\partial }{\partial t}& 0& \frac{1}{\left(R\text{d}\varphi \right)}\frac{\partial }{\partial t}\\ -\frac{1}{\left({\text{e}}^b\text{d}r\right)}\frac{\partial }{\partial t}& 0& 0& \frac{{\text{e}}^{-b}}{\left(R\text{d}\varphi \right)}\frac{\partial }{\partial r}\\ 0& 0& 0& 0\\ -\frac{1}{\left(R\text{d}\varphi \right)}\frac{\partial }{\partial t}& -\frac{{\text{e}}^{-b}}{\left(R\text{d}\varphi \right)}\frac{\partial }{\partial r}& 0& 0\end{array}\right)$

Then,

$D_{\sigma \gamma }=\eta \text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{1}{{\text{e}}^b\text{d}r}\frac{\partial }{\partial t}& 0& -\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}\\ -\frac{1}{{\text{e}}^b\text{d}r}\frac{\partial }{\partial t}& 0& 0& \frac{{\text{e}}^{-b}}{R\text{d}\varphi }\frac{\partial }{\partial r}\\ 0& 0& 0& 0\\ -\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}& -\frac{e^{-b}}{R\text{d}\varphi }\frac{\partial }{\partial r}& 0& 0\end{array}\right)$

We now proceed to find the Ricci tensor ${\hat{R}}_{\sigma \gamma }$ , there are three steps:

Step1: The components of the kinetic matrix

$\begin{array}{l}\left(k\right)=D_{\sigma \gamma }{F}_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{1}{{\text{e}}^b\text{d}r}\frac{\partial }{\partial t}& 0& -\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}\\ -\frac{1}{{\text{e}}^b\text{d}r}\frac{\partial }{\partial t}& 0& 0& \frac{{\text{e}}^{-b}}{R\text{d}\varphi }\frac{\partial }{\partial r}\\ 0& 0& 0& 0\\ -\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}& -\frac{{\text{e}}^{-b}}{R\text{d}\varphi }\frac{\partial }{\partial r}& 0& 0\end{array}\right)\\ \times \left(\begin{array}{cccc}0& -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& -\frac{\partial R}{\partial t}\text{d}\varphi \\ -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& 0& {\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi \\ 0& 0& 0& 0\\ -\frac{\partial R}{\partial t}\text{d}\varphi & -{\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi & 0& 0\end{array}\right)\end{array}$

The nonzero components of kinetic matrix $k_{\sigma \gamma }$ are

$k_{tt}=\frac{1}{{\text{e}}^b\text{d}r}\frac{\partial }{\partial t}\left(\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r\right)+\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}\left(\frac{\partial R}{\partial t}\text{d}\varphi \right)$

$k_{tt}=\frac{1}{{\text{e}}^b}\frac{{\partial }^{2}b}{\partial t^2}{\text{e}}^b+\frac{1}{{\text{e}}^b}{\left(\frac{\partial b}{\partial t}\right)}^2{\text{e}}^b+\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}=\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2+\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}$

$k_{tr}=\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}\left({\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi \right)=-\frac{{\text{e}}^{-b}}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial r}+\frac{{\text{e}}^{-b}}{R}\frac{{\partial }^{2}R}{\partial t\partial r}$

$k_{rr}=\frac{1}{{\text{e}}^b\text{d}r}\frac{\partial }{\partial t}\left(\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r\right)-\frac{{\text{e}}^{-b}}{R\text{d}\varphi }\frac{\partial }{\partial r}\left({\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi \right)$

$k_{rr}=\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}$

$k_{\varphi \varphi }=\frac{1}{R\text{d}\varphi }\frac{\partial }{\partial t}\left(\frac{\partial R}{\partial t}\text{d}\varphi \right)-\frac{{\text{e}}^{-b}}{R\text{d}\varphi }\frac{\partial }{\partial r}\left({\text{e}}^{-b}\frac{\partial R}{\partial r}\text{d}\varphi \right)$

$k_{\varphi \varphi }=\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}$

The next step is to find the components of the potential filed matrix:

Using option 2 as shown in (3.12), we can calculate the square of vector $\stackrel{\to }{F}$ form as,

$F_{\sigma \gamma }^2=F_{\sigma \gamma }{F}_{\sigma \gamma }=\left(\begin{array}{cccc}0& 0& 0& -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r\frac{\partial R}{\partial r}{\text{e}}^{-b}\text{d}\varphi \\ 0& 0& 0& \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r\frac{\partial R}{\partial t}\text{d}\varphi \\ 0& 0& 0& 0\\ \frac{\partial R}{\partial r}{\text{e}}^{-b}\text{d}\varphi \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& \frac{\partial R}{\partial t}\text{d}\varphi \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& 0\end{array}\right)$

Then we obtain

$\begin{array}{c}{P}_{\sigma \gamma }={\tau }_{\sigma \gamma }{F}_{\sigma \gamma }^2\\ =\left(\begin{array}{cccc}0& -\frac{1}{{\text{e}}^b\text{d}t\text{d}r}& 0& -\frac{1}{R\text{d}t\text{d}\varphi }\\ \frac{1}{{\text{e}}^b\text{d}t\text{d}r}& 0& 0& \frac{1}{{\text{e}}^{b}R\text{d}r\text{d}\varphi }\\ 0& 0& 0& 0\\ \frac{1}{R\text{d}t\text{d}\varphi }& \frac{1}{{\text{e}}^{b}R\text{d}r\text{d}\varphi }& 0& 0\end{array}\right)\\ \times \left(\begin{array}{cccc}0& 0& 0& -\frac{\partial b}{\partial t}{\text{e}}^b\text{d}r\frac{\partial R}{\partial r}{\text{e}}^{-b}\text{d}\varphi \\ 0& 0& 0& \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r\frac{\partial R}{\partial t}\text{d}\varphi \\ 0& 0& 0& 0\\ \frac{\partial R}{\partial r}{\text{e}}^{-b}\text{d}\varphi \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& \frac{\partial R}{\partial t}\text{d}\varphi \frac{\partial b}{\partial t}{\text{e}}^b\text{d}r& 0& 0\end{array}\right)\end{array}$

And the matrix $P_{\sigma \gamma }$ is

$P_{\sigma \gamma }=\left(\begin{array}{cccc}0& 0& 0& 0\\ 0& \frac{1}{R}\frac{\partial R}{\partial t}\frac{\partial b}{\partial t}& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& \frac{1}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial t}\end{array}\right)$

Step3: Summation kinetic matrix and potential matrix

By summation (K) and (P) we get Ricci tensor ${\hat{R}}_{\sigma \gamma }$ as shown in (3.6)

Thus

${\hat{R}}_{tt}=\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2+\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}$

${\hat{R}}_{tr}=-\frac{{\text{e}}^{-b}}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial r}+\frac{{\text{e}}^{-b}}{R}\frac{{\partial }^{2}R}{\partial t\partial r}$

${\hat{R}}_{rr}=\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}+\frac{1}{R}\frac{\partial R}{\partial t}\frac{\partial b}{\partial t}$

${\hat{R}}_{\varphi \varphi }=\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}+\frac{1}{R}\frac{\partial R}{\partial t}\frac{\partial b}{\partial t}$

4.2. Calculate the Components of the Stress-Energy Tensor T

EXAMPLE 2

Using the results of Example above, find the components of the stress-energy tensor T

SOLUTION

First, we need to find the Ricci scalar R

$\begin{array}{c}R=tr\left({\hat{R}}_{\sigma \gamma }\right)=\left[\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2+\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}\right]\\ +\left[\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}+\frac{1}{R}\frac{\partial R}{\partial t}\frac{\partial b}{\partial t}\right]\\ +\left[\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}+\frac{1}{R}\frac{\partial R}{\partial t}\frac{\partial b}{\partial t}\right]\end{array}$

$R=2\left[\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2+\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}+\frac{1}{R}\frac{\partial R}{\partial t}\frac{\partial b}{\partial t}\right]$

Using (3.9), we can find the components of Einstein’s tensor as the following

${\hat{G}}_{tt}=\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}-\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}-\frac{1}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial t}$

${\hat{G}}_{tr}=\frac{{\text{e}}^{-b}}{R}\frac{{\partial }^{2}R}{\partial t\partial r}-\frac{{\text{e}}^{-b}}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial r}$

${\hat{G}}_{rr}=-\frac{1}{R}\frac{{\partial }^{2}R}{\partial t^2}$

${\hat{G}}_{rr}=-\left[\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2\right]$

Now, using transformation

$G=g_{\sigma \gamma }\hat{G}$

we obtain

$\frac{8\text{π}G}{c^4}T=\left(\begin{array}{cccc}-\frac{{\text{e}}^{-2b}}{R}\frac{{\partial }^{2}R}{\partial r^2}+\frac{{\text{e}}^{-2b}}{R}\frac{\partial b}{\partial r}\frac{\partial R}{\partial r}+\frac{1}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial t}& -\frac{1}{R}\frac{{\partial }^{2}R}{\partial t\partial r}+\frac{1}{R}\frac{\partial b}{\partial t}\frac{\partial R}{\partial r}& 0& 0\\ 0& -\frac{{\text{e}}^{2b}}{R}\frac{{\partial }^{2}R}{\partial t^2}& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& -R^2\left[\frac{{\partial }^{2}b}{\partial t^2}+{\left(\frac{\partial b}{\partial t}\right)}^2\right]\end{array}\right)$

4.3. Solution of the Robertson-Walker Metric

EXAMPLE 3

Consider the Robertson-Walker metric

$\text{d}{s}^2=-\text{d}{t}^2+\frac{a^2}{1-k{r}^2}\text{d}{r}^2+a^2\left(t\right)r^2\text{d}{\theta }^2+a^2\left(t\right)r^2{\text{sin}}^2\theta \text{d}{\varphi }^2$

and use Abou Layla’s structure equations to find the components of the Ricci tensor and the Ricci scalar

SOLUTION

For the given metric, we see that we can define the orthonormal vector z as the following

$z=-\text{d}tL+\frac{a}{\sqrt{1-k{r}^2}}\text{d}ri+a\left(t\right)r\text{d}\theta j+a\left(t\right)r\sin \theta \text{d}\varphi k$

Then we have

$\nabla =\frac{\partial }{\partial t}L+\frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}i+\frac{1}{ar}\frac{\partial }{\partial \theta }j+\frac{1}{ar\sin \theta }\frac{\partial }{\partial \varphi }k$

Now we calculate the equation of 4-interval angular $\text{d}\stackrel{\to }{\phi }$

$\nabla \times z=\left|\begin{array}{cccc}L& i& j& k\\ \frac{\partial }{\partial t}& \frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}& \frac{1}{ar}\frac{\partial }{\partial \theta }& \frac{1}{ar\sin \theta }\frac{\partial }{\partial \varphi }\\ \text{d}t& \frac{a}{\sqrt{1-k{r}^2}}\text{d}r& ar\text{d}\theta & ar\sin \theta \text{d}\varphi \end{array}\right|$ , ( $\text{d}{x}_i$ as a constant value)

then we get

$\begin{array}{c}\text{d}\stackrel{\to }{\phi }=\nabla \times z\\ =\left(\frac{\partial }{\partial t}\frac{a}{\sqrt{1-k{r}^2}}\text{d}r-\frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}\text{d}t\right)L\times i\\ +\left(\frac{\partial }{\partial t}ar\text{d}\theta -\frac{1}{ar}\frac{\partial }{\partial \theta }\text{d}t\right)L\times j\\ +\left(\frac{\partial }{\partial t}ar\sin \theta \text{d}\varphi -\frac{1}{ar\sin \theta }\frac{\partial }{\partial \varphi }\text{d}t\right)L\times k\end{array}$

$\begin{array}{l}+\left(\frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}ar\text{d}\theta -\frac{1}{ar}\frac{\partial }{\partial \theta }\frac{a}{\sqrt{1-k{r}^2}}\text{d}r\right)i\times j\\ +\left(\frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}ar\sin \theta \text{d}\varphi -\frac{1}{ar\sin \theta }\frac{\partial }{\partial \varphi }\frac{a}{\sqrt{1-k{r}^2}}\text{d}r\right)i\times k\\ +\left(\frac{1}{ar}\frac{\partial }{\partial \theta }ar\sin \theta \text{d}\varphi -\frac{1}{ar\sin \theta }\frac{\partial }{\partial \varphi }ar\text{d}\theta \right)j\times k\end{array}$

$\begin{array}{c}\text{d}\stackrel{\to }{\phi }=\left(\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r\right)v^{01}+\left(\stackrel{\dot{}}{a}r\text{d}\theta \right)v^{02}+\left(\stackrel{\dot{}}{a}r\sin \theta \text{d}\varphi \right)v^{03}\\ +\left(\sqrt{1-k{r}^2}\text{d}\theta \right)v^{12}+\left(\sqrt{1-k{r}^2}\sin \theta \text{d}\varphi \right)v^{13}+\left(\cos \theta \text{d}\varphi \right)v^{23}\end{array}$

$\text{d}{\phi }_{\sigma \gamma }=\left(\begin{array}{cccc}0& \frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r& \stackrel{\dot{}}{a}r\text{d}\theta & r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi \\ -\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r& 0& \sqrt{1-k{r}^2}\text{d}\theta & \sin \theta \sqrt{1-k{r}^2}\text{d}\varphi \\ -r\stackrel{\dot{}}{a}\text{d}\theta & -\sqrt{1-k{r}^2}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ -r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi & -\sin \theta \sqrt{1-k{r}^2}\text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)$

We can find the components of matrix $\text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }$ directly from matrix $\text{d}{\stackrel{\to }{\phi }}_{\sigma \gamma }$ as the following

$\text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }=\left(\begin{array}{cccc}0& \frac{\sqrt{1-k{r}^2}}{a\text{d}r}\frac{\partial }{\partial t}& \frac{1}{ar\text{d}\theta }\frac{\partial }{\partial t}& \frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial t}\\ -\frac{\sqrt{1-k{r}^2}}{a\text{d}r}\frac{\partial }{\partial t}& 0& \frac{\sqrt{1-k{r}^2}}{a}\frac{1}{ar\text{d}\theta }\frac{\partial }{\partial r}& \frac{\sqrt{1-k{r}^2}}{a}\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\\ -\frac{1}{ar\text{d}\theta }\frac{\partial }{\partial t}& -\frac{1}{ar\text{d}\theta }\frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}& 0& \frac{1}{ar}\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\\ -\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial t}& -\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\sqrt{1-k{r}^2}}{a}\frac{\partial }{\partial r}& -\frac{1}{ar\sin \theta }\frac{1}{ar\text{d}\varphi }\frac{\partial }{\partial \theta }& 0\end{array}\right)$

We now proceed to find the Ricci tensor ${\hat{R}}_{\sigma \gamma }$

Step1: The components of the kinetic matrix

$D_{\sigma \gamma }=\eta \text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{\sqrt{1-k{r}^2}}{a\text{d}r}\frac{\partial }{\partial t}& -\frac{1}{ar\text{d}\theta }\frac{\partial }{\partial t}& -\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial t}\\ -\frac{\sqrt{1-k{r}^2}}{a\text{d}r}\frac{\partial }{\partial t}& 0& \frac{\sqrt{1-k{r}^2}}{a^{2}r\text{d}\theta }\frac{\partial }{\partial r}& \frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\\ -\frac{1}{ar\text{d}\theta }\frac{\partial }{\partial t}& -\frac{\sqrt{1-k{r}^2}}{a^{2}r\text{d}\theta }\frac{\partial }{\partial r}& 0& \frac{1}{a^2{r}^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\\ -\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial t}& -\frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}& -\frac{1}{a^2{r}^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }& 0\end{array}\right)$ $F_{\sigma \gamma }=\eta \text{d}{\stackrel{\to }{\phi }}_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r& -\stackrel{\dot{}}{a}r\text{d}\theta & -r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi \\ -\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r& 0& \sqrt{1-k{r}^2}\text{d}\theta & \sin \theta \sqrt{1-k{r}^2}\text{d}\varphi \\ -r\stackrel{\dot{}}{a}\text{d}\theta & -\sqrt{1-k{r}^2}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ -r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi & -\sin \theta \sqrt{1-k{r}^2}\text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)$

By product of two matrix $D_{\sigma \gamma }$ and $F_{\sigma \gamma }$ we can find the components of the kinetic field matrix $k_{\sigma \gamma }$ as the following

$k_{00}=\frac{\sqrt{1-k{r}^2}}{a\text{d}r}\frac{\partial }{\partial t}\left(\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\right)+\frac{1}{ar\text{d}\theta }\frac{\partial }{\partial t}\left(r\stackrel{\dot{}}{a}\text{d}\theta \right)+\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial t}\left(r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi \right)$

$k_{00}=\frac{\ddot{a}}{a}+\frac{\ddot{a}}{a}+\frac{\ddot{a}}{a}=3\frac{\ddot{a}}{a}$

$\begin{array}{c}{k}_{11}=\frac{\sqrt{1-k{r}^2}}{adr}\frac{\partial }{\partial t}\left(\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r\right)+\frac{\sqrt{1-k{r}^2}}{a^{2}r\text{d}\theta }\frac{\partial }{\partial r}\left(-\sqrt{1-k{r}^2}\text{d}\theta \right)\\ +\frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\left(-\sin \theta \sqrt{1-k{r}^2}\text{d}\varphi \right)\end{array}$

$k_{11}=\frac{\ddot{a}}{a}+\frac{k}{a^2}+\frac{k}{a^2}=\frac{a\ddot{a}+2k}{a^2}$

$\begin{array}{c}{k}_{22}=\frac{1}{ar\text{d}\theta }\frac{\partial }{\partial t}\left(\stackrel{\dot{}}{a}r\text{d}\theta \right)-\frac{\sqrt{1-k{r}^2}}{a^{2}r\text{d}\theta }\frac{\partial }{\partial r}\left(\sqrt{1-k{r}^2}\text{d}\theta \right)\\ +\frac{1}{a^2{r}^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\left(-\cos \theta \text{d}\varphi \right)\end{array}$

$k_{22}=\frac{\ddot{a}}{a}+\frac{k}{a^2}+\frac{1}{a^2{r}^2}$

$\begin{array}{c}{k}_{33}=\frac{1}{ar\sin \theta \text{d}\varphi }\frac{\partial }{\partial t}\left(r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi \right)-\frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\left(\sin \theta \sqrt{1-k{r}^2}\text{d}\varphi \right)\\ -\frac{1}{a^2{r}^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\left(\cos \theta \text{d}\varphi \right)\end{array}$

$k_{33}=\frac{\ddot{a}}{a}+\frac{k}{a^2}+\frac{1}{a^2{r}^2}$

Step 2: The components of the potential filed matrix:

By using (3.5) we can find matrix ( $\tau$ ) as follows

${\tau }_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{\sqrt{1-k{r}^2}}{a\text{d}t\text{d}r}& -\frac{1}{ar\text{d}t\text{d}\theta }& -\frac{1}{ar\sin \theta \text{d}t\text{d}\varphi }\\ \frac{\sqrt{1-k{r}^2}}{a\text{d}r\text{d}t}& 0& \frac{\sqrt{1-k{r}^2}}{r{a}^2\text{d}r\text{d}\theta }& \frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}r\text{d}\varphi }\\ \frac{1}{ar\text{d}\theta \text{d}t}& \frac{\sqrt{1-k{r}^2}}{a^{2}r\text{d}\theta \text{d}r}& 0& \frac{1}{a^2{r}^2\sin \theta \text{d}\theta \text{d}\varphi }\\ \frac{1}{ar\sin \theta \text{d}\varphi \text{d}t}& \frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}\varphi \text{d}r}& \frac{1}{a^2{r}^2\sin \theta \text{d}\varphi \text{d}\theta }& 0\end{array}\right)$

Using option 1 as shown in (3.10) and (3.11), we can calculate the product of two matrix ${\tau }_{\sigma \gamma }$ and $\text{d}{\phi }_{\sigma \gamma }$ as follows,

$\begin{array}{l}\left(\begin{array}{cccc}0& -\frac{\sqrt{1-k{r}^2}}{a\text{d}t\text{d}r}& -\frac{1}{ar\text{d}t\text{d}\theta }& -\frac{1}{ar\sin \theta \text{d}t\text{d}\varphi }\\ \frac{\sqrt{1-k{r}^2}}{a\text{d}r\text{d}t}& 0& \frac{\sqrt{1-k{r}^2}}{r{a}^2\text{d}r\text{d}\theta }& \frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}r\text{d}\varphi }\\ \frac{1}{ar\text{d}\theta \text{d}t}& \frac{\sqrt{1-k{r}^2}}{a^{2}r\text{d}\theta \text{d}r}& 0& \frac{1}{a^2{r}^2\sin \theta \text{d}\theta \text{d}\varphi }\\ \frac{1}{ar\sin \theta \text{d}\varphi \text{d}t}& \frac{\sqrt{1-k{r}^2}}{a^{2}r\sin \theta \text{d}\varphi \text{d}r}& \frac{1}{a^2{r}^2\sin \theta \text{d}\varphi \text{d}\theta }& 0\end{array}\right)\\ \times \left(\begin{array}{cccc}0& -\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r& -\stackrel{\dot{}}{a}r\text{d}\theta & -r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi \\ -\frac{\stackrel{\dot{}}{a}}{\sqrt{1-k{r}^2}}\text{d}r& 0& \sqrt{1-k{r}^2}\text{d}\theta & \sin \theta \sqrt{1-k{r}^2}\text{d}\varphi \\ -r\stackrel{\dot{}}{a}\text{d}\theta & -\sqrt{1-k{r}^2}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ -r\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi & -\sin \theta \sqrt{1-k{r}^2}\text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)\end{array}$

Thus we have

$\left(\begin{array}{cccc}3\frac{\stackrel{\dot{}}{a}}{a\text{d}t}& 2\frac{\sqrt{1-k{r}^2}}{ar\text{d}t}& -\frac{\left(1-k{r}^2\right)}{a\text{d}t\text{d}r}\text{d}\theta +\frac{\cos \theta }{ar\sin \theta \text{d}t}& -\frac{\left(1-k{r}^2\right)\sin \theta }{a\text{d}t\text{d}r}\text{d}\varphi -\frac{\cos \theta }{ar\text{d}t\text{d}\theta }\text{d}\varphi \\ -\frac{2\stackrel{\dot{}}{a}\sqrt{1-k{r}^2}}{a^2\text{d}r}& -\frac{\stackrel{\dot{}}{a}}{a\text{d}t}-2\frac{\left(1-k{r}^2\right)}{r{a}^2\text{d}r}& -\frac{\stackrel{\dot{}}{a}r\sqrt{1-k{r}^2}}{a\text{d}r\text{d}t}\text{d}\theta -\frac{\cos \theta \sqrt{1-k{r}^2}}{r{a}^2\sin \theta \text{d}r}& -\frac{r\stackrel{\dot{}}{a}\sin \theta \sqrt{1-k{r}^2}}{a\text{d}r\text{d}t}\text{d}\varphi +\frac{\cos \theta \text{d}\varphi \sqrt{1-k{r}^2}}{a^{2}r\text{d}r\text{d}\theta }\\ -\frac{2\stackrel{\dot{}}{a}}{a^{2}r\text{d}\theta }& -\frac{\stackrel{\dot{}}{a}}{ar\sqrt{1-k{r}^2}\text{d}\theta \text{d}t}\text{d}r-\frac{\sqrt{1-k{r}^2}}{a^2{r}^2\text{d}\theta }& -\frac{\stackrel{\dot{}}{a}}{a\text{d}t}+\frac{\left(1-k{r}^2\right)}{a^{2}r\text{d}r}-\frac{\cos \theta }{a^2{r}^2\sin \theta \text{d}\theta }& -\frac{\stackrel{\dot{}}{a}\sin \theta \text{d}\varphi }{a\text{d}\theta \text{d}t}+\frac{\left(1-k{r}^2\right)\sin \theta }{a^{2}r\text{d}\theta \text{d}r}\text{d}\varphi \\ -2\frac{\stackrel{\dot{}}{a}}{a^{2}r\sin \theta \text{d}\varphi }& -\frac{\stackrel{\dot{}}{a}}{ar\sin \theta \sqrt{1-k{r}^2}\text{d}\varphi \text{d}t}\text{d}r-\frac{\sqrt{1-k{r}^2}}{a^2{r}^2\sin \theta \text{d}\varphi }& -\frac{\stackrel{\dot{}}{a}\text{d}\theta }{a\sin \theta \text{d}\varphi \text{d}t}+\frac{\left(1-k{r}^2\right)}{a^{2}r\sin \theta \text{d}\varphi \text{d}r}\text{d}\theta & -\frac{\stackrel{\dot{}}{a}}{a\text{d}t}+\frac{\left(1-k{r}^2\right)}{a^{2}r\text{d}r}+\frac{\cos \theta }{a^2{r}^2\sin \theta \text{d}\theta }\end{array}\right)$

Using the fact $\frac{1}{\text{d}{x}_j\text{d}{x}_k}\text{d}{x}_i\equiv 0$ together with (3.11), we obtain

$S_{\sigma \gamma }=\left(\begin{array}{cccc}3\frac{\stackrel{\dot{}}{a}}{a\text{d}t}& 2\frac{\sqrt{1-k{r}^2}}{ar\text{d}t}& +\frac{\cos \theta }{ar\sin \theta \text{d}t}& 0\\ -\frac{2\stackrel{\dot{}}{a}\sqrt{1-k{r}^2}}{a^2\text{d}r}& -\frac{\stackrel{\dot{}}{a}}{a\text{d}t}-2\frac{\left(1-k{r}^2\right)}{r{a}^2\text{d}r}& -\frac{\cos \theta \sqrt{1-k{r}^2}}{r{a}^2\sin \theta \text{d}r}& 0\\ -\frac{2\stackrel{\dot{}}{a}}{a^{2}r\text{d}\theta }& -\frac{\sqrt{1-k{r}^2}}{a^2{r}^2\text{d}\theta }& -\frac{\stackrel{\dot{}}{a}}{a\text{d}t}+\frac{\left(1-k{r}^2\right)}{a^{2}r\text{d}r}-\frac{\cos \theta }{a^2{r}^2\sin \theta \text{d}\theta }& 0\\ -2\frac{\stackrel{\dot{}}{a}}{a^{2}r\sin \theta \text{d}\varphi }& -\frac{\sqrt{1-k{r}^2}}{a^2{r}^2\sin \theta \text{d}\varphi }& 0& -\frac{\stackrel{\dot{}}{a}}{a\text{d}t}+\frac{\left(1-k{r}^2\right)}{a^{2}r\text{d}r}+\frac{\cos \theta }{a^2{r}^2\sin \theta \text{d}\theta }\end{array}\right)$

potential matrix is

$P_{\sigma \gamma }=S_{\sigma \gamma }\text{d}{\phi }_{\sigma \gamma }=\left(\begin{array}{cccc}0& 0& 0& 0\\ -\frac{2\stackrel{\dot{}}{a}}{r\sqrt{1-k{r}^2}}& \frac{2{\stackrel{\dot{}}{a}}^2}{a^2}& 0& 0\\ \frac{r\stackrel{\dot{}}{a}\cos \theta }{\sin \theta }& \frac{\cos \theta }{\sin \theta }\sqrt{1-k{r}^2}& 2\frac{{\stackrel{\dot{}}{a}}^2}{a^2}-\frac{\left(1-k{r}^2\right)}{a^2{r}^2}& 0\\ 0& 0& 0& 2\frac{{\stackrel{\dot{}}{a}}^2}{a^2}-\frac{\left(1-k{r}^2\right)}{a^2{r}^2}\end{array}\right)$

Step3: The components of the Ricci tensor ${\hat{R}}_{\sigma \gamma }$ in the local frame

${\hat{R}}_{00}=3\frac{\ddot{a}}{a}$

${\hat{R}}_{11}=\frac{\ddot{a}}{a}+2\frac{k}{a^2}+\frac{2{\stackrel{\dot{}}{a}}^2}{a^2}=\frac{a\ddot{a}+2k+2{\stackrel{\dot{}}{a}}^2}{a^2}$

${\hat{R}}_{22}=\frac{\ddot{a}}{a}+\frac{k}{a^2}+\frac{1}{a^2{r}^2}+2\frac{{\stackrel{\dot{}}{a}}^2}{a^2}-\frac{\left(1-k{r}^2\right)}{a^2{r}^2}=\frac{\left(a\ddot{a}+2k+2{\stackrel{\dot{}}{a}}^2\right)}{a^2}$

${\hat{R}}_{33}=\frac{\ddot{a}}{a}+\frac{k}{a^2}+\frac{1}{a^2{r}^2}+2\frac{{\stackrel{\dot{}}{a}}^2}{a^2}-\frac{\left(1-k{r}^2\right)}{a^2{r}^2}=\frac{\left(a\ddot{a}+2k+2{\stackrel{\dot{}}{a}}^2\right)}{a^2}$

Now, we can find the Ricci tensor $R_{\sigma \gamma }$ as follows:

$R_{\sigma \gamma }=g_{\sigma \gamma }{\hat{R}}_{\sigma \gamma }$

Then,

$R_{\sigma \gamma }=\left(\begin{array}{cccc}-3\frac{\ddot{a}}{a}& 0& 0& 0\\ 0& \frac{a\ddot{a}+2k+2{\stackrel{\dot{}}{a}}^2}{1-k{r}^2}& 0& 0\\ 0& 0& \left(a\ddot{a}+2k+2{\stackrel{\dot{}}{a}}^2\right)r^2& 0\\ 0& 0& 0& \left(a\ddot{a}+2k+2{\stackrel{\dot{}}{a}}^2\right)r^2{\sin }^2\theta \end{array}\right)$

we can find the Ricci Scalar R as follows:

$R=tr\left({\hat{R}}_{\sigma \gamma }\right)=6\left(\frac{a\ddot{a}+2k+{\stackrel{\dot{}}{a}}^2}{a^2}\right)$

4.4. Solution of the Schwarzschild Metric

EXAMPLE 4

For the Schwarzschild metric

$\text{d}{s}^2=-{\text{e}}^{2a}\text{d}{t}^2+{\text{e}}^{2b}\text{d}{r}^2+r^2\text{d}{\theta }^2+r^2{\text{sin}}^2\theta \text{d}{\varphi }^2$

calculate the components of the Ricci tensor and the Ricci scalar using Abou Layla’s methods

SOLUTION 3

we can define the orthonormal vector z as the following

$z={\text{e}}^a\text{d}tL+{\text{e}}^b\text{d}ri+r\text{d}\theta j+r\sin \theta \text{d}\varphi k$

Now we can calculate the equation of 4-interval angular $\text{d}\stackrel{\to }{\phi }$ by

$\nabla \times z=\left|\begin{array}{cccc}L& i& j& k\\ {\text{e}}^{-a}\frac{\partial }{\partial t}& {\text{e}}^{-b}\frac{\partial }{\partial r}& \frac{1}{r}\frac{\partial }{\partial \theta }& \frac{1}{r\sin \theta }\frac{\partial }{\partial \varphi }\\ {\text{e}}^a\text{d}t& {\text{e}}^b\text{d}r& r\text{d}\theta & r\sin \theta \text{d}\varphi \end{array}\right|$

then we get

$\begin{array}{c}\text{d}\stackrel{\to }{\phi }=\nabla \times z\\ =\left({\text{e}}^{-a}\frac{\partial }{\partial t}{\text{e}}^b\text{d}r-{\text{e}}^{-b}\frac{\partial }{\partial r}{\text{e}}^a\text{d}t\right)v^{01}\\ +\left({\text{e}}^{-a}\frac{\partial }{\partial t}r\text{d}\theta -\frac{1}{r}\frac{\partial }{\partial \theta }{\text{e}}^a\text{d}t\right)v^{02}\\ +\left({\text{e}}^{-a}\frac{\partial }{\partial t}r\sin \theta \text{d}\varphi -\frac{1}{r\sin \theta }\frac{\partial }{\partial \varphi }{\text{e}}^a\text{d}t\right)v^{03}\end{array}$

$\begin{array}{l}+\left({\text{e}}^{-b}\frac{\partial }{\partial r}r\text{d}\theta -\frac{1}{r}\frac{\partial }{\partial \theta }{\text{e}}^b\text{d}r\right)v^{12}\\ +\left({\text{e}}^{-b}\frac{\partial }{\partial r}r\sin \theta \text{d}\varphi -\frac{1}{r\sin \theta }\frac{\partial }{\partial \varphi }{\text{e}}^b\text{d}r\right)v^{13}\\ +\left(\frac{1}{r}\frac{\partial }{\partial \theta }r\sin \theta \text{d}\varphi -\frac{1}{r\sin \theta }\frac{\partial }{\partial \varphi }r\text{d}\theta \right)v^{23}\end{array}$

Therefore, we have

$\begin{array}{c}\text{d}{\phi }_{\sigma \gamma }=\left(\begin{array}{cccc}0& -{\text{e}}^{-b}\frac{\partial }{\partial r}{\text{e}}^a\text{d}t& 0& 0\\ \stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}r& 0& {\text{e}}^{-b}\frac{\partial }{\partial r}r\text{d}\theta & {\text{e}}^{-b}\frac{\partial }{\partial r}r\sin \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\text{d}\theta & 0& \frac{1}{r}\frac{\partial }{\partial \theta }r\sin \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\sin \theta \text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)\\ =\left(\begin{array}{cccc}0& -\stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& 0\\ \stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& {\text{e}}^{-b}\text{d}\theta & {\text{e}}^{-b}\sin \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\sin \theta \text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)\end{array}$

And the components of matrix $\text{d}{\stackrel{\leftarrow }{\phi }}_{\sigma \gamma }$ will be as the following

$\text{d}\stackrel{\leftarrow }{\phi }=\left(\begin{array}{cccc}0& -\frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial t}& 0& 0\\ \frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial t}& 0& \frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}& \frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\\ 0& -\frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}& 0& \frac{1}{r\left(r\sin \theta \text{d}\varphi \right)}\frac{\partial }{\partial \theta }\\ 0& -\frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}& -\frac{1}{r\left(r\sin \theta \text{d}\varphi \right)}\frac{\partial }{\partial \theta }& 0\end{array}\right)$

Step1: The components of the kinetic field matrix K

In our case, we have

$F=\text{d}\stackrel{\to }{\phi }$

$D=\text{d}\stackrel{\leftarrow }{\phi }$

${\tau }_{\sigma \gamma }=\frac{1}{z_{\sigma }{z}_{\gamma }}$

Then we get

$D_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial t}& 0& 0\\ \frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial t}& 0& \frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}& \frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\\ 0& -\frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}& 0& \frac{1}{r^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\\ 0& -\frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}& -\frac{1}{r^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }& 0\end{array}\right)$

we can find the kinetic field matrix in the coordinate basis (K) as

$\begin{array}{l}{k}_{\sigma \gamma }=\left(\begin{array}{cccc}0& -\frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial t}& 0& 0\\ \frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial t}& 0& \frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}& \frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\\ 0& -\frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}& 0& \frac{1}{r^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\\ 0& -\frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}& -\frac{1}{r^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }& 0\end{array}\right)\\ \times \left(\begin{array}{cccc}0& -\stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& 0\\ \stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& {\text{e}}^{-b}\text{d}\theta & {\text{e}}^{-b}\sin \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\sin \theta \text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)\end{array}$

The nonzero components of $k_{\sigma \gamma }$ matrix are

$k_{00}=-\frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial r}\left(\stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t\right)=-\ddot{a}{\text{e}}^{-2b}-{\stackrel{\dot{}}{a}}^2{\text{e}}^{-2b}+\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}{\text{e}}^{-2b}$

$k_{11}=\frac{{\text{e}}^{-b}}{{\text{e}}^a\text{d}t}\frac{\partial }{\partial r}\left(-\stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}r\right)+\frac{{\text{e}}^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}\left(-{\text{e}}^{-b}\text{d}\theta \right)+\frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\left(-{\text{e}}^{-b}\sin \theta \text{d}\varphi \right)$

$k_{11}=-\ddot{a}{e}^{-2b}-{\stackrel{\dot{}}{a}}^2{\text{e}}^{-2b}+\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}{\text{e}}^{-2b}+2\frac{\stackrel{\dot{}}{b}}{r}{\text{e}}^{-2b}$

$k_{22}=-\frac{e^{-b}}{r\text{d}\theta }\frac{\partial }{\partial r}\left({\text{e}}^{-b}\text{d}\theta \right)+\frac{1}{r^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\left(-\cos \theta \text{d}\varphi \right)$

$k_{22}=\frac{\stackrel{\dot{}}{b}}{r}{\text{e}}^{-2b}+\frac{1}{r^2}$

$k_{33}=-\frac{{\text{e}}^{-b}}{r\sin \theta \text{d}\varphi }\frac{\partial }{\partial r}\left({\text{e}}^{-b}\sin \theta \text{d}\varphi \right)-\frac{1}{r^2\sin \theta \text{d}\varphi }\frac{\partial }{\partial \theta }\left(\cos \theta \text{d}\varphi \right)$

$k_{33}=\frac{\stackrel{\dot{}}{b}}{r}{\text{e}}^{-2b}+\frac{1}{r^2}$

Step2: The components of the potential field matrix P

${\tau }_{\sigma \gamma }=\frac{1}{z_{\sigma }{z}_{\gamma }}=\left(\begin{array}{cccc}0& \frac{1}{{\text{e}}^a{\text{e}}^b\text{d}t\text{d}r}& \frac{1}{r{\text{e}}^a\text{d}t\text{d}\theta }& \frac{1}{r{\text{e}}^a\sin \theta \text{d}t\text{d}\varphi }\\ \frac{1}{{\text{e}}^a{\text{e}}^b\text{d}t\text{d}r}& 0& \frac{1}{r{\text{e}}^b\text{d}r\text{d}\theta }& \frac{1}{r{\text{e}}^b\sin \theta \text{d}r\text{d}\varphi }\\ \frac{1}{r{\text{e}}^a\text{d}\theta \text{d}t}& \frac{1}{r{\text{e}}^b\text{d}\theta \text{d}r}& 0& \frac{1}{r^2\sin \theta \text{d}\theta \text{d}\varphi }\\ \frac{1}{r{\text{e}}^a\sin \theta \text{d}\varphi \text{d}t}& \frac{1}{r{\text{e}}^b\sin \theta \text{d}\varphi \text{d}r}& \frac{\sin \theta }{r^2\sin \theta \text{d}\varphi \text{d}\theta }& 0\end{array}\right)$

The square of vector $\text{d}\stackrel{\to }{\phi }$ form as,

$\begin{array}{c}{\left(\text{d}\phi \right)}^2=\left(\begin{array}{cccc}0& -\stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& 0\\ \stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& {\text{e}}^{-b}\text{d}\theta & {\text{e}}^{-b}\sin \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\sin \theta \text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)\\ \times \left(\begin{array}{cccc}0& -\stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& 0\\ \stackrel{\dot{}}{a}{\text{e}}^{a-b}\text{d}t& 0& {\text{e}}^{-b}\text{d}\theta & {\text{e}}^{-b}\sin \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\text{d}\theta & 0& \cos \theta \text{d}\varphi \\ 0& -{\text{e}}^{-b}\sin \theta \text{d}\varphi & -\cos \theta \text{d}\varphi & 0\end{array}\right)\end{array}$

thus

$\begin{array}{l}{F}^2={\left(\text{d}\phi \right)}^2\\ =\left(\begin{array}{cccc}0& 0& -\stackrel{\dot{}}{a}{\text{e}}^{a-2b}\text{d}t\text{d}\theta & -\stackrel{\dot{}}{a}{\text{e}}^{a-2b}\sin \theta \text{d}t\text{d}\varphi \\ 0& 0& 0& {\text{e}}^{-b}\cos \theta \text{d}\theta \text{d}\varphi \\ -\stackrel{\dot{}}{a}{\text{e}}^{a-2b}\text{d}\theta \text{d}t& 0& 0& -{\text{e}}^{-2b}\sin \theta \text{d}\theta \text{d}\varphi \\ -\stackrel{\dot{}}{a}{\text{e}}^{a-2b}\sin \theta \text{d}\varphi \text{d}t& \cos \theta {\text{e}}^{-b}\text{d}\varphi \text{d}\theta & -{\text{e}}^{-2b}\sin \theta \text{d}\varphi \text{d}\theta & 0\end{array}\right)\end{array}$

Using (3.12) the matrix of potential field matrix $P_{\sigma \gamma }$ will be

$P_{\sigma \gamma }=\left(\begin{array}{cccc}-2\frac{\stackrel{\dot{}}{a}}{r}{\text{e}}^{-2b}& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& -\frac{\stackrel{\dot{}}{a}}{r}{\text{e}}^{-2b}-\frac{1}{r^2}{\text{e}}^{-2b}& 0\\ 0& 0& 0& -\frac{\stackrel{\dot{}}{a}}{r}{\text{e}}^{-2b}-\frac{1}{r^2}{\text{e}}^{-2b}\end{array}\right)$

Step3: The components of the Ricci tensor

${\hat{R}}_{\sigma \gamma }=k_{\sigma \gamma }+P_{\sigma \gamma }$

Then we get

${\hat{R}}_{00}=-\ddot{a}{\text{e}}^{-2b}-{\stackrel{\dot{}}{a}}^2{\text{e}}^{-2b}+\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}{\text{e}}^{-2b}-2\frac{\stackrel{\dot{}}{a}}{r}{\text{e}}^{-2b}$

${\hat{R}}_{11}=-\ddot{a}{\text{e}}^{-2b}-{\stackrel{\dot{}}{a}}^2{\text{e}}^{-2b}+\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}{\text{e}}^{-2b}+2\frac{\stackrel{\dot{}}{b}}{r}{\text{e}}^{-2b}$

${\hat{R}}_{22}=\frac{\stackrel{\dot{}}{b}}{r}{\text{e}}^{-2b}+\frac{1}{r^2}-\frac{\stackrel{\dot{}}{a}}{r}{\text{e}}^{-2b}-\frac{1}{r^2}{\text{e}}^{-2b}$

${\hat{R}}_{33}=\frac{\stackrel{\dot{}}{b}}{r}{\text{e}}^{-2b}+\frac{1}{r^2}-\frac{\stackrel{\dot{}}{a}}{r}{\text{e}}^{-2b}-\frac{1}{r^2}{\text{e}}^{-2b}$

Now, we can find the Ricci tensor $R_{\sigma \gamma }$ as follows:

$R_{\sigma \gamma }=g_{\sigma \gamma }{\hat{R}}_{\sigma \gamma }$

Then,

$R_{00}=\ddot{a}{\text{e}}^{2a-2b}+{\stackrel{\dot{}}{a}}^2{\text{e}}^{2a-2b}-\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}{\text{e}}^{2a-2b}+\frac{2}{r}\stackrel{\dot{}}{a}{\text{e}}^{2a-2b}$

$R_{11}=-\ddot{a}-{\stackrel{\dot{}}{a}}^2+\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}+\frac{2}{r}\stackrel{\dot{}}{b}$

$R_{22}=r\stackrel{\dot{}}{b}{\text{e}}^{-2b}+1-r\stackrel{\dot{}}{a}{\text{e}}^{-2b}-{\text{e}}^{-2b}$

$R_{33}=r{\sin }^2\theta \stackrel{\dot{}}{b}{\text{e}}^{-2b}+{\sin }^2\theta -r\stackrel{\dot{}}{a}{\text{e}}^{-2b}{\sin }^2\theta -{\text{e}}^{-2b}{\sin }^2\theta$

The Ricci Scalar R can be found by calculating $tr\left({\hat{R}}_{\sigma \gamma }\right)$ :

$R=tr\left({\hat{R}}_{\sigma \gamma }\right)=-2\frac{\ddot{a}+{\stackrel{\dot{}}{a}}^2-\stackrel{\dot{}}{a}\stackrel{\dot{}}{b}}{{\text{e}}^{2b}}-4\frac{\stackrel{\dot{}}{a}-\stackrel{\dot{}}{b}}{r{\text{e}}^{2b}}+2\frac{{\text{e}}^{2b}-1}{r^2{\text{e}}^{2b}}$

5. Conclusion

Dealing with general relativity theory in a matrix form will ease the way to find the Ricci tensor, Ricci scalar, and Einstein Field Equation. The most important, it will open the doors to unify this theory with electromagnetic theory in the near future.

Acknowledgements

I would like have the opportunity to obtain the Canadian citizenship in the near future due to the current conditions in my city.

I’ll be thankful to Canadian authorities if they help me in this issue, especially it will help in publishing of more of my scientific research.

Conflicts of Interest

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

References

[1]	Abou Layla, A.K. (2018) Old Mechanics, Gravity, Electromagnetics and Relativity in One Theory: Part I. Journal of High Energy Physics, Gravitation and Cosmology, 4, 529-540. https://doi.org/10.4236/jhepgc.2018.43031
[2]	Wikipedia, General Relativity. https://en.wikipedia.org/wiki/General_relativity