Spherically symmetric brane spacetime with bulk f(R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathcal {R})$$\end{document} gravity

Introducing f(R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathcal {R})$$\end{document} term in the five-dimensional bulk action we derive effective Einstein’s equation on the brane using Gauss–Codazzi equation. This effective equation is then solved for different conditions on dark radiation and dark pressure to obtain various spherically symmetric solutions. Some of these static spherically symmetric solutions correspond to black hole solutions, with parameters induced from the bulk. Specially, the dark pressure and dark radiation terms (electric part of Weyl curvature) affect the brane spherically symmetric solutions significantly. We have solved for one parameter group of conformal motions where the dark radiation and dark pressure terms are exactly obtained exploiting the corresponding Lie symmetry. Various thermodynamic features of these spherically symmetric space-times are studied, showing existence of second order phase transition. This phenomenon has its origin in the higher curvature term with f(R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathcal {R})$$\end{document} gravity in the bulk.


Introduction
Our four dimensional world might be embedded in a five dimensional space-time was proposed in [1,2] in order to explain the observed hierarchy between Electroweak and Planck scale. Such extra dimensional models also have their origin in some suitable compactifications of ten dimensional E 8 × E 8 heterotic string theory [3].
This scenario has attracted considerable attraction due to its elegant nature and simplicity. In this brane world scenario the standard model fields are confined on a 3brane, while gravity can propagate both in the brane and the bulk. A single 3-brane, which is embedded in a five dimensional bulk has the five dimensional line element, ds 2 = e −A(y) η μν dx μ dx ν + dy 2 . The warp factor e −A(y) can be tuned properly to induce Einstein gravity on the brane a e-mail: sumantac.physics@gmail.com; sumanta@iucaa.ernet.in b e-mail: tpssg@iacs.res.in as a leading order term. We could have also considered a two brane system, which comes with an additional field known as radion, representing separation between the branes, with interesting features [4,5]. However we will restrict ourselves only to the single brane system for the rest of the discussion.
However due to the presence of extra dimensions, we should expect deviation from Einstein theory, which play a significant role at high energies [6,7]. Gravity sector also gets modified at electroweak scale ∼1 TeV, changing the cosmological implications, which have been extensively studied in Ref. [8][9][10][11][12]. The effect of extra dimension on formation of black hole has been studied in References [13][14][15]. Also these models have very interesting properties from the point of view of particle phenomenology [16][17][18][19][20].
In General Relativity the exterior space-time of a spherically symmetric black hole or a compact object is standard Schwarzschild geometry. However due to the presence of an extra dimension in the brane world scenario the Schwarzschild solution gets modified non-trivially. This originates due to high energy corrections, Weyl stress on gravitons propagating in the bulk. One such solution was obtained in References [21], in the form of Reissner-Nördstrom solution. The interior solution can be matched to a brane world star having constant energy density [22][23][24]. A non singular solution for black holes in these models can be obtained by relaxing the condition of zero scalar curvature while retaining null energy condition [25,26]. Also the Gauss-Codazzi equations can be solved in Randall-Sundrum type II model to get exterior solution for spherically symmetric star [27]. The various classes of vacuum solutions has been obtained in Reference [28] by solving the vacuum field equations on the brane obtained from Gauss Codazzi equation. The results of various such calculations suggest that brane world black hole horizons has the peculiar structure of a "pancake".
In recent years, there has been a new concept in General Relativity suggesting modifications of Einstein-Hilbert action in order to explain the late time cosmic accelera-tion to inflation. This is achieved by introducing higher curvature terms in the action, and a very promising candidate among such modifications is f (R) gravity theories (for recent reviews see [29][30][31]). The main difficulties with these modifications are, they become infected with ghost modes. However f (R) theory on a constant curvature hyper surface is shown to be ghost free [32][33][34][35]. The modification due to introduction of f (R) term in the Lagrangian can address variety of problems e.g. four cosmological phases [36,37], late time cosmic acceleration [38,39], initial power law inflation [40], rotation curves of spiral galaxies [41,42], detection of gravitational waves [43][44][45] and many others. This theory has also the potential to pass through all known tests of general relativity.
Motivated by such striking properties of f (R) gravity it is also introduced in brane world models, where the five dimensional action is modified by introduction of f (R) term in the bulk, with R being the Ricci scalar of the five dimensional theory. In particular for bulk geometry with high curvature ∼ Planck scale, such higher order corrections to gravity are expected to become extremely relevant. Effective gravitational equations on the brane have been obtained in References [46][47][48][49] while perturbations on the scalar and tensor modes on the brane has been studied in Reference [50]. Cosmology on these brane world models Reference [51] along with brane world sum rules have also been discussed in these f (R) gravity models [52]. The nature of warped geometric models in this f (R) gravity theory with constant bulk curvature has been obtained in Reference [53] and the graviton KK mode masses in these models have been examined in the light of recent ATLAS data in LHC.
Ever since the pioneering works of Regge and Wheeler [54][55][56], the stability of a four dimensional black hole under linear perturbation has been investigated extensively. The importance of linear stability of a black hole can be understood as follows: the black hole solutions should describe the final state of gravitational collapse and thus they should be stable against small fluctuations. Also technically, this implies that at the order of linear perturbation, Einstein equation reduces to a simple set of wave equations. For the static situation, these equations resemble Schrödinger equation with time dependent Hamiltonian. Thus the stability analysis becomes equivalent to a simple, quantum mechanical problem. We also mention that there are solutions which describe naked singularity, and stability of a naked singularity is an important issue from the viewpoint of cosmic censorship conjecture. In this work we will use the wave equations to study the stability [57][58][59][60].
An important aspect of black hole physics, pioneered by Bekenstein, shows a remarkable similarity between black hole and a thermodynamic system. The similarity arises from the fact that just like a thermodynamic system one can attribute temperature to a black hole (known as Hawking temperature) which is proportional to the surface gravity and also an entropy proportional to the horizon area [61][62][63][64][65][66]. Any arbitrary black hole can be characterized by three parameters, its mass, charge and angular momentum. The thermodynamic stability of such a system can be determined by the sign of heat capacity just like any normal thermodynamic system. For a black hole the criteria c v < 0 makes the system thermodynamically unstable. However if the specific heat changes sign as well as diverges in its parameter space, then it indicates a second order phase transition [67,68]. Phase transitions in various black hole solutions have been studied extensively in Einstein gravity as well as in alternative gravity theories [69][70][71][72][73][74][75].
The purpose of this work is to consider various spherically symmetric vacuum space-times on the brane obtained from f (R) action on the bulk. In order to achieve this we consider the decomposition of electric part of the Weyl tensor into dark radiation and dark pressure terms. It turns out that these determine the space-time geometry we are considering. Moreover some simple integrability conditions lead to different classes of vacuum solutions. These issues are addressed in Sects. 2 and 3. Then we have discussed stability of black holes and naked singularities in these spacetime in Sect. 4.
Next we consider vacuum space-time related to Lie groups of transformation. As a simple situation we consider spherically symmetric and static solutions with the metric tensor admitting one parameter group of conformal motion. With proper integrability condition an exact solution corresponding to a brane with one parameter group of motions can be obtained (see Sect. 5).
Finally we consider the thermodynamics of these black hole solutions. As these solutions are induced on the brane due to bulk action, the thermodynamic properties are related to the dark pressure and radiation terms coming from the electric part of Weyl tensor and thus the thermodynamic properties of the brane black holes are directly related to those of bulk space-time (see Sect. 6). We finally conclude with a discussion on our results.

