Regular black holes and its thermodynamics in Lovelock gravity

In this work two new families of non-singular or regular black hole solutions are displayed. These black holes behave as de Sitter space near its center and have a well defined AdS asymptotic region for negative cosmological constant. These solutions are constructed on a general ground through the introduction of a finite density of mass/energy. This removes the usual singularity of a black hole and also introduces a new internal geometry. The thermodynamic properties of these solutions are discussed as well.


Introduction
One of the most relevant predictions of General Relativity was the existence of black holes and nowadays there is substantial evidence that this is a usual phenomenon in nature. Now, despite the Schwarzschild black hole solution has been known for over a century it was not until the end of 60's that was shown that, under general classical conditions, that the formation of a black holes is unavoidable provided the energy (density) within a region of the space surpasses a certain limit. Indeed, the matter in that region collapses until an (event) horizon is formed. However, under the same argument Penrose determined that the final stage of that gravitational collapse gives rise to a singularity as predicted by the existence of Schwarzschild and Kerr solutions.
One unwanted consequence of the presence of singularities is that they break predictability. Fortunately, it is certain that classical general relativity cannot be valid at all scales [1] and thus at Planck scales the description of nature must change drastically. It is precisely this new description which is expected to provide a tamed version of the singularities once quantum effects are considered. This scenario is supported by results in either String Theory or Loop Quantum a e-mail: milko.estrada@ua.cl b e-mail: raros@unab.cl Gravity. Results in LQG, for instance, determine that before matter can reach the Planck density, quantum (gravity) fluctuations actually generate enough pressure to counterbalance weight. For the physics of a black hole, this implies that the gravitational collapse stops before a singularity can be formed. Furthermore, this can be understood as the formation of a dense central core whose density is of the order of magnitude of the Planck density. These objects are called Planck stars [2]. Once one has embraced the idea that inside a black hole, instead of a singularity, a dense core exists, one has to propose model for it. The first approximation to do this is to treat the problem as a classical gravitational problem with an energy-momentum density which condenses the quantum effects, in particular the existence of a pseudo repulsive force at the origin. In practice Planck stars can be studied as a geometry which far away from the core recovers a standard black hole solution, says Schwarzschild for instance, but whose center, although contains a dense core, can still be treated as a manifold. Moreover, the core of the geometry must approach, in a first approximation, a de Sitter space as such the geodesics diverge mimicking the repulsing force mentioned above. This kind of solutions are called non-singular or regular black holes and in the context of this work can be considered synonyms to Planck star.
Historically one of the first regular black hole solutions was found by Bardeen [3]. This corresponds to the spherically symmetric space described by where f (r ) = 1 − 2m(r )/r with m(r ) = Mr 3 /(r 2 + e 2 ) 3/2 . Evidently e is a regulator, but can also be understood, see [3], as due to an electric charge density whose electrostatic repulsion prevents the singularity to occur. One can check that f (r ) has a zero for r = r + > 0, showing the existence of a horizon. By the same token, it is direct to check the absence of singularities. Furthermore, with K > 0 showing that this space behaves as a de Sitter space near r = 0. After [3] others regular black holes have been also studied. See for instance [4][5][6][7][8][9][10][11][12][13][14][15][16][17] and for the higher dimensional case see [18][19][20]. As mentioned one must expect an energy density concentrated around r ≈ 0 which rapidly decays as r grows. This simple idea led, for instance, in [4] to propose an Gaussian energy density profile, i.e., ρ ∝ exp(−r 2 ) and to Dymnikova in [5,6] to propose that ρ ∝ exp(−r 3 ). For another approach to obtain a regular black hole from string theory, see for instance [21].
Lovelock gravity One of the fundamental aspects of GR, which almost singles it out in three and four dimensions, is that its equations of motion are of second order and thus causality is guarantied. In higher dimensions, d > 4, however having second order equations of motion is a property of much larger families of theories of gravity. Among them Lovelock gravities have a predominant rôle [22,23].
The Lagrangian of a Lovelock gravity in d dimensions is the addition, with arbitrary coefficients, of the lower dimensional topological densities [22,23], 1 i.e., where N = d 2 − 1 for even d or N = d−1 2 for odd d and where R αβ μν is the Riemann tensor and δ μ 1 ...μ n ν 1 ...ν n is the generalized n-antisymmetric Kronecker delta [24]. {α p } is a set of arbitrary coupling constants. The normalization in Eq. (4) is merely a convention.
The first two terms in this series are L 1 ∝ R, the Ricci scalar, and μν , the Gauss Bonnet density. In addition L 0 ∝ 1 is introduced to represent a cosmological constant term. In this case the equations of motion are the generalization of the Einstein equations given by where G μν (L L) = 1 √ g δ δg μν (L √ g) and T μν is the energy momentum tensor of the matter fields. Notice that ∇ μ G μν (L L) ≡ 0 is an identity.
To analyze the potential asymptotic behaviors of the solutions one needs to do a small digression. Let us consider that α p = 0 for p > I . It is direct to check that the equations of motion can be rewritten as This seems to indicate, as expected, that any Lovelock gravity should have for ground states constant curvature manifolds.
To analyze those backgrounds one can introduce the ansatz R Now, it is direct to demonstrate that, in general, the κ i can be complex numbers, even though ∀α p ∈ R. This does not only restrict the possible constant curvature solutions, and so the potential ground states of the theory, but also severely constraints the space of solutions with a well defined asymptotic region for arbitrary {α p }. Indeed, the only allowed behaviors of a ground sate as those which match a constant curvature and, in turn, are related with the zeros of P l (see Eq. (7)) in the real numbers. 2 For positive or null κ p , which correspond to locally AdS or flat, is possible to define an asymptotic region. That region also corresponds to the allowed asymptotic behavior of the solutions of that branch. The case κ p negative stands apart as in this case the ground state is locally dS and there is no asymptotic region. This can be called a dynamical selection of ground sates and simultaneously of asymptotic behaviors. It is worth mentioning that, as noticed in [25], in certain cases the definition of a ground state can be extended to non-constant curvature spaces. Those cases will be ignored in this work, however. There are several known black hole solutions of Lovelock gravity in vacuum (T μν = 0). See for instance [23,[26][27][28][29][30][31] and reference therein. However, there are not many known solutions in presence of matter fields, see for instance [32]. This is mostly due to the non-linearities of any theory of gravity, which makes difficult, if not impossible, to solve analytically its equations of motion for an arbitrary matter field configuration. Indeed, only highly symmetric configuration can be studied analytically. For the case of our interest, classically regular black holes have been studied in [19,20] within Einstein Gauss Bonnet theories.
In this work two new families of regular black holes will be displayed. These solutions share to belong to families of solutions which have a single locally AdS ground state. Moreover, these families have a single well defined asymptotically locally AdS region which approaches the ground state. These two new families correspond to the generalization of the Pure Lovelock solutions [28,30,32] and those discussed in [27] which have a n-fold degenerated ground state.
During the next sections, first the general conditions to be satisfied by the mass density will be discussed. Next, it will be obtained the two families of solutions and analyzed their behavior. Finally their thermodynamics will be displayed.

A well posed mass definition
In order to be able to solve analytically a black star, one can impose a highly symmetric geometry. With this in mind one can consider, as a first step, to neglect the presence of angular momenta. In principle, one could consider non spherical transverse sections as long as they are compact constant curvature (d − 2)−(-sub-)spaces, but for now this only would complicate the analysis. As it is well known in Schwarzschild coordinates a static spherical symmetric geometry can be described by It is worth to recall that the existence of event horizons is merely determined, due to the geometry, by the zeros of f (r ).
On the other hand, the energy momentum tensor of a fluid living in this geometry, given the symmetries of the space, must have the form Moreover, it must be satisfied that ρ = −p r as the lapse function is unitary. Finally, due to ∇ μ T μν = 0, In general this fluid is usually called an anisotropic fluid. Now we can proceed to analyze the general behavior of the mass density ρ. In the next section will be shown that it is convenient, not only to simplify the notation, to define the effective mass function In order to have a well posed physical situation it must be satisfied the following conditions; 1. ρ must be a positive due the weak energy condition and a continuous differentiable function to avoid singularities. This implies that m(r ) is a positive monotonically increasing function (m(r ) > 0 ∀r and m(r 1 ) > m(r 2 ) if r 1 > r 2 ) which vanishes at r = 0. 2. ρ must have a finite single maximum at r = 0, the core, (ρ(0) > ρ(r ) ∀r > 0) and to rapidly decrease away from the core. This yields the condition with K > 0 proportional to ρ(0). The finiteness and snootness of ρ(0) forbid the presence of a curvature singularity at r = 0 [4]. However, it must be noted that this is not enough to ensure a dS behavior near the center of the geometry, and thus additional conditions will be imposed in the next sections. 3. For the space to have a well defined asymptotic region, such as those to be studied, and to describe a physical object, ρ(r ) must be such that m(r ) be bounded for 0 < r < ∞ , i.e., with a well defined limit for r → ∞. Therefore, for M some constant. Later it will be shown that M is proportional to the total mass of the geometry. This implies that : 4. As mentioned above the idea of a regular black hole is to mimic the exterior of black hole. For this to happen the density ρ must be such that there is a radius r = r * where is satisfied m(r * ) ≈ M and d dr m(r * ) ≈ 0. In general, one can also expected that for large masses that P r * r + be satisfied. This condition, however, is not satisfied for masses within the range of Plank scales but still the thermodynamics can be studied [5,6].

First family of solutions: regular black holes in pure Lovelock theory
As a first step it will be considered a gravitational theory in d dimensions whose Lagrangian is a single term in Eq. (3) plus a cosmological constant, i.e., L = α n L n + α 0 L 0 . From now on α 0 , which can be understood as the cosmological constant, will be normalized such that with l 2 > 0. The symmetries of the ansatz considered (see Eq. (8)) leave just one (relevant) equation of motion to be solved. This is given by, The direct integration of Eq. (16) defines where with G n a constant of units L d−2n , where L represents a unit of length. Following the definition above, is the mass function defined above. For simplicity, in this work, we consider arbitrarily that the constant G n has a magnitude equal to 1 (see Appendix A) . Finally, by manipulating Eq. (17),

Global analysis
As mentioned above, any zero of f (r ) defines an event horizon in the geometry. This fact significantly simplifies the analysis. Now, to proceed, the cases Λ > 0 and Λ < 0 will be discussed separately.

Λ < 0 or negative cosmological constant
The first to notice is that for even n = 2k, with k ∈ N, and Λ < 0 (− sign in Eq. (20) ) f (r ) can take imaginary values.
To avoid that it is necessary that which occurs only for a certain ranges of r . This rules out the existence of an asymptotic region, defined by r → ∞, and even an interior region. Notice that there is not an (locally) AdS solution in this case and thus there is no ground state in the spectrum of the solutions. By observing Eq. (7), this can traced back to the fact the equation x 2k = −1 has no roots in R.
On the other hand, for odd n = 2k + 1 there are no constraints to be satisfied, and since m(r ) > 0, then there could be a well defined asymptotically AdS region in this case. f (r ) can be rewritten as Before to proceed a comment is to be made about the space of solutions. In principle one could have expected to have more than a potential asymptotic behavior in the spectrum of the solutions. However, from Eq. (7), one can notice that the allowed asymptotic behaviors are determined by x 2k+1 = −1, whose only real root is x = −1. This implies that an asymptotically (locally) AdS behavior of radius l is the only allowed asymptotic behavior in this case. As a consequence, in this work, it will be considered the case where n is odd when Λ < 0 in the Pure Lovelock solution as a representative of the general features. Now, assuming n = 2k + 1 one can analyze the structure of the horizons. The zeros of f (r ) can be studied qualitatively from the equation By recalling that m(r ) > 0, and in general, one can demonstrate the existence of up to two zeros of f (r ) and therefore the presence of up to two horizons. Those radii will be called r − < r + . Later it will be identified r + with the external radius of the black hole horizon. As usual there is an extreme case when both horizons merge (r + → r − ) and thus the space has zero temperature. To address analysis of the different cases requires of numerical analysis which be will carried out later in this work. For now, it is noteworthy that the absence of a horizon, unlike for the usual black hole solution, does not rule out the solution, as no singularities are presented. It only rules out a direct thermodynamic interpretation.

Λ > 0 or positive cosmological constant
For a positive cosmological constant (Λ > 0) the analogous of Eq. (21) in Eq. (20) is always positive and thus the family of solution contains a well defined ground state. Once again, this could have been foreseen from Eq. (7) which in this case becomes x n = 1 and has always zeros in R. Indeed, for n = 2k + 1, x = 1 is the only solution in R and determines that dS spaces is a solution and a ground state in this case.
Remarkably for n = 2k, x = ±1 are both solutions, and therefore, for Λ > 0 it is also possible to construct solutions with an AdS asymptotia. This produces a solution whose behavior mimics the previous case with Λ < 0 and n = 2k+1 mentioned above, thus no further attention is necessary in this case. Before to proceed, let us have a digression about the positive root which is connected with a dS ground state. The horizons can be analyzed by studying graphically the relation Since m(r ) > 0 strictly, it is direct to check the existence of up to three horizons in this case. These can be called r − < r + < r ++ , respectively. One utmost relevant consideration to be made is that r ++ defines a cosmological horizon but also the outer spatial boundary of the space, ruling out the existence of an asymptotic region, in the sense of a r → ∞ region, for the space. This condition separates this case from the asymptotically AdS case fundamentally. Finally, in the in-between region r − < r < r + can be interpreted as Kantowski-Sachs cosmology.

Internal geometry
The previous construction was performed without much concern of the details of ρ(r ). One can, in fact, use any of the densities mentioned above [1,[3][4][5]. Fortunately, the same can done for small r . As mentioned above, m(r ) must be such that with K > 0 a constant, and thus Eq. (20) for r → 0 behaves as: Moreover, one can check that the curvature invariants are finite and thus the geometry is smooth everywhere. Now, in order to have a non-collapsing region near r = 0, i.e., to describe a regular black hole [1], it is necessary to impose that at r ≈ 0, the geometry approximates a locally dS space [19]. This last condition determines that K > 1 l 2n , which in turn corresponds to a constraint on the mass density ρ(r )| r =0 . It is worth to mention that regular black holes with an effective AdS behavior near their origins have been studied in [19].

The non spherical symmetric solutions
As mentioned above it is possible to explore solutions with null and negative constant curvature transverse sections, instead of only the spherical symmetry. The geometry to be considered in this case is given by Here it was replaced the spherical transverse geometry in Eq. (8) by a closed (d − 2-) manifold of constant curvature γ . See [33] for a discussion, In this case the solution is given by It can be noticed that the analysis of the solution is analogous to the spherical symmetry. This is mostly due to the presence of the n-root of the same expression, whose analysis is independent of the value of γ . Therefore, the families of solutions with asymptotically locally AdS region are defined accordingly. For instance, the solution exists only for n = 2k + 1 and Λ < 0 or for one of the roots of Λ > 0, but the analysis of the horizon is the same for both cases. After recalling that m(r ) > 0 must be strictly satisfied, one can notice that for γ = 0 f (r ) has always two roots, r = 0 and r + = 0. One can also notice that r = 0 defines an irrelevant horizon for the discussion in this work and thus only r = r + defines a relevant horizon. For γ = −1 the situation is slightly different as there is always a single horizon. This also will be r + , as usual.

Second family of solutions: regular black holes with n-fold degenerated ground state
As mentioned above in general the Lovelock gravity might have more than a single constant curvature ground state, which makes those ground states unstable under dynamical evolution. One way to avoid this, observe Eq. (7), is by choosing the α p 3 such that P l (x) becomes P l (x) = (x ± l −2 ) n 3 Roughly speaking one can avoid this to occurs if the equation of motions have the form creating n−fold degenerated ground state of constant curvature ±l −2 . From now on only the negative cosmological constant case will be discussed as this the case where a genuine asymptotic region exists. Still some comment on the positive cosmological constant case will be done when they are straightforward. The static black hole solutions of these theories with negative cosmological constant have been studied in [27,33]. Remarkably, their generalization to the case under discussion is straightforward. The geometry to be consider in this case is given by Eq. (27) where the spherical transverse geometry in Eq. (8) has been replaced by closed (d − 1-) manifold of constant curvature γ , As for the previous solution the symmetries are strong enough to just left one equation to be solved. This is given by After some manipulations, where m(r ) can be in principle any function that satisfies the criteria introduced in Sect. 2 . It can be noticed, from Eq. (34), that, besides d > 2n, there is no restriction on the power n or the dimension d.

Horizons
As previously, the structure of zeros of f (r ), and thus of the horizons, can be studied qualitatively from analyzing graphically One first must recall that m(r ) > 0 which is a difference with case discussed in [33]. It is straightforward that the roots This can be obtained provided the parameters α p of Lovelock Lagrangian equation (3) are given by [27] where α n is a global coupling constant.
of f (r ) depend on γ = 0, ±1. For γ = 1, f (r ) has up to two different roots, to be called r − < r + , which can coalesce in an extreme case, As before, due to the absence of singularities, the solutions with no horizons cannot be ruled out. For γ = 0, r = 0 is a root and there is always a second root which will be called r + , For γ = −1 there is always a single root, to be called r + as well.

Limits of this solution
Recalling Sect. 2 one has the conditions to be satisfied by m(r ), which in this case imply and thus these solutions approach the locally AdS ground state for large r , as expected. In the same fashion, as mentioned in Sect. 2, as r → 0 is satisfied that m(r ) ≈ K r d−1 . This determines It can be noticed that this constrains K > l −2n in order to have a nearly (locally) dS geometry at the origin as required to model a regular black hole.

The thermodynamics before the thermodynamics
After studying the properties of these solutions in a general framework one can proceed to analyze the thermodynamics of the black holes in the spectrum of solutions. As mentioned previously, the central concern in this work are only spaces with a well defined asymptotically locally AdS region, but still some general comments will be made as well. For instance, although the number of horizons, and their properties, depends on m(r ), there is always a particular horizon which can be unambiguously cast as the horizon of the black hole. The associate zero will be denoted r + . This is even independent of the existence, for Λ > 0, of an even outer (cosmological) horizon on the geometry. It is utmost relevant to stress that the thermodynamics of black holes whose ground states are locally AdS or flat spaces differs of the thermodynamics of those whose ground state is a locally dS space. This is due to the obvious fact, mentioned already, that while for the latter there exists a well defined asymptotic region, for later there is no asymptotic region at all. It is well known that temperature and entropy can be defined independently of the existence of asymptotic region. However, the existence of conserved charges cannot, in particular the mass. In this way, a standard first law of thermodynamics, i.e. T d S = d E + . . ., can only be developed for spaces with a well defined asymptotic region. The lack of an asymptotic region, or the presence of a cosmological horizon, forbids such a kind of relations to exist. Moreover, in [34] was established that the presence of a black hole in geometry which also has a cosmological horizon can be interpreted as non-equilibrium thermodynamic system. That system, through the evaporation processes of both horizons, evolves into a de Sitter space. One can conjecture, for the solutions above, that the same can happen. That analysis, however, will be carried out elsewhere.
In what follows it will be analyzed the thermodynamics for asymptotically locally AdS spaces whose have a single well defined locally AdS ground state.

Thermodynamics of asymptotically AdS solutions
In this case the thermodynamics can be obtained following standard Wald's prescriptions based on the adiabatic change of the Noether charges in the region ]r + , ∞[ [35]. First, the mass of a solution can be determined in terms of the Noether charge associated with a timelike Killing vector in the asymptotic region. In this case above ∂ t in Eq. (8) is that Killing vector. In principle a second contribution to the Noether charge coming from action principle of the matters fields should be considered. However, and also in general, that second part of Noether charge usually becomes negligible at the asymptotic region. In the case at hand those additional parts of Noether charge will be assumed that do not contribute to the asymptotic value of the Noether charge.

Charges
In this particular case, after a straightforward computation, the Noether charge is given in terms of f (r ) and n = 2k + 1, is given by Here Reg is a regulator that removes the divergences associated with the AdS volume. 4 It can be noticed that to evaluate the total mass is only necessary to know the asymptotic value of m(r ). Indeed, if this value is fixed such that lim r →∞ m(r ) = M, then the final result, after regularization following [39], is wherẽ Here E 0 = 0 for even dimensions and the AdS vacuum energy for odd dimensions [40]. 5 This confirms that the mass parameter M indeed is proportional to the mass/energy of the solution.

Temperature
Given the form of the metric in Schwarzschild coordinates, one can follow the standard approach to obtain the expression of the temperature. This is given by One can notice that the last term, which is always finite, is responsible for the vanishing of the temperature for certain r + . However, until a particular ρ(r ) is given there is no form to determine analytically any further feature of T (r + ). In general, one can predict the existence of two inflection point, as T (r + ) can have a local maximum and minimum for certain values of r + . If so, then there is a range in [r + ] where three different r + yield the same temperature. As these two inflection points merge the temperature becomes a strictly increasing function. See below for a particular model.

Entropy
The form of the space together with the action principle for gravity allows to define a canonical ensemble for the thermodynamics of these solutions. By following the standard procedure defined by Wald, the first law of thermodynamics will emerged in its usual form with E = M the mass of the black hole and q i are the additional charges due to the presence of the matter fields and Φ i their corresponding potentials. Therefore, to address beyond the mass and the temperature is necessary to fix additional thermodynamic variables, in this case at the horizon, in order to obtain the first law of thermodynamics. Unfortunately to do that analysis one needs to have an action principle for the matter fields. Nonetheless, due its pure geometrical original, the entropy can be computed following the approach proposed by Wald [41]. The final result in d dimensions is given by where S 0 is a finite term independent of the value of r + and due to the topological pedigree of the regulator mentioned above [37]. It must be stressed that this expression is totally generic and functional independent of m(r ). The value of r + is not independent of the function m(r ), though. It can be noticed as well that this entropy does not follow an area law by a power r 2n−2 + but still defines an increasing function of r + . It is worth to stress that for the particular case n = 1 the area law is well defined.
To pursue any further into an analytic analysis of the thermodynamics properties, for instance the evaporation of these black holes, would require to define the temperature as an function of the mass and additional charges of the solution, which in turn needs at least to consider a ρ(r ) in particular. Unfortunately, this is not enough in general to obtain close expressions. Because of that in the next sections will be considered a particular case which will be analyzed through numerical methods.
6.2 Analysis of the second family of solutions

Charges
As before, mass can be computed by computing the Noether charge associated with ∂ t . After regularization (one can follow [37]), is given where E 0 = 0 for even dimensions and corresponds to the vacuum energy for odd dimensions. This is in complete analogy with [27,33] but for the Chern Simons case. The difference with results in [27,33] arise due to the location of the horizons.

Temperature
In complete analogy with the previous case the temperature can be computed. For the outer horizon this yields Notice that the last term is a correction to the expression obtained in [27] due to the finite density of mass. This correction is responsible for the vanishing of the temperature for a value of r + .

Entropy
In general the entropy can be computed as well in this case. This yields where Σ γ is the area of the transverse section. It is matter of fact that equation (47) is same expression obtained for the black hole solutions in [33] and the differences are due to the value of r + for a given value ofM. S 0 is due to topological term [37] added to regularize the action principle [42]. For large r + this entropy always approach an area law. It can be noticed as well that for n = 1 the usual area law is recovered. This also occurs for γ = 0 for any value of n or d. It is direct to check that for γ = 1 and γ = 0 these functions are monotonically increasing functions of r + . Far more interesting is the fact that S γ =−1 (r + ) can be negative which imposes severe constrains in the space of parameters. As previously, r + cannot be determined analytically and so to proceed a numerical approach would be required. This will be also considered for a next work.

Horizons
We see the behavior of the function f (r ) for the first family of solutions in Fig. 1 and for the second family of solutions in Fig. 2. The most relevant observation corresponds to the appearance of mass value where the solution evolves from a black hole into a regular star in the sense of a very dense object without a horizon. In this transition the temperature vanishes as mentioned above and can be seen in Figs. 1 and 2. As mentioned previously, that transition for non-regular black hole is meaningless, or nonphysical, as would produce a naked singularity. Therefore, for both solutions it can be distinguished three cases: -for M < M cri f (r ) has non zeros in the real numbers and thus the space is not a black hole. Moreover, due to the absence of singularities the space has a causal structure similar to a Minkowski space. -for M = M cri f (r ) has a double zero. This represents a regular extreme black hole with degenerate Killing horizon.
(b) f (r) for n = 3, d = 8. In general, one can notice that the temperature function has both a local maximum and a local minimum, says at r + = r * and r + = r * * respectively. As dT dr + vanishes in those points, the heat capacity C Q = d M dr + / dT dr + diverges. See bellow. It must be noticed, however, that for GR (n = 1) in four dimensions (d = 4) and for the second family of solutions with n = 2 and d = 6 the temperature is always an increasing function of r + and thus r * and r * * disappear.   The heat capacity can be defined for Q constant as This definition, following the rules of thermodynamics, allows to determine that black holes are locally stable/unstable under thermal fluctuation provided C Q > 0 or C Q < 0 respectively. With this in mind one can proceed with the analysis of the solutions. Figures 5 and 6 display examples of the behavior of the heat capacity. Notice that Figure  5a represents both families of solutions as n = 1.
In general, it can be observed a phase transitions at both r * and r * * . Going from C Q > 0 for r + > r * * to C Q < 0 for r * < r + < r * * . Finally C Q > 0 for r + < r * . Moreover, one can notice as well that C Q vanishes as T → 0.  It is worth to mention that for GR (n = 1) in four dimensions (d = 4) and for the second family of solutions with n = 2 and d = 6 there is no phase transitions. Moreover, although cannot be observed at the figures, the specific heat in this case is a positive function of r + which vanishes as T → 0 as well.

Conclusions and discussions
In this work two families of regular black hole solutions have been discussed. Although each family is a solution of a different Lovelock gravity in d dimensions both share to have a single ground state, which is approached asymptotically by the solutions and is defined by a single cosmological constant. These theories correspond to the Pure Lovelock theory [28,30,32] and to the theory discussed in [27] which have a n-fold degenerated ground state.
First it must be noticed that solutions have a minimum value of the parameter M, called M cri above, to represent a black hole geometry. As expected, these solutions asymptot-   ically, says for large r , are indistinguishable from the previously known solutions [27,30] in vacuum. On the other hand, both families of solutions near the origin, for r ∼ 0, behaves as maximally symmetric spaces which can be fixed, unambiguously, to be of positive curvature, i.e. the parameters can be fixed such as the solutions approach de Sitter spaces, as required to model a regular black hole, at their origins.
Concerning the thermodynamics of the solutions it was shown that, although both families of solutions differ, in general terms their thermodynamics presents the same features. The associate temperatures of the horizons have, in general, a local maximum and a local minimum at r * and r * * respectively. Indeed, one can notice that there can be three different values of r + defining the same temperature. Because of this last, the heat capacity changes from positive for r + < r * to negative at r * < r + < r * * and finally to positive for r * * < r + , signing out two phase transitions of the system, from thermodynamically stable (C Q > 0) into unstable (C Q < 0) and viceversa. Therefore, in general, there is a range of the black hole radii, r * < r + < r * * , or equivalently of the masses, where the black holes become thermodynamically unstable.
The existence of a phase transition for the vacuum solutions in [27] was known, and thus the phase transition at r * * could have been anticipated. However the existence of a second phase transition at r + = r * , and thus the existence a second stable range for r + < r * , is a new feature proper of these regular black holes.
In order to analyze the thermodynamics in detail it was considered as model the generalized Hayward energy density [7] defined in Eq,(48). Now, as mentioned in Sect. 2, the conditions to be satisfied by the mass density are quite general, and thus it can be quite interesting to explore other options in order to determine to what extent the arise of phase transitions is a model dependent feature. For instance, it would be quite interesting to study a generalization of the proposal in [5,6] A Units The unit system used throughout the text is slightly different from the ones presented in other articles. Because of that the unit system will be reviewed here in some detail.
To begin with it has been fixedh = 1 and c = 1, the light speed. Therefore, the mass/energy units is actually L −1 , where L is some distant unit. Moreover, it must be noticed that sinceh = 1 the action principle must be dimensionless.
As the α n coefficients must be fixed such that α n L n , with L n the n term in the Lovelock Lagrangian, be dimensionless, and as expected for consistency.