# The Deng Algorithm in Higher Dimensions

Authors: Y. Nyonyi; S. D. Maharaj; K. S. Govinder
Journal: Advances in Mathematical Physics
 (2014)

---

## Abstract

We extend an algorithm of Deng in spherically symmetric spacetimes to higher dimensions. We show that it is possible to integrate the generalised condition of pressure isotropy and generate exact solutions to the Einstein field equations for a shear-free cosmological model with heat flow in higher dimensions. Three new metrics are identified which contain results of four dimensions as special cases. We show graphically that the matter variables are well behaved and the speed of sound is causal.


---

## Body

## 1. Introduction

Spherically symmetric gravitating models with heat flow, in the absence of shear, are important in the study of various cosmological processes and the evolution of relativistic astrophysical bodies. For a variety of applications in the presence of inhomogeneity, see Krasinski [1]. Heat flow models are also important in analysing gravitational collapse and relativistic stellar processes. Astrophysical studies in which heat flow is important include the shear-free models of Wagh et al. [2], Maharaj and Govender [3], Misthry et al. [4], and Herrera et al. [5]. By studying shear-free models, we avail ourselves with a rather simpler avenue where we only need to provide solutions to the generalised condition of pressure isotropy containing two metric functions. A complete study of shear-free heat conducting fluids with charge was completed by Nyonyi et al. [6] using Lie’s group theoretic approach applied to differential equations. Shearing models where heat flow is significant have been recently studied by Thirukkanesh et al. [7] for radiating spherically symmetric spheres. It turns out that the resulting nonlinear equations with shear are much more difficult to analyse.A generic method of obtaining new solutions to the Einstein field equations with heat flow was provided by Deng [8]. Using this general method, we can regain existing results and obtain new classes of solutions. Nyonyi et al. [6], Ivanov [9], and Msomi et al. [10] have obtained new solutions using the Lie group theoretic approach and other methods, by solving the underlying pressure isotropy condition. These investigations are applicable to four dimensions. Extensions to higher dimensions have also been considered by many authors because of physical requirements; for example, Bhui et al. [11] showed the absence of horizons in nonadiabatic gravitational collapse. Studies of this type motivated the Lie symmetry analysis of heat conducting fluids by Msomi et al. [12] in dimensions greater than four. In the present treatment, we extend the Deng [8] algorithm to higher dimensions and show that new results are possible.


## 2. The Model

We consider the line element of a shear-free, spherically symmetric(n+2)-dimensional manifold in the form
(1)ds2=-D2dt2+1V2(dr2+r2dXn2),
where n≥2. The gravitational potential components D and V are functions of r and t with
(2)Xn2=dθ12+sin2θ1dθ22+⋯+sin2θ1sin2θ2⋯sin2θn-1dθn2.
For a heat conducting fluid, the energy momentum tensor is given by
(3)Tab=(ρ+p)UaUb+pgab+qaUb+qbUa,
where ρ is the energy density, p is the kinetic pressure, q is the heat flux tensor, and U is a timelike (n+2)-velocity vector. For a comoving observer, we have Ua=(1/D,0,0,…,0) and qa=(0,q,0,…,0).Utilizing (1)–(3), we obtain the Einstein field equations(4a)ρ=n(n+1)Vt22D2V2-n(n+1)VVr22+nVVrr+n2VVrr,(4b)p=-nDrVVrD+nDrV2rD+n(n-1)Vr22-n(n-1)VVrr+nVttD2V-n(n+3)Vt22D2V2-nDtVtD3V,(4c)p=DrrV2D-(n-1)VVrr+n(n-1)Vr22+(n-1)DrV2rD-(n-1)2VVrr-(n-2)DrVVrD+nVttD2V-n(n+3)Vt22D2V2-nDtVtD3V,(4d)q=-nVVtrD+nVrVtD+nDrVVtD2,which are consistent with the derivation of Bhui et al. [11]. Equations (4b) and (4c), together with the transformation u=r2, give the pressure isotropy condition
(5)VDuu+2DuVu-(n-1)DVuu=0
which is the master equation for the gravitating fluid in (n+2)-dimensions.Deng [8] provided a general recipe for generating a series of solutions of the isotropy condition (5) for n=2. This technique may be extended to the master equation (5). Note that the isotropy condition is an ordinary differential equation in u (since no time derivatives appear) which may be reduced to a simpler differential equation in D if V is known and vice versa. In this technique, the elementary forms of either V or D are chosen that are in turn substituted into the equation to be solved so as to obtain the form of the remaining term. Below we give a brief outline of the method.(1)
Take a simple form ofV, say V=V1, and substitute it into the master equation to find the most general solution of D, say D=D1. The pair V=V1 and D=D1 provides the first class of solutions to (5).(2)
TakeD=D1 and substitute it into the master equation. This gives an equation in V with V=V1 already a solution. We are now in a position to obtain a second solution V=V2 linearly independent of V1. The linear combination V3=aV2+bV1 gives the general solution that satisfies the master equation. The pair V=V3 and D=D1 is the second class of solutions to (5).(3)
TakeV=V3 and substitute it into the master equation. We obtain an equation in D with D=D1 already a solution. We are then in a position to obtain D=D2 in the same way we obtained V2. The pair V=V3 and D=cD1+dD2 is the third class of solutions to (5).(4)
Repeat the above process to obtain an infinite sequence of solutions.It is important to note that, in principle, this is a nonterminating process for obtaining solutions, and an infinite number of solutions can be listed. The difficulty arises in obtaining subsequent solutions in the process because the integration may become more complicated. However, the algorithm proves to be a powerful mechanism for generating new solutions.


