Regular black holes with $\Lambda>0$ and its evolution in Lovelock gravity

In this work it is shown that the thermodynamics of regular black holes with a cosmological horizon, which are solutions of Lovelock gravity, determines that they must evolve either into a state where the black hole and cosmological horizons have reached thermal equilibrium or into an extreme black hole geometry where the black hole and cosmological horizons have merged. This differs from the behavior of Schwarzschild de Sitter geometry which evolves into a de Sitter space, the ground state of the space of solutions. This occurs due to a phase transition of the heat capacity of the black hole horizon. To perform that analysis it is shown that at each horizon a local first law of thermodynamics can be obtained from the gravitational equations.


I. INTRODUCTION
Motivated by the many known analogies between the physics of black holes and thermodynamics, there has been a great interest to analyze how deep these relations hold and how much they can be extended [1][2][3][4]. Indeed, nowadays there are several research programs exploring the idea that gravity may be an emergent theory with thermodynamic roots. In reference [5] was shown that the Friedmann equations can also be interpreted as thermodynamic identities. Certainly, one can check that this is well studied topic in the literature [6][7][8][9][10][11]. For this work, however, it is more relevant to mention that the connection between the gravitational field equations at the horizons and the first law of thermodynamics. For instance, in [12] was shown that the Einstein equations in the presence of an anisotropic fluid (with a mass density ρ and a radial pressure P = −ρ) on any horizon can be shaped up into the form of a first law of thermodynamics P dV = T dS − dU. The temperature T is defined as the inverse of the Euclidean period 1 at the horizon. In [13] this was generalized to the Kerr-Newman solution and a time dependent horizons. Two different questions arise after this. First, what does happen when the geometry has more than a single horizon and thus there is more than a single temperature?. Considering thermodynamics the answer is that this must represent a non-equilibrium system and therefore a system that must evolve consequently [14]. The second question is, how are connected the local thermodynamics relations defined in [12], with the standard black hole thermodynamic relations which are essentially non-local? This work aims to shed some light on those problems. .

Schwarzschild de Sitter space and thermodynamics
Originally, the first law of thermodynamics of black holes was discovered by analyzing black hole geometries with a well defined asymptotic region. Roughly speaking, this law can be understood as the relation between the changes of the mass, angular momenta and electric charges, which are to be computed at the asymptotic region, and the changes at the horizon. In fact, the first law was re-derived in terms of the analysis of the variation of Noether's charges in [15].
From the above, it seems that spaces with cosmological horizons may challenge the exis- 1 This coincides with the thermal definition in [4] tence of a first law of thermodynamics since they have no asymptotic region and consequently conserved charges cannot be defined in the usual way. Fortunately, this is not the case. Following [16], for instance, the first law of the thermodynamics turns into a sort of balance equation of the heat flows between the two horizons. This is complete consistent with the picture of a non-equilibrium system mentioned above. This can be visualized by considering the Schwarzschild-dS solution in d-dimensions where the derivation à la Wald [15], with the black hole and cosmological horizons defined by r + and r ++ respectively, yields [14] T where T + and T ++ are the temperatures at the black hole and cosmological horizons respectively. dS + and dS ++ are the variations of the entropies/areas of both horizons. Indeed, the conservation can be understood as dQ + = −dQ ++ .
In this moment one may question what is the role of the local thermodynamic relations at each horizon of the proposal in [12]. Indeed, the existence of a local and non local first laws of thermodynamics may seem contradictory but actually it is not. For Schwarzschild-dS solution, with a mass parameter M, the local relations are respectively dM = T + dS + and dM = −T ++ dS ++ , and thus Eq.(1) is nothing but the combination both local thermodynamic relations. It must be stressed that the scenario in the presence of an (anisotropic) fluid is not that simple. It is worth to mention [17] where the local approach was applied to multiple horizon geometries.

