# Higher Dimensional Charged Black Hole Solutions inf(R) Gravitational Theories

Authors: G. G. L. Nashed
Journal: Advances in High Energy Physics
 (2018)

---

## Abstract

We present, without any assumption, a class of electric and magnetic flat horizonD-dimension solutions for a specific class of f(R)=R+αR2, all of which behave asymptotically as Anti-de-Sitter spacetime. The most interesting property of these solutions is that the higher dimensions black holes, D>4, always have constant electric and magnetic charges in contrast to what is known in the literature. For D=4, we show that the magnetic field participates in the metric on equal foot as the electric field participates. Another interesting result is the fact that the Cauchy horizon is not identical with the event horizon. We use Komar formula to calculate the conserved quantities. We study the singularities and calculate the Hawking temperature and entropy and show that the first law of thermodynamics is always satisfied.


---

## Body

## 1. Introduction

The most effective gravitational theory in the last century is the theory of general relativity (GR). This theory is a fully accepted one that depicts the macroscopic geometrical properties of spacetime. Using an isotropic and homogeneous symmetry, the field equations of GR give Friedmann equations which depict the evolution of the universe with radiation and then matter dominated epochs. However, recent observations indicate that our universe goes through a phase of accelerated expansion [1–3]. This fact cannot be explained in the frame of GR using ordinary matter as a source. Another issue that GR cannot explain is the cosmological era which is known as inflation [4]. This phase of the universe is believed to have occurred before the radiation era which could relax some issues of cosmology like horizon, flatness singularities, and so on [5]. Moreover, using baryonic matter, GR is not able to discuss the observed density limited by the fitting of the standard Λ cold dark matter (ΛCDM) of the Wilkinson Microwave Anisotropy Probe results for 7 years of observations data (WMAP7) [6], the recent measurements of Baryon Acoustic Oscillations (BAO) [7], and the Hubble constant H0 [8]. Therefore, GR needs to impose an extra component known as the dark matter (DM) which constitutes about 23% of the energy content of the universe [7]. In spite of the fact that there are many possible roots of such component [9–18], DM is assumed in a form of thermal relics which naturally freeze-out with the right abundance in the extensions of the standard model of particles [19–31]. Coming experiments enable us to distinguish between large number of candidates and model by direct and indirect detection prepared for their search [32–34], or even at large hadron collider (LHC) where they could be produced [35–42].Another puzzling issue is the one of the accelerated expansion of our universe. Many explanations have been setup to demonstrate such phenomena. Among these explanations is the one which assumes the validity of GR and suggests the presence of extra fluid called dark energy (DE). The equation of state of DE takes the formp=ωρ (where p, ρ, and ω are pressure, density, and a dimensionless parameter, resp.) with ω<-1/3 to create an accelerated cosmic expansion era [43–45]. There is another model which could explain the DE which includes the cosmological constant in the field equations of GR and assumes the equation of state to have the form ω=-1. However, such model suffers from the discrepancy which comes from the fact that if we postulate the cosmological constant to represent the quantum vacuum energy then its value has higher orders of magnitude than those of observations [46]. It has been shown that, in the Palatini formalism of f(R)=R-α2/3R+R2/3β, α and β are dimension parameters, the R2 term cannot lead to an early time inflation [47].To overcome the problem of acceleration a modification of GR has been considered by modifying Einstein-Hilbert action [48–62]. Some examples of these modifications are as follows:(i)
Brans-Dicke theory whose interaction is considered by GR tensor and scalar field [63](ii)
String theory which include Gauss-Bonnet term [64](iii)
Lovelock theories which are of second derivatives at most [65].(iv)
Thef(T) gravitational theory whose field equations are of second order [66–80](v)
The one known asf(R) theories which we focused the present study on [81]It is shown thatf(R) gravitational theories are able to describe the whole story of cosmology starting from inflation to the present accelerated expansion epoch [82]. Many applications of f(R) gravitational theories have been carried out [83, 84]. Also local tests on f(R) have been achieved to constrain f(R) theories [85–88]. To study modified theories of gravity one requires to assure or reject their validity by deriving solutions that could investigate the evolution of the universe [89, 90] and the occurrence of GR astrophysical prediction. There are many BH derived in the frame of f(R) by assuming a constant scalar curvature, R=R0 [91–93]. It is the aim of the present study to abandon this condition R=R0 in a class of f(R)=R+αR2 and try to derive D-dimension solutions for a flat horizon metric spacetime.The arrangements of this study are as follows: In Section2, gradients elements of Maxwell-f(R) gravitational theory are presented. In Section 3, a metric spacetime with one unknown and two unknown functions is presented and applied to the charged field equation of f(R). Exact classes of charged black holes are derived in Section 3. In Section 4, the relevant physics of these classes is discussed by calculating the singularities. In Section 5, we calculate the conserved charges related to each class by using Komar method. In Section 6, we calculate the thermodynamical quantities like Hawking temperature and entropy, and so on. Also in this section, we have shown that the first law of thermodynamics is always satisfied for all the solutions derived in this study. The main results are discussed in the final section.


## 2. Fundamentals of Maxwell-f(R) Gravitational Theories

