Journal of High Energy Physics, Gravitation and Cosmology
Vol.05 No.03(2019), Article ID:92328,25 pages

On the Thermodynamical Black Hole Stability in the Space-Time of a Global Monopole in f(R)-Gravity

Francisco Bento Lustosa1*, Maria Emilia Xavier Guimarães2, Cristine Nunes Ferreira3, Joaquim Lopes Neto4, José Abdalla Helayël-Neto5

1Departamento de Fsica Teórica, Instituto de Fsica, Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil

2Instituto de Fsica, Universidade Federal Fluminense, Campus da Praia Vermelha, Niterói, Brazil

3Núcleo de Estudos em Fsica, Instituto Federal de Educação Ciência e Tecnologia Fluminense, Rio de Janeiro, Brazil

4Instituto de Fsica, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

5Centro Brasileiro de Pesquisas Fsicas, Rio de Janeiro, Brazil

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: March 6, 2019; Accepted: May 7, 2019; Published: May 10, 2019


In this work, we re-assess a class of black hole solutions in a global monopole spacetime in the framework of an f ( R ) -gravity model. Our main line of investigation consists in considering a region close enough to the black hole, but such that the weak field approximation is still valid. The stability of the black hole is studied in terms of its thermodynamical properties, with the radial coordinate written as a power-law function with the status of the main factor underneath the stability of the model. We obtain the explicit expressions for the thermodynamical quantities of the black hole as functions of the event horizon, by considering both the Hawking and the local temperatures. The phase transitions that may occur in this system, including the Hawking-Page phase transition, are inspected with particular attention. We work out and contemplate a solution of special interest in which one of the parameters is related to the cosmological constant. Our main result sets out to establish a comparison between both the Hawking and the local formalisms for the black hole in the framework of the f ( R ) -gravity in the particular space-time adopted here.


Black Hole, f(R)-Gravity, Hawking Temperature, Global Monopole, Cosmological Constant

1. Introduction

There are still some important issues for which current cosmological standard theories in Physics have not yet found complete answers. General Relativity (GR) has not been able to clarify questions such as the non-renormalization of the gravity theory, the singularity problems in black hole physics and the physics of the early Universe, leading us to the necessity of finding alternative approaches. In this context, one of the most intriguing facts is the accelerating Universe. Despite of the existence of some alternative approaches to explain this behavior, one of those without adding dark matter or dark energy [1] [2] [3] , and the other with f ( R ) gravities [4] [5] that has received much attention as one of the strongest candidates to explain the current accelerating universe [6] . The f ( R ) gravity is constructed by replacing the Ricci scalar in the Einstein-Hilbert action by an arbitrary function of the Ricci scalar f ( R ) . Such theory is well known to lead to an extra scalar degree of freedom and it has been shown, for example in [7] [8] [9] , that observationally acceptable models may be built up. In some new contributions, f ( R ) -gravity can exhibit dependence on higher powers of the curvature scalar [10] .

In 1972, Bekenstein published his first article demonstrating the relation between thermodynamic quantities and gravitational properties of a black hole. His works [11] [12] [13] [14] were followed by the model of particle creation around BH proposed by Hawking [14] . Since then, the thermodynamical behavior of quasi-classical systems has been widely explored in other gravitational frameworks. Most recent contributions were made in the study of the thermodynamics of the black hole in modified Schwarszchild [15] , Born-Infeld-anti-de Sitter [16] , Horava-Lifshitz [17] , torsion [18] and f ( R ) theory [19] .

It has been shown that topological defects such as cosmic strings, monopoles and domain walls could be formed as a result of spontaneously broken symmetry in a vacuum phase transition of the early universe [20] [21] . In the context of the GR, Barriola and Vilenkin [22] studied the gravitational effects of a global monopole as a spherically symmetric topological defect. In the context of f ( R ) -theories, global monopoles heve been recently studied by [23] [24] [25] [26] , and references therein. In this solution, there appears a term which corresponds to a black hole solution. It is natural, then, to imagine that there is a region with a global monopole that could have been swallowed by a black-hole. Another interesting situation where the monopole appears as an important ingredient involves boson stars. When these compact objects are present along with a non-minimally coupled global monopole, a black-hole-type solution can be mimicked. This approach is important to unveil the non-linear gravitational effects and the gravitational back reaction [24] .

On the other hand, the analysis of the thermodynamical properties of that case, in the context of GR, was carried out in [27] and, more recently, in the context of an f ( R ) -theory in [28] , where the authors adopt the weak field approximation solution presented in [23] with a specific ansatz for

f ( R ) = F ( R ) = d f ( R ) d R

as a linear function of the radial coordinate.

In the present paper, we consider thermodynamical aspects of black holes in the space-time of a global monopole in the framework of f ( R ) -gravity. We study a general case in which f ( R ) is a power law function of the radial coordinates. We anticipate that we obtain the explicit expressions for the local thermodynamical quantities of the black hole as a function of the event horizon, the parameter describing the monopole and the measurable corrections on the usual GR gravity due to the f ( R ) -extension.

Our work is organized according to the following outline: In Section 2, motivated by the analysis of thermodynamic aspects, we revise the solution for a black hole in the region containing a global monopole with the use of a f ( R ) -model; the investigation is pursued by adopting the weak field approximation [26] . The new result in this section is the ψ n < 0 -case. Section 3 is split into subsections. In the first part, we analyze the entropy expression taking into account the corrections inherent to the case of an f ( R ) -theory. The stability of the black hole is also contemplated―in the Hawking’s approach―by taking into account the positivity of temperature and specific heat. In Part B, the black hole stability is re-assessed by following an alternative treatment. In part C, one studies the phase transitions of the system and a comparison is established between the Hawking’s and the local case approaches. Finally, in Section 4, we cast our Closing Remarks and present, more carefully, a discussion on the particular n = 2-case; in this special case, the ψ 2 -parameter is related with both situations of a positive and a negative cosmological constant that, in the presence of the global monopole, exhibits a solid angle deficit.