Lovelock Gravity
In higher dimensions (d > 4) Lovelock gravities arise as a simple generalization of the Einstein gravity. For instance, Lovelock gravity's equations of motion are of second order.
The Lovelock gravity Lagrangian is given by where {α p } is a set of arbitrary coupling constants and L p is conventionally normalized such as where R αβ µν is the Riemann tensor and δ µ 1 ...µn ν 1 ...νn is the generalized n-antisymmetric Kronecker delta.
In principle, the Lovelock gravity equations of motion may have as many different constant curvature solutions as the highest power of the Riemann tensor presented in the Lagrangian [18,19]. This defines a set of effective cosmological constants [Λ 1 , . . . Λ n ]. These constant curvature solutions are to be identified with the ground states of the Lovelock gravity considered.
In [20] was studied the asymptocally locally AdS regular black hole geometries with a unique AdS ground state. In this work it will be studied the case with a unique de Sitter ground state. There are two cases of interest. The first one is when the equation of motion have, roughly speaking, the form (R − Λ) n = 0. In this case there is a single, but n-fold degenerated, ground state of curvature Λ > 0. This can be considered a kinetic selection of the ground state 2 . The second case is when only one of the effective cosmological constant is a real positive number. For simplicity, the case of interest in this work is when the Lagrangian is just a single term in the Lovelock series plus a cosmological constant, i.e., L = α n L n + α 0 L 0 . This theory is known as Pure Lovelock. See for instance [22][23][24][25]. This theory, for Λ > 0, has a unique de Sitter ground state ∀n. For n even, however, there is an additional anti de Sitter ground state which is not relevant for this work, and thus it will be ignored.
Finally, it is worth to mention references [26][27][28][29] where the connection between generic Lovelock gravity and the local thermodynamic relation at the horizons were studied.

Thermodynamic evolution
The thermodynamics of a solution with a de Sitter ground state is expected to differ from those which are asymptotically flat or AdS in many aspects. The most remarkable difference is due to an effective positive cosmological constant which usually implies the presence of a cosmological horizon, and thus the lack of an asymptotic region. Now, as mentioned above, if a geometry has a black hole and a cosmological horizon, this can naturally be interpreted as non equilibrium thermodynamic system [14,30] and 2 This case is analogous to the negative cosmological constant case analyzed in [21] therefore a system that must evolve. The evolution of the vacuum solutions of the nfold Lovelock gravity with a de Sitter ground state was studied in [14]. In that case the thermodynamics determines that geometry evolves from a two-horizons geometry into the de Sitter space, the ground state of the theory. This happens due to the temperature of the black hole horizon is always larger than the cosmological horizon one, while the heat capacities, defined as are always negative for the black hole horizon and positive for the cosmological horizon. By treating this system as any ordinary thermodynamic system, this implies that while the cosmological horizon aims to increase its temperature due to the hotter radiation it absorbs. The black hole horizon, on the other hand, aims to increase its temperature as well due to the negative heat capacity. The only compatible scenario is therefore that M must decrease during the process. This process will continue as long M = 0. The whole evolution can be interpreted as a process in which black horizon emits more and more radiation until eventually M is entirely radiated away into the cosmological horizon. Conversely, the analysis of the Pure Lovelock vacuum solutions shows that the temperature of the cosmological horizon in fact can be larger than the black hole one (within a region of the space of parameters). This has a profound impact in the evolution of the Pure Lovelock solutions as will be discussed below. The analysis of a a generic Lovelock theory in this respect will be presented in a future work.