Static, spherically symmetric field equations on the brane
To obtain the vacuum solution we start from the bulk action with f (R) term as, where L m is the matter Lagrangian, g AB is the bulk metric and R is the bulk Ricci scalar. The bulk indices A, B runs through 0 . . . 4 i.e. over all the space-time dimensions. The variation of the action S with respect to bulk metric g AB leads to, Here the negative vacuum energy density on the bulk and the brane energy-momentum tensor are the sources of the gravitational field. Eq. (2) can be put into the form, where τ μν is the brane energy-momentum tensor and λ T is the corresponding brane tension. Also the quantity h μν is the induced metric on y = constant hypersurfaces. The effective four-dimensional gravitational equations on the brane are, where, with, Note that for f (R) = R, we retrieve the usual Gauss-Codazzi equation for a pure Einstein gravity in the bulk. We now proceed to simplify the expression for Q μν . The normal to y = constant hypersurface being n A = ∂ A y, we have n μ = 0. In addition if we assume that ∂ μ R = 0 then using the relations: Now the scalar curvature for the bulk must be a well behaved quantity, and we can expand it in a Taylor series around y = 0 hypersurface, as, R = R 0 + R 1 y + R 2 y 2 /2 + O(y 3 ). Since bulk curvature depends only on the extra dimension y, all the coefficients are constants. Thus all the derivatives calculated at y = 0 yield a constant contribution which does not depend on any of the brane coordinates.
The electric part of the Weyl tensor E μν has its origin in the nonlocal effect from free bulk gravitational field. This is the projection of bulk Weyl tensor such that, E AB = C ABC D n C n D along with E AB = E μν δ μ A δ ν B on the brane (y → 0). From the Gauss-Codazzi equation we also have conservation of energy momentum tensor as, D μ T μν = 0, where D μ is the brane covariant derivative. This also imposes restrictions on projected Weyl tensor from Bianchi identities. Following Reference [6] the projected Weyl tensor can be expanded as, with k = k 5 / √ 8π G N and ξ μν = h μν +u μ u ν . This decomposition is with respect to the four velocity field u μ . The respective terms in the above expression are, the "Dark Radiation" term, is a spatial, trace free, symmetric tensor. For static solutions, Q μ = 0, while the constraint becomes dependent on dark radiation U (r ), vector A μ = A(r )r μ and a ten- Here r μ is unit radial vector.
In order to obtain solution in a source free region on the brane, brane energy momentum tensor appearing on the right hand side of effective Einstein's equation is taken to be zero. Thus we readily obtain τ μν = 0 = π μν . Also from the previous discussion it is evident that R is dependent only on y and on the brane (at y = 0) all its derivatives with respect to coordinates become constants. Then the Einstein equation becomes, Now we choose an ansatz for spherically symmetric solution in the form, ds 2 = −e ν(r ) dt 2 + e λ(r ) dr 2 + r 2 d 2 (13) For this choice the effective Einstein's equation and energy-momentum conservation equation on the brane become, where we have denoted a ≡ da/dr . Now Eq. (14) can be solved for e −λ to yield, where C 1 is an arbitrary constant of integration. The quantity Q(r ) is defined as, We can interpret the term Q as equivalent to gravitational mass originating from dark radiation and henceforth will be referred as dark mass. In the limit f (R) → R, 4 → 0 as well as U → 0 we retrieve the standard Schwarzschild solution. This helps us to identify the arbitrary constant as C 1 = 2G M, M being the constant mass of the gravitating body. Also we can obtain the differential equations that are satisfied by dark radiation U (r ) and dark pressure P(r ) in static spherically symmetric space-time. Eliminating ν from Eqs. (17) and (15) and using e −λ from Eq. (18) we obtain: where we introduce two extra parameters, α = (1/4π Gλ T ) and χ = F(R)− 4 . Now we define the following quantities in order to transform the above differential equation into a more convenient form which will be used extensively later, In terms of these variables the differential equations satisfied by the dark radiation and dark pressure are, Thus the Eqs. (14)- (17) are the effective field equations, on the brane, while the Eqs. (23)-(24) represent equations for the source terms in the bulk i.e. dark pressure and dark radiation.

