Generation of spherically symmetric metrics in $f\left( R\right) $ gravity

In $D-$dimensional spherically symmetric $f\left( R\right) $ gravity there are three unknown functions to be determined from the fourth order differential equations. It is shown that the system remarkably integrates to relate two functions through the third one to provide reduction to second order equations accompanied with a large class of potential solutions. The third function which acts as the generator of the process is $F\left( R\right) =\frac{df\left( R\right) }{dR}.$ We recall that our generating function has been employed as a scalar field with an accompanying self-interacting potential previously which is entirely different from our approach. Reduction of $f\left( R\right) $ theory into system of equations seems to be efficient enough to generate a solution corresponding to each generating function. As particular examples, besides known ones, we obtain new black hole solutions in any dimension $D$. We further extend our analysis to cover non-zero energy-momentum tensors. Global monopole and Maxwell sources are given as examples.


I. INTRODUCTION
f (R) gravity is one of the modified theories of Einstein's general relativity that attracted much attention in recent times [1][2][3][4]. In [5] f (R) = R + αR 2 with α > 0 has been introduced as the model of inflated universe while f (R) = R − α/R n (α > 0, n > 0) was considered as a candidate for the dark energy model [6][7][8][9][10][11]. This model, however, is not a viable model for dark energy and instead f (R) = R − αR n with α > 0 and 0 < n < 1 emerged as alternative which has been proposed in [12,13]. Later on more viable models were studied in [14][15][16][17][18]. A detailed review of these models are given in [19] (Other review papers are given in [20][21][22][23]). Some recent works on solutions in f (R) gravity are given in . For f (R) = R, in D = 4 dimensional spacetime it coincides with the standard general relativity but otherwise it stands with an action that is arbitrarily dependent on the Ricci scalar. Finding exact solutions in this theory with fourth order derivatives of the metric tensor is both important and challenging. Apart from exact analytic solutions there are f (R) models that can only be expressed implicitly in non-polynomial expressions. Each particular model has advantages / disadvantages as far as experimental tests are concerned [47][48][49][50][51][52][53]. There are even models that lack the Einstein's R−gravity limit. Among other expectations the UV / IR behaviors at near / far distances, quantum renormalizability with power counting of the counter terms are prominent. At any cost preferring to abide by the classical regime we confine ourselves first to sourceless (vacuum) f (R) models that admits exact integrals. In the last section we extend our discussion to cover external sources such as global monopole [54] and electromagnetic field. Let us add that the equivalence of f (R) gravity to Brans-Dicke (BD) theory (ω = 0) with a potential has also been high- * zahra.amirabi@emu.edu.tr † mustafa.halilsoy@emu.edu.tr ‡ habib.mazhari@emu.edu.tr lighted extensively in the past as a transition between Jordan and Einstein frames. In this approach the exact solutions can be generated by adopting scalar field ansatzes which in general brings into the Lagrangian intricate potentials. Our method will be confined entirely to the Jordan frame without reference to the BD field or any scalar potential.
The vacuum of f (R) gravity is known to carry its own curvature sources. By vacuum in this theory it is meant the absence of an external energy-momentum tensor T µν , of any physical source [55]. Carames and de Mello in latter work have considered the spherically symmetric vacuum solutions of f (R) gravity in higher dimensions. We shall rederive most of their results anew together with some additional extensions which we are interested in to develop further. Addition of external T µν = 0, no doubt makes the problem technically more complicated, but following the lesson learned from the vacuum / empty solutions of f (R) = R theory of gravity we will attempt to derive the most general equations and in some specific cases the solutions as well. The D dimensional spherically symmetric line element that we shall consider will be in which A (r) and B (r) are metric functions to be determined while dΩ 2 D−2 represents the (D − 2) −dimensional unit spherical line element. In particular integrals we have the restrictive case A (r) = B (r) included, but more general cases with A (r) = B (r) must be interesting as well. Beside the metric functions A (r) and B (r) we shall employ a third function denoted by F (R) = df (R) dR , which characterizes the type of the f (R) gravity. Let us add that not in all cases of f (R) models the explicit form of f (R) can be expressed analytically in terms of the variable R, the Ricci scalar. Instead, it involves a transcendental part that can't be inverted in the form of r = r (R) , these are hybrid forms. We add that even these hybrid forms don't prevent us from calculating df dR > 0 and d 2 f dR 2 > 0 which are crucial terms to deter-mine the absence of ghosts and thermodynamic stability, respectively.
In brief what has been achieved in this paper is to show that the metric functions A (r) and B (r) are related through an integral expression for the function F (R) (or F (r)). This amounts to the fact that once F (r) is given it acts as a generator to generate a new set of (A (r), B (r)) pair. The function A (r) is expressed in terms of B (r) and F (r) and the remaining equation is reduced into a master equation satisfied by B (r). Once we give an ansatz for F (r) our master equation can be integrated in principle to obtain B (r) . In this manner we can obtain an infinite class of metrics in f (R) gravity generated from an infinite set of F (r) . No doubt the dimensionality of spacetime D(= d + 1) also plays role in the derivation. In particular, we present examples of new black hole solutions in D ≥ 3, by the method described above. We wish to add also that in the reduction process the system of differential equations in f (R) gravity reduce naturally from the fourth order to the second order.
The paper is organized as follows. In Section II we rederive the f (R) field equations in D−dimensions which is comparable in some sense with [55]. A number of examples to justify the effectiveness of our method are given. Generalization to T µν = 0 is analyzed in Section III. We end our discussion with Conclusion in Section IV.