2. Field Equations Solution for the f(R) Gravity in the Metric Formalism: A Review

In this section, for the sake of understanding, we make a small review of the solution of the field equations for an f ( R ) -theory with a spherically symmetric space-time obtained in [26] . In this reference it is shown that it is possible to find a black hole solution in a global monopole region by using an f ( R ) modification of the GR gravity within the weak field approximation. For an f ( R ) theory, the action associated to the matter field coupled with gravity is given by:

S = 1 2 κ d 4 x g f ( R ) + S m ( g μ ν , ψ ) , (1)

where κ = 8 π G , G is the Newton constant, g is the determinant of the metric g μ ν with μ , ν = 0 , 1 , 2 , 3 , R is the curvature scalar, S m is the action associated with the matter fields and f ( R ) is an analytic function of the Ricci scalar. In this model, the Ricci scalar in Einstein-Hilbert action is replaced by f ( R ) . The monopole which introduces an angular deficit in the space-time metric, gives us some interesting effects that we shall discuss in the sequel. Here, we assume that the Christoffel symbol is a function of the metric, its derivatives and its inverse. The general form of the time-independent metric with spherical symmetry in (1 + 3) dimensions is given by

d s 2 = B ( r ) d t 2 A ( r ) d r 2 r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) (2)

The g μ ν field equations read as below:

R μ ν F ( R ) 1 2 f ( R ) g μ ν ( μ ν g μ ν ) F ( R ) = κ T μ ν (3)

with F ( R ) = d f ( R ) d R , and is the usual notation for covariant D’Alembert

operator μ μ . The only term associated to the matter action is related to the global monopole described by the Lagrangian density

L = 1 2 ( μ ϕ a ) ( μ ϕ a ) 1 4 λ ( ϕ a ϕ a η 2 ) 2 , (4)

where λ and η are the monopole field parameters and the triplet field that will result in a monopole configuration that can be described by ϕ a = η h ( r ) r 2 x a , with a = 1 , 2 , 3 and x a x a = r 2 . The function h ( r ) is dimentionless and is constrained by the conditions h ( 0 ) = 0 and h ( r > η ) 1 [22] . The energy-momentum tensor associated with that field configuration is

T t t = T r r η 2 r 2 T θ θ = T ϕ ϕ = 0 (5)

We can now focus on the gravitational effects of the global monopole in the f ( R ) approach to gravity, by taking the trace of the field Equation (3) in the presence of the monopole [23] [25]

f ( R ) = 1 2 R F ( R ) + 3 2 F ( R ) κ 2 T (6)

where F ( R ) = d f ( R ) d R . By replacing (6) in (3), we have

F ( R ) R μ ν μ ν F ( R ) T μ ν = g μ ν C (7)

where the quantity C = 1 4 ( F ( R ) R F ( R ) κ T ) is a scalar quantity.

In a spherically symmetric space-time, with the metric given by (2) and with the energy-momentum tensor expressed by (5), we find, by virtue of (7), that the equations may be written as

2 r F β ( r F + 2 F ) = 0 , (8)

4 B ( A 1 ) + r ( 2 F F β ) ( r B 2 B ) + 2 r 2 B 4 α A B F = 0 (9)

where α and β are:

α = 2 κ η 2 , with κ = 8 π G (10)

β = ( A B ) A B . (11)

We make the assumption of time-independent solutions, i.e. B = B ( r ) and A = A ( r ) , which yields the metric as given in (2). With these assumptions, we consider F ( R ) = F ( r ) and F and F are the first and the second derivatives with respect to r, respectively. The solution to these equations is the exact description of the global monopole in f ( R ) theories. This prescription has however an analytical solution only in the weak field approximation. In this approximation, the field equations read as follows:

F ( R ) = F ( r ) = 1 + ψ (R)

B ( r ) = 1 + b ( r ) , A ( r ) = 1 + a ( r ) (12)

with | b ( r ) | , | a ( r ) | and | ψ ( r ) | much smaller than one. These redefinitions are used in Equation (8) and Equation (9), and with the help of

F F ψ , F F ψ (13)

A A a , B B b (14)

the field Equation (8) and Equation (9) can then be written in a linear form as:

β r ψ ( r ) , (15)

2 a ( r ) 2 r ψ ( r ) + r β + r 2 b ( r ) 2 κ η 2 = 0, (16)

where β has been defined in (11). This is the solution to our equations in the weak field approximation for a metric with spherical symmetry in an f (R)

-theory with d f ( R ) d R = 1 + ψ ( r ) where R is the function of the radial coordinate, r. Following in [26] , we assume that the function ψ ( r ) is a power law-like function of the radial coordinate, namely ψ ( r ) = ψ ( n ) r n , where ψ n is a constant parameter in r. Notice, however, that in our work it shows a dependence on n. With this ansatz, we can find solutions for the equations above. Noting that Equation (11) with Equation (14) can be rewritten as:

a ( r ) + b ( r ) = n ( n 1 ) ψ n r ( n 1 ) (17)

we can solve (16) for b ( r ) , obtaining:

b ( r ) = c 1 r ψ n r n κ η 2 . (18)

We define the integration constants c 1 = 2 G M and, for convenience, we take c 2 = 0 . This term corresponds to the cosmological constant term. Equation (11) and Equation (15) yield the relation:

A ( r ) B ( r ) = a 0 e ( n 1 ) ψ n r n . (19)

Defining the integration constant a 0 = 1 , we then have the full form of the metric as

d s 2 = B ( r ) d t 2 e ( n 1 ) ψ n r n B ( r ) 1 d r 2 r 2 d Ω 2 , (20)

B ( r ) = 1 κ η 2 2 G M r ψ n r n (21)

In Figure 1, we plot the function B ( r , n ) for n varying from 1 to 4. It can be seen, in the case of the figure on the left, with ψ n > 0 , that the function B ( r , n ) grows to a value and then decreases for all orders of n. In the case ψ n < 0 , the function transits to ever increasing values.