Various classes of solutions on the brane
Equations (20) and (21) can not be solved for dark radiation U and dark pressure P simultaneously unless we have a relation connecting them. We therefore choose some possible relations between the dark radiation U and dark pressure P which essentially define different equations of state. For different such choices we get different solutions. In this section we impose certain conditions on dark radiation U and dark pressure P, to obtain the corresponding solution. It turns out that the solutions are very distinct for different choices.

Case-I: U = 0
This condition comes with vanishing dark radiation, which imply readily Q = 0. In this scenario, one of the metric elements can be given by, The differential equation satisfied by the dark pressure P(r ) is given by, while the differential equation satisfied by ν is given by, Solution for these two differential equations give the pressure and metric for this case. Note that in this situation the metric element e ν is solely determined from the pressure, which can be seen directly from Eqs. (27) and (26) as, ν = −2P /P − 6/r . This equation can be integrated to yield, exp(ν) = C 2 /r 6 P 2 , where C 2 is an arbitrary constant of integration. Thus once pressure equation is solved, the metric element is also known.
In order to obtain the pressure two quantities r 1 and d would be important with the following expressions: With these variables the solution for the pressure is obtained as: From the above expression it is evident that at r = r 1 the metric element e ν vanishes. Thus the space-time has an event horizon located at r = r 1 with its characteristic thermodynamic features.
3.2 Case-II: P = 0 In this situation Eqs. (23) and (24) reduces to the following form, These two equations can be combined to yield a single differential equation such that, The transformations dq/dθ = 1/v and v = w(6 − 6q + ) −2/3 lead to the following differential equation, The above differential equation has a particular solution, w = − 1 q (6 − 6q + ) 2/3 . However for a wider class of solutions we define a new variable η = (6 − 6q + ) −1/3 . This leads to the differential equation, It is hard to find an exact solution of this differential equation. Therefore we resort to approximated methods. For that purpose we choose the differential equation (33) and making Laplace transform of this equation we get, Then using the convolution theorem in the form, we readily obtain the following integral solution, where we have the following functions, q 0 (θ ) = A 1 e −θ + A 2 e 2θ (40) Having obtained an integral solution we now move forward to determine the metric. However the solution is usually obtained by successive approximation methods, which invokes iterations. At zeroth order we get the solution by using only the linear part of the differential equation (33) and will be denoted by q 0 . Then we can write our full solution as a limiting process, such that q(θ ) = lim m→∞ q m (θ ). In this situation for m ∈ N , we have the iterative solution at mth order connected to (m − 1)th order by the following integral equation, Then following Ref. [28] the zeroth order static and spherically symmetric solution to the field equations turn out to be, where we have C 0 as an arbitrary integration constant. After using one more iteration i.e. up to first order approximation the metric components are obtained as, Note that the dependence on f (R) gravity appears through the A 1 factor. However the dependence is quiet complicated and affects both the metric elements.
where P 0 is an arbitrary integration constant. Also the dark mass can be calculated from Eq. (21) as, where again Q 0 is an integration constant. For this particular choice we have from Eqs. (14) and (15) ν = −λ . Hence the metric elements are given by, This solution has several interesting features which we discuss now. Firstly this solution is asymptotically dS (AdS) or flat depending on the sign of (F(R) − 4 ) being negative (positive) or zero. Then there is an analogous charge term which is the coefficient of 1/r 2 term and is given by −3α P 0 /2. Finally we have a mass term given by, 2G M + Q 0 . Thus we note that the charge term is coming solely from the dark pressure term and thus has its origin in the bulk geometry. Similar argument hold true for the mass term also. However the effect of f (R) gravity on the bulk actually induces a dS (AdS) nature to the vacuum solutions.
3.4 Case-IV: U + 2P = 0 Here we consider a different condition on the dark radiation and dark pressure terms. In this case Eq. (20) leads to the expression for the dark mass Q as, along with the the solution for dark radiation term and dark pressure term as, The metric elements in this case can be evaluated as, e ν = C 0 r 2 (56) Note that this solution actually represents a naked singularity since the event horizon is determined by the equation, e ν = 0. Thus though the f (R) model modifies the e λ term however it yields a naked singularity solution. Moreover e −λ = 0 determines the null surface, however in this situation the null surface exists only if 4 > F(R) and is located at, r h = √ 4 − F(R). Hence by imposing appropriate conditions we obtain either black hole solution with event horizon or solution with naked singularity.
In this context we should mention that naked singularities are just not some artifact, these can be used to probe structures as well. For example we can use naked singularity to take part in gravitational lensing and time delay, with centroid deformation of astrophysical objects [76,77].