II. THE FIELD EQUATIONS IN D−DIMENSIONS
The D−dimensional vacuum f (R) −gravity is represented by the action in which f (R) is a function of Ricci scalar R and D ≥ 3. Variation of the action I with respect to g µν provides the field equations (in metric formalism) in which F = df dR and is the covariant Laplacian. The general spherically symmetric line element is given by (1) and the field equations (3) are explicitly given by and and Herein and in the rest of the paper, D k = D−k, a prime stands for the derivative with respect to r and F = df dR . A particular combination of these three equations leads to two equations which are independent of f . The first equation (4) may be written as which integrates to where The second equation (5) upon considering (11) becomes independent of A too. The closed form of the second equation reduces to a linear equation for B (r) , which is given by where and Finally, the explicit form of f in terms of r is given by To complete our analysis we give the explicit form of the Ricci scalar in terms of r which reads In summary, the only equation to be solved is Eq. (13) which is second order and linear for B. Therefore the procedure is reduced to set a F (r) -which eventually represents the form of f (R) -and solve the only equation (13) to find B (r) and consequently A (r) .

A. Applications of the method
Before we give certain applications for our formalism we would like to compare our approach with the work of Carames and Bezerra de Mello [55]. The main difference can be seen from the fact that in [55] there are two generating functions (so to say) which are F and Y . In other words Eqs. (16) and (17) of [55] are coupled and one must consider them together to find a solution to the field equations. In our formalism we have only one generating function which is F and Eq. (13) is the only equation to be solved. In the following three cases we shall show that for simple cases (F = 1 and F = 1 + αr) our results overlap with [55] but for more complicated case (F = αr a ) our solution is the general one while the solution given in [55] is a restricted one (look at Eq. (46) in [55] and (35) and (36) in this work).

F (r) = 1
We start with the simplest case, with F = 1 or equiv- in which −2Λ is an integration constant to be interpreted as the cosmological constant. The main equations (11) and (13) admit for integration constants C 1 and C 2 and consequently Finally the solution becomes in which we set where M is the ADM mass. The solution for D > 3 is Schwarzschild de / anti de-Sitter black hole solution and for D = 3 it is the BTZ black hole. Let's add that F (r) = ξ =constant, does not change the nature of the solution for the metric with constant scalar curvature. We note that our results in this section expectedly is the same as section 3.1 in [55].
2. F (r) = 1 + αr, with α =constant As of Sec. 3.2 of Ref. [55] we consider F ′′ (r) = 0 in our general formalism but instead of going through a general D−dimensional solution, we investigate the cases in closed form for f (R) . From Eq. (10) with F ′′ = 0 one obtains A = B (up to a constant which one can set it unity via a redefinition of time). For an arbitrary D the solution for B may not be possible in a closed form but for specific dimensions we may find.
a. D = 3 In three dimensional spacetime the solution for B becomes where C 1 and C 2 are integration constants with the curvature scalar The hybrid relation between f (r (R) , R) and R is given by Note also that from the foregoing expressions we can identify the cosmological constant by C 2 = −Λ. Furthermore, setting α = 0, one recovers the previous example with C 1 scaled which suggests that it is related to the mass of the central object. b. D = 4 In four dimensions the solution is given by and R = −12C 2 − 72C 1 α 3 ln 1 + 1 αr + 6αC 1 −1 + 6α 2 r 2 + 6αr (1 + 2αr) We see clearly the role of C 2 = − Λ 3 and we wish to proceed with C 1 = 0 which implies , α = 0 (30) so that This is a black hole solution with a singularity at r = 0 such that If we set Λ = 0, with α < 0 we may introduce a horizon for the solution located at while for α > 0 the solution possesses a naked singularity at r = 0. To complete our investigation let's determine the absence of ghosts and thermodynamic stability of the explicit f (R) found in (31). We see that df Clearly both conditions can not be satisfied simultaneously.