Considering the mass term as cosmologically relevant is equivalent to adopting the scenario with a black hole in the space-time of a global monopole [22] . To confirm that our solution is actually a black-hole-type solution for all degrees n of the ψ ( r ) -function, we have to find if all of them have an event horizon. To do this, we have to keep in mind the approximations that we have used to solve the field equations, namely h 1 , which we have employed to define the energy-momentum tensor outside the core of the global monopole, and | ψ ( r ) | 1 . These approximations require that we restrict our analysis to a region defined by:

δ < r < 1 | ψ n r n | (22)

where δ ( λ η 1 / 2 ) is the order of magnitude of the monopole’s core.

In the following sections, we shall concentrate on the cases where ψ n < 0 and ψ n > 0 in the limit given by (22). These solutions do not contemplate the asymptotic case when n > 0 ; however, its analysis is important if we wish to consider regions closer to the black hole.

3. The Thermodynamics of the Black Hole in an f(R) Global Monopole for Generic n

In this section, the thermodynamic behavior is reported for two cases. In one of the cases, with ψ < 0 , the stability of the BH is ensured for all powers; the other

Figure 1. Graphs for B(r, n) in function of position for each n-orders going from 1 to 4, where G M = 1 and κ η 2 = 0.3 × 10 5 . In the figure (a) ψ n = 0.2 × 10 2 n and (b) ψ n = 0.2 × 10 2 n .

case, where ψ > 0 , the technique of local thermodynamics has been adopted to analyze the stability of the BH.

3.1. The Hawking Thermodynamics and the Power Law Series for a Generic n

In this sub-section, we shall be considering the thermodynamical behavior of a black hole in a region with a topological defect, by following the Hawking procedure. We can analytically with prove that the metric has an event horizon for the solution of an ( n + 1 ) -degree equation:

1 κ η 2 2 G M r H ψ n r H n = 0. (23)

We have written the parameter ψ n as a power law-like function of radial coordinate for a generic n. This is essential to respect the condition (22) without narrowing the region that we are going to analyze. Considering only the region where | ψ ( r ) | 1 implies that the other roots of Equation (23) are out of the region where our solution is valid. In addition to the limit for ψ n and κ η 2 , used in the previous graphs, the black hole (BH) stability was analyzed in detail considering the functions of temperature, heat capacity and its phase.

From (23), we get the dependence of the energy (GM) on the horizon r H

G M = 1 2 [ ( 1 κ η 2 ) r H ψ n r H n + 1 ] (24)

The Hawking temperature can be readily obtained from the GM expression, which corresponds to the energy of the black hole. The Hawking temperature is, therefore, the derivative of the energy with respect to the entropy of the black hole [29]

S = F ( R ) A / 4 1 4 ( 1 + ψ n r H n ) π r H 2 , (25)

where A is the area of the event horizon.

In Figure 2, which shows graphically the behavior of the entropy with the radius of the horizon, it indicates that there are critical points for different values of n given by the expression to r H :

r H max = 2 1 / n ( ψ n ( n + 2 ) ) 1 / n (26)

In this expression it can be easily verified that r H max applies only when ψ < 0 , as can be seen by the graphs of Figure 2. This fact has given us an important result which was related to the critical range between allowed and forbidden region, but throughout this work will be analyzed if this transition is within the limit of validity of the solution found.

In Figure 3, we plot the energy in terms of the entropy of the black hole. It is interesting to investigate how the analogous relations for gravity and thermodynamics become; with the development of a quasi-classical theory [12] [30] [31] , in a widely accepted relation in which the area of the black hole corresponds

Figure 2. Graph for S ( r H ) for n going from 2 to 4 and Schwarzschild black hole with global monopole, κ η 2 = 0.3 × 10 5 where in left panel ψ n = 0.4 × 10 2 n and in the right panel ψ n = 0.4 × 10 2 n .

Figure 3. Graphs for G M H W ( S , n ) for n going from 2 to 4 and the Schwarzschild black hole case, κ η 2 = 0.3 × 10 5 and ψ n = ± 0.4 × 10 2 n .

in fact to its entropy and consequently, its superficial gravity is responsible for measuring the temperature that the radiation emitted by the black hole will have. Although we are far from testing these results experimentally, the apparent lack of inconsistencies in the theoretical model and its application to several cases of interest [4] make it widely accepted.

The prescription for the Hawking’s temperature is given by the second low T = d G M d S given by

T H W = d G M d r H ( d S d r H ) 1 (27)


d S d r H = π r H 2 [ 1 + ( 1 + n 2 ) ψ n r H n ] (28)

d G M H W d r H = 1 2 ( 1 κ η 2 ( n + 1 ) ψ n r H n ) (29)

T H W ( η , ψ n ) = 1 π r H [ 1 κ η 2 ( n + 1 ) ψ n r H n ] [ 1 + ( 1 + n 2 ) ψ n r H n ] , (30)

which, for n = 1 , recovers the result obtained in [28] , and, for η = 0 , ψ ( r ) = 0 , the Schwarzschild result is reproduced. In this work, we concentrate in the case n > 1 with a Schwarzschild-like solution in presence of a global monopole.

Figure 4 shows that, for all values of n, there are two types of the black hole temperature: one with the temperature decreasing with the horizon, r H , and the other with the temperature growing with the horizon. The minimum point between these two regimes is calculated in this model by the derivative of the Hawking temperature as a function of the horizon radius d T H W r H = 0 . The expression resulting from this calculation is a simple polynomial of degree n, for n 1 . In this graph it can also be seen that the temperature behavior in this region has the same shape and transition. This fact has not been applied to the whole region of space.

Despite the positivity of the temperature, to analyze the stability of the black hole, it is necessary to verify the prescription for the heat capacity. Using the