The Lagrangian off(R) theory has the form (1)L≔Lg+Lem, with Lg representing the gravitational Lagrangian given by (2)Lg≔12κ∫dDx-gfR-Λ, with Λ being the cosmological constant and κ being the D-dimensional gravitational constant defined as κ=2(D-3)ΩD-1GD, where GD is the Newton’s constant in D-dimensions (the units for the D-dimensional gravitational constant are GD=G4LD-4, where G4 is the gravitational constant in 4-spacetime dimensions and L is a unit of length). Here ΩD-1 refers to the volume of (D-1)-dimensional unit sphere that is defined as(3)ΩD-1=2πD-1/2ΓD-1/2,where Γ is the gamma function of the argument (for D=4, one gets 2(D-3)ΩD-1=8π),  R is the Ricci scalar of the spacetime, g is the determinant of the metric, and f(R) is the analytic function of the considered theory. In this study Lem refers to the action of Maxwell field that is defined as (4)Lem≔-12F∧⋆F,whereF=dA, with A=Aμdxμ being the 1-form electromagnetic potential [94].By carrying out variations of (2) with respect to the metric tensor gμν and the vector potential Aμ one can obtain the following field equations of f(R) gravitational theory [95, 96]:(5)RμνfR-12gμνfR-2gμνΛ+gμν□fR-∇μ∇νfR=2κTμν,(6)∂ν-gFμν=0, with Rμν being the Ricci tensor defined by(7)Rμν=Rρμρν=∂ρΓρμν-∂νΓρμρ+ΓρρβΓβμν-ΓρνβΓβμρ=2Γρμν,ρ+2ΓρβρΓβνμ,where Γρμν is the Christoffel symbols second kind and the square brackets mean(8)Aμ,ν=12Aμ,ν-Aν,μ. The D’Alembert operator □ is defined as □=∇α∇α, where ∇αVβ is the covariant derivatives of the vector Vβ and fR=df(R)/dR. In this study Tμν is defined as (9)Tμν≔gρσFνρFμσ-14δμνgλρgϵσFλϵFρσ, which is the energy momentum tensor of the electromagnetic field.The trace of (5) yields (10)RfR-2fR-8Λ+3□fR=0. Now we are going to discuss some basic property of the above f(R) theories.


### 2.1.f(R)-Gravitational Theories and Their Viable Conditions

The most important conditions and restrictions [81] that are usually put on f(R) gravitational theories to give consistency on both of gravitational and cosmological evolutions are as follows.(a) The first condition is given by(11)f′′R≥0,whenR≥f′′R, which represents the stability condition for curvature [97]. Condition (11) represents the existence of a matter dominated era in cosmological evolution. The relevant physics of (11) is that if the constant, (12)Geff=Gf′R,has a defined value, then dGeff/dR is fixed by the sign of f′′(R).(b) The conditionf′(R)>0 guarantees the positivity of the effective gravitational constant. This condition from quantum point of view avoids the graviton from becoming a ghost [98, 99].(c) The conditionf′(R)<0 ensures the recovery of GR behavior at early times.


## 2.1.f(R)-Gravitational Theories and Their Viable Conditions

The most important conditions and restrictions [81] that are usually put on f(R) gravitational theories to give consistency on both of gravitational and cosmological evolutions are as follows.(a) The first condition is given by(11)f′′R≥0,whenR≥f′′R, which represents the stability condition for curvature [97]. Condition (11) represents the existence of a matter dominated era in cosmological evolution. The relevant physics of (11) is that if the constant, (12)Geff=Gf′R,has a defined value, then dGeff/dR is fixed by the sign of f′′(R).(b) The conditionf′(R)>0 guarantees the positivity of the effective gravitational constant. This condition from quantum point of view avoids the graviton from becoming a ghost [98, 99].(c) The conditionf′(R)<0 ensures the recovery of GR behavior at early times.


## 3. Analytic Solution inf(R)=R+αR2

In this section, we are going to apply the field equations of Maxwell-f(R) with cosmological constant to two different metric spacetimes having flat horizon:


### 3.1. Flat Horizon Metric with One Variable

Let us assume the first metric spacetime possessing one unknown function has the form(13)ds2=-ardt2+1ardr2+r2∑i=1D-2dθ2i,where a(r) is an unknown function of the radial coordinate, r. Using (13) we get the Ricci scalar in the form(14)R=-r2a′′+2D-2ra′+D-2D-3ar2,where a′=da(r)/dr and a′′=d2a(r)/dr2. The nonvanishing components of the Maxwell-f(R) field equations, (5) and (6), when f(R)=R+αR2 take the form (the detailed calculations of the Ricci curvature tensor are given in Appendix B)(15)ℵtt=12r4αr4a′′2-4r4Λ+2r4q′2+2sθD-42-D-2r22rαa′a′′-ra′-aD-3-4αr2aa′′D-23D-11-4αraa′D3-10D2+24D-4-4αa24D-11-2αr3a′′′ra′+6aD-2+2raa′′′′-2αr2a′2D-2D-5=0,ℵrr=12r4αr4a′′2-4r4Λ+2r4q′2+2sθD-42-D-2r22rαa′a′′-ra′-aD-3-8αr2aa′′D-22-4αraa′D3-7D2+9D+8-4αa22D-12D-7-2αr3a′′′ra′+2aD-2-2αr2a′2D-2D-5=0,ℵθ1θ1=Sθ2θ2=⋯=SθD-2θD-2=12r4r4αa′′-1a′′+4r4Λ+2r4q′2+2sθD-42+4αr3a′a′′3D-7-2D-3r3a′-r2αa′2D-22D-9-4αr3a′′′ra′+a3D-7+raa′′′′-2αa8D2-13D+34ra′-a2D-1D-3D-6=0,where sθD-4=ds(θD-4)/dθD-4 and q′=dq(r)/dr with q(r) and s(θ) being two unknown functions related to the electric and magnetic charges of the system and defined from the gauge potential as(16)A=def.qrdt+sθD-4dθD-4.The solution of the differential equations (15) has the form(17)qr=c1+δ44c2r,sθD-4=c3+δ44c4θD-4,ar=2r2Λ3+c5r+c22+c42r21-16αΛδ44+δhhD-2rD-11±1-16DD-4αΛ/D-222DD-1D-4αrD-3+c6rD-3,where h≥5 and ci, i=1⋯6 are constants. Solution (17) is decomposed into two parts one for D=4 and the other for D≥0. The reason for this decomposition is the electromagnetic field. In the four dimensions there is a charged solution; however, for D≥5 there is no charged solution. In fact this is in conflict with the spherically symmetric case [91]. It is important to mention here that solution (17) when α=0 will reduce to the well known AdS/dS solution in the case of D=4. However, in the case of D>4 this solution does not allow the parameter α to be vanishing; therefore, it has no analogue in GR. In the noncharged case and when D=4 solution (17) coincides with what is know in GR. However, when D>4 and as long as α and Λ have no relation between them solution (17) has no analogue in GR but when Λ=1/16α it will coincide with what is known in GR; that is, it behaves as Ads/dS.The metric of solutions (17) has the form(18)ds2=-2r2Λ3+c5r+c22+c42r21-16αΛδ44+δhhD-2rD-11±1-16DD-4αΛ/D-222DD-1D-4αrD-3+c6rD-3dt2+dr22r2Λ/3+c5/r+c22+c42/r21-16αΛδ44+δhhD-2rD-11±1-16DD-4αΛ/D-22/2DD-1D-4αrD-3+c6/rD-3+r2∑i=1D-2dθ2i.