II. EVOLUTION OF REGULAR BLACK HOLES
To begin with the analysis it is worth to recall that regular black hole solutions emerge from the idea that corrections, due to would-be quantum gravity, should give rise to the removal of the gravitational singularities presented in standard black hole geometries. The point is that those corrections can be mimicked by the introduction an anisotropic fluid which must satisfy among other conditions, see [20] for a discussion, to be strongly concentrated at the origin. For simplicity in this work it will be described only the static case. In Schwarzschild coordinates, the energy momentum tensor of any fluid must be given by where ρ = −p r . Moreover, due to the conservation law ∇ µ T µν = 0, This determines that this is an anisotropic fluid. Before to proceed it is convenient to define the mass function where Ω d−2 is the volume of the d − 2 dimensional transverse section. In order to mimic the quantum effects the mass density ρ must have a singles maximum at r = 0 and rapidly to go to zero with r. This implies that [20] m(r)| r≈0 ≈ C 1 r d−1 with C 1 a positive constant. This also implies that the geometry around r = 0 is a de Sitter space.
For spaces with an asymptotic region, i.e. for spaces where r → ∞ is a proper limit, m(r) must tend to a constant M which is proportional to the mass of the solution [20]. Now, although for spaces with a cosmological horizon the mass of a solution is not well defined, still it is necessary that the solution have a mass parameter M, whose vanishing defines the ground state, i.e., a de Sitter space. By the same token, m(r) must be also promoted to a function m(r, M) such as m(r, M)| M =0 = 0. Moreover, the mass function m(r, M) must increase with the mass parameter. Therefore, for any value of r. M will be called from now on the mass parameter of the solution In the next sections will be shown that the general case has three horizons, defined by the radii r − < r + < r ++ . As mentioned above, the construction of a non-local thermodynamics from the local thermodynamic is not direct in the presence of fluid, but still the analysis can be carried out to a certain degree thanks to the presence of the mass parameter M. First, it must be noticed that to perform the analysis f (r) must be formally promoted to f (r, M).
Once this is established one can proceed. Let r = a be a zero of f (r, M). Now, in order to map a solution whose horizon is defined by r = a with a mass parameter M into a solution with a horizon at r = a + da then one must consider the transformation in the space of solutions a → a + da and M → M + dM such as f (a + da, M + dM) = 0 or equivalently This defines, in turn, the relation a long the curve of solutions, It must be noticed that the definition of the temperature of a cosmological horizon requires certain considerations due to the orientation of the normal vectors. While the radial derivative is outward to the black hole horizon, it is inward to the cosmological horizon.
This defines the direction of the period clockwise for the cosmological horizon. Because of that it is necessary to define the temperatures as With this convention the both temperatures are positive. See for instance [17,31]. Under this convention one obtains the non-local relation This relation will be made explicit for the two Lovelock theories considered in the next sections.
In order to study the thermodynamics as for any other standard thermodynamic system one can define a heat capacity. In this case, however, there are two potential generalizations.
First, one can define [32,33] where sign = 1 for r = r + the black hole horizon. Analogously, sign = −1 for r = r ++ , cosmological horizon. This is a direct generalization of the definition of the heat capacity for vacuum (black hole) solutions. It must be noticed that although this is a local definition for each horizon, due to the mass parameter the changes are not independent and thus this provides with an element to analyze the evolution. For the vacuum solutions the mass parameter coincides with a local definition of energy M = U a=r + = −U a=r ++ [30,32].
The second defintion of a heat capacity comes from the local thermodynamic relations proposed in [12], in the form of Unlike the conventional definition of heat capacity, the equation (16) can not represent the usual (∂Q/∂T ) V or (∂Q/∂T ) S as the change of volume and entropy are both functions of horizon's radii, and thus not mutually independent. However, as we will see below, Eq. (16) will be indeed useful to analyze the evolution of the solution. On the other hand, it is straightforward to notice that Eq.(16) differs from Eq.(15) in some formal aspects. In the the next sections, however, when equations (15) and (16) will be analyzed for some particular examples, it will be noticed that both provide essentially the same information about the evolution.

III. PURE LOVELOCK REGULAR BLACK HOLES
Pure Lovelock theory, as mentioned above, is defined by the Lagrangian L n + L 0 and its black hole solutions in vacuum have been studied. See [22,34,35] for instance. In [22] are discussed the thermodynamics features of the vacuum solution with a de Sitter ground state.
For Λ > 0 the equations of motion for a regular black hole were discussed originally in [20]. The gravitational equations are actually just reduced to solve where Ω d−2 stands for volume of the transverse section. By direct integration, where d − 2n − 1 > 0 and m(r) is the mass function defined in Eq. (7).

A. Structure of the horizons
Although, no analytic expression for m(r, M) has been given yet, still it is possible to analyze the structure of zeros of f (r, M). This can be done by drawing the functions whose intersection defines the zeros of f (r, M). In the Fig Although, for this picture one particular m(r, M) has been chosen, it is straightforward to realize that the behavior must be generic, as along as ρ satisfies the criteria mentioned in [20]. Therefore, there are up to three zeros of f (r, M) for an arbitrary function m(r, M).
These are to be called r − < r + < r ++ . The outer zero defines the presence of a cosmological horizon r ++ , and thus defines the outer spatial boundary of the space. r − and r + , represent an internal and black hole horizon, respectively.