Stability of the solutions
Stability of black holes under linearized perturbation is considered as an important problem in black hole physics. Here we consider gravitational perturbation in a static spherically symmetric background. Gauge invariant formalisms were developed in an arbitrary static background metric having the form −g tt = g rr = f (r ). It turns out that for certain ranges of the parameter space the Hamiltonian is positive guarantying the self-adjoint extension of it under suitable boundary condition.
The perturbation can be grouped into three types: scalar, vector and tensor perturbations. Expansion of each of these perturbations in harmonic functions leads to a set of equations expressed in terms of gauge covariant variables. Further reduction of these equations then reduces them to a set of decoupled wave equation in the form: where as usual, represents the d'Alembertian operator with respect to the two dimensional metric. Also = S , V and T represent scalar, vector and tensor perturbations respectively. The potential function for each of these perturbation modes corresponds to [57]: where we have used the following expressions: It should be noted that the total number of independent components of the scalar, vector and tensor modes adds up to 2, the number of independent degrees of freedom for graviton in the brane. Since the tensor mode has no degrees of freedom we need to concentrate only on the vector and scalar modes.
Let us now consider the black hole and naked singularity solutions obtained in the previous section using effective gravitational field equation on the brane. Most of these solutions are quiet complex and we shall focus into some appropriate limiting cases. • Next we discuss the case of vanishing dark pressure. In this case the solutions are not exact and even the zeroth order solution for e −λ looks like Schwarzschild de-Sitter. However the other one is merely a constant. Thus from the expressions for the potential it turns out they depend on the e −λ at the outside and thus will represent stable solution for the range of parameter space where f (r ) > 0. From large r limit we observe that stability requires the condition A 2 > 0, which is acceptable since this in turn implies that dark radiation to be positive from Eq. (46). Thus positivity of the dark radiation term ensures stability of this solution at zeroth order. Since we have higher order solutions in a perturbative form, the stability of the full solution is expected to be dominated by the zeroth order term. • The most important case in our hand is the situation where dark pressure and dark radiation satisfies the constraint relation 2U + P = 0. In this case we can determine stability exactly. For this solution it turns out that V T and V S are positive for all choices of F(R) − 4 . However though V V is positive for F(R) > 4 it becomes negative for the other choice. Hence All these modes are positive ensuring stability of the solution for the parameter space: F(R) > 4 . Otherwise, the solution is though stable under the tensor and scalar perturbations, is not so under vector perturbation.
• Another important aspect of this solution comes into picture when P 0 = 0. Then the solution represents a Schwarzschild (A)de-Sitter spacetime, which under proper limit leads to the Nariai spacetime. This has the peculiar property that a black hole in Nariai spacetime has increasing surface area due to quantum corrections as shown by Bousso and Hawking [78][79][80]. This phenomenon of antievaporation was then generalized for Nariai black holes in f (R) gravity [81], with f (R) gravity playing the role of anomaly induced effective action leading to anti evaporation. In our case as well with P 0 = 0, we have Nariai black hole as one limit and thus our solutions will also exhibit antievaporation. However for P 0 = 0, our solution cannot be reduced to the Nariai form and thus in general the solution presented here will not exhibit antievaporation. • Finally we consider the solution which corresponds to the other constraint relation with U + 2P = 0. In this case we have both black hole and naked singularity depending on 4 > F(R) or otherwise. In this case at large r limit both the solutions can be taken as 1 + Cr 2 . It turns out that, V T and V V are positive for all choices between F(R) and 4 , however V V ensures stability for the black hole solution not for the naked singularity. Thus the black hole solution is stable under all these perturbation, while the global naked singularity is stable only under tensor and scalar perturbation, but not under vector perturbation.
Thus we observe that the solutions present here are mostly stable under perturbations, except in some specific cases where the vector mode of the perturbation shows instability. Also we have pointed out that our solution reduces to the Nariai form and thus exhibits anti-evaporation in f (R) gravity, similar to previously obtained results.