### 3.2. Metric with Two Variables

The metric of a flat horizon with two unknown functions has the form(19)ds2=-ardt2+1brdr2+r2∑i=1D-2dθ2i,where a(r) and b(r) are two unknown functions of the radial coordinate, r. Using (19) we get the Ricci scalar in the form(20)R=-a′b′2a+ba′′a-ba′22a2+D-2ba′ar+D-2b′r+D-2D-3br2. Using (19) we get the nonvanishing components of the field equations, (5) and (6), when f(R)=R+αR2 in the form (21)ℵtt=18r4a436αr4a2b2a′′2-2αr3baa′329rb′+26D-2b+8rba3a′′′′+8a3asθD-42+r4bq′2+49αr4b2a′4+αr2a′2a22bD-225D-18b+30rb′+9r2b′2-8αa2br3ab′′′2D-2a+ra′-2a′′′3rba′-a2D-2b+3rb′+4r2αaa′′rba′a27rb′+22D-2b-29rba′-a216D-2rbb′+3b′2+4D-2D-4b2-4r2a2αb′′2D-2a222D-7b+rb′+raa′6D-2b+rb′-6r2ba′2-8D-2rαa3a′rb′+2D-4brb′-b-4a3D-2rab′3D-10αrb′-r2-23D-11D-4+43D-10D-5αb+a4r4Λ-D-2D-3r2b-44D-11αb2=0,ℵrr=18r4a44αr4a2b2a′′2+2αr3baa′33rb′-4D-2b+8a3asθD-42+r4bq′2-7αr4b2a′4+αr2a′2a2r2b′2+4bD-23D-4b+rb′-8αa2b2r3a′′′2aD-2+ra′-4r2αaba′′ra′2arb′-3D-2b-3rba′+4D-2a2rb′+D-2b+ab′′2D-2a+ra′2-4D-2ra3a′23D-8αrbb′-r2b-4αb2-αr2b′2+4a3D-22αr2ab′2-a4r4Λ+D-24rbαD-2D-4-D-3r2b-44D-9αb2=0,ℵθ1θ1=ℵθ2θ2=⋯=ℵθD-2θD-2=18r4a428αr4a2b2a′′2-2αr3baa′323rb′+63D-7b+8rba3a′′′′-8a3asθD-42+r4bq′2+39αr4b2a′4+αr2a′2a22bD-22D-2D-5+2D-5D-4b+43D-7rb′+7r2b′2-2r2bα-8αa2br3ab′′′2D-2a+ra′-a′′′5rba′-a22D-5b+6rb′+4r2αaa′′rba′a28D-19b+22rb′-23rba′-a229D-22rbb′+3r2b′2-42D-3D-5-D-4b2-r2bα-8r2a2bb′′-4r2a2αb′′2D-2a24D-4b+rb′+raa′24D-9b+rb′-5r2ba′2+2ra3a′r3b′+2D-3r2b-42D-1D-5+D-43D-11αrbb′-23D-7r2αb′2+85D-19αb2-4a33D-4D-2αr2ab′2-D-3rab′r2+42D-1αb+a4r4Λ-D-4D-3r2b-2D-6D-32D-1-DD-4D-5αb2=0.The above system of differential equations (21) has the following solution:(22)1qr=c7,sθD-4=c8+δ44c9θD-4,ar=c10r2r3Λ+3c1116αΛ-1-3c122r2δ44+δhhc20D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3,br=r2r3Λ+3c1116αΛ-1-3c1223r216αΛ-1δ44+δhhD-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3,2qr=c14+δ44c15r,sθD-4=c16+δ44c17θD-4,ar=c1522Λr4+3c18+rc193r2c172+c1816αΛ-1δ44+δhhc20D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3,br=3c18+rc19+2r2Λ3r2δ44+δhhD-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3,where h≥5. We must mention here that solution (22) has the same property of solution (17); that is, when α=0 it will reduce to the well known AdS/dS solution in the case of D=4. However, in the case of D>4 this solution does not allow the parameter α to be vanishing; therefore, this solution has no analogue in GR. In the noncharged case and when D=4 and α≠0 solution (22) coincides with what is known in GR. However, when D>4 and as long as α and Λ have no relation between them solution (22) has no analogue in GR but when Λ=1/16α it will coincide with GR solution that behaves as Ads/dS. The metric spacetimes of solutions (22) have the form(23)1ds2=-c10232r4αΛ2+48c11αΛr-2r4Λ+3c122+3c11rr2δ44+δhhc20D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3dt2+32r4αΛ2+48c11αΛr-2r4Λ-3c122-3c11r3r216αΛ-1δ44+δhhD-2rD-11±1-16DD-4αΛ/D-222DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i,2ds2=-c1522Λr4+3c18+rc193r2c172+c1816αΛ-1δ44+δhhc20D-2rD-11±1-16DD-4αΛ/D-22±c13rD-3dt2+3c18+rc19+2r2Λ3r2δ44+δhhD-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i,where ci, i=6⋯19 are constants.


## 3.1. Flat Horizon Metric with One Variable