Figure 4. Graph for the Hawkings temperature, T in function of the horizn radio r H , with κ η 2 = 0.3 × 10 5 and n going from 2 to 4, with Schwarzschild black hole with global monopole ψ n = 0.4 × 10 2 n .

laws of thermodynamics for black holes, it has been found that it is possible to obtain the heat capacity of the black hole by considering the entropy of the system, i.e. d E = T d S . In the case where the volume is constant, the energy, E, can be related to the mass of the black hole.

We shall now look for the heat capacity, which can be calculated from the energy. The expression is obtained as follows:

C H W = ( d G M d T H ) r = [ d G M H d r H ( d T H d r H ) 1 ] r (31)

u = d G M H d r H = 1 2 ( 1 κ η 2 ( n + 1 ) ψ r H n )

v = d T H d r H = 1 4 π [ 1 + ( 1 + n 2 ) ψ n r H n ] 2 r H 2 × { 1 κ η 2 + ( n 2 1 ) [ 1 + ( 1 + n 2 ) ( n 1 ) ( 1 κ η 2 ψ n r H n ) ] ψ n r H n }

C H = u v (32)

By considering Equation (30) and Equation (24), we obtain the expression:

C HW =2π r H 2 [ 1κ η 2 ( n+1 ) ψ n r H n ] [ 1+( 1+ n 2 ) ψ n r H 2 ] 2 1κ η 2 +( n 2 1 )[ 1+ ( 1+ n 2 ) ( n1 ) ( 1κ η 2 ψ n r H n ) ] ψ n r H n (33)

In Figure 5, we see the heat capacity plotted for the classical case and the n = 2 to 4 and compared with the Schwarzschild black hole.

By inspecting Figure 4 and Figure 5, one can discuss the validity regions of our solution of the BH in terms of the temperature and heat capacity. It can be seen that, for the region δ r H < r H * , the temperature is positive but the specific heat is negative. This region is therefore prohibited because it has a negative heat capacity because it is not in accordance with the laws of thermodynamics. There is a discontinuity in the heat capacity equation at the point r H = r H * , which can be seen in the graph of Figure 5. This discontinuity is the point where the specific heat changes behavior, becoming positive. In this region both temperature and heat capacity are positive, with the possibility of obtaining an allowable solution for the black hole.

Due to the maximum point given by the Equation (26) and which can be seen graphically in Figure 2, we have another phase transition in this system. For r H > r max the black hole becomes unstable again, although the heat capacity is positive. The heat capacity at this point goes from decreasing to the radius of the

Figure 5. Graph for heat capacity, C, with κ η 2 = 0.3 × 10 5 with ψ n = 0.4 × 10 2 n and panel right.

horizon, to increasing with the radius of the horizon, and this destabilizes the BH presenting negative temperatures, in this way this region becomes forbidden as we can saw in Figure 6. In the next section we will discuss the local case, where we will compare the local temperature with the Hawking temperature, this comparison gives us new results.

In the graph of Figure 6 it is easy to see that the transition is located exactly at the maximum point of the entropy curve given by the expression (26).

It may be verified that, whenever ψ n > 0 , both regions are prohibited, because, for a positive temperature, the heat capacity is negative, and when the heat capacity is positive, the temperature is negative. Then, in the case of ψ n > 0 it is necessary to work with another prescription that, in this paper, we take as the local prescription.

3.2. The Local Thermodynamics of the Black Hole for a Generic n

In the previous section, we have treated the thermodynamical properties by considering the Hawking formalism, where it is possible to investigate the properties in the case ψ n < 0 . In the present section, we go through a formulation where it is possible to discuss the stability conditions independently of the signal ψ n . We have nevertheless to recall the limitations imposed in the metric analysis (20). The solution that led us to this configuration was found by using weak field approximations, which consists in disregarding any cross-terms. In addition, the metric (20) is only valid locally, for a region δ < r < 1 / | ψ n | , also we analyzed in Hawking’s case, which forces us to work with local thermodynamical quantities. The Schwarzschild metric does not present such limitations, but for the purposes of comparison, in this section the local thermodynamic quantities are going to

Figure 6. Graph for temperature and heat specific, C in function of the horizon ratio for n = 2, with κ η 2 = 0.3 × 10 5 and ψ n = 0.4 × 10 2 n .

be derived. We use the prescription given in [32] to obtain the local temperature of the black body radiation:

T l o c = T g 00 (34)

with this relation the temperature was obtained as a function of the position as

T l o c = 1 π [ 1 κ η 2 r H ( n + 1 ) ψ n r H n 1 ] [ 1 + ( 1 + n 2 ) ψ n r H n ] × r ψ n r H n + 1 r H ( 1 8 π G η 2 ) + r ( 1 8 π G η 2 ) ψ n r n + 1 (35)

The local temperature for n = 2 to 4 and the Schwarzschild black hole is plotted in Figure 7 as a function of the event horizon where we fixed the position at r = 10 . The graph is zoomed to show the slight variation difference between the minimum values of the temperature for different values of the degree n.

An important characteristic to observe at this temperature is the existence of a minimum in its variation with respect to the event horizon. This minimum can be called the critical temperature [33] , which would be the minimum temperature in the region around the black hole. In Figure 7 the local temperature curves for n = 2, 3, 4 and the Schwarzschild were calculated numerically. Displays that, for all values of n, there are two possible black holes: one for which the temperature diminishes with the horizon, to r H > r H * , and the other with the temperature growing with the horizon, to r H > r min . These phase transitions were analyzed in more detail in the local case in connection with the heat capacity formalism.

Figure 7. Graph for T l o c ( r H ) for n = 2 , 3 , 4 and Schwarzschild, for r = 10 .

It is important to notice, from this and the following graphs, that we fix a position for the measurement of the temperature and we are analyzing how the temperature in that position varies as the event horizon increases. That means, we are finding what size the black hole has at the moment in which the temperature in the position r = 10 is minimal. In this sense, we see that as the horizon approaches the chosen position the temperature decreases rapidly to the minimum and then increases asymptotically when it approaches r = 10.