Static spherically symmetric brane with conformal motion
We can use symmetries to explore the connection between geometry and matter through Einstein's equation. The most important of such symmetries can be realized through the use of conformal Killing vectors. The symmetry under which the space-time manifold admits conformal Killing vectors are known as, conformal motion. In this section we derive a particular metric which admits conformal motions. For the spherically symmetric and static solutions on the brane if one requires to have one-parameter group of conformal motion, the following condition results, L ξ h μν = ξ μ;ν + ξ ν;μ = φ(r )g μν (63) In the above relation ξ is the conformal Killing vector and φ(r ) is the conformal factor, while the above symmetry of the metric is known as conformal motion. The above relation should hold for all the individual metric components. In this relation h μν is the metric determining the vacuum space-time configuration, ξ μ is a vector field in this spacetime with respect to which the Lie variation has been taken and φ(r ) is an arbitrary function of the radial coordinate. Then following the procedure adopted in Reference [82] to determine interior structure of stellar objects, here also we can impose some symmetry requirement like, ξ μ u μ = 0. This symmetry enables one to determine all the unknowns exactly using the effective Einstein's equation. Thus using the metric ansatz given by Eq. (13), the above equation is shown to be equivalent to [82], where A, C and D are arbitrary constants. With the above results the Einstein equations (14), (15) and (16) reduce to, 1 From Eqs. (66) and (67) we obtain the dark radiation and dark pressure in terms of the unknown function φ as, Then from Eq. (65) and the expression for dark radiation, the differential equation satisfied by φ(r ) turns out to be, This can be solved with little effort to yield the general solution as, where, B is an integration constant. Thus full solution corresponding to this one parameter symmetry group of conformal motion leads to, There exists another important properties of the field equations. Having obtained a single solution we can make a transformation such that, r →r (r ), U →Ū (U ), P →P(P) and Q →Q(Q) [83], called homology transformations. The homology properties of the equations determining dark radiation and dark pressure can be simplified by assuming γ = P(U )/U = constant and c s = d P/dU = constant. The above transformations are being generated with the infinitesimal generator as,L = ζ(r )∂/∂r + ψ 1 (U )∂/∂U + ψ 2 (Q)∂/∂ Q. Then in order to have consistent solutions we must have, ζ = 0, ψ 1 = U and ψ 2 = Q + 2G M. Thus with inclusion of f (R) gravity the infinitesimal generator for the homologous transformation becomes restricted compared to that in Einstein gravity.