## 3. Results

We start with a simple case(6)D1=1.
Equation (5) reduces to
(7)Vuu=0.
This equation can be solved directly to obtain
(8)V1=au+b,
where a and b are arbitrary functions of t. The pair of equations (6) and (8) gives the first class of solutions
(9)ds2=-dt2+1(au+b)2(dr2+r2dXn2).
Observe that the dimension n is absent in metric (9). In this class of solutions the potentials are independent of the dimensionality; however, the spatial volume depends on n and the matter variables are affected through the field equations. When n=2, we regain the results of Bergmann [13].On substituting (8) back into (5) we obtain
(10)(au+b)Duu+2aDu=0.
The general solution to (10) is
(11)D2=cu+dau+b,
with c and d arbitrary functions of t. The pair of equations (8) and (11) gives the second class of solutions
(12)ds2=-(cu+dau+b)2dt2+1(au+b)2(dr2+r2dXn2).
Again, we observe that the dimension n does not appear explicitly in (12). This means that the potentials in (12) are independent of the dimension. When n=2, we regain the solutions obtained by Maiti [14] and later generalised by Modak [15] and Sanyal and Ray [16]. We make the general point that, for a linear form of V, the parameter n does not appear in (5). Thus all solutions with a linear form of V do not contain the dimension n, thereby leading to the metrics (9) and (12).Now, substituting (11) into (5), we obtain
(13)Vuu-2n-1((bc-ad)/(au+b)2(cu+d)/(au+b))Vuhh+2n-1(a(bc-ad)/(au+b)3(cu+d)/(au+b))V=0.
We require two independent solutions V1 and V2 of (13). Note that V1 in (8) is a solution to (13). We propose the second solution to (13) to be given by
(14)V2=α(u,t)V1,
where the function α(u,t) has to be found explicitly. On substituting (14) into (13), we obtain
(15)αuu+2[aau+b-1n-1((bc-ad)/(au+b)2(cu+d)/(au+b))]αu=0.
On integrating (15), we obtain α expressed as
(16)α=∫ue(cs+d(as+b)n)2/(n-1)ds.
Consequently, the second solution V2 will depend on the dimension n. To evaluate integral (16), we need to consider two cases: ad=bc and ad≠bc.


### 3.1. Case I

Whenad=bc we have
(17)α=e(dbn)2/(n-1)∫‍(1+abu)-2du=a2gu+(abg-kb2)a(au+b),
where g is an arbitrary function of t. Therefore V2 becomes
(18)V2=agu+(bg-kab2).
This implies that V2 is proportional to V1 and is therefore not a second linearly independent solution. The case ad=bc is degenerate.


### 3.2. Case II

Whenad≠bc, we have
(19)α=ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g,
where a, b, c, d, e, and g are arbitrary functions of t. Therefore the second solution V2 becomes
(20)V2=(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)(au+b).
And since the general solution to (5) is a linear combination of V1 and V2, we obtain
(21)V3=((ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)h(t)hh+j(t)(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g))h×(au+b),
where we have introduced for convenience h(t) and j(t). The third class of solutions is therefore given by (11) and (21) with metric
(22)ds2=-(cu+dau+b)2dt2+(+j(t)(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)][(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)h(t)+j(t)hh×(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)]×(au+b)+j(t)(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)])-2×(dr2+r2dXn2).
This is a new class of solutions and it is evident that it certainly depends on the dimension n≥2. Therefore we can conclude that the dimensionality of the problem does indeed affect the dynamics of the gravitational field with heat flow. The next class of solutions can be obtained by substituting V3 into (5) and then solving the resulting equation for D3. This may be continued to obtain further new solutions. The integration process gets more complicated for further iterations.We now consider the special case of four dimensions. Whene=1 and n=2, the line element (22) becomes
(23)ds2=-(cu+dau+b)2dt2+([h-j(13(ad-bc)(cu+dau+b)3+g)]×(au+b)[h-j(13(ad-bc)(cu+dau+b)3+g)])-2×(dr2+r2(dθ2+sin2θ2dϕ2)).
We can rewrite (23) in the equivalent form
(24)ds2=-(cu+dau+b)2dt2+(-j3a(c2a2+caau+bcu+d+(cu+dau+b)2)(h+κ)(au+b)-j3a(c2a2+caau+bcu+d+(cu+dau+b)2))-2×(dr2+r2(dθ2+sin2θ2dϕ2)),
where the function κ is given by
(25)κ=(g-c33a3(ad-bc)).
When we set κ=0 in (24), we regain the result of Deng [8]. We interpret (22) as the higher dimensional generalisation of the Deng model with heat flow.


## 3.1. Case I

Whenad=bc we have
(17)α=e(dbn)2/(n-1)∫‍(1+abu)-2du=a2gu+(abg-kb2)a(au+b),
where g is an arbitrary function of t. Therefore V2 becomes
(18)V2=agu+(bg-kab2).
This implies that V2 is proportional to V1 and is therefore not a second linearly independent solution. The case ad=bc is degenerate.


## 3.2. Case II

Whenad≠bc, we have
(19)α=ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g,
where a, b, c, d, e, and g are arbitrary functions of t. Therefore the second solution V2 becomes
(20)V2=(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)(au+b).
And since the general solution to (5) is a linear combination of V1 and V2, we obtain
(21)V3=((ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)h(t)hh+j(t)(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g))h×(au+b),
where we have introduced for convenience h(t) and j(t). The third class of solutions is therefore given by (11) and (21) with metric
(22)ds2=-(cu+dau+b)2dt2+(+j(t)(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)][(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)h(t)+j(t)hh×(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)]×(au+b)+j(t)(ead-bc(1-nn+1)(cu+dau+b)(n+1)/(n-1)+g)])-2×(dr2+r2dXn2).
This is a new class of solutions and it is evident that it certainly depends on the dimension n≥2. Therefore we can conclude that the dimensionality of the problem does indeed affect the dynamics of the gravitational field with heat flow. The next class of solutions can be obtained by substituting V3 into (5) and then solving the resulting equation for D3. This may be continued to obtain further new solutions. The integration process gets more complicated for further iterations.We now consider the special case of four dimensions. Whene=1 and n=2, the line element (22) becomes
(23)ds2=-(cu+dau+b)2dt2+([h-j(13(ad-bc)(cu+dau+b)3+g)]×(au+b)[h-j(13(ad-bc)(cu+dau+b)3+g)])-2×(dr2+r2(dθ2+sin2θ2dϕ2)).
We can rewrite (23) in the equivalent form
(24)ds2=-(cu+dau+b)2dt2+(-j3a(c2a2+caau+bcu+d+(cu+dau+b)2)(h+κ)(au+b)-j3a(c2a2+caau+bcu+d+(cu+dau+b)2))-2×(dr2+r2(dθ2+sin2θ2dϕ2)),
where the function κ is given by
(25)κ=(g-c33a3(ad-bc)).
When we set κ=0 in (24), we regain the result of Deng [8]. We interpret (22) as the higher dimensional generalisation of the Deng model with heat flow.


## 4. Example

We illustrate the validity of our solutions by considering a simple example with physically viable conditions. For the line element (22), we make the simple choice: a=d=0, b=c=1, h+jg=1, and e=j=1. This gives the simplified forms of the potentials
(26)D=r2,V=1+(1-nn+1)tr(2(n+1)/(n-1)).
Even with this simple example, a qualitative analysis of the matter variables and energy conditions for the interior matter distribution is arduous. We therefore generate graphical plots on a constant timelike hypersurface to illustrate the validity of our solutions using this example. Figures 1, 2, and 3 are the plots for the energy density ρ, the pressure p, and the heat flow q for three different dimensions: n=2 (dashed line), n=3 (solid line), and n=4 (dotted line). It is clearly evident that the matter variables are positive and they decrease with the increase in dimension. This is due to the fact that an increase in dimension translates to an increase in the number of degrees of freedom leading to a decrease in the mass per unit volume of the gravitating fluid. In Figure 4, we have plotted the speed of sound. From Figure 4, we observe that causality is not violated for the dimensions n=2, 3, and 4. In Figures 5, 6, and 7, we have plotted the quantities A=ρ-p+Δ, B=ρ-3p+Δ, and C=2p+Δ, where Δ=(p+q)2-4q2 for n=3. We observe that A, B, and C are positive; hence, the weak, dominant, and strong energy conditions are satisfied. Therefore the matter distribution for this example is physically reasonable.Figure 1
Energy densityρ.Figure 2
Pressure.Figure 3
Heat flow.Figure 4
Speed of sound.Figure 5
Weak energy condition.Figure 6
Dominant energy condition.Figure 7
Strong energy condition.


## 5. Discussion

We obtained new generalised classes of exact solutions to the Einstein field equations for a neutral relativistic fluid in the presence of heat flow in a higher dimensional manifold. We found new solutions to the coupled Einstein system by solving the higher dimensional pressure isotropy condition which is a second order nonlinear differential equation. We solved the master equation by making use of the Deng algorithm [8] and obtained three new metrics. The first metric (9) generalises the Bergmann [13] line element. The second metric (12) generalises the Maiti [14], Modak [15], and Sanyal and Ray [16] line elements. It is remarkable that the potentials in (9) and (12) are independent of the dimension. The third metric (22) depends on the dimension n and generalises the Deng [8] line element. We conclude that the dimension of the spacetime affects the dynamics of the heat conducting gravitating fluid. We briefly studied the physical features by graphically plotting the matter variables. The energy conditions are found to be positive and causality is not violated.



---

*Source: 290459-2014-01-09.xml*