# Generalised spheroidal spacetimes in 5-D Einstein–Maxwell–Gauss–Bonnet gravity

###### Abstract

The field equations for static EGBM gravity are obtained and transformed to an equivalent form through a coordinate redefinition. A form for one of the metric potentials that generalises the spheroidal ansatz of Vaidya–Tikekar superdense stars and additionally prescribing the electric field intensity yields viable solutions. Some special cases of the general solution are considered and analogous classes in the Einstein framework are studied. In particular the Finch–Skea ansatz is examined in detail and found to satisfy the elementary physical requirements. These include positivity of pressure and density, the existence of a pressure free hypersurface marking the boundary, continuity with the exterior metric, a subluminal sound speed as well as the energy conditions. Moreover, the solution possesses no coordinate singularities. It is found that the impact of the Gauss–Bonnet term is to correct undesirable features in the pressure profile and sound speed index when compared to the equivalent Einstein gravity model. Furthermore graphical analyses suggest that higher densities are achievable for the same radial values when compared to the 5–dimensional Einstein case. The case of a constant gravitational potential, isothermal distribution as well as an incompressible fluid are studied. All exact solutions derived exhibit an equation of state explicitly.

###### pacs:

## I Introduction

Recently exact interior metrics were reported hm (); mh (); chil-hans () for spherically symmetric perfect fluid distributions in the Einstein–Gauss–Bonnet (EGB) gravity theory. The exterior metric was derived by Boulware and Deser boulware () for neutral spheres and by Wiltshire wilt () for the charged counterpart a few decades ago. Since then various configurations of entities were studied in EGB gravity. The collapse of inhomogeneous dust was discussed by Jhingan and Ghosh jhingan () and a distinction was drawn between shell focussing and shell crossing hyperspheres. Ghosh and Jhingan ghosh () examined the case of quasispherical gravitational collapse in 5D EGB framework and showed that naked singularities result and the cosmic censorship hypothesis is violated. Kleihaus et al blackstring () studied spinning black strings in the context of 5D EGB theory. However for some three decades compact objects were not investigated in this theory on account of the highly complicated nonlinear field equations involved. For example, the Kerr-Schild ansatz known to linearise the Einstein tensor was attempted in EGB theory anabalon () and it was found that the solution of the trace of the EGB equations did not solve all the field equations. Success in solving the EGB equations was achieved through a coordinate transformation hm (); mh (); chil-hans (). This transformation was customarily used in the standard Einstein theory to convert the equation of pressure isotropy to a linear differential equation in any of the two gravitational potentials. While the defining master equation is not rendered as a linear differential equation, it nevertheless is possible to isolate a number of exact solutions. Examination of the resulting models revealed that they satisfy physically reasonable requirements demanded of stellar models. These include satisfying the energy conditions, the existence of a pressure free hypersurface determining the boundary of the sphere as well as the condition of causality. A similar programme proved successful in constructing compact pure Lovelock stars dad-hans ().

One of the key elements in developing a model corresponding to realistic matter is the existence of an equation of state. The equation of state for a star is known to be directly linked to the mass-radius ratio for observed objects. For example, in the work of Seager et al seager () mass-radius relationships for solid exoplanets were studied in an effort to understand the constituents of such planets by analysing the possible equations of state. A key finding reported in the aforesaid work is that the mass-radius relationships for cold terrestrial mass planets did not obey a simple power law form although the equations of state for solid planets are often taken to be modelled by a polytropic equation of state. In examining data on the planet HD 149026b, the mass-radius relations suggest that almost two thirds of the mass is contained within the core. sato (). The planet GJ 436b displayed mass-radius behaviour akin to Neptune butler (); gillon (). Commencing with a prescribed equation of state imposes severe restrictions on the system of nonlinear differential equations and impedes the finding of exact solutions nils1 (); nils2 (). What has been found is that equations of state may more readily be found when other mathematical prescriptions are made to complete the model. For example see the Finch–Skea fs () model (and its charged counterpart hansraj ()) which has been shown to be compatible with astrophysical predictions in the theory of Walecka wal (). In the current problem, the field equations are more formidable yet we find that the equation of state may be determined. In fact, this strongly suggests that in constructing exact models of stars, the route of specifying a geometrical component and/ or another aspects of the physics (such as the electric field in our case), is a more productive means of determining the equation of state.