Some thermodynamic features
In this section we will discuss thermodynamics associated with these spherically symmetric vacuum spacetime. Our main motive is to observe if there exists any thermodynamic interpretation which is induced solely by the bulk. We focus on the line element obtained for the condition 2U + P = 0 which has the following expression, The horizon is determined by setting coefficient of g tt to zero, which in turn leads to the equation, Then the mass term equivalent to internal energy of a thermodynamic system can be obtained in terms of the horizon radius as, The surface area of the event horizon is given by, A = πr 2 h , while the entropy for the black hole is given by, S = k B A/4h = k B πr 2 h /4hr. Choosingh = 1 and Boltzmann constant appropriately we readily obtain, Thus the mass of the black hole in terms of the entropy becomes, This leads to the first law of black hole mechanics as, from which the black hole temperature turns out to be: while the chemical potential has the following expression: From the expression of temperature as a function of entropy it turns out that the specific heat has the following behavior: (84) Figure 1 shows that while temperature T and potential φ are continuous with both the entropy and the quantity F(R)− In order to understand the physics behind these results, it is always illuminating to discuss some limiting cases. For example, if we assume pure Einstein gravity, where F(R) = 0, then with the assumption of 4 ∼ 0, we arrive at the Reissner-Nordström solution. From Eq. (84) the specific heat turns out to be C v = −S(6α P 0 + 4S)/(2S + 9α P 0 ). This can also be divergent, provided the entropy satisfies the criteria: S = −(9α P 0 /2). In general P 0 is taken to be positive and thus the above relation cannot be satisfied in general. Hence the bulk term with positive dark pressure cannot lead to second order phase transition. However the other limit is interesting. For P 0 = 0, we get the divergence of specific heat to correspond to the condition: S = 1/(F(R) − 4 ). Thus for F(R) > 4 we have second order phase transition. Our calculations therefore confirm that Schwarzschild anti-de Sitter solution shows second order phase transition.
The case for which the dark radiation vanishes, i.e. U = 0 also exhibits the appearance of black hole horizon. The actual calculations are quiet complex, and we have presented them in Appendix A. However here we consider some limiting cases and discuss the corresponding thermodynamic features. The first case corresponds to F(R) − 4 = 0. In this situation the solution for the metric elements resemble Schwarzschild solution with no associated phase transition. We cannot take P = 0 as in that case the metric elements would diverge. Thus another obvious choice is M = 0.
Then also horizon appears and the specific heat diverges for 4 > F(R). Thus this configuration exhibits an opposite effect in respect to 2U + P = 0 case.
We therefore observe that in both the black hole solutions the specific heat diverge showing second order phase transition, due to the presence of F(R) gravity in the bulk. Thus bulk F(R) gravity plays a crucial role in determining the thermodynamic feature of the brane world black holes.

Discussion
In this work we have considered a bulk action with a f (R) term, where R is the bulk curvature. Starting from the bulk action we have derived the full effective Einstein's equation on the brane located at y = 0, which under f (R) → R limit goes to the usual Gauss Codazzi equation in Einstein gravity. In order to get spherically symmetric solutions we have assumed that in the region of interest there is no matter field present on the brane and also the four dimensional scalar curvature is constant. Under these conditions the Einstein equation simplifies considerably, however the Weyl tensor on bulk has non trivial decomposition on the brane leading to the appearance of dark pressure and dark radiation in the effective Einstein's equation. Also the induced four dimensional cosmological constant and contribution from f (R) term have significant effects on the solutions of the effective Einstein's equation on the brane.
Due to the presence of f (R) gravity in the bulk, Einstein's equation on the brane picks up an extra contribution which acts as an effective cosmological constant having expression: F(R) − 4 . Thus though the four dimensional parameter 4 is not small, an effective small cosmological constant can be generated by fine tuning 4 and F(R). Hence we can argue that the observed smallness of four dimensional cosmological constant is due to a fine tuning of induced cosmological constant on the brane with the f (R) term in the bulk.
From the effective Einstein's equation we can solve for the metric elements as well as for dark radiation and dark pressure term provided a relation between dark pressure and dark radiation term is assumed. For four such choices the equations get sufficiently simplified such that analytic solutions can be obtained. We have derived all the metric elements for these four choices. Among the four solutions two of them show the presence of event horizon and thus is important from thermodynamic point of view. On the other hand the other two solutions lead to naked singularity and thus does not have much astrophysical importance. The important features of these solutions are the asymptotic non-flatness due to presence of f (R) term. This might be of some relevance in the context of AdS-CFT correspondence.
After obtaining various solutions leading to either a black hole or a naked singularity, we have performed a stability analysis of our solutions in some appropriate limit. It turns out that the solutions are stable under tensor and scalar perturbations, while under certain choices of parameters the vector mode leads to instability. Also some solutions can be reduced to Nariai form, where the well known anti-evaporation in f (R) gravity takes place leading to an increase in the area of the event horizon. However we have argued that in general the solutions are stable under perturbations.
In order to get some idea about solutions representing stellar interior, a symmetry transformation, known as conformal motion is invoked. For this particular symmetry class we can solve the field equations exactly. This leads to direct evaluation of dark pressure and radiation using these symmetries. Also there exists another class of transformations known as homology transformations. For this class of solutions the homology operator has been evaluated and it turns out that f (R) term makes the homology class restricted compared to that in Einstein gravity.
Finally we consider thermodynamical behavior of these spherically symmetric space-times. Since thermodynamics is intimately connected to existence of a horizon, we consider only the two relevant cases. Here also the f (R) term plays a dominant role in determining the thermodynamic behavior. In both the cases, the temperature and chemical potentials are found to be continuous, while the specific heat turns out to be discontinuous along a surface indicating a second order phase transition. Such features of these spherically symmetric solutions have their origin in the f (R) term in the bulk action and only because of the presence of higher curvature terms in the action, the black hole solutions exhibit a phase transition, which, is second order in nature.
Then by the previous conditions:h = 1 and an appropriate choice of Boltzmann constant we get entropy to be S = r 2 h . From the first law of black hole mechanics as presented in Eq. (81) the temperature turns out to be, This again shows that the black hole solution presented by the condition of vanishing dark radiation has a divergent behavior on the above surface which in turn indicates that the black hole undergoes a second order phase transition on this surface.