Let us assume the first metric spacetime possessing one unknown function has the form(13)ds2=-ardt2+1ardr2+r2∑i=1D-2dθ2i,where a(r) is an unknown function of the radial coordinate, r. Using (13) we get the Ricci scalar in the form(14)R=-r2a′′+2D-2ra′+D-2D-3ar2,where a′=da(r)/dr and a′′=d2a(r)/dr2. The nonvanishing components of the Maxwell-f(R) field equations, (5) and (6), when f(R)=R+αR2 take the form (the detailed calculations of the Ricci curvature tensor are given in Appendix B)(15)ℵtt=12r4αr4a′′2-4r4Λ+2r4q′2+2sθD-42-D-2r22rαa′a′′-ra′-aD-3-4αr2aa′′D-23D-11-4αraa′D3-10D2+24D-4-4αa24D-11-2αr3a′′′ra′+6aD-2+2raa′′′′-2αr2a′2D-2D-5=0,ℵrr=12r4αr4a′′2-4r4Λ+2r4q′2+2sθD-42-D-2r22rαa′a′′-ra′-aD-3-8αr2aa′′D-22-4αraa′D3-7D2+9D+8-4αa22D-12D-7-2αr3a′′′ra′+2aD-2-2αr2a′2D-2D-5=0,ℵθ1θ1=Sθ2θ2=⋯=SθD-2θD-2=12r4r4αa′′-1a′′+4r4Λ+2r4q′2+2sθD-42+4αr3a′a′′3D-7-2D-3r3a′-r2αa′2D-22D-9-4αr3a′′′ra′+a3D-7+raa′′′′-2αa8D2-13D+34ra′-a2D-1D-3D-6=0,where sθD-4=ds(θD-4)/dθD-4 and q′=dq(r)/dr with q(r) and s(θ) being two unknown functions related to the electric and magnetic charges of the system and defined from the gauge potential as(16)A=def.qrdt+sθD-4dθD-4.The solution of the differential equations (15) has the form(17)qr=c1+δ44c2r,sθD-4=c3+δ44c4θD-4,ar=2r2Λ3+c5r+c22+c42r21-16αΛδ44+δhhD-2rD-11±1-16DD-4αΛ/D-222DD-1D-4αrD-3+c6rD-3,where h≥5 and ci, i=1⋯6 are constants. Solution (17) is decomposed into two parts one for D=4 and the other for D≥0. The reason for this decomposition is the electromagnetic field. In the four dimensions there is a charged solution; however, for D≥5 there is no charged solution. In fact this is in conflict with the spherically symmetric case [91]. It is important to mention here that solution (17) when α=0 will reduce to the well known AdS/dS solution in the case of D=4. However, in the case of D>4 this solution does not allow the parameter α to be vanishing; therefore, it has no analogue in GR. In the noncharged case and when D=4 solution (17) coincides with what is know in GR. However, when D>4 and as long as α and Λ have no relation between them solution (17) has no analogue in GR but when Λ=1/16α it will coincide with what is known in GR; that is, it behaves as Ads/dS.The metric of solutions (17) has the form(18)ds2=-2r2Λ3+c5r+c22+c42r21-16αΛδ44+δhhD-2rD-11±1-16DD-4αΛ/D-222DD-1D-4αrD-3+c6rD-3dt2+dr22r2Λ/3+c5/r+c22+c42/r21-16αΛδ44+δhhD-2rD-11±1-16DD-4αΛ/D-22/2DD-1D-4αrD-3+c6/rD-3+r2∑i=1D-2dθ2i.


## 3.2. Metric with Two Variables

The metric of a flat horizon with two unknown functions has the form(19)ds2=-ardt2+1brdr2+r2∑i=1D-2dθ2i,where a(r) and b(r) are two unknown functions of the radial coordinate, r. Using (19) we get the Ricci scalar in the form(20)R=-a′b′2a+ba′′a-ba′22a2+D-2ba′ar+D-2b′r+D-2D-3br2. Using (19) we get the nonvanishing components of the field equations, (5) and (6), when f(R)=R+αR2 in the form (21)ℵtt=18r4a436αr4a2b2a′′2-2αr3baa′329rb′+26D-2b+8rba3a′′′′+8a3asθD-42+r4bq′2+49αr4b2a′4+αr2a′2a22bD-225D-18b+30rb′+9r2b′2-8αa2br3ab′′′2D-2a+ra′-2a′′′3rba′-a2D-2b+3rb′+4r2αaa′′rba′a27rb′+22D-2b-29rba′-a216D-2rbb′+3b′2+4D-2D-4b2-4r2a2αb′′2D-2a222D-7b+rb′+raa′6D-2b+rb′-6r2ba′2-8D-2rαa3a′rb′+2D-4brb′-b-4a3D-2rab′3D-10αrb′-r2-23D-11D-4+43D-10D-5αb+a4r4Λ-D-2D-3r2b-44D-11αb2=0,ℵrr=18r4a44αr4a2b2a′′2+2αr3baa′33rb′-4D-2b+8a3asθD-42+r4bq′2-7αr4b2a′4+αr2a′2a2r2b′2+4bD-23D-4b+rb′-8αa2b2r3a′′′2aD-2+ra′-4r2αaba′′ra′2arb′-3D-2b-3rba′+4D-2a2rb′+D-2b+ab′′2D-2a+ra′2-4D-2ra3a′23D-8αrbb′-r2b-4αb2-αr2b′2+4a3D-22αr2ab′2-a4r4Λ+D-24rbαD-2D-4-D-3r2b-44D-9αb2=0,ℵθ1θ1=ℵθ2θ2=⋯=ℵθD-2θD-2=18r4a428αr4a2b2a′′2-2αr3baa′323rb′+63D-7b+8rba3a′′′′-8a3asθD-42+r4bq′2+39αr4b2a′4+αr2a′2a22bD-22D-2D-5+2D-5D-4b+43D-7rb′+7r2b′2-2r2bα-8αa2br3ab′′′2D-2a+ra′-a′′′5rba′-a22D-5b+6rb′+4r2αaa′′rba′a28D-19b+22rb′-23rba′-a229D-22rbb′+3r2b′2-42D-3D-5-D-4b2-r2bα-8r2a2bb′′-4r2a2αb′′2D-2a24D-4b+rb′+raa′24D-9b+rb′-5r2ba′2+2ra3a′r3b′+2D-3r2b-42D-1D-5+D-43D-11αrbb′-23D-7r2αb′2+85D-19αb2-4a33D-4D-2αr2ab′2-D-3rab′r2+42D-1αb+a4r4Λ-D-4D-3r2b-2D-6D-32D-1-DD-4D-5αb2=0.The above system of differential equations (21) has the following solution:(22)1qr=c7,sθD-4=c8+δ44c9θD-4,ar=c10r2r3Λ+3c1116αΛ-1-3c122r2δ44+δhhc20D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3,br=r2r3Λ+3c1116αΛ-1-3c1223r216αΛ-1δ44+δhhD-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3,2qr=c14+δ44c15r,sθD-4=c16+δ44c17θD-4,ar=c1522Λr4+3c18+rc193r2c172+c1816αΛ-1δ44+δhhc20D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3,br=3c18+rc19+2r2Λ3r2δ44+δhhD-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3,where h≥5. We must mention here that solution (22) has the same property of solution (17); that is, when α=0 it will reduce to the well known AdS/dS solution in the case of D=4. However, in the case of D>4 this solution does not allow the parameter α to be vanishing; therefore, this solution has no analogue in GR. In the noncharged case and when D=4 and α≠0 solution (22) coincides with what is known in GR. However, when D>4 and as long as α and Λ have no relation between them solution (22) has no analogue in GR but when Λ=1/16α it will coincide with GR solution that behaves as Ads/dS. The metric spacetimes of solutions (22) have the form(23)1ds2=-c10232r4αΛ2+48c11αΛr-2r4Λ+3c122+3c11rr2δ44+δhhc20D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3dt2+32r4αΛ2+48c11αΛr-2r4Λ-3c122-3c11r3r216αΛ-1δ44+δhhD-2rD-11±1-16DD-4αΛ/D-222DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i,2ds2=-c1522Λr4+3c18+rc193r2c172+c1816αΛ-1δ44+δhhc20D-2rD-11±1-16DD-4αΛ/D-22±c13rD-3dt2+3c18+rc19+2r2Λ3r2δ44+δhhD-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i,where ci, i=6⋯19 are constants.