F (r) = αr a
Our next example is a power-law form for F (r) , i.e., with constants α and a which upon (11) yields Substituting into (13) one finds in which C 1 and C 2 are the integration constants. Following A and B we also find and We also note that although α and a are two arbitrary constants a must satisfy a = −D 2 , 1 ± √ D 1 . This is remarkable to observe that in Sec. 3.3.2 of [55] the same ansatz for F has been considered but the solutions to the field equations (see Eq. (46)-(49) of [55]) are not the same as what we found here is more general. As a matter of fact our solutions (35)- (38) with C 2 = 0 reduce to their solutions. This shows that reducing the field equations into a master equation with a single generating function makes advances in finding exact solutions in f (R) gravity.
In D = 3 dimensions the solution becomes rather specific since the last term vanishes for all values of a. The functions then read as and This is nothing but the solution found by Zhang, Liu and Li in [56] with their parameters and Note that which yields f ∼ R 2 for the specific choice a = 1 3 . For D ≥ 4 one may set C 2 = 0 and therefore with an analytic relation for f (R) given by in whichᾱ is a constant which can be set to unity (by a fine choice of α). In the case of a = 0 the theory gives R−gravity. The solution is a black hole with singularity at r = 0. Let us also add that with C 1 = 0 and a = 1 the solution reduces to and This represents a global monopole-type solution with a deficit angle. Another interesting setting is for a = −2 and C 2 = 0 which upon a proper choice of α one finds and the solution becomes demanding AB = 0 which is not acceptable. For D = 10 the particular solution can be obtained as and with the Ricci scalar The metric function A (r) follows accordingly from (11) which is given in (54).
We comment that f =ᾱR 1− a 2 does not satisfy df dR > 0 and d 2 f dR 2 > 0 simultaneously unless we set a to be negative (Note thatᾱ = 1 is needed to have the Einstein gravity recovered.). For instance withᾱ = 1 and a = −2 which implies f = R 2 , and both conditions are satisfied.

A new black hole solution in D = 3
In [56] where f (R) = R d+1 in three dimensional spacetime with d = const. has been studied, the solution does not cover the case d = − 1 2 which makes f (R) = √ R. This can be seen from the Eq. (12) of [56] and Eq. (45) (note that 3−a 2(1−a) = 1 2 has no answer) and in that paper both do not cover the case f (R) = √ R. However, in what follows we wish to show that the solution f (R) = √ R can easily be obtained in three dimensions. To do so let's consider with α and β two real constants. This yields and Note that, C 1 and C 2 are two integration constants such that C 2 effectively plays the role of a cosmological constant. The solution is a black hole with a horizon located with the condition that −C1 C2 > 0. Also from R we see that the solution is singular if and only if C 1 = 0. In the sequel we are interested in C 2 = 0 which makes the solution to be rather simple but singular. Accordingly the form of f and R are given by which upon tuning the free parameter β by 4αβ 2 C 1 = 1 the form of f becomes Here α and C 1 are positive constants and the line element finally reads in which r 0 = 1 4αβ 2 is also a positive constant.

A general class of solutions in 2 + 1−dimensions
In three dimensional spacetime in addition to what we found by now, we wish to show that there exists an important class of solutions yet to be discovered. This specific class, however is a characteristic feature of only three dimensions. To see this solution let's set D = 3 in (13) which yields where we substituted B (r) by A (r) using (11) i.e., Eq. (67) possesses a trivial solution for A (r) irrespective of the form of F (r) which is given by in which C 0 is an integration constant. In this situation the field equations are all satisfied provided B (r) and A (r) satisfy the condition (68) i.e., One can see easily that with F = ξ =constant we obtain The solution given by (69), (71) and with η = const., constitutes a particular class. Note that since the choice of F (r) in (70) is arbitrary this can be used to generate an infinite class of solutions. This will not be searched any further here.
Previously we found an exact solution for the model of gravity in the form with α = 0. The solution to the field equations is given by The solution is a black hole solution with α < 0 and therefore one may write in which the horizon is shown as r + . A change of variables of the form t = √ 2T, r = ρ √ 2 reduces (76) to which is the Schwarzschild black hole with a deficit angle caused by a cosmic string.