Charged fluid distributions have been traditionally studied despite the notion that celestial bodies are generally understood to be neutral. In earlier times efforts to apply solutions of the Einstein–Maxwell system to model the electron failed on account of an unrealistic mass-radius ratios. Nevertheless, it is mathematically worthwhile to examine the possibilities for charged compact objects to exist in modified theories of gravity. We attempt to construct such models in the present work. Since the problem in EMGB is similar to the Einstein theory in that there are 4 equations governing six unknowns which means that effectively there are 15 different two element sets of these six variables that may be chosen a priori to close the system. A listing of the various two element choices that have been made historically is contained in ivan (). As in the simplified Einstein case, it is a trivial exercise to find exact solutions if the two metric potentials are prescribed at the outset. However, it is by no means trivial if other two choices are made besides both metric potentials. For example, the following problems are under study. Models of charged dust by setting the pressure to zero then require one further prescription and the process is by no means trivial. Such a model is interesting as Coulombic forces oppose the pressure to prevent the collapse of the star to a point singularity. Then if a linear barotropic equation of state is imposed, one further choice needs to be made to close the system. Again, this is a nontrivial problem and will be discussed elsewhere.

The EGB action principle is a quadratic modification of the Einstein action which is linear in the Ricci tensor. Here quadratic forms of the Riemann tensor, Ricci tensor and Ricci scalar are used to construct an action principle that yields at most second order equations of motion. EGB belongs to a larger class of order polynomials discovered by Lovelock lovelock () that represent the most general Lagrangians that generate up to second order differential equations governing the evolution of stars. EGB is the case of the Lovelock polynomial. Vacuum solutions have been determined for Lovelock gravity in general whit (); whee () as well as for dimensionally continued Ban1 (); Ban2 (); Myers (); wilt (), and also for pure Lovelock black holes probes3 (); cai (). Sharma Sharm () showed that the limiting case of the polytropic fluid may be regarded as isothermal and is important in the study of clusters of stars. It is demonstrated by Dadhich and coworkers dkm () that the Schwarzschild interior solution is universal in describing a uniform density sphere in all dimensions in higher dimensional Einstein or Lovelock theory. In fact Dadhich dad () has strongly argued that the pure Lovelock theory in which only the order term (and not the sum of terms up to ) constitute the more accurate theory of gravity. Indeed the case is general relativity so the results in general relativity still hold for that special case.

Higher derivative gravity theories have been proposed in order to explain the cosmic accelerated expansion without resorting to exotic matter. Strong support for EGB theory lies in the fact that the effective action in heterotic string theory involves a Gauss-Bonnet term gross (). Of course string theory relies upon the existence of higher dimensions than 4 - something which is generally unpalatable in a theory of gravitation. However, if a grand unified theory exists then higher dimensional effects must exist in gravity theory as they do in the quantum regime - whether they are accessible is an open question. Generally, higher dimensions are explained away by claiming that they are hidden topologically. These issues were also considered in dad ().

There are other theories of the gravitational field gaining popularity.The trace-free Einstein equations (also known as unimodular gravity) have been proposed by Ellis ellis1 (); ellis2 () as the true theory of gravity. This idea, first proposed by Weinberg weinberg (), meant to address the great difference in the value of the vacuum energy density predicted by quantum theory in comparison to the actual value of the cosmological constant from astronomical observations. Starobinsky’s theory staro () is also well studied, however, unlike Lovelock theory fourth order equations of motion emerge. The theory has been shown to be conformally equivalent to scalar tensor theory plus general relativity. Of late theory harko () has been proposed where the Lagrangian is functionally dependent on the Ricci scalar as well as the trace of the energy momentum tensor . These investigations show that modified theories of gravity are indeed worthwhile areas of study and corrections to the standard theory may result. It is also important to note that studies of data collected by observations of celestial phenomena are grounded in the results of general relativity. If such computer codes were written using the modifications to the theory (such as the EGB corrections herein) the results that emerge may well be different and could convey a different portrait of the universe, stellar structures and galaxy formation and evolution. In this spirit, it is believed that the implications of extended theories are important.

In this paper we develop the governing equations for a spherically symmetric distribution of matter coupled with an electric field in the EGB framework. Exact solutions for the resulting Einstein–Maxwell-Gauss–Bonnet (EMGB) field equations are then sought. The equation of pressure isotropy is recast into an equivalent form with the help of a coordinate transformation. The analysis amounts to solving a system of four partial differential equations in six unknowns: energy density , pressure , electric field intensity , proper charge density and two metric potentials and . As is characteristic of the EGB framework the isotropy equation is second order in one variable and first order in the other. Unlike the Einstein theory, nonlinearity is inherent in both formulations of the isotropy equation.