## 4. Physical Properties of the Analytic Solutions

### 4.1. Metric with One Variable

From (18) we can deduce the following properties.(i) In case of 4 dimensions we get(24)ds2=-2r2Λ3+c5r+c22+c42r21-16αΛdt2+dr22r2Λ/3+c5/r+c22+c42/r21-16αΛ+r2dθ21+dθ22,from which it is clear that the metric behaves asymptotically as dS/AdS. Equation (24) shows that the effect of the higher dimension curvature is related to the electric field as well as the magnetic field and also (25)Λ≠116α. (ii) In case of more than 4 dimensions we get(26)ds2=-D-2rD-11±1-16DD-4αΛ/D-22+2DD-1D-4αc62DD-1D-4αrD-3dt2+2DD-1D-4αrD-3dr2D-2rD-11±1-16DD-4αΛ/D-22+2DD-1D-4αc6+r2∑i=1D-2dθ2i.Equation (26) shows that the dimensional parameter α must not be equal to zero; otherwise, we will have a singular metric. Also the spacetime of metric (26) behaves as dS/AdS and when the cosmological constant Λ takes the form(27)Λ=D-2216DD-4α, then, (26) reduces to(28)ds2=-D-2rD-1+2DD-1D-4αc62DD-1D-4αrD-3dt2+2DD-1D-4αrD-3dr2D-2rD-1+2DD-1D-4αc6+r2∑i=1D-2dθ2i, which is asymptotically dS/AdS and cannot reduce to GR.


### 4.2. Metric with Two Variables

From (23), we can show the following.(iii) In case of 4 dimensions we get(29)1ds2=-c1022r4Λ16αΛ-1+48c11αΛr+3c122+3c11rr2dt2+3r216αΛ-12r4Λ16αΛ-1+48c11αΛr-3c122-3c11rdr2+r2dθ21+dθ22.2ds2=-c1522Λr4+3c18+rc193r2c162+c1816αΛ-1dt2+3r23c18+rc19+2r4Λdr2+r2dθ21+dθ22.First set of equations (29) shows that Λ≠1/16α and when the constant c10=1/3(16αΛ-1) then the metric behaves as dS/AdS. For the second set, it is allowed to put Λ=1/16α and when the constant c15=c16 then the second equation of (29) behaves as dS/AdS.(iv) In case of more than 4 dimensions we get from (23) the two sets(30)ds2=-D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3dt2+D-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i.Equation (30) shows that the dimensional parameter α must not be equal to zero; otherwise, we will have a singular metric. The asymptotes of (30) behave as dS/AdS. It is important to stress that metric (30) cannot reduce to that of GR and hence we can say that solution (22) is a new solution. Using condition given by (27) in (30) we get(31)1ds2=-±D-2rD-1±c13rD-3dt2+±D-2rD-12DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i. We can conclude from the above discussion of the metrics given by (30) or (31) that there is no charged solution for the form of f(R)=R+αR2. Also spacetime metrics of (30) or (31) instruct us that the dimension parameter α must not be equal zero.


### 4.3. Singularities