From the first law of thermodynamics d G M l o c = T l o c d S , the thermodynamical local energy, G M l o c , can be derived. Following the prescription discussed in last section, we can assume that S = F ( R ) A / 4 1 4 ( 1 + ψ n r H n ) π r H 2 , where A is the area of the event horizon.

G M loc = E 0 + S 0 S T loc dS= E 0 + π 2 r H ( M=0 ) r H ( M ) T loc ( ξ )[ 1+( 1+ n 2 ) ψ n ξ n ]ξdξ (36)

where we choose E 0 = 0 for simplicity. Using the expression for the local temperature in the integral using

r H ( M = 0 ) = [ ( 1 κ η 2 ) / ψ n ] 1 / n (37)

with the integration constants conveniently, we get

G M l o c = r ( 1 8 π G η 2 ψ n r n ) r r ( 1 8 π G η 2 ψ n r n ) r H ( 1 8 π G η 2 ψ n r H n ) , (38)

In Figure 8, we plot the Energy as a function of the entropy. It can be seen that, in the local framework, the energy and the temperature T are both positive.

We are now going to look for the local heat capacity, which can be calculated from the energy. The expression is obtained as follows:

Figure 8. Graph for G M ( S , n ) l o c for n going from 2 to 4 and Schwarzschild r = 15, κ η 2 = 10 5 and ψ n = 0.4 × 10 2 n .

C l o c = ( d G M l o c d T l o c ) r = [ d G M l o c d r H ( d T l o c d r H ) 1 ] r (39)

The subscript is there to show that we are calculating this quantity in a fixed position. We work with Equation (39) and give us

d G M l o c d r H = r 2 ( 1 8 π G η 2 ( n + 1 ) ψ n r H n ) × [ r ( 1 8 π G η 2 ψ n r n ) r H ( 1 8 π G η 2 ψ n r H n ) ] 1 / 2 , (40)

d T loc d r H =v r × [ r( 18πG η 2 ψ n r n ) r H ( 18πG η 2 ψ n r H n ) ] 1/2 + T H 2 r ( 18πG η 2 ( n+1 ) ψ n r H n+2 ) × [ r( 18πG η 2 ψ n r n ) r H ( 18πG η 2 ψ n r H n ) ] 3/2 (41)

In Figure 9, we see the local heat capacity plotted for the classical case and the n = 2 to 4 and the local form for the Schwarzschild black hole. It can be observed that for the local case, it is possible to obtain an analysis for the thermodynamics of black holes for both ψ < 0 and ψ > 0 to r = 10. It can be seen that the curves for heat capacity are similar for the case n > 2 . In fact in both cases the r limit of truncates the function.

The transition from negative to positive values happens as the horizon is closer to the chosen position, r also vanishes present only the positive part. We can notice this feature more clearly by analyzing the heat capacity as a function of the position, with a fixed horizon. The expected negative heat capacity will be

Figure 9. Graph for C l o c ( r H , n ) for the n = 2 , n = 3 and n = 4 and the Schwarzschild for black hole cases, for fixed r = 10 , κ η 2 = 10 5 and ψ n = 0.4 × 10 2 n .

observed as r grows and the transition to positive values happens even for bigger black hole. This transition is yet to be explained and it occurs even for the classical Schwarzschild case in the local analysis, but not in the most usual global calculation where it remains negative for every value of r H . It can be observed that for the local case, it is possible to obtain case for ψ>0 , that to n > 2 has been similar to case ψ < 0 , in the case of small r.

With Figure 10 it can be seen more easily that the function of the specific heat is stopped precisely at the point r l i m i t . Now we can analyze the case n = 2 , where another phase transition appears where the specific heat changes direction and direction in case ψ > 0 . It is found that when r is large this would happen for all n. However, when r is very large, the specific heat becomes totally negative and BH will not even form.

Turning to the transition of the entropy phase, due to the correction introduced by the f ( R ) theory discussed in [29] , we can have two cases the first for ψ < 0 , which is analyzed in a shown in Figure 11 on the left to n = 2. The Hawking temperature introduces a maximum limit for the horizon radius r H , which divides a permitted region with, T H W > 0 , of a forbidden region T H W < 0 . This maximum limit in the radius of the horizon r H can be calculated using the Equation (26) and corresponds to the solid line of the graph.

An interesting analysis is that the limit of r max that appears in the temperature of Hawking is also present in the local case. In the local case, this boundary also divides a permitted region of a forbidden region, and this happens when the position r = r H max . When r > r H max , the minimum temperature decreases, but the allowed region always has the same size as r H max .

This behavior can be seen in the graph of Figure 11 on the right, showing the value of r H max for n = 2, considering different values of r. It can be shown that for r < r H max there is no forbidden region.

Figure 10. Local heat capacity for n = 2 as a function of the position, κ η 2 = 10 5 and ψ n = 0.4 × 10 2 n .

Figure 11. Local temperature in the vicinity of r H max for n = 2 in panel left with ψ n = 0.4 × 10 2 n compared with Hawking temperature where r H max = 111.8033989 and in panel right the heat capacity.

The existence of a r H limit for black holes that has been studied in the case ψ < 0 repeats to ψ > 0 existing a value of r H that can not be passed r H max . The Hawking temperature for ψ > 0 has no minimum point becoming negative, the r H max can be found when that temperature is zero by doing the zero expression (30) that gives us:

r H max = ( 1 κ η 2 ψ n ( n + 1 ) ) 1 / 2 (42)

In the local case we have a behavior similar to ψ < 0 the maximum value for r is r H max , however in this case there will always be a transition from the allowed region to the forbidden region, even when r < r H max in the vicinity of r max . If we analyze the heat capacity we can observe that the validity limit to r is very small and far from r H max because the heat capacity became negative and the black hole is not stable.

3.3. The Thermodynamical Phase Transitions for the Black Hole Systems