Ordinarily an equation of state is selected to help close the system however one further assumption may be made to generate a unique model. This has generally proved a difficult route in the standard Einstein gravity. An alternative viable route is to propose functional forms for two of the six variables and endeavor to solve the system. Note that the four dynamical variables may all be expressed in terms of the metric potentials and and so specifying the geometry will trivially produce a unique model without needing to perform any integrations. This approach has been followed by krori () to generate a singularity free model but the drawback is that all control over the pressure and density is sacrificed. We elect to specify one of the metric potentials and the electric field intensity in order to solve the system of field equations. This approach has been shown to yield physically viable charged star models in Einstein theory and importantly consisting of a barotropic equation of state hansraj (). Therefore this route is being pursued in this work. The form of the potential selected corresponds to a generalised version of a spheroidal spacetime and contains a number of well studied ansatze in the Einstein theory and we investigate their behavior in this modified EMGB theory.

The paper is arranged as follows: We recall basic elements of EGB gravity and then extend the theory to incorporate the electric field by supplementing the field equations with the Maxwell’s equations. The field equations for a static spherically symmetric perfect fluid in a charged field are written and converted to an equivalent form. We then postulate a general form for the metric potential and the electric field that permits the complete integration of the field equations. Special cases of the potential are then considered and the exact model is exhibited in these cases.

## Ii Einstein–Gauss–Bonnet Gravity

The action principle in the standard Einstein theory of relativity is the Einstein–Hilbert action given by

(1) |

where is the determinant of the metric tensor , is the Ricci scalar and where is the Newton’s gravitational constant and is the speed of light in vacuum. The Einstein field equations

(2) |

consequently arise where is the energy momentum tensor for the matter configuration. If the cosmological constant is included then the Lagrangian has the form

(3) |

and the associated Einstein field equations are given by

(4) |

The Lovelock lovelock () Lagrangian is written as

(5) |

where and is the Riemann or curvature tensor. Also is the required Kronecker delta.

The Lovelock action (5) may be expanded as

from which we define the Gauss–Bonnet (GB) term as

denoted as . This term arises in the low energy effective action of heterotic string theory gross (). As a result the Einstein–Gauss–Bonnet field equations are given by

(6) |

where

The Gauss–Bonnet action is written as

(7) |

where is the GB coupling constant. The constant is linked with the string tension in string theory gross (). The remarkable feature of the GB action lies in the fact that despite the Lagrangian being quadratic in the Ricci tensor, Ricci Scalar and the Riemann tensor, the equations of motion turn out to be second order quasilinear. The GB term has no effect for but is generally non–zero for .

## Iii Field equations for charged spheres

The generic 5-D line element for static spherically symmetric spacetimes may be expressed as

(8) |

where and are the gravitational potentials. The Einstein–Maxwell-Gauss-Bonnet (EMGB) system of field equations is given by

(11) |

where is the total energy–momentum tensor, is the energy–momentum tensor for neutral matter, is the contribution of the electromagnetic field, is the electromagnetic field tensor and is the five–current density.

The electromagnetic contribution E to the total energy–momentum tensor is given by

(12) |

where is skew–symmetric. The five–current density for a non–conducting fluid can be written as

(13) |

where is the proper charge density. The electromagnetic field tensor is defined in terms of the five–potential by

(14) |

Gauge freedom allows us to choose the five–potential as

so that the effects of the magnetic field are suppressed and only the electric field contributes. Only one non-zero component of the Faraday tensor survives and is given by

(15) |

where we have utilised (14). The corresponding contravariant component has the form

where we have put

(16) |

following Herrera and Ponce de Leon herrera () for the 4 dimensional Einstein–Maxwell equations. The quantity is interpreted as the electrostatic field intensity. The components of the electromagnetic field energy tensor (12) are given by

(17) |

where we have used (15). The energy–momentum tensor for the comoving fluid velocity vector has the form

(18) |

for uncharged matter.

The field equation (11) gives the relationship

(19) |

for the value . The conservation laws generate the equation

(20) |

which can be substituted for one of the field equations.

For this line element the EMGB field equations are given by

(24) |

in the canonical spherical coordinates.