Now, let us explain the singularities and the horizons of solutions (17) and (22). For this reason, we have to find at which values of r do the functions a(r) and b(r) turn out to be zero or infinity, due to the fact that singularities may be coordinate ones which are physical singularities. Usually, to study singularities one calculates all the invariants constructed from Riemann tensor and its contractions. The curvature invariants that arise from solution (17), in case of 4 dimensions, have the form(32)RμνλρRμνλρ=4r21-16αΛ28Λ2r6+9c52+24c22+c4261-16αΛ+7c22+c423r81-16αΛ2,RμνRμν=44r8Λ21-16αΛ2+c22+c422r81-16αΛ2,R=-8Λ.And in case of more than 4 dimensions we get(33)RμνλρRμνλρ=r4a′′+2D-2r2a′2+2D-2D-3a2r4,RμνRμν=DD-2r2a′2+2D-2r3a′a′′+r4a′′+4D-2D-3raa′+D-2D-32a22r4,R=-2D-2ra′+r2a′′+D-2D-3ar4.Equation (32), for the 4 dimensions case, shows the following.(a)Λ≠1/16α; otherwise, we will have a singularity for both invariants RμνλρRμνλρ and RμνRμν.(b) Also (32) tells us that there is a singularity at r=0 which can not be removed for the invariants RμνλρRμνλρ and RμνRμν.WhenD>4 and by using (17), we can show the following.(c) We have a true singularity atr=0.(d) From (33), after using (17), we can show that the dimension parameter α must not be equal to zero; otherwise, we get a singularity.(e) Finally, (33) shows that there is a singularity at (34)c6=-D-2rD-12DD-1D-4α1±1-16DD-4αΛD-22, and when constraint (27) is used, (34) takes the form(35)c6=-D-2rD-12DD-1D-4α.Repeating the same calculations for solution (22) we get for D=4(36)1RμνλρRμνλρ=4r21-16αΛ28Λ2r6+9c112+9c11c122r1-16αΛ+168c1143r81-16αΛ2,2RμνλρRμνλρ=36r2c19+144rc19c18+168c184+32r8Λ23r8,1RμνRμν=4c124-128r8Λ3α1-8αΛ+4r8Λ2r81-16αΛ2,2RμνRμν=4c182+4r8Λ2r8,1R=-8Λ2R=-8Λ.We can apply the same discussion applied for solution (17). In case of D>4 we get(37)RμνλρRμνλρ=frr2D-1α2D-2rD-11±1-16DD-4αΛ/D-22±2DD-1D-4αc132,RμνRμν=f1rα2D-2rD-11±1-16DD-4αΛ/D-22±2DD-1D-4αc132,R=f2rαD-2rD-11±1-16DD-4αΛ/D-22±2DD-1D-4αc13,where f(r), f1(r), and f2(r) are lengthy functions of r. Using (27) in (37) we get(38)RμνλρRμνλρ=D-22r2D-1+2DD-12D-3D-42α2c1322DD-1D-42α2,RμνRμν=D-224DD-42α2,R=-D-22D-4α,which indicates that the parameter α must not be equal to zero.The horizons of solutions (17) and (22) are the zeros of the metric grr; therefore, all the above singularities are far from these horizons. The study of perturbations of solutions (17) and (22) is an important issue to study their stability and then discuss the formulation of weak cosmic censorship. This problem will be discussed elsewhere.


## 4.1. Metric with One Variable

From (18) we can deduce the following properties.(i) In case of 4 dimensions we get(24)ds2=-2r2Λ3+c5r+c22+c42r21-16αΛdt2+dr22r2Λ/3+c5/r+c22+c42/r21-16αΛ+r2dθ21+dθ22,from which it is clear that the metric behaves asymptotically as dS/AdS. Equation (24) shows that the effect of the higher dimension curvature is related to the electric field as well as the magnetic field and also (25)Λ≠116α. (ii) In case of more than 4 dimensions we get(26)ds2=-D-2rD-11±1-16DD-4αΛ/D-22+2DD-1D-4αc62DD-1D-4αrD-3dt2+2DD-1D-4αrD-3dr2D-2rD-11±1-16DD-4αΛ/D-22+2DD-1D-4αc6+r2∑i=1D-2dθ2i.Equation (26) shows that the dimensional parameter α must not be equal to zero; otherwise, we will have a singular metric. Also the spacetime of metric (26) behaves as dS/AdS and when the cosmological constant Λ takes the form(27)Λ=D-2216DD-4α, then, (26) reduces to(28)ds2=-D-2rD-1+2DD-1D-4αc62DD-1D-4αrD-3dt2+2DD-1D-4αrD-3dr2D-2rD-1+2DD-1D-4αc6+r2∑i=1D-2dθ2i, which is asymptotically dS/AdS and cannot reduce to GR.


## 4.2. Metric with Two Variables

From (23), we can show the following.(iii) In case of 4 dimensions we get(29)1ds2=-c1022r4Λ16αΛ-1+48c11αΛr+3c122+3c11rr2dt2+3r216αΛ-12r4Λ16αΛ-1+48c11αΛr-3c122-3c11rdr2+r2dθ21+dθ22.2ds2=-c1522Λr4+3c18+rc193r2c162+c1816αΛ-1dt2+3r23c18+rc19+2r4Λdr2+r2dθ21+dθ22.First set of equations (29) shows that Λ≠1/16α and when the constant c10=1/3(16αΛ-1) then the metric behaves as dS/AdS. For the second set, it is allowed to put Λ=1/16α and when the constant c15=c16 then the second equation of (29) behaves as dS/AdS.(iv) In case of more than 4 dimensions we get from (23) the two sets(30)ds2=-D-2rD-11-16DD-4αΛ/D-22±1±c13rD-3dt2+D-2rD-11-16DD-4αΛ/D-22±12DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i.Equation (30) shows that the dimensional parameter α must not be equal to zero; otherwise, we will have a singular metric. The asymptotes of (30) behave as dS/AdS. It is important to stress that metric (30) cannot reduce to that of GR and hence we can say that solution (22) is a new solution. Using condition given by (27) in (30) we get(31)1ds2=-±D-2rD-1±c13rD-3dt2+±D-2rD-12DD-1D-4αrD-3±c13rD-3-1dr2+r2∑i=1D-2dθ2i. We can conclude from the above discussion of the metrics given by (30) or (31) that there is no charged solution for the form of f(R)=R+αR2. Also spacetime metrics of (30) or (31) instruct us that the dimension parameter α must not be equal zero.


## 4.3. Singularities