We shall take care in this section of the possible phase transitions of the system we are considering. The first one is the Hawking-Page phase transition. It occurs whenever the function that gives the behavior of the Helmholtz free energy with the radius of the horizon has a root which is a minimum of this function.

F H P H W | r H = r H P = G M T H P S = 1 2 [ ( 1 κ η 2 ) ψ n r H P n ] r H P 1 4 ( 1 + ψ n r H P n ) π T H P r H P 2 = 0

( F H P H W r H ) r H = r H P = 1 2 [ 1 κ η 2 ( n + 1 ) ψ n r H P n ] π 2 T H P [ 1 + ( n 2 + 1 ) ψ n r H P n ] r H P = 0 (43)

where T H P is the Hawking-Page temperature. We find the radius and the Hawking-Page temperature as given below:

T H P H W = 1 π r H P [ 1 κ η 2 ( n + 1 ) ψ n r H P n ] [ 1 + ( 1 + n 2 ) ψ n r H P n ] (44)

r H P = [ n 1 2 α ( n + 1 ) ± 1 2 Δ 1 / 2 ψ n ] 1 / n (45)

Δ = α 2 ( n + 1 ) 2 4 α ( n 2 + n + 1 ) + 4 ( n 2 + 1 ) (46)

In Figure 12, the free energy question can be analyzed by considering the Hawking-Page phase transition showing its phase change for each value of n using Hawking’s formalism. In the local case we can have this transition for both ψ < 0 and ψ > 0 . Repeating what happens in the case of the Hawking-Page radius, these values are very close, so we chose n = 2 to do the analysis as shown in Figure 13.

In Figure 13, we present the possible regimes for the free energy in the case n = 2 . The same analysis applies for the others orders of n. The Hawking case that corresponds to ψ 2 < 0 was plotted in the panel of the left. In this case, we have a negative cosmological constant. In the right plot, we have ψ 2 > 0 that have similar results to ψ 2 < 0 in the local case, which can represents when ψ > 0 the case with a positive cosmological constant.

The region of temperatures higher than the Hawking-page temperature of the Black Hole will always be stable. The minimum temperature puts into evidence a black hole’s unstable region. In both cases, for T < T min , where T min is represented

Figure 12. Graph of F H W for the n-order going from 2 to 5, where G M = 1 , κ η 2 = 0.3 × 10 5 and ψ n = 0.4 × 10 2 n .

Figure 13. Graph for F for each n = 2, κ η 2 = 0.3 × 10 5 . The solid line is the Hawking-Page free energy, the dashed line is the minimal energy the dashed dotted line is the energy when T > T H P and the dotted line is the energy when T < T min and in (a) the Hawking framework ψ n = 0.4 × 10 2 n ; (b) The local framework ψ n = 0.4 × 10 2 n .

by the dashed curve, we have a region of pure radiation. The solid curve represents the Hawking-Page temperature that gives us the Hawking-Page Phase transition. The expression, for the minimum of the temperature in Hawking case, can be readily worked out and gives us

T min = 1 ( 1 κ η 2 r min ( n + 1 ) ψ 2 r min ( n 1 ) ) (47)

r min = ( ( 1 κ η 2 ) ( n 2 1 ) ψ 2 ) 1 / 2 (48)

As previously discussed, this work considers the boundary of the weak field approximation in a region in the vicinity of the BH with the defect. For this reason, we do not contemplate the asymptotic limits. However, whenever n = 2 and ψ is negative, the behavior of BH in the f(R)-model considering the Hawking framework resembles the one which occurs in the case of the negative cosmological factor corresponding to the Anti-de-Sitter case. This case is extremely important for the study of the quantum behavior of gravitational systems.

So far, the case analyzed considers ψ n < 0 , which can be treated by the Hawking’s formalism. In the case of the local formalism, it is possible to verify that independently of ψ n , it was possible to study both behaviors.

In the graph of Figure 14, the study of the compatibility between the two formalisms can be found. It can be verified that both present the same critical behavior if we choose the value in the local quantities of the fixed point where local quantities are calculated.

In Figure 14, it can be seen that both the case of ψ n < 0 to local and Hawking cases, can present the same set of minima and Hawking-Page temperature, showing compatibility between the two theories. It is always possible to find the same set of values if we regulate the values of r in local formalism. In this Figure, this happens when, for a given Hawking-Page temperature, in local formalism, we set r so that the Hawking-Page temperature value is consistent with Hawking’s formalism.

In Figure 15, a comparison was made between the heat capacity calculated in

Figure 14. Graph for T for each n = 2, κ η 2 = 0.3 × 10 5 and ψ n = 0.4 × 10 2 n with r = 105.57.

Figure 15. Graph for the comparation between C l o c and C H for each n = 2, κ η 2 = 0.3 × 10 5 and ψ n = 0.4 × 10 2 n with T H P = 1.3408 × 10 2 for both cases.

Hawking’s formalism ( C H W ) and Local ( C l o c ) in the case of psi < 0. In the figure to the left it can be seen that in the case of Hawking formalism, as indicated in the graph, there are three distinct regions, the first comprises a negative C H W for r H < r * . The other region corresponds to a positive C H W for r * < r H < r c and the latter comprises the region given by r H > r c . In the figure to the right, upper part, this behavior was further detailed. Analyzing the graph to the right, it can be clearly seen that r * coincides with the forbidden/allowed transition region. The region r H < r * is a forbidden layer because it has specific heat negative. The region r * < r H < r c is permissive where we find a stable black hole. In this region the temperature starts from the minimum and grows with the radius of the horizon while the heat capacity decreases. In this region the black hole is stable. When r H > r c we have again a forbidden region.