B. Thermodynamics redone
In order to analyze the thermodynamics of these spaces, r = a will denote a generic horizon subjected, for now, only to satisfy f (a) = 0. The direct manipulation of equation (18) yields One can noticed that although m(M, r) is a function of M and r, once evaluated on a horizon becomes just a function of the horizon radius, in this case a. Following ref. [13], by evaluating the equation of motion (17) at r = a, identifying ρ = −P , yields, or equivalently This relation can be understood as a constraint in the evolution of the solution, within the space of parameters, along the change a and a + da. As mentioned above, the same relation can be obtained by requiring that f (a + da) = 0 along the transformation a → a + da.
The interpretation of Eq.(22) as a local thermodynamics at each horizon arises directly [12,22,32]. For a = r + one identifies, Notice that the f ′ | r=r + > 0 is satisfied strictly. In addition, m(r + ) can be interpreted as a local notion of energy. With this in mind, the equation takes the thermodynamic form Analogously, If a = r ++ one needs to identify Temperature cosmological horizon ( , Energy cosmological horizon (U ++ ) = −m(r ++ ), to obtain P (−dV ) = T dS − dU. Notice that f ′ | r=r ++ < 0. On the other hand, P (−dV ) has the correct sign since the inaccessible region (where f < 0) is now outside the cosmological horizon and thus the volume of the outer region change as −dV if r ++ → r ++ + dr ++ with dr ++ > 0 [12].
Finally, the heat capacities C + and C ++ mentioned above in equation (16) in this case One can notice that the second factor of in Eqs. (33) and (34) is expression of the temperature of the vacuum solution. It is straightforward to demonstrate that the second factor in both expression, for r + and r ++ , is positive. Therefore, the signs of the specific heats C + and C ++ depend only on the sign of ∂T /∂a.
On other hand, from equation (15), following the convention in [33], the K heat capacity of the black hole horizon (r + ) is given by The expression for r ++ is analogous. Obviously, an explicit expression for m(r, M) is required to obtain the exact expressions, however, it is direct to see that the three first factors are always positive for a = r + and a = r ++ . Therefore, the sign of the specific heat depends only on the sign of the factor ∂T /∂a as previously for Eq. (33,34).
Before to continue it is worth to notice the relation between the variation of the mass parameter and the the horizon's radii defined by Eq. (14). In this case these are given by and therefore one can argue, as there is a single mass parameter, that as M increases then r + must increase but r ++ must decrease, and viceverse. This will be very useful in the next sections to connect the thermodynamic evolution with the evolution of the horizons r + and r ++ .

IV. LOVELOCK THEORY WITH n-FOLD DEGENERATED GROUND STATE AS A THERMODYNAMIC IDENTITY
The asymptotically locally AdS vacuum solutions of the n-degenerated ground state Lovelock have been studied in references [21,37]. The regular solution for negative and positive cosmological constant were obtained in [20], to our knowledge for first time. The corresponding equation of motion reduces to solve whose solution is Here d − 2n − 1 > 0. The mass function m(r, M) is defined by the equation (7) as previously.
If the mass function fulfills the criteria listed in reference [20] the solution (38) also can be a regular black hole. Moreover, the structure zeros of f (r, M) in this case can be obtained by drawing the pair of functions whose intersection defines the zero. One can notice that this shares the generic behavior seen in Fig.(1). Therefore, there are at most three horizons, r − < r + < r ++ for a range of

A. Local Thermodynamics
To deduce the local thermodynamics one must consider a horizon defined by r = a.
Therefore, from Eq.(38) one can obtain It is worth to stress that m(M, r), once evaluated on a horizon, becomes just a function of the horizon radius. Now, by evaluating the equation of motion (37) at r = a and ρ = −P and multiplying by da yields In this expression one can recognize the definition of variation of entropy obtained in [37].
The heat capacities are given by Here one can recognize the third factor as the previously known definition of the temperature of the vacuum solutions. This term is always positive. Therefore the sign of the heat capacity of the each horizon depends only on the sign of the factor ∂T /∂a.
Again, from equation (15), and following the convention of reference [33], the K-heat capacity is given by and analogously for r ++ . One can notice that the four first factors are positive for both a = r + and a = r ++ . Therefore, the sign of the K + and K ++ depends only on the sign of ∂T /∂a.
Before to continue it is worth to notice the relation between the variation of the mass parameter and the the horizon's radii defined by Eq. (14). In this case these are given by As for the Pure Lovelock solution, in this case one can argue that as M increases then r + must increase but r ++ must decrease, and viceverse.

V. A PARTICULAR MASS DENSITY
Although the analysis above was carried out without the use of an explicit mass density ρ, still one can consider to analyse one in particular to notice how the results above become manifest. Because of that, in this section it is analyzed a new d dimensional generalization of the mass density proposed in [17,38]: where M is the mass parameter. R is a parameter of units of length L. V is given by The direct integration of the mass density in Eq.(7) yields the mass function The mass function Eq.(49) behaves as Eq. (8) for r ≪ R ≪ l, and therefore this defines proper Lovelock regular black holes. which implies that as the mass decreases r + must decrease but r ++ must increase.

C. Energy at the horizons
The next analysis to be done comprehends the change of the local definition of energy mentioned in Eqs. (27,32). In Figs.(8, 9, 10) are depicted, for EH, Pure Lovelock and n-fold degenerated GS Lovelock solutions respectively, the change of energy as a function of the horizon radii. It is direct to notice that Although the exact form of the figures above depends on the precise model considered, it is straightforward to show that this behaviour is generic for any other ρ(r), as along as the conditions discussed in [20] are satisfied. One can notice that the curve has a local maximum. This corresponds to r + = r ++ and therefore when the geometry has a external zero temperature horizon. In this case C + = C ++ vanish, which is in agreement with the usual notion of a zero temperature object.       the existence of a particular value r + = r * such that Therefore, temperature of black hole horizon reaches a local maximum at r + = r * . As dT dr * vanishes at r * then the heat capacity C + = dU + dr + / dT + dr + diverges in this point. See below.
It must be noticed that, at cosmological horizon, the temperature is always an increasing function of r ++ , therefore:       Black hole horizon temperature and cosmological horizon temperature: In order to simplify the discussion it is convenient to define M = M * as the mass parameter associated with black hole horizon r + = r * mentioned above. In figure (14),

E. Heat Capacity
First of all, it is worth stressing that at the cosmological horizon, by using equations (51) and (53), the heat capacity is always positive, i.e., therefore during the interaction the cosmological horizon would always aim to reach thermal equilibrium with the radiation it is absorbing [14]. The heat capacity of the black hole horizon is displayed for some particular values of d and n in figures (16) for the EH solution,  whereas at the right side the heat capacity is negative.
For values of the M parameter M ∈]M cri1 , M cri2 [, where there are three horizons, we observe that, for each ordered pair (r + , r ++ ), such that r + > r * , C + < 0 and C ++ > 0.
Conversely, for r + < r * , the local heat capacity becomes C + > 0. This is due to the phase transition mentioned above. The evolution will be discussed in the next section.

VI. EVOLUTION OF THREE HORIZON SOLUTION
It is worth to stress that the local definitions of heat capacity K and C share the same behavior. Both have changes of sign and phase transitions in the same points, due that • Above the phase transition: In region of M * < M < M cri2 let us consider the case where the temperatures satisfy T ++ < T + 3 . Recall that that K + < 0 and K ++ > 0. This implies that the cosmological horizon increases its temperature as is absorbing hotter radiation coming from the black hole horizon. The black hole horizon also further increases its temperature as it is absorbing, but with a K + < 0, cooler radiation from the cosmological horizon. The only consistent scenario is that the mass parameter M decreases. This can directly translate, see Eqs. (36,46) into a decrease of r + and an increase of r ++ .
• Below the phase transition: For M just below M * the temperatures satisfy T ++ < T + , but now both K + > 0 and K ++ > 0. As previously, the cosmological horizon increases its temperature as it is absorbing hotter radiation coming from the black hole horizon. Now, due to K + > 0, the black hole horizon also aims to thermalize with the cooler radiation coming form the cosmological horizon and thus it must reduce its temperature. Since the mass parameter M is unique (for both horizons) the only consistent scenario is that M must further decrease. As previously, this implies that a further decrease of r + and a further increase of r ++ . This process is expected to continue until both horizons reach the same temperature or M = M eq .
The evolution process can be sketched from figure (19).
If instead of M > M eq the starting point would have been a point where M cri1 < M < M eq , or equivalently T + < T ++ , then the evolution must be in the opposite direction. This is due to K + > 0 and K ++ > 0 and therefore the black hole horizon must increase its temperature simultaneously with a decrease of the temperature of the cosmological horizon.
In this case M must increase its value until thermal equilibrium between both horizons is reached. This implies that the black hole horizon radius r + increases while the cosmological horizon radius decreases of r ++ . This evolution has be sketched in figure (20). The previous line of thought determines that any three horizon regular black hole must evolves until both black and cosmological horizons reach thermal equilibrium. This pictures is a complete different scenario compared with non-regular black hole geometries, such as Schw-dS, which evolve into a de Sitter space, the ground state [14]. M CR < M cri2 . Now, although this a almost singular behavior still it is very relevant due to its consequences in the evolution.
In this case the evolution is complete different than the previous case. The heat capacity of the black horizon is always negative and the heat capacity of the cosmological horizon is always positive. As the black hole radiation has a lower temperature than the temperature of the cosmological horizon, this later must aim to reduce its temperature. Conversely, as the temperature of the black hole horizon is lower than the temperature of the cosmological horizon, then the black hole horizon is absorbing radiation and aims to reduce its temperature due to its negative heat capacity. In this way, both horizon aims to reduce their temperatures. The only compatible scenario is the mass parameter must increase and therefore the geometry must evolve into the merging of the black hole and cosmological horizons. Beyond the merging point the evolution cannot be studied by thermodynamics.
This behavior was sketched in figure (21).
On the other hand, if M ∈]M cr1 , M CR [ the evolution is into the thermal equilibrium at M = M eq . The analysis is analogous to the previous subsection (VI A).
It is remarkable feature that in this case M CR actually defines a transition in the evolution of the system from evolving into thermal equilibrium between the horizon or into a extreme black hole geometry. Moreover, it is also remarkable that the mass value M CR doesn't define a phase transition in the (local) thermodynamics of any of the horizons, but define a property of the system/geometry as a whole. M ∈]M * , M cr2 [ is contained in this case as will be shown. The first thing to notice is that the cosmological horizon decreases its temperature, since K ++ > 0, and thus, the mass parameter must increases. By other hand, the black hole horizons has specific heat K + > 0, and thus, increases its temperature up to its maximum value in M = M * . Next, for the range in M ∈ [M * , M cri2 ], the cosmological horizon continue to decrease its temperature due to K ++ > 0, and the black hole horizon also aims to further reduce its temperature due to K + < 0. Thus, both horizons reduce their temperatures until both reach T + = T ++ = 0, or equivalently r + = r ++ and M = M cri2 . This process has be sketched in figure (22).

VII. CONCLUSIONS AND DISCUSSION
In this work it has been studied the problem of black solution in the presence of a anisotropic fluid that mimics the some of quantum effect that could prevent the formation of singularities. This is done for two different Lovelock theories that share to have a single de Sitter ground state. Similar to the vacuum scenario these solution present a cosmological and black hole horizons. However, the presence of the fluid introduces a third, and inner, horizon. These solutions can have two and one horizons as well, being the two solution is what can usually refer to as an extreme or zero temperature solution. There is always at least one horizon, but this can happen in two very different scenarios. In fact, it is shown that there is a range of the mass parameter of solution M ∈ [M cr1 , M cr2 ] where the solution present three horizons.
The thermodynamics of these solutions is explored by a local method that provides with a local first law of thermodynamics at each horizon. In this respect in this work is shown that the local method can be implemented for Lovelock theories in the presence of anisotropic fluids. Although, the more traditional, and non-local, thermodynamic relations cannot be obtained directly, in this work some non-local relations are obtained as well.
Analysis of the solutions by means of the local and non local thermodynamic relations mentioned above share a similar thermodynamic behavior. This, on the other hand, precludes that any solution must evolve until black and cosmological horizons reach thermal equilibrium or in some particular cases evolve into an extreme (black hole) geometry. The first scenario differs from the vacuum solutions, which under a similar analysis, evolve into the ground state, a de Sitter space. The evolution into the extreme black hole geometry is shared by the vacuum solutions in the corresponding cases.
We provide a new toy model of regular black hole, whose energy density is a d dimensional generalization of [17]. This one confirms the results just mentioned explicitly.