Now, let us explain the singularities and the horizons of solutions (17) and (22). For this reason, we have to find at which values of r do the functions a(r) and b(r) turn out to be zero or infinity, due to the fact that singularities may be coordinate ones which are physical singularities. Usually, to study singularities one calculates all the invariants constructed from Riemann tensor and its contractions. The curvature invariants that arise from solution (17), in case of 4 dimensions, have the form(32)RμνλρRμνλρ=4r21-16αΛ28Λ2r6+9c52+24c22+c4261-16αΛ+7c22+c423r81-16αΛ2,RμνRμν=44r8Λ21-16αΛ2+c22+c422r81-16αΛ2,R=-8Λ.And in case of more than 4 dimensions we get(33)RμνλρRμνλρ=r4a′′+2D-2r2a′2+2D-2D-3a2r4,RμνRμν=DD-2r2a′2+2D-2r3a′a′′+r4a′′+4D-2D-3raa′+D-2D-32a22r4,R=-2D-2ra′+r2a′′+D-2D-3ar4.Equation (32), for the 4 dimensions case, shows the following.(a)Λ≠1/16α; otherwise, we will have a singularity for both invariants RμνλρRμνλρ and RμνRμν.(b) Also (32) tells us that there is a singularity at r=0 which can not be removed for the invariants RμνλρRμνλρ and RμνRμν.WhenD>4 and by using (17), we can show the following.(c) We have a true singularity atr=0.(d) From (33), after using (17), we can show that the dimension parameter α must not be equal to zero; otherwise, we get a singularity.(e) Finally, (33) shows that there is a singularity at (34)c6=-D-2rD-12DD-1D-4α1±1-16DD-4αΛD-22, and when constraint (27) is used, (34) takes the form(35)c6=-D-2rD-12DD-1D-4α.Repeating the same calculations for solution (22) we get for D=4(36)1RμνλρRμνλρ=4r21-16αΛ28Λ2r6+9c112+9c11c122r1-16αΛ+168c1143r81-16αΛ2,2RμνλρRμνλρ=36r2c19+144rc19c18+168c184+32r8Λ23r8,1RμνRμν=4c124-128r8Λ3α1-8αΛ+4r8Λ2r81-16αΛ2,2RμνRμν=4c182+4r8Λ2r8,1R=-8Λ2R=-8Λ.We can apply the same discussion applied for solution (17). In case of D>4 we get(37)RμνλρRμνλρ=frr2D-1α2D-2rD-11±1-16DD-4αΛ/D-22±2DD-1D-4αc132,RμνRμν=f1rα2D-2rD-11±1-16DD-4αΛ/D-22±2DD-1D-4αc132,R=f2rαD-2rD-11±1-16DD-4αΛ/D-22±2DD-1D-4αc13,where f(r), f1(r), and f2(r) are lengthy functions of r. Using (27) in (37) we get(38)RμνλρRμνλρ=D-22r2D-1+2DD-12D-3D-42α2c1322DD-1D-42α2,RμνRμν=D-224DD-42α2,R=-D-22D-4α,which indicates that the parameter α must not be equal to zero.The horizons of solutions (17) and (22) are the zeros of the metric grr; therefore, all the above singularities are far from these horizons. The study of perturbations of solutions (17) and (22) is an important issue to study their stability and then discuss the formulation of weak cosmic censorship. This problem will be discussed elsewhere.


## 5. Total Conserved Charge

In this section, we are going to studythe conserved quantities of the solutions derived in Section3. For this purpose we are going to make a brief review of the geometry used in this calculations. The Lagrangian of Einstein-Cartan theory is defined by [100] (the fundamental entities of this theory Appendix A) (39)Sϑi,Rjk≔-12κRij∧ηij-2Λη,where ϑi is the one-form coframe and Rij is the two-form curvature tensor. Carrying out variation of (39) with respect to the coframes ϑi and Rij lead to [100, 101] (40)Ei≔-∂S∂ϑi=-12κRjk∧ηijk-2Ληi,Bij≔-∂SRij=12κλij,which are the canonical energy-momentum and rotational gauge field momentum, respectively. The translational momentum and the spin 2-forms are defined as (41)Hi≔-∂S∂Ti=0,Eij≔-ϑi∧Hj=0. The conserved quantity of the gravitational field has the form [100](42)ȷξ=12κd​∗dk+ξϑi∧Ti,wherek=ξiϑi,ξi=ξϑi, with ∗ being the Hodge duality, ξ being a vector field ξ=ξi∂i, and ξi being D parameters ξ0, ξ1,…,ξD-1. For the solutions having spinless matter or vacuum ones, the torsion is vanishing, Ti, and therefore the total charge of (42) takes the form(43)Qξ=12κ∫∂Sd∗k. This invariant conserved quantity Q[ξ] was given before in [102–106].Now let us apply (43) to solution (17) and calculate the necessary components; we finally get (44)ϑ0=ardt,ϑ1=drar,ϑ2=rdθ1,ϑ3=rdθ2,⋮ϑD-1=rdθD-2.Using (44) in (42) we get(45)k=a2rξ0dt-ξ1dr-r2arξ2dθ2+ξ3dθ2+⋯+ξD-2dθD-2ar. The total derivative of (45) gives(46)dk=ξ0a′rdr∧dt-2rξ2dr∧dθ1+ξ3dr∧dθ2+⋯+ξD-1dr∧dθD-2.Calculating the inverse of (46) and using it in (46) and (43) after applying the Hodge-dual to the output of dk, we get conservation of the charge in the form (47)Qξt=D-2ΩD-2M8π,Qξr=Qξθ1=⋯=QξθD-2=0, where ΩD-2 is the volume of the unit (D-2)-sphere; we have substitute c5=-2M and have used the regularization method, by putting the physical quantities vanishing, and subtract the resulting expression from the original one to remove any divergence term [100] (physical quantities mean the constants c2, c4, c5, and c6 for solution (17); constants 10, 11, 12, and 13 of first solution of (22); and constants 13, 15, 18, and c18 of the second solution of (22)).Applying the same above techniques we get for the first solution of (22) in case of D=4(48)Qξt=M1-16αΛ,Qξr=Qξθ1=⋯=QξθD-2=0, where c11=-2M/c10 and for the second one we get (49)Qξt=Mc181-16αΛ-c172,Qξr=Qξθ1=⋯=QξθD-2=0,where c19=-2M/c15.In caseD>4 we get (50)Qξt=D-2ΩD-2Mα8π,Qξr=Qξθ1=⋯=QξθD-2=0, where (51)c13=-2MαDD-3D-4c15.


## 6. Thermodynamics of Black Holes