The other curve of the specific heat of Figure 15 left corresponds to the local case. For the comparison the parameters of the local model were adjusted to r = 105.57 . With these conditions, the Hawking-Page temperature can be obtained. In Figure 15 the right panel down can analyze the detail of the local case. It can be seen that in this case the solution to the black hole presents two possibilities corresponding to both heat capacity and temperature greater than zero. It can be seen that there is a region between r * < r H < r min , where a black hole solution is possible, but this is unstable because the temperature decreases with the radius of the horizon going to a more stable location that is r min . This region is called a small black hole. The other region presents, as in the case of the Hawking formalism, the temperature increasing with the increase of r h , showing that for the increase of the black hole this is giving more energy, for that reason warmer. This region is called the Great Black Hole.

In Figure 16, for n = 2, the comparison of the black hole in the local formalism

Figure 16. Graph for the comparation between C l o c for ψ n < 0 and C l o c for ψ n > 0 each n = 2, κ η 2 = 0.3 × 10 5 and ψ n = 0.4 × 10 2 n r = 10 and T H P = 0.1078 .

with ψ < 0 and ψ > 0 for the same Hawking-Page temperature was analyzed. As in the case of Figure 15, where the ψ < 0 case was analyzed, it can be seen that the limits for large r H depend on the value of r. Comparing these two cases, it can be analyzed that when ψ > 0 , has only one allowed region for r * < r H < r c , where we have both temperature and calorific capacity positive. In this region there may be several possibilities, the temperature decreases as the heat capacity decreases, for increasing r H , small black hole. A region where temperature increases and specific heat decreases and a region where both specific heat and temperature increase with increasing horizon radius r H this region is the big black hole.

4. Closing Remarks

In this work, the solution of a black hole in a region containing a global monopole in a class of f(R) gravity models has been investigated using a power law approach. To pursue our investigation, we have considered the weak field approximation for small radial coordinate values, assuming the validity of the condition of the Equation (22). In this his region, it is possible to analyze the behavior of BH solutions for n > 0 and ψ n positive and negative.

The very goal of this work is the analysis of the thermodynamical behavior of the black hole, considering two possible formalisms: Hawking’s, on the one hand, and local formulation on the other hand.

In the case of Hawking’s approach, it was possible to study the case where ψ n < 0 as being stable, by studying the positivity of the temperature with the positivity of the heat capacity. This black hole was called a large black hole, located for the radius of the horizon, r H > r min . For r H < r min , we have an unstable black hole despite having a region between r H * < r H < r min , where we the heat capacity is positive and the temperature is positive.

The case with ψ n > 0 has been tackled in the local formalism; this framework is actually valid for both ψ n positive and negative. The local case was studied in more details for ψ n > 0 and a BH with similar characteristics as in the Hawking’s treatment was obtained, if we work in the region of validity. In this work, we also compared these two regimes, by analyzing not only the relationship between heat capacity and temperature, but also the phase transitions that the system may undergo. It has been found that, for small values of r H , the two formalisms are compatible. When r H grows, it reaches a value where the BH undergoes another transition, now related to a change of behavior in the heat capacity. In the case of the thermodynamics with the Hawking formalism, the heat capacity passes from a region where it decreases with increasing temperature to a region where it increases with temperature. In the local case the heat capacity for r H > r m i n always decreases with temperature going to negative values.

Based on the assumption that heat capacity and temperature are not sufficient to guarantee stability, the Hawking-Page phase transition, which applies to both Hawking’s formalism and local formalism, was also considered. It has been checked that, in both formalisms, a Hawking-Page phase transition exists as another attribute in favor of the stability of BH.

The case with n = 2 has also been contemplated. The relevance of this case lies in the fact that the power law parameter, ψ 2 , can act as a cosmological constant, that can be positive or negative in this region.

The idea now is to inspect the asymptotic region that occurs for the case n < 0 . With this analysis, it will be possible to study the limits of-Sitter and Anti-de Sitter and the behavior on the borders. Besides the asymptotic behaviors, it is also interesting to study the potential generators of the power law presage that can be clarified when discussing the problem all over the space. In future works, we intend to investigate this case, that might be important for the AdS/CFT context for the n = 2 situation. These aspects are under study and shall be the subject of a forthcoming work. It is possible to solve the stability problem for a negative ψ 2 and to use the same thermodynamical procedure that we have used here. The importance, in the AdS context, is the fact that this black hole presents an angular deficit, and it is important to study the holographic principle of this object. We have shown in [34] that the structure of the angular deficit in a theory in D = ( 1 + 3 ) is preserved on the ( 1 + 2 ) boundary and, in our case, the defect is preserved in D = ( 1 + 2 ) . The AdS4/CFT3 is very important in studying planar materials like graphene and topological insulators. The introduction of the defect in this type of systems gives us current properties in these materials and it may become a relevant aspect of our work, opening up new application possibilities for our model. Another issue that we can pick up to study is the behavior of the objects near these systems [35] with both regimes for ψ 2 . We can intensively study quark/anti-quark interactions and discuss the confining and non-confining transitions as well as the chiral symmetry breaking critical temperature [36] .