III. GENERALIZATION TO f (R) GRAVITY COUPLED TO MATTER SOURCES
In this Section we extend our vacuum analysis to the presence of matter coupled with gravity. Therefore the action becomes in which L m is the matter Lagrangian density. The field equations become with T ν µ = diag (−ρ, p, q, q) the energy momentum tensor of the matter source. The line element is going to be a spherically symmetric as (1) and without going through the details of the field equations, we give the changes in the field equations (11) and (13). The corresponding field equation to Eq. (11) reads while the main equation (13) i.e. the master equation for B (r) takes the form (81) in which P , Q and H are given in (14), (15) and (12), respectively and Note that, unlike the source-free case, here the master B equation is not a linear equation. It can also be observed that the particular choice of ρ + p = 0 removes the nonlinearity in the B equation (81). Further choice of p = q leaves us with the same equation for B as in the sourceless case. Yet with q = 0, the source shows itself in the Ψ function as given in the sequel.
The closed form of f is given by The foregoing expressions are not impressive much unless we provide concrete examples. This is our aim in the section that follows.
A. Applications One immediate application of Eq. (81) is the extension of the f = R + 2α √ R in 3 + 1−dimensional vacuum to the gravity coupled to the global monopole [54] whose energy momentum at very large distance is given by (we assume A = B in the line element) in which η represents the global monopole charge. The solution for the metric is given by and F (r) = 1 + αr.
Again we must have α = 0 and in the limit η = 0 one recovers the vacuum solution.

2 + 1−dimensional Maxwell electric field coupled to f (R) gravity
Let's consider now Maxwell electric field to be coupled to f (R) gravity in 2 + 1−dimensions. Using the line element (1) and the standard Maxwell Lagrangian together with F = 1 + αr one finds in which Q is the electric charge. The solution with small α up to first order simply reads and Clearly by setting α = 0 one recovers the charged BTZ black hole with cosmological constant Λ = C 2 . We add that the solution when the energy momentum tensor of the matter source is of the form of a fluid, one should consider it as an interior solution. Therefore one has to make sure that the Israel junction conditions are satisfied at the interface between exterior and interior solutions [57].
In two examples we studied here, the energy momentum tensors are long range fields which allow us to consider our solutions to be exterior.

IV. CONCLUSION
The integrability of f (R) vacuum gravity with spherical / circular symmetry is reduced first to a set of master relations, i.e., Eq. (11) and a master equation for B (r) , i.e., Eq. (13). Given any ansatz generating function F = df dR in terms of the coordinate r, our method generates a solution pair of A (r) , B (r) and a solution for the function f (R). Most of our solutions encountered are of hybrid nature, that is f (R) can't be expressed explicitly in terms of R. The power-law form for instance, of the form f (R) ∼ R k , with a rational number k, is obtained easily in our method. Some of the solutions presented as applications are already known. Yet, new and rare type of solutions can also be obtained easily. In the second part of the paper, we extend the integrability of vacuum case to the non-vacuum f (R) theories. For this case we found also a master relation, i.e., Eq. (80) and a master equation, i.e., Eq. (81). Considering F (R) to be the generating function for our formalism and T µ ν = diag [−ρ, p, q, q] to be our energy momentum tensor one finds a solution to the master equation (81). We presented two examples. In the first example we set F = 1 + αr with a global monopole coupled to the gravity in 3 + 1−dimensions. The second example considers the Maxwell electric field coupled to the gravity in 2 + 1−dimensions with the same generating function. In this case we found the solutions approximately for small α. Let's also add that among the few examples we have studied we found some closed form of f (R) such as f (R) = R + 2α √ R − 2Λ − 2Λ and f (R) =ᾱR 1− a 2 . We found that the first one can not satisfy the conditions for the absence of ghosts and thermodynamic stability simultaneously while the second one for specific a does satisfy. Finally it will be in order to state that our method of reduction for the spherically symmetric f (R) gravity has a large scope as far as solutions are concerned. Physical implications of the solutions obtained such as dark matter connection is not considered in the present study.
Extension of our formalism to problems with different symmetries, such as stationary axial symmetry requires a separate investigation.