Thin shells surrounding black holes in F(R) gravity

In this article, we consider spherical thin shells of matter surrounding black holes in F(R) theories of gravity. We study the stability of the static configurations under perturbations that conserve the symmetry. In particular, we analyze the case of charged shells outside the horizon of non--charged black holes. We obtain that stable static thin shells are possible if the values of the parameters of the model are properly selected.


Introduction
The observations concerning the accelerated expansion of the Universe, the rotation curves of galaxies, and the anisotropy of the microwave background radiation can be explained within General Relativity by the presence of dark matter (∼ 25 %) and dark energy (∼ 70 %), besides the ordinary barionic matter (∼ 5 %). In the concordance or ΛCDM model, the dark energy contribution comes in the form of a cosmological constant Λ and the cold dark matter in the form of non-relativistic fluid, supplemented by an inflationary scenario driven by a scalar field called the inflaton. Although successful, this model is not free of difficulties, such as the extremely small observed value of Λ compared to the expected one if thought as originated from a vacuum energy in particle physics, or the unclear nature of dark matter (although several candidates exist). Other approaches can be adopted, such as modified gravity theories, in order to try to avoid these problems and explain the observed features of the Universe without dark matter and dark energy. Quantum gravity also provides motivation for modified gravity. One well known theory is F (R) gravity [1][2][3] in which the Einstein-Hilbert Lagrangian is replaced by a function F (R) of the Ricci scalar R. In recent years, several solutions of the field equations in F (R) gravity have been found, including static and spherically symmetric black holes [1,[4][5][6][7][8], traversable wormholes [9], and branes [10].
The Darmois-Israel junction conditions [11] provide the tools for matching two solutions across a hypersurface in General Relativity. These conditions allow for the study of thin shells of matter, by relating the energy-momentum tensor at the joining hypersurface with the space-times at both sides of it. The formalism has been broadly adopted in many different scenarios due to its flexibility and simplicity; among them it is used to model vacuum bubbles and thin layers around black holes [12][13][14], gravastars [15,16], and wormholes [17][18][19][20]. In the case of highly symmetric configurations, the stability analysis is usually easy to perform, at least for perturbations preserving the symmetry.
The junction conditions in F (R) theories [21,22] are more restrictive than in General Relativity. For non-linear F (R), the trace of the second fundamental form should always be continuous at the matching hypersurface [22]. Except in the quadratic case, the curvature scalar R should also be continuous there [22]. In quadratic F (R) gravity, the hypersurface has in general, in addition to the standard energy-momentum tensor, an external energy flux vector, an external scalar pressure (or tension), and another energy-momentum contribution resembling classical dipole distributions. All these contributions have to be present [22,23] in order to have a divergence-free energy-momentum tensor, which guarantees local conservation. It was recently shown that these features are shared by any theory with a quadratic lagrangian [24]. The junction conditions in F (R) were recently applied to the construction of thin-shell wormholes [25,26] and bubbles [27,28]. A particularly interesting example of a pure double layer in quadratic F (R) was found [27].
In this work, we construct spherical thin shells surrounding non-charged black holes by using the junction conditions in F (R) gravity and we analyze the stability of the static configurations under perturbations that preserve the symmetry. In Sec. 2, we review the general formalism for spherical geometries with constant curvature scalars at both sides of the shell. In Sec. 3, we perform the construction and we study the stability of charged thin shells outside the black hole event horizon. Finally, in Sec. 4, we present the conclusions of the paper. We adopt a system of units in which c = G = 1, with c the speed of light and G the gravitational constant.

Spherical thin shells: construction and stability
We begin by reviewing the formalism for spherical thin shells in four dimensional F (R) gravity introduced in Ref. [28]. We consider a manifold composed of two regions with a constant curvature scalar in each, separated by a thin shell of matter. For this purpose, we take two different spherically symmetric solutions in F (R) gravity, with metrics where r 1,2 > 0 are the radial coordinates corresponding to each geometry, and 0 ≤ θ ≤ π and 0 ≤ ϕ < 2π are the angular coordinates. We proceed with the construction of a new manifold by selecting a radius a ≡ a 1 = a 2 and cutting two regions M 1 and M 2 defined as the inner 0 ≤ r 1 ≤ a and the outer r 2 ≥ a parts of the geometries 1 and 2, respectively. These regions are pasted to one another at the surface Σ with radius a. This construction results in the spacetime M = M 1 ∪ M 2 , with the inner zone corresponding to M 1 and the exterior one to M 2 . The angular coordinates have been naturally identified everywhere from the beginning. We define a new global radial coordinate r ∈ [0, +∞) by identifying r with r 1 in M 1 and with r 2 in M 2 , respectively. Then, the global coordinates are X α 1,2 = (t 1,2 , r, θ, ϕ), while on the surface Σ, corresponding to G(r) ≡ r − a = 0, we adopt the coordinates ξ i = (τ, θ, ϕ), with τ the proper time. In what follows, we take the radius of the surface Σ as a function a(τ ) of the proper time. The equality of the proper time at the sides of the shell requires that in which the free signs were fixed by choosing all times t 1,2 and τ to run to the future, andȧ is the proper time derivative of a. We denote the first fundamental form by h µν , the second fundamental form (or extrinsic curvature) by K µν , and the unit normals at the surface Σ by n 1,2 γ (pointing from M 1 to M 2 ). The first fundamental form associated with the two sides of the shell is given by and the second fundamental form has components with the unit normals (n γ n γ = 1) determined by We adopt at the surface Σ the orthonormal basis {eτ = e τ , eθ = a −1 e θ , eφ = (a sin θ) −1 e ϕ }. Then, for the geometry given by (1), the first fundamental form results h 1,2 ı = diag(−1, 1, 1), the unit normals read and the second fundamental form non-null components are and where the prime on A 1,2 (r) represents the derivative with respect to r. From now on, we denote with a prime on F (R) the derivative with respect to the curvature scalar R and the jump of any quantity Υ across Σ by [Υ] ≡ (Υ 2 − Υ 1 )| Σ . The junction formalism in F (R) gravity theories provides the conditions that should be fulfilled at Σ. One of them is the continuity of the first fundamental form i.e. [h µν ] = 0. It is straightforward to verify that this condition is satisfied by our construction. Another one is the continuity of the trace of the second fundamental form, i.e. [K µ µ ] = 0, which by using Eqs. (6) and (7), takes the form When F ′′′ (R) = 0 the continuity of R across the Σ is also required i.e.
[R] = 0. The field equations at Σ read [22] with κ = 8π and S µν the energy-momentum tensor at the shell. If F ′′′ (R) = 0 (quadratic case), the curvature scalar can be discontinuous at Σ, and the field equations adopt the form [22] κS supplemented by three other contributions: an external energy flux vector with∇ the intrinsic covariant derivative on Σ, an external scalar pressure or tension and a two-covariant symmetric tensor distribution with δ Σ the Dirac delta on Σ. This last expression admits an equivalent form for any test tensor field Ψ µν . In quadratic F (R), the shell has, in addition to the standard energymomentum tensor S µν , an external energy flux vector T µ , an external scalar tension/pressure T , and a double layer tensor distribution T µν of Dirac "delta prime" type, having a resemblance with classical dipole distributions. All these contributions are necessary in order to ensure the energymomentum tensor to be divergence-free, a condition that is required for conservation locally [22].

The same constant curvature scalar R 0
We firstly study the case with a constant curvature scalar R 0 at both sides of Σ. Then, the condition [R] = 0 is automatically fulfilled when required, and Eqs. (9) and (10) both simplify to give In quadratic F (R) the contributions T , T µ and T µν are all zero due to their proportionality to [R]. The energy-momentum tensor in the orthonormal basis takes the form Sî = diag(σ, pθ, pφ), with σ the surface energy density and pθ = pφ = p the transverse pressures, so from Eq. (15) we find that and In F (R) gravity, the inequality F ′ (R) > 0 implies that the effective Newton constant [5], so preventing, from a quantum point of view, the graviton to be a ghost. An interesting discussion about this issue, within a wormhole scenario, is presented in Ref. [29]. We assume the absence of ghosts, so we demand that F ′ (R 0 ) > 0 from now on. Normal matter satisfies the weak energy condition, which in the orthonormal frame requires that σ ≥ 0 and σ + p ≥ 0; if it does not, the matter is exotic. From Eqs. (8), (16), and (17) we obtain the equation of state σ − 2p = 0.
By taking the time derivative of the equation of state and using Eqs. (16) and (17), we obtain the conservation equation where A = 4πa 2 is the area of the shell. The first term in this equation represents the internal energy change while the second one is the work done by the internal forces at the shell. For static configurations with a constant radius a 0 , Eq. (8) reduces to The surface energy density σ 0 and the pressure p 0 take the form and respectively; they fulfill the equation of state σ 0 − 2p 0 = 0. We analyze the stability of static solutions under perturbations preserving the spherical symmetry. By using thatä = (1/2)d(ȧ 2 )/da and with the definition z = A 2 (a) +ȧ 2 − A 1 (a) +ȧ 2 , we can rewrite Eq. (8) to give az ′ (a) + 2z(a) = 0; so by solving this equation we find an expression forȧ 2 in terms of an effective potentialȧ where It is easy to see that V (a 0 ) = 0 and by using Eq. (20) that V ′ (a 0 ) = 0. The second derivative of the potential at a 0 reads A static configuration having a radius a 0 is stable if and only if V ′′ (a 0 ) > 0, corresponding to a minimum of the potential.

Different and constant curvature scalars R 1 and R 2
Now we take two different curvature scalars R 1 and R 2 at the sides of the shell Σ. The jump [R] = 0 restricts our study to the quadratic case F (R) = R − 2Λ + αR 2 case. Consequently, we only demand the continuity of the first fundamental form and of the trace of the second fundamental form, i.e. [h µν ] = 0 and [K µ µ ] = 0. We proceed as above, but now with constant R 1 = R 2 . The shell radius a has to satisfy Eq. (8). From Eq. (10), the field equations take the form then, using that Sî = diag(σ, p, p) in the orthonormal basis, we find that the energy density and the transverse pressure at the shell are respectively. As it was explained before, we assume that F ′ (R 1 ) = 1 + 2αR 1 > 0 and F ′ (R 2 ) = 1 + 2αR 2 > 0, in order to avoid the presence of ghosts. The matter is normal at Σ if it satisfies the weak energy condition. From Eq. (11) we obtain that T µ = 0 and from to Eq. (12) the external scalar tension/pressure T is given by or, using Eq. (8), by the expression From Eqs. (27), (28), and (30) we obtain the equation of state relating σ, p, and T By taking the time derivative of Eq. (31) and with the help of Eqs. (27) and (28), we readily find the generalized conservation equation with the area A defined above. In the left hand side of this equation, the first term is thought as the change in the total energy of the shell, the second one as the work done by the internal pressure, while the right hand side represents an external flux. The double layer distribution T µν , should satisfy Eq. (14), which in our case adopts the form for any test tensor field Ψ µν . The components in the orthonormal basis of the double layer distribution strength are which only depend on α and [R]; as a consequence, the dependence of Tî with the metric comes from the unit normal and the covariant derivative.
In the static configurations, the radius a 0 should fulfill Eq. (20), and from Eqs. (27), (28), and (30), the surface energy density σ 0 , the pressure p 0 , and the external tension/pressure T 0 read and Following the same procedure of the previous sub-section, the stability of the static configurations is found from Eq. (23) in terms of the potential of Eq. (24), with its second derivative given by Eq. (25); again V ′′ (a 0 ) > 0 correspond to the stable ones.

Charged thin shells
We start from the action corresponding to F (R) gravity coupled to Maxwell electrodynamics where g = det(g µν ) is the determinant of the metric tensor and F µν = ∂ µ A ν − ∂ ν A µ is the electromagnetic tensor. In the metric formalism, the field equations obtained from this action, considering an electromagnetic potential A µ = (V(r), 0, 0, 0) and a constant curvature scalar R, admit a spherically symmetric solution given by Eq. (1), in which the metric function [5,6] has the form with Q the charge and M the mass. The electromagnetic potential is V(r) = −Q/r and the cosmological constant is related with the curvature scalar by 4Λ = R.

Curvature scalar R 0 at both sides
For the construction of the thin shell Σ, we take the metric function given by Eq. (39), with mass M 1 = 0 and null charge for the internal region M 1 , and mass M 2 = 0 and charge Q for the external region M 2 . At both sides of Σ we adopt the same value R 0 for the curvature scalar. Therefore, the metric functions read for the internal zone, and for the external one. The possible horizons result from the zeros of the expressions A 1 (r) and A 2 (r). Both metrics are singular at r = 0. When A 1 (r) = 0 we get a polynomial of third degree if R 0 = 0, its real and positive roots correspond to the radii of the different horizons. For R 0 0 there is only the presence of an event horizon. When 0 < R 0 < 4/(9M 2 1 ) there is a cosmological horizon in addition to the event horizon. For the metric function A 2 (r), the horizons are determined by the solutions of the quadratic equation when R 0 = 0, while for R 0 = 0 they are given by the real and positive roots of a fourth degree polynomial. The critical charge value Q c has an important role in the study of the solutions, since it determines the number of horizons of the geometry. If R 0 < 0 and 0 < |Q| < Q c there are two horizons, the internal and the event ones. When |Q| = Q c , they fuse into one, and if |Q| > Q c only a naked singularity is left. For R 0 > 0 the metric has a cosmological horizon. Besides it, when 0 < |Q| < Q c there exist an internal and an event horizons, if |Q| = Q c both merged into one, and they finally disappear when |Q| > Q c , so that there is a naked singularity at the origin.
In order to start the construction of the shell we choose a radius a satisfying Eq. (8), larger than the event horizon radius in M 1 , so the black hole is always present, and when R 0 > 0, also smaller than the cosmological horizon radius coming from the original geometry for this region. On the other hand, this radius a should be large enough to avoid the presence of the event horizon and the singularity of the geometry used for the region M 2 . When R 0 > 0, it also has to be smaller than the cosmological horizon of this outer region. As we mentioned in the previous section, it is necessary that F ′ (R 0 ) > 0 to avoid ghosts. It is also preferable that the matter on the shell satisfies the weak energy condition, in order to guarantee the presence of normal matter on Σ. The energy density and the pressure are given by Eqs. (16) and (17), respectively, which fulfill the equation of state (18).
In the static case, the radius a 0 should satisfy Eq. (20), while σ 0 and p 0 are given by Eqs. (21) and (22). The potential (24) allows for the stability analysis of the solutions, which is determined from the sign of V ′′ (a 0 ), given by Eq. (25); the stable ones correspond to V ′′ (a 0 ) > 0. The results are presented graphically in Fig. 1, in which the most representative ones are shown. All quantities are adimensionalized with the mass M 2 of the outer region; the relations M 1 /M 2 = 0.5 and α/M 2 2 = 0.1 are adopted in all plots. The meshed zones represent those that satisfy the weak energy condition, and the gray areas have no physical meaning. Solid lines represent stable solutions, while unstable solutions are drawn with dotted lines. The behavior of the solutions does not vary significantly with the modulus of the value of the curvature scalar R 0 , but mainly with its sign, obtaining that • When R 0 < 0, for |Q| < Q c , there is only an unstable solution composed by normal matter.
For larger values of |Q|, there are two solutions made of exotic matter, one of them is stable while the other, close to the event horizon of the black hole, is unstable.
• For R 0 > 0 and |Q| < Q c there is an unstable solution constituted by normal matter. For values |Q| > Q c and a restricted range of charge, there are three solutions composed by exotic matter, one of them is stable. For values of |Q| much larger than Q c , there is only an unstable solution close to the event horizon of the black hole.
The function F (R 0 ), which is present through its derivative, does not produce significant changes in the qualitative behavior of the solutions, it only affects them by modifying their scale. The quotient |Q|/ F ′ (R 0 ) can be interpreted as an effective charge.

Curvature scalars
Analogously to the previous sub-section, we construct the shell by taking the mass M 1 = 0 and a null charge for the internal region M 1 , and the mass M 2 = 0 and the charge Q for the external one M 2 . But now, we adopt different curvature scalars at the sides of the shell Σ, so that [R] = 0. Then, the metric functions are and The possible horizons are found as explained in the previous sub-section, with the replacement of R 0 by R 1 or R 2 as appropriate, retaining the same characteristics described therein.
The radius a of the shell is properly chosen, in the same way as done in the previous sub-section, and it should satisfy Eq. (8). The surface energy density is obtained by Eq. (27), the pressure by Eq. (28) and the external tension/pressure by Eq. (30). These three equations determine, together with Eq. (8), the equation of state at the shell (31). We also have that T µ = 0 and, because [R] = 0, the non-null tensor T µν , with a dipolar density P µν given by Eq. (34). In the static case, the construction of the shell is done by choosing a radius a 0 that satisfies the Eq. (20). Besides, we assume that F ′ (R 1 ) > 0 and F ′ (R 2 ) > 0 in order to avoid the presence of ghosts. Again, we prefer matter fulfilling the weak energy condition. The surface energy density, pressure, and external tension/pressure are obtained from Eqs. (35), (36), and (37), respectively. As it is mentioned above, a non-null dipolar distribution given by Eq. (34) is also present. The stability of the solutions is determined by using the Eq. (25), through the study of the sign of V ′′ (a 0 ), which guarantees the stability of the solution when V ′′ (a 0 ) > 0. The results are displayed in Fig. 2  Again, the quotient |Q|/ F ′ (R 2 ) can be understood as an effective charge. In particular, with a suitable choice of parameters, i.e. R 1 = 0.39 and R 2 = 0.3, we have found that a stable solution made of normal matter and without charge is possible, as shown in Fig. 3.

Conclusions
In this work, we have studied spherically symmetric thin shells of matter around black holes within F (R) theories of gravity. We have adopted constant curvature scalars at both sides of the shell and we have analyzed two scenarios: one in which the curvature scalars are equal to the same value R 0 , and the other in which the values of the curvature scalars R 1 and R 2 are different. The case with both curvature scalars equal to R 0 does not impose any limitations on the function F (R). The matter at the shell should fulfill the equation of state σ − 2p = 0 that relates the surface energy density σ and the pressure p. When the curvature scalars are different, it is necessary to restrict the analysis to quadratic F (R). Then, for R 1 = R 2 the shell is composed by matter that satisfies the equation of state σ − 2p = T , which depends on the external tension/pressure T ; it also requires the presence of the vector T µ = 0 and the tensor T µν = 0 contributions. In particular, we have constructed a thin shell of matter surrounding a static non-charged black hole with mass M 1 . The geometry outside the shell corresponds to a solution with mass M 2 and charge Q.
In the case with the same value of the curvature scalar at both sides of the shell, the behavior of the solutions is determined by the sign of R 0 , and it is always possible to find unstable solutions with normal matter. On the other hand, stable solutions constituted by normal matter are not found; those that are stable are present, but always built with exotic matter, within a small range of |Q| when R 0 > 0 or for large values of |Q| when R 0 < 0. This result is similar to the one obtained in Ref. [28] for bubbles, with the main difference being that for the shells around black holes there is an extra solution, unstable and constituted by exotic matter, with values of the shell radius close to the event horizon one.
For different values of the curvature scalar at the regions separated by the shell, the solutions have a distinct behavior depending on the relation between R 1 and R 2 . Only in the case with R 1 > R 2 it is possible to find stable solutions constituted by normal matter for |Q| close to Q c . The rest of the solutions are unstable and made of exotic matter, with the exception of a small range of |Q| for which the shell is stable but, once again, it is composed by exotic matter. When R 1 < R 2 , we have only found unstable solutions constituted by normal matter for |Q| < Q c , or by exotic matter for |Q| > Q c . The only stable solution for this case is made of exotic matter for a broad range of values of the charge. These results are similar to those found in Ref. [28] in the case of R 1 = R 2 , with the difference that, in shells surrounding black holes, an additional unstable solution emerges close to the event horizon. Compared with the case of the Ref. [28], it is more complicated to find an appropriate set of parameters that allow the construction of a stable shell around the black hole, made of normal matter and without charge. However, we have found that it is possible to build such a case, for a limited range of values of R 1 and R 2 , and we have shown an example with R 1 = 0.39 and R 2 = 0.3.
It is well known that there is an equivalence between any F (R) gravity theory and a properly taken scalar-tensor theory [1,2]; in particular, quadratic F (R) is equivalent to Brans-Dicke theory with a parameter ω = 0, with a potential V (φ) = 2Λ + (φ 2 − 2φ − 3)/(4α), where the scalar field φ and the curvature scalar are related by φ = 2αR − 1. Then, it is worthy to note that the results obtained here can be translated to the corresponding scalar-tensor theory.