In this section, we are going to study the thermodynamical quantities of solutions (17) and (22). The temperature of Hawking of any solution can be derived by requiring the singularity at the horizon to be vanishing in the Euclidean continuation of the black hole solutions. One can obtain the temperature of the outer event horizon at r=rh, for solution (17) in case of D=4 in the form (52)T=14πdgttrdrrh=14π2rhΛ+c22+c4216αΛ-1r3h, and in case of D>4 we get (53)T=14πdgttrdrrh=rh4πD-21±1-16DD-4αΛ/D-222DD-4α.For the first solution of (22) in case D=4, we get the Hawking temperature as (54)T=14πdgttrdrrh=14π3c102Λr4h1-16αΛ+c112r3h, and for the second solution of (22) in case of D=4, we get (55)T=14πdgttrdrrh=14πc152c18-2Λr4hr3hc1816αΛ-1+c162. In the case of D>4 we get (56)T=14πdgttrdrrh=rhc13D-1D-24π1-16DD-4αΛD-22±1. We give a brief discussion of the entropy of black hole in f(R) gravity. For this purpose, we use the arguments presented in [107]. From the Noether method used to calculate the entropy associated with black holes in f(R) theory that have constant Ricci scalar, one finds [95] (57)S=A4Gf′Rr=rh, where A is the area of the event horizon. Using solution (17) we get for D=4(58)S=1-16Λαπrh2, and for D>4 we get (59)S=rhD-2ΩD-241-D-2D-41+1-16DD-4αΛD-22. Using (17) and (27) in (59) we get the entropy in the form (60)S=rhD-2ΩD-22D-4.For solution (22) and in case of D=4 we get the same expression of entropy as given by (58) and in case D>4 we get(61)S=rhD-2ΩD-241+2αRrh. When (27) is used we get the entropy in the form of (60).Utilizing (52) and (58), the Smarr relation in the extended phase space can be obtained in the case D=4 of solution (17) as (62)M=2TS+ΦQ,whereΦ=c22+c421-Λαrh,Q=c22+c421-Λα. For the first solution of (22) we get the the Smarr relation in the extended phase as (63)M=2TS+ΦQ,whereΦ=c11r1-2Λα,Q=c111-2Λα, where we have put c10=1/3(1-2Λα)2. Finally, for the second solution of (22) we get the Smarr relation as (64)M=2TS+ΦQ,whereΦ=2c18rh,Q=2c18, where we have put c15=c18(1-16αΛ)-c162/16αΛ-1.In the caseD>4, we introduce the extended phase space where the cosmological constant is identified as the thermodynamic pressure while the conjugate quantity is regarded as the thermodynamic volume. We adopt the following definition of pressure that is commonly used in the literatures of extended phase space [108] (65)P=-D-1D-2Λ48π. Using solution (17) in the case D>4 we get the Smarr relation in the extended phase as [108] (66)M=D-2D-3TS-2D-3VP,wherePis  given  by  Eq.65andVhas  the  formV=3D-2rhD-1ΩD-22+D-21-16DD-4αΛ/D-221+1-16DD-4αΛ/D-224αΛDD-1D-42.It is important to note that Smarr relation given by (66) has no charge term because the higher dimension solution given by Eq. (17) has constant electric and magnetic charges. Using constrains (27) in (66) we get (67)V=24rhD-1ΩD-2D-1D-2D-4.Using the same above procedure we get for solution (22) in case of D>4 the the Smarr relation in the extended phase as given by (66) where the volume has the form (68)V=2rhD-1ΩD-2D-4.


## 7. Main Results and Discussion

In this study, we have presented Maxwell-f(R) gravity in D-dimensions and have checked the flat horizon solutions. We have applied two metrics, the first with one unknown and the second with two unknown functions to Maxwell-f(R) field equations using the special case f(R)=R+αR2. The resulting differential equations are solved analytically without any assumption and general solutions containing three classes have been obtained. These classes are classified as follows.(i)
For the metric with one unknown, the electric and magnetic charges of this class are not constant forD=4 and are constants for D>4.(ii)
For the metric with two unknowns, we have two sets. The first one has a constant charge and a nontrivial magnetic field forD=4 and for D>4 we have constant electric and magnetic charges.(iii)
For the second set of solutions (22) and for D=4 the electric and magnetic charges are not trivial and are constants for D>4.To understand the physics of these solutions we wrote the metric of each solution and have shown the following.(a)
We have shown that the metric behaves asymptotically as dS/AdS for the first solution forD=4 and it is not allowed to put Λ=1/16α. When D>4 the metric also behaves asymptotically as dS/AdS and have shown that the following relation between the cosmological constant and the dimension parameter α holds Λ=(D-2)2/16D(D-4)α.(b)
For the metric with two unknowns, whenD=4, we have shown that the metric behaves as dS/Ads for both sets. However, for the first solution Λ≠1/16α but for the second one Λ=1/16α.(c)
WhenD>4, for solution (22), we have shown that the parameter α can not be equal to zero and therefore this solution is new.We have calculated the singularities of each class by calculating the invariants of curvature. For the first solution we have shown that the invariants of curvature have singularity atr=0 and Λ=1/16α. However, when D>4 we have shown that there is singularities at r=0 and α=0. For the first set of (22) the same discussion of solution (17) can be applied for the invariants of curvature. However, for the second set the invariants of curvature have a singularity at r=0 which represents horizon. For D>4 the same discussion of solution (17) can be applied for the two sets.To understand the meaning of the constants of each solution we have calculated the conserved quantities. In the case ofD=4 we have shown that the only nonvanishing physical quantity is the energy. We have shown that the value of energy of solution (17) is not affected by the dimension parameter α while it is affected for solution (22) and reduces to the well known formula of ADM when α=0. When D>4 the energy of solution (17) does not depend on the parameter α while it depends on it for solution (22) and can not be equal to zero.Finally, we have calculated the thermodynamical quantities like Hawking temperature, entropy, and so on. For solution (17) and when D=4 we have shown that the temperature depends on the electric and magnetic charges and for D>4 the parameter α is not allowed to be zero; otherwise, the temperature will be indefinite. We have shown that the entropy, for solution (17), depends on α and when it is vanishing we return to GR. For this solution, we have shown that the first law of thermodynamics is satisfied for general D≥4 using the extended phase space [108].For solution (22) we have shown that temperature will be finite if α=0 and will return to GR after rescaling the constants of integration. For the entropy and when D=4 the same discussion of solution (17) can be applied. The entropy will have a finite value only when the constraint (27) is used. Finally we have shown that the first law of thermodynamics is satisfied for general D≥4.



---

*Source: 7323574-2018-01-21.xml*