F. B. Lustosa would like to thank CNPq-Brasil for his MSc Fellowship at the Instituto de Fsica of Universidade Federal Fluminense. The authors express their gratitude to T.R.P. Caramês for clarifying discussions.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Lustosa, F.B., Guimarães, M.E.X., Ferreira, C.N., Neto, J.L. and Helayël-Neto, J.A. (2019) On the Thermodynamical Black Hole Stability in the Space-Time of a Global Monopole in f(R)-Gravity. Journal of High Energy Physics, Gravitation and Cosmology, 5, 587-611.


  1. 1. Buchdahl, H.A. (1970) Non-Linear Lagrangians and Cosmological Theory. Monthly Notices of the Royal Astronomical Society, 150, 1-8.

  2. 2. Morais, J., Bouhmadi-Lpez, M. and Capozziello, S. (2015) Can f(R) Gravity Contribute to (Dark) Radiation? JCAP, 1509, 041.

  3. 3. Takahashi, K. and Yokoyama, J. (2015) Equation of State of Dark Energy in f(R) Gravity. Physical Review D, 91, Article ID: 084060.

  4. 4. Nojiri, S. and Odintsov, S.D. (2011) Unified Cosmic History in Modified Gravity: From F(R) Theory to Lorentz Non-Invariant Models. Physics Reports, 505, 59-144.

  5. 5. Sotiriou, T.P. and Faraoni, V. (2010) f(R) Theories of Gravity. Reviews of Modern Physics, 82, 451.

  6. 6. Perlmutter, S., et al. (1999) Measurements of Ω and Λ from 42 High-Redshift Supernovae. The Astrophysical Journal, 517, 565-586.

  7. 7. Hu, W. and Sawicki, I. (2007) Models of f(R) Cosmic Acceleration That Evade Solar System Tests. Physical Review D, 76, Article ID: 064004.

  8. 8. Appleby, S.A. and Battye, R.A. (2007) Do Consistent F(R) Models Mimic General Relativity plus Lambda? Physics Letters B, 654, 7-12.

  9. 9. Starobinsky, A.A. (2007) Disappearing Cosmological Constant in f(R) Gravity. JETP Letters, 86, 157-163.

  10. 10. Ohta, N., Percacci, R. and Vacca, G.P. (2015) Flow Equation for f(R) Gravity and Some of Its Exact Solutions. Physical Review D, 92, Article ID: 061501.

  11. 11. Corda, C. (2009) Interferometric Detection of Gravitational Waves: The Definitive Test for General Relativity. International Journal of Modern Physics D, 18, 2275-2282.

  12. 12. Bekenstein, J.D. (1972) Black Holes and the Second Law. Lettere al Nuovo Cimento, 4, 737-740.

  13. 13. Bekenstein, J.D. (1973) Black Holes and Entropy. Physical Review D, 7, 2333.

  14. 14. Bekenstein, J.D. (1974) Generalized Second Law of Thermodynamics in Black-Hole Physics. Physical Review D, 9, 3292.

  15. 15. Hawking, S.W. (1975) Particle Creation by Black Holes. Communications in Mathematical Physics, 43, 199-220.

  16. 16. Cai, R.-G., Cao, L.-M. and Ohta, N. (2010) Black Holes in Gravity with Conformal Anomaly and Logarithmic Term in Black Hole Entropy. JHEP, 1004, 082.

  17. 17. Myung, Y.S., Kim, Y.-W. and Park, Y.-J. (2008) Thermodynamics and Phase Transitions in the Born-Infeld-Anti-De Sitter Black Holes. Physical Review D, 78, Article ID: 084002.

  18. 18. Biswas, R. and Chakraborty, S. (2011) Geometry of the Thermodynamics of the Black Holes in Hořava-Lifshitz Gravity. General Relativity and Gravitation, 43, 41-50.

  19. 19. de la Cruz-Dombriz, A., Dobado, A. and Maroto, A.L. (2009) Black Holes in f(R) Theories. Physical Review D, 80, Article ID: 124011.

  20. 20. Kibble, T.W.B. (1976) Topology of Cosmic Domains and Strings. Journal of Physics A, 9, 1387.

  21. 21. Vilenkin, A. (1985) Cosmic Strings and Domain Walls. Physics Reports, 121, 263-315.

  22. 22. Barriola, M. and Vilenkin, A. (1989) Gravitational Field of a Global Monopole. Physical Review Letters, 63, 341.

  23. 23. Caramês, T.R.P., Bezerra de Mello, E.R. and Guimarães, M.E.X. (2011) Gravitational Field of a Global Monopole in a Modified Gravity. International Journal of Modern Physics: Conference Series, 3, 446.

  24. 24. Marunovic, A. and Murkovic, M. (2014) A Novel Black Hole Mimicker: A Boson Star and a Global Monopole Nonminimally Coupled to Gravity. Classical and Quantum Gravity, 31, Article ID: 045010.

  25. 25. Multamaki, T. and Vilja, I. (2006) Spherically Symmetric Solutions of Modified Field Equations in f(R) Theories of Gravity. Physical Review D, 74, Article ID: 064022.

  26. 26. Caramês, T.R.P., Bezerra de Mello, E.R. and Guimarães, M.E.X. (2012) On the Motion of a Test Particle around a Global Monopole in a Modified Gravity. Modern Physics Letters A, 27, Article ID: 1250177.

  27. 27. Yu, H.-W. (1994) Black Hole Thermodynamics and Global Monopoles. Nuclear Physics B, 430, 427-440.

  28. 28. Man, J. and Cheng, H. (2013) Thermodynamic Quantities of a Black Hole with an f(R) Global Monopole. Physical Review D, 87, Article ID: 044002.

  29. 29. Faraoni, V. (2010) Black Hole Entropy in Scalar-Tensor and ƒ(R) Gravity: An Overview. Entropy, 12, 1246-1263.

  30. 30. Bardeen, J.M., Carter, B. and Hawking, S.W. (1973) The Four Laws of Black Hole Mechanics. Communications in Mathematical Physics, 31, 161-170.

  31. 31. Hawking, S.W. (1974) Black Hole Explosions? Nature, 248, 30-31.

  32. 32. Tolman, R.C. (1930) On the Weight of Heat and Thermal Equilibrium in General Relativity. Physical Review, 35, 904.

  33. 33. Davies, P.C.W. (1978) Thermodynamics of Black Holes. Reports on Progress in Physics, 41, 1313-1355.

  34. 34. Bayona, C.A.B., Ferreira, C.N. and Otoya, V.J.V. (2011) A Conical Deficit in the AdS4/CFT3 Correspondence. Classical and Quantum Gravity, 28, Article ID: 015011.

  35. 35. Ballon-Bayona, A., Ferreira, C.N. and Otoya, V.J.V. (2013) DBI Equations and Holographic DC Conductivity. Physical Review D, 87, Article ID: 106007.

  36. 36. Ballon Bayona, C.A. and Ferreira, C.N. (2008) Static Strings in Global AdS Space and Quark-Antiquark Potential. Physical Review D, 78, Article ID: 026004.