Employing the transformations , and ( a constant) the field equations (III – 24) may be rewritten as

(28) |

where (III) is the equation of pressure isotropy. Equation (III) has been arranged as a second order differential equation in . However, it should be noted that (III) may also be regarded as a first order ordinary differential equation in . Essentially the system (III) – (28) comprise an under-determined system of four coupled partial differential equations in six unknowns. To complete the system two choices for the geometric or matter variables may be made or alternatively a functional dependence of one quantity on another, such as an equation of state, may be selected and then a second choice made. Clearly all the matter variables may be written explicitly in terms of the metric functions and . This means that any metric can satisfy the EMGB field equations. This is also true for the simpler Einstein–Maxwell field equations. However, arbitrary choices of metrics will not easily yield dynamical variables that may be considered physically reasonable. It would hardly be expected that an equation of state may exist if random metrics are selected. Prescribing an equation of state then creates the problem of an intractable system - nevertheless this is a viable direction that is being pursued by the author in a different work. For the purposes of this investigation we elect to nominate an electric field intensity that behaves essentially as the reciprocal of the radius to conform with the Newtonian case. The metric ansatz chosen, namely the Vaiyda–Tikekar one, has been extensively studied and shown to generate models that do satisfy the elementary requirements for physical acceptability. Hence, this is a route we pursue in order to study the effect of the higher curvature terms on the gravitational behaviour of the charged fluid when compared to general relativity (GR).

Note that on setting in the system (III) – (28) these equations take the form

(32) |

which constitute the 5-dimensional Einstein equations for static charged perfect fluid spheres. In what follows we report new exact models for the EMGB field equations and then analyse the role of the Gauss–Bonnet coupling by setting it to zero to facilitate a comparison with the 5–D Einstein charged spacetime.

The following conditions are usually imposed on models as elementary physical requirements. It is required that the energy density () and pressure () are positive. The pressure should vanish for some radial value . The sound speed should be less than the speed of light, that is, . Across the boundary the interior metric should match with the exterior charged Boulware–Deser spacetime wilt ()

(33) |

where in the absence of the cosmological constant. Recall that is a coupling constant related to the string tension in string theory and is an arbitrary constant. Note that and represent the gravitational mass and charge of the fluid as measured by an observer at spatial infinity. It should be disclosed that these boundary conditions are extrapolated from the standard Einstein case. The junction conditions for the EGB framework was discussed by Davis davis () however these results have yet to be transformed to an explicit useable form for modelling purposes.

## Iv Generalised Vaidya–Tikekar ansatz

As the EMGB field equations constitute a system of four field equations in six unknowns, functional forms for any two of the matter or geometrical quantities may be postulated and by integration the remaining variables may be determined. The alternative approach is to specify an equation of state relating the isotropic particle pressure and the energy density, that is . This approach has has given rise to difficulties in the simpler version of the Einstein’s equations therefore is not pursued in this framework at present. For example, see the work of Nillson and Uggla nils1 (); nils2 () where the analysis of both the barotropic gamma law equation of state as well as the polytropic distribution was discussed numerically in the context of general relativity. Exact solutions were not located. A more productive approach appears to be selecting a form of one of the gravitational potentials or and also prescribing how the electric field behaves. This is the route followed here.

In our investigation we study the gravitational potential in the form

(34) |

where and are two real parameters. This may be recognised as a generalisation of the Vaidya–Tikekar ansatz vt (). Setting then the parameter is understood as the spheroidal parameter as discussed in vt ()and which has been used in the construction of models of superdense stars. Setting and regains the Finch–Skea metric fs (). The case and corresponds to the Schwarzschild interior metric in Einstein gravity and which has been shown to be a persistent solution in other theories of gravity such as Lovelock gravity lovelock (). This case deserves special consideration and will be considered in a different article. Putting reduces the metric to a Minkowski spacetime for . More general constant gravitational potentials deserve detailed treatment in their own right as they may generate isothermal fluid spheres with the Newtonian behaviour density falling off according to the inverse square law and a linear barotropic equation of state. It has been shown that isothermal fluids are universal in Lovelock gravity isothermal () and that the necessary and sufficient condition for isothermal behaviour is a constant gravitational potential. The Einstein version sas () was regained as a special case.

With the general ansatz (34), the master field equation (III) assumes the form

(35) |

where we have put . It now remains to select suitable forms for the electrostatic function to allow for the complete integration of (35). The choice

(36) |

for the electric field is motivated by the fact that the last term on the left of (35) vanishes thus effectively reducing the order of the differential equation to first order. Moreover it is a physically viable choice and has the property that the electric field disappears at the stellar centre . Also the right hand side is positive provided . The electric field has a maximum or minimum value (depending on whether or not, of when .

With this prescription, the field equation (35) is solved exactly by

(38) | |||||

where and are constants of integration. Note immediately that the choices and are excluded from the solution (38). These cases must be treated separately. The energy density, pressure and proper charge density may now be obtained with the help of (III), (III) and (28), however, given the lengthy form of the solution (38), it will be instructive to examine the physical behaviour for specific cases with a known Einstein analogue.

### iv.1 Vaidya–Tikekar case ,

The case was investigated thoroughly by Vaidya and Tikekar vt () for a four dimensional perfect fluid source without charge. The electric field in our case simplifies to

(39) |

while the gravitational potential has the form

(40) | |||||

for . The energy density is given by

(41) |

while the pressure has the form

(42) |

where we have made the redefinitions

for simplicity. The expressions for the sound speed is easily achievable but is omitted for brevity. Similarly it is straightforward to write the constants and in terms of , and using the zero boundary pressure condition as well as the continuity of across since this amounts to solving two algebraic equations linear in and however we do not display the lengthy expressions. An important observation in this model is that an equation of state exists and is easily calculated. From equation (41) may be expressed in terms of and the resulting expression for may be substituted in (42) to give the functional dependence of pressure on density as desired.

### iv.2 Finch–Skea case

The Finch–Skea fs () case is obtained for and that is . In this case the other gravitational potential is given by

(43) |

for integration constants and . This case is of particular importance as it has been demonstrated to correspond to realistic stellar distributions in Einstein gravity according to the theory of Walecka wal (). The electric field intensity is given by

(44) |

and correspondingly the proper charge density assumes the form

(45) |

The energy density and pressure have the forms

(47) | |||||

respectively. Note that solving the cubic equation (IV.2 for and plugging into (47) gives a barotropic equation of state for this model. The equation of state is explicitly given by

(48) | |||||

where we have substituted

(49) | |||||

to shorten the lengthy expression for . Although cumbersome and complicated the equation of state does indeed exist and this underscores the methodology used that of specifying a gravitational potential and the electric field intensity. If an equation of state were imposed early in the model the mathematical complexity may not have permitted a complete discovery of the model.

The sound speed may be expressed as

(50) | |||||

The mass function is computed from and is given by

(52) |

for the EGB and GR frameworks respectively. Solving the vanishing boundary condition and the continuity of metric potentials across we obtain

for the integration constants and where we have put and

. Note that a boundary certainly exists in this model as setting amounts to solving a cubic algebraic equation and the existence of at least one real-valued solution is guaranteed.
In order to examine the physical properties of our model graphically we select the following parameter values: . Plots of the dynamical quantities have been generated with the help of Mathematica 11.

The plots display some important information about the viability of our model to represent realistic distributions of charged fluids. The pressure plot Fig. 2 reveals that the pressure is zero for the radial value units (EGB) and units (GR). This suggests that spheres of larger radius are admitted through the introduction of higher derivative terms in the lagrangian density. Additionally note that the pressure decreases monotonically outwards from the centre for both EGB and GR. Within these radii, Fig. 1 shows that the energy density is always positive and also has a negative gradient. Moreover for the same radius the sphere in EGB has a greater density in general. Fig. 3 demonstrates that the sound speed remains lower than the light speed in that within the boundary for the EGB case. On the other hand, the graphs suggest that causality is violated in the radial interval for the GR framework. The energy conditions are plotted in Fig 4 and Fig 5 and we can observe that the weak, strong and dominant energy conditions are satisfied everywhere within the fluid for both gravity models: and . From Fig. 6 and Fig. 7 it can be noted that the electric field increases for a while and then slowly drops off towards the boundary and the proper charge density is always positive and a smoothly decreasing function outwardly. Since and are independent of , the plots for EGB and GR are the same. Finally Fig. 8 demonstrates a higher compactification values for the function in the EGB framework when compared to GR. Interestingly both compactification functions have the same . All of these are reasonable conditions for a perfect fluid sphere in an Einstein–Maxwell field. Clearly the higher curvature terms have a marked impact on the gravitational behaviour of the static charged fluid sphere and the indications are that undesirable features such as causality violation may be eliminated through the introduction of higher curvature terms.

### iv.3 Schwarzschild metric ansatz

Setting and generates the Schwarzschild potential . The solution (35) simplifies to

(55) |

which is also identical to the remaining Schwarzschild potential for a constant density perfect fluid in Einstein gravity. The reason for this is that for this choice of the electric field vanishes. This means that a different functional form for should be chosen to obtain a model of a charged fluid. This will be pursued in a different article.

### iv.4 Isothermal type potential

The case of a constant gravitational potential for some constant, generates the solution

(56) |

The complete solution for the energy density, pressure, charge density and speed of sound respectively is given by

(61) |

Demanding a positive energy density requires from (IV.4) while ensuring the right hand side of (IV.4) and (IV.4) remains positive restricts as . Therefore this model will only be physically reasonable if . The vanishing of the pressure is also possible although this requires solving a non-algebraic equation - however, such a radial value exists demarcating the boundary of the fluid. Observe that the energy density obeys the inverse square law which is characteristic of isothermal fluid spheres in Newtonian gravity. However, the linear barotropic equation of state () is not in effect here. An equation of state does indeed exist for this solution as may easily be found in terms of via (IV.4) and this may then be substituted in (IV.4) to give the equation of state explicitly. It has also been established by Dadhich et al isothermal () that a constant gravitational potential is a necessary and sufficient condition for isothermal behaviour for pure Lovelock gravity where the action is constructed using the order term in the Lovelock polynomial. To ensure a subluminal sound speed, it is required that (61) lie between 0 and 1. The forms of the expressions do not allow for an analytic treatment and the use of graphical representations may prove useful. The vanishing of the surface pressure at the boundary together with the matching of gravitational potentials of the charged Boulware–Deser (33) metric fix the values of the integration constants as

(63) |

in terms of the mass and charge of the fluid and where we have put

for simplicity.

## V Discussion

We have written the EMGB field equations governing the behaviour of a static ball of perfect fluid matter in the presence of an electric field. Exact solutions of the field equations were obtained by prescribing one of the metric potentials as well as by prescribing the electric field intensity to behave dimensionally as the reciprocal of the radius. The metric ansatz utilised generalised the spheroidal geometry discussed by Vaidya and Tikekar for their superdense star models. Special cases included the Vaidya–Tikekar vt () geometry, the Finch–Skea fs () metric, the Schwarzschild interior metric as well as the constant gravitational potential fluid sphere. The Schwarzschild choice coupled with the electric field prescribed failed to produce a charged distribution. A different electric field intensity must be selected to construct an EMGB model. In the remaining cases, it was possible to explicitly obtain exact solutions yielding the metric potentials, the energy density, pressure, electric field intensity and proper charge density. It was also demonstrated that in each case a hypersurface of zero pressure existed identifying the fluids boundary. Across this boundary the matching of gravitational potentials allowed to the settling of all integration constants. In the case of the Finch–Skea ansatz it was demonstrated graphically that suitable parameter values existed in order to generate a model that satisfied elementary physical requirements. In particular, it was shown that the sound speed always remained lower than the speed of light so that causality was maintained. The higher curvature terms support spheres of greater radius and density in the EGB regime when compared to the 5-dimensional Einstein case. Moreover, the GB term ensured that the model was causal while the Einstein case suffered the defect of the fluid being superluminal from the centre to a radial value well inside the boundary. This illustrates the potential of the Gauss–Bonnet higher curvature terms in correcting the physical behaviour of realistic objects which is a clear improvement on general relativity. Importantly, it has been shown in each case that a barotropic equation of state exists in the EGB framework. This investigation therefore leads us to conclude that compact star models indeed do exist in the EMGB framework and that GB higher curvature terms improve the likelihood of models conforming to realistic distributions. The important open questions under investigation include the behaviour of simple charged dust as well as the consequences of imposing a linear equation of state at the outset on the nonlinear system. In addition it will be useful to check the impact of the higher curvature terms on the properties of compact objects such as neutron stars and fluid planets. These studies are presently being undertaken.

## References

- (1) S Hansraj, B Chilambwe and S D Maharaj, Eur. Phys. J. C 27 277 (2015)
- (2) S D Maharaj, B. Chilambwe and S Hansraj, Phys. Rev. D 91, 084049 (2015)
- (3) B Chilambwe, S. Hansraj and S D Maharaj, Int. J. Mod. Phys. D 24 1550051 (2015)
- (4) D G Boulware and S Deser,Phys. Rev. Lett. 55 2656 (1985)
- (5) D L Wiltshire, Phys. Rev. D 38, 2445 (1988)
- (6) S. G. Ghosh and S. Jhingan Phys. Rev. D 82, 024017 (2010)
- (7) S. Jinghan and S. G. Ghosh Phys.Rev. D 81 024010 (2010) arXiv:1002.3245
- (8) B Kleihaus, J Kunz and E Radu Phys. Lett. B 713 110 (2012)
- (9) A. Anabalón, N. Deruelle, Y. Morisawa, J. Oliva, M Sasaki, D. Tempo and R. Troncoso Class. Quant. Grav. 26 065002 (2009)
- (10) N. Dadhich, S. Hansraj and B Chilambwe arXiv:1607.07095
- (11) S Seager, M Kuchner, C A Hier–Majumder and B Militzer The Astrophysical Journal669 1279 (2007)
- (12) B Sato, et al., ApJ, 633, 465 (2005)
- (13) R P Butler, S S Vogt, G W Marcy, D A Fischer, J T Wright, G W Henry, G Laughlin, and J J Lissauer, ApJ, 617, 580 (2004)
- (14) M Gillon, F Pont, B O Demory, F Mallmann, M Mayor, T Mazeh, D Queloz, A Shporer, S Udry and C. Vuissoz Astronomy and Astrophysics 472 L13 (2007)
- (15) B V Ivanov, Phys. Rev. D 65 104001 (2002)
- (16) D Lovelock, J. Math. Phys. 12 498 (1971).
- (17) B Whitt, Phys. Rev. D 38, 3000 (1988).
- (18) J T Wheeler, Nucl. Phys. B 268, 737 (1986); 27, 732 (1986)
- (19) M Bañados, C Teitelboim and J Zanelli, Phys. Rev. D, 49 975 (1994)
- (20) M Bañados, C Teitelboim and J Zanelli, Phys. Rev. Lett., 69 1849 (1992).
- (21) R C Myers and J Simon, Phys. Rev. D 38, 2434 (1988)
- (22) N Dadhich, K Prabhu, J Pons, Gen. Relativ. Gravit. 45 1131 (2013)
- (23) R G Cai and N Ohta, Phys. Rev. D 74 064001 (2006)
- (24) J P Sharma, Astophysics and Space Science, 163 109 (1990)
- (25) N Dadhich, A Molina, A Khugaev, Phys. Rev. D 81, 104026 (2010); arxiv:1001.3952
- (26) N. Dadhich Eur. Phys. J. C 76 104 (2016)
- (27) G F R Ellis, H van Elst, J. Murugan and J-P Uzan Class. Quantum Grav. 28 225007 (2011)
- (28) G F R Ellis, Gen. Relativ. Gravit. 46 1619 (2014)
- (29) S Weinberg Rev. Mod. Phys. 61 1 (1989)
- (30) A. A. Starobinsky, Phys. Lett. B 91 99 (1980)
- (31) Harko, T., Lobo, F.S.N., Nojiri, S. and Odintsov, S.D. Phys. Rev. D 84 024020 (2011)
- (32) K D Krori and J Barua J. Phys. A: Mathematical and General 8 508 (1975)
- (33) S Hansraj and S D Maharaj Int. J. Mod. Phys. D 15 1311 (2006)
- (34) D Gross Nucl. Phys. Proc. Suppl. 74 426 (1999)
- (35) L. Herrera and J. J. Ponce de Leon J. Math. Phys., 26 2302 (1985)
- (36) S C Davis, Phys. Rev. D 67, 024030 (2002)
- (37) Nilsson U S and Uggla C Annals of Physics 286 278 (2000)
- (38) Nilsson U S and Uggla C Annals of Physics 286 292 (2000)
- (39) PC Vaidya and R Tikekar, J. Astrophys. 3, 325 (1982)
- (40) M R Finch and J E F Skea Class. Quant. Grav. 6 467 (1989)
- (41) N. Dadhich, S Hansraj, S D Maharaj, Phys. Rev. D 93 044072 (2016)
- (42) W. C. Saslaw, S. D. Maharaj and N. Dadhich The Astrophys. J. 471 571 (1996)
- (43) J D Walecka , Phys. Lett. B 59 (1975) 109