Charged de Sitter-like black holes: quintessence-dependent enthalpy and new extreme solutions

We consider Reissner-Nordstr\"{o}m black holes surrounded by quintessence where both a non-extremal event horizon and a cosmological horizon exist besides an inner horizon ($-1\leq \omega<-1/3$). We determine new extreme black hole solutions that generalize the Nariai horizon to asymptotically de Sitter-like solutions for any order relation between the squares of the charge $q^2$ and the mass parameter $M^2$ provided $q^2$ remains smaller than some limit, which is larger than $M^2$. In the limit case $q^2=9\omega^2 M^2/(9\omega^2-1)$, we derive the general expression of the extreme cosmo-black-hole, where the three horizons merge, and discuss some of its properties. We also show that the endpoint of the evaporation process is independent of any order relation between $q^2$ and $M^2$. The Teitelboim's energy and Padmanabhan's energy are related by a nonlinear expression and are shown to correspond to different ensembles. We also determine the enthalpy $H$ of the event horizon, as well as the effective thermodynamic volume which is the conjugate variable of the negative quintessential pressure, and show that in general the mass parameter and the Teitelboim's energy are different from the enthalpy and internal energy; only in the cosmological case, that is, for Reissner-Nordstr\"{o}m-de Sitter black hole we have $H=M$. Generalized Smarr formulas are also derived. It is concluded that the internal energy has a universal expression for all static charged black holes, with possibly a variable mass parameter, but it is not a suitable thermodynamic potential for static-black-hole thermodynamics if $M$ is constant. It is also shown that the Reverse Isoperimetric Inequality holds. We generalize the results to the case of the Reissner-Nordstr\"{o}m-de Sitter black hole surrounded by quintessence with two physical constants yielding two thermodynamic volumes.


I. INTRODUCTION
The inclusion of the P -V term in the first law of thermodynamics or in its familiar equivalent laws [1]- [11] has led to the notion of the effective thermodynamic volume, which is in general different from the geometric volume excluded by, say, the event horizon. The thermodynamic volume is the conjugate variable, with respect to some appropriate thermodynamic potential, of the pressure exerted on the horizon attributable to the presence of a constant cosmological density, or a quintessence, or both.
From this point of view, much more progress has been made for anti-de Sitter black holes [12]- [19] thanks to the AdS/CFT correspondence, the applicability of which has ever been extended [20]- [25]. The dS/CFT emerged as a possible dual relation relating quantum gravity on a de Sitter space to a Euclidean conformal field theory on a boundary of the same space [26]- [29]. Both these correspondences have motivated the classical and quantum investigations of the de Sitter-like and anti-de Sitter spacetimes.
The inclusion of the P -V yields, on the one hand, a generalized Smarr formula preserving a scaling law between thermodynamic variables and, on the other hand, an identification of the mass parameter with the enthalpy of the event horizon. These properties apply to both static and rotating black holes. In the static (nonrotating) case, however, a potential problem exists as noticed by Dolan [8]: the thermodynamic volume V is a function of the entropy S, and conversely, so one of the two variables, S or V , is redundant. This implies that the internal energy is not a suitable thermodynamic potential for the thermodynamic description of the static de Sitter and anti-de Sitter black holes. When rotation is included, the volume no longer depends on the entropy only, and so it is an independent thermodynamic variable.
In this work we consider Reissner-Nordström and Reissner-Nordström-de Sitter black holes surrounded by quintessence where both a non-extremal event horizon and a cosmological horizon exist besides an inner horizon. These are the asymptotically de Sitter solutions that correspond to −1 ≤ ω < −1/3. The case of asymptotically flat solutions corresponding to −1/3 ≤ ω < 0, where only a non-extremal event horizon and an inner horizon exist, was treated in Ref. [30], so we will not consider it here. As we shall see, some conclusions drawn, and results derived, in this work apply to the case of asymptotically flat solutions too.
In Sec. II we briefly review the Reissner-Nordström black holes surrounded by quintessence derived in Ref. [31]; some other of their properties are discussed in Ref. [30].
As is well known, extreme black holes, while instable, are important ingredients in the theory of quantum gravity where one can find a microscopic explanation of the Bekenstein-Hawking entropy [32]. Another type of extreme black holes, also instable, known as Nariaitype solutions [33]- [36] are also used in quantum gravity, where some singularities may be replaced by a Nariai type universes [37], besides their use for generating new solutions [38,39]. Some special Nariai black holes with quintessence have been discussed in [47]. In Sec. III we will determine explicitly new extreme black hole solutions that generalize the Nariai horizon [33,34] to all asymptotically de Sitter-like solutions (−1 ≤ ω < −1/3) for any order relation between the squares of the electric charge q 2 and the mass parameter M 2 provided q 2 remains smaller than some limit, which depends on ω and remains proportional to, but larger than, M 2 . In the limit case q 2 = 9ω 2 M 2 /(9ω 2 − 1), the three horizons merge and we derive the general expression of the extreme cosmo-black-hole and discuss some of its properties.
In the first part of Sec. IV, we consider the thermodynamics of the event horizon and investigate the evaporation process and its endpoint by providing the final values of the mass parameter and the radius of the event horizon.
The second part of Sec. IV is devoted to a discussion of the conserved charges, mainly, the energy. Because of the nonexistence of global timelike Killing vector for the de Sitter-like spacetimes, there is no notion of spatially asymptotically conserved charges which is similar to that for asymptotically flat or anti-de Sitter spacetimes. Other notions of conserved charges, however, exist but generally lead to different values of the charges. Using different approaches, some authors [40]- [46] were led to simple prescriptions when applied to asymptotically de Sitter-like black holes, among which we will discuss the Teitelboim's energy [42,43] and the Padmanabhan's energy [44][45][46]. We will relate these two notions of energy and show that they correspond to different ensembles. This will be clarified noticing, beforehand, that the notion of ensembles for the de Sitter-like spacetimes is larger than that of classical thermodynamics.
In Sec. V we will determine the enthalpy H of the event horizon, as well as the effective thermodynamic volume, and show that in general the mass parameter and the Teitelboim's energy are different from the enthalpy; only in the cosmological case ω = −1, that is, for Reissner-Nordström-de Sitter black hole we have H = M . Generalized Smarr formulas are also derived. It is concluded that the internal energy is not a suitable thermodynamic potential for the thermodynamic description of the static de Sitter-like black holes.
In Sec. VI we generalize the results concerning the thermodynamics to the case of the Reissner-Nordström-de Sitter black hole surrounded by quintessence with two physical constants. We conclude in Sec. VII.

II. REISSNER-NORDSTRÖM BLACK HOLES SURROUNDED BY QUINTESSENCE
In 4-dimensional spacetime, a spherical symmetric Reissner-Nordström black hole plunged into the field of a spherical symmetric quintessence has the metric [30,31] with f (r) = 1− 2M r + q 2 r 2 − 2c r 3ω+1 , −1 ≤ ω < 0 and c > 0. (2) With this notation and the convention G = = 1, the density of energy and isotropic pressure of quintessence are Here and in Ref. [30] c is half the opposite of its original value [31]. The convention used in Ref. [31] is such that 4πG = 1 where the expressions of (ρ q , p q ) have different constant coefficients. This same convention was used in Refs. [47]- [49] and partly in Ref. [30]. The black holes described by (1) and (2) are classified according to their asymptotic behavior [30] −1 3 ≤ ω < 0 : asymptotically flat solutions (4) − 1 ≤ ω < −1 3 : asymptotically de Sitter solutions, (5) with different physical properties depending on the sign of 3ω + 1. We worked out the case of asymptotically flat solutions in [30]. In this work we shall consider the asymptotically de Sitter solutions. This corresponds to −2 ≤ 3ω + 1 < 0 (−1 ≤ ω < −1/3). This will extend the work done in [30], which necessitated a special treatment different from the one we are aiming to perform here, to asymptotically de Sitter solutions.
As is well known, the thermodynamics of singular horizons is a subtle issue. In the case of metric (1), we assert on exploring its physical properties that all the scalar invariants diverge only at the singularity r = 0, as the density of energy and isotropic pressure do [Eq. (3)]. Particularly the curvature scalar takes the form This implies that all the horizons r h > 0, which are solutions to f (r) = 0, are regular. This point is important for the thermodynamic treatment we will present in Sec's. IV, V, and VI where no singular horizon is present.

III. NARIAI TYPE HORIZONS -EXTREME BLACK HOLES INSIDE COSMOLOGICAL HORIZONS -EXTREME COSMO-BLACK-HOLES
From now on we restrict ourselves to asymptotically de Sitter solutions where 3ω + 1 < 0. For fixed (M 2 , q 2 , ω), the number and nature of the horizons depend on the quintessence charge c. Setting u = 1/r, f (r) = 0 implies M 2 , respectively. We consider these cases separately and we determine, fixing (M 2 , q 2 , ω), the constraint(s) on c for which two or three horizons merge, the corresponding values of all the horizons of the solution (1), and the metric f .
In plot (a) of figure 1 the solution (1) has three horizons, a cosmological horizon (u ch < M/q 2 with M/q 2 being the value of u where the parabola has a minimum), an event horizon (u eh < M/q 2 ), and an inner horizon (u ah > M/q 2 ).
In plot (b) of figure 1, the two curves have a common tangent at u 1 < M/q 2 where u ch and u eh merge: this is a generalization of Nariai horizon [33,34] to all asymptotically de Sitter solutions −1 ≤ ω < −1/3. These solutions (1) possess another horizon u ah . This happens when c = c 1 (see similar discussion following Eqs. (2.11) and (2.12) of Ref. [30]) with and u ch = u eh = u 1 .

The cosmological constant case ω = −1
In this case the metric (1) reads Here u 1 is the common value of the horizons u ch and u eh . In the case q 2 ≤ M 2 < (9/8)M 2 of this section, 3M > m 2 > 0 and it is straightforward to show that u n < 0. Thus, there are only three positive roots to f = 0: u ah > u ch = u eh > 0. The Nariai type solution (10) generalizes the known neutral solution [33,34] to charged one. This is shown as follows. In the limit q → 0, we have lim q→0 Λ = 1/(9M 2 ), lim q→0 (1/u 1 ) = 3M , and the inner horizon disappears in the limit q → 0 as expected [33,34]. We also find: lim q→0 (1/u n ) = −6M and lim q→0 q 2 u ah = 2M yielding as in [33,34].

The case ω = −2/3
This case yields another new charged Nariai type solution: Here again u 1 > 0 is the common value of the horizons u ch and u eh . In the limit q → 0, we have lim q→0 c 1 = 1/(16M ), lim q→0 (1/u 1 ) = 4M , lim q→0 q 2 u ah = 2M (u ah disappears in the limit q → 0), and The plot (a) of figure 1 corresponds to c < c 1 and the plot (c) of the same figure, where the solution (1) has only a cosmological horizon u ch > M/q 2 , corresponds to c > c 1 .
The two curves will have two common tangents (for two different values of c), as shown in plots (b) and (c) of figure 2, provided (where the non-asymptotic flat condition −1 Constraints (14) being satisfied, the common tangents occur at (a) u 1 = u ch = u eh < M/q 2 if c = c 1 . In the cases ω = −1 and ω = −2/3, the new charged Nariai type solutions are still given by (10) and (12). Since the leftmost hand side of (14) is 9/8 and 4/3 for ω = −1 and ω = −2/3, respectively, these two new charged Nariai type solutions generalize the previous cases (10) and (12) where the event horizon merges with the inner horizon yielding an extreme black hole inside a cosmological horizon [plot (c) of figure 2]. c 2 and u 2 are given by For this case b), we again consider separately the cosmological constant case ω = −1 and the case ω = −2/3. The solution reads This is an extreme black hole inside a cosmological horizon where u 2 is the common value of u eh and u ah . In the limit q → M + we have lim q→M + c 2 = 0, lim q→M + (1/u 2 ) = M , and the cosmological horizon disappears (lim q→M + u ch = lim q→M + u ′ n = 0): this is the extreme Reissner-Nordström black hole, as expected, This is another extreme black hole inside a cosmological horizon: In the limit q → M + we have lim q→M + c 2 = 0, lim q→M + (1/u 2 ) = M , and the cosmological horizon disappears (lim q→M + u ch = 0): this is again the extreme Reissner-Nordström black hole where f = (r − M ) 2 /r 2 . If the constraints (14) are still satisfied but c 2 < c < c 1 , the three horizons exist as in plot (a) of figure 2. Otherwise, if q 2 /M 2 ≥ 1 + 1/(9ω 2 − 1), only the cosmological horizon survives as in the plot (d) of figure 2. However, when the equality holds, q 2 /M 2 = 1 + 1/(9ω 2 − 1), the three horizons merge, as in figure 3, and c 1 = c 2 . This case deserves a special treatment.
For this case we have simple expressions for u ch and c 1 derived as follows. By Eqs. (8), (9), (15) and (16) we Note that for the asymptotically de Sitter solutions (−1 ≤ ω < −1/3) all factors ω, 3ω + 1, and 3ω − 1 in (19) are negative, resulting in positive factors ω/(3ω + 1) and ω(3ω −1). We term this type of solutions where the three horizons merge by the extreme cosmo-black-holes. Eliminating M in (19), we write the radius r H = 1/u H of the extreme cosmo-black-hole as This is the most general relation expressing r H in terms of (c H , ω), (M, ω), or (|q|, ω) for asymptotically de Sitter solutions.
For ω held constant, r H appears as increasing linear function of M and of |q| with slopes 3ω/(3ω + 1) and [(3ω−1)/(3ω+1)] 1/2 , respectively. The slopes themselves are increasing functions of ω varying from 3/2 and √ 2, For M held constant, r H takes its minimum value 3M/2 at ω = −1 and increases monotonically to indefinitely large values as ω approaches −1/3 from the left. However, for fixed M , c H does not always vary monotonically. For instance, for M = 1, c H increases monotonically to (1/2) − , and for M = 0.3, c H first decreases to some minimum value then approaches (1/2) − as ω → (−1/3) − . Using (19), we obtain for M held constant have thus a huge (r H → ∞) extreme Reissner-Nordström cosmo-black hole with a huge electric charge but finite mass surrounded by a finite quintessence charge c → (1/2) − . From a physical point of view it is rather easier to figure out configurations where q is held constant than configuration where the mass parameter M is. It is straightforward to establish that, for q held constant, r H takes its minimum value √ 2|q| at ω = −1 and increases monotonically to indefinitely large values as ω approaches −1/3 from the left. Similarly to the previous case, c H does vary monotonically with ω and, using (19), we obtain for q held constant This also shows that lim ω→(−1/3) − c H = (1/2) − is independent of the value of q (held constant) as is lim ω→(−1/3) − M = 0. In the limit ω → (−1/3) − , we have thus a huge (r H → ∞) but massless extreme Reissner-Nordström cosmo-black hole with a finite electric charge surrounded by a finite quintessence charge c → (1/2) − . Finally, let us discuss the case where M is taken proportional to 3ω +1: M = −α(3ω +1) and α > 0. Eq. (19) leads to q 2 ∝ (3ω + 1) and 1 In the limit ω → (−1/3) − , M → 0, q → 0, and r H → α. There remains a pure quintessence state with a finite cosmological horizon r H → α and a finite quintessence charge c → 1/2. Now, we consider the special cases ω = −1 and ω = −2/3. For the cosmological constant case ω = −1, on applying Eqs. (19) and (20), we obtain the extreme cosmo-black-hole The case ω = −2/3 yields another simple extreme cosmoblack-hole The case ω = −2/3 was treated in details in [47] where more or less equivalent formulas to (12) and (18) were given, but no general formulas as (8), (9), (15), and (16), which are valid for the whole range of ω, were derived. Similarly, the general formulas (19) and (20) were not derived in [47] but only the relation r H = 1/(6c H ) was given (Eq. (30) of [47]) along with expressions of M and q in terms of c. In our notation [30], c is half its value in [47] and half the opposite of its original value [31].

IV. CONSERVED CHARGES AND THERMODYNAMICS
In this work, we only consider the thermodynamics of the event horizon and, from now on, we restrict ourselves to non-extremal solutions, that is, to cases where the three horizons do not merge with each other. For the asymptotically de Sitter solutions (−1 ≤ ω < −1/3), these are the solutions satisfying either one of the two following constraints: 1 In this case c H behaves as: c H = 1 2 − 3 4 (3 + 2 ln α)ǫ + 9 8 (3 + 6 ln α + 2 ln 2 α)ǫ 2 + O(ǫ 3 ) and ǫ ≡ −ω − 1 3 > 0.
Solutions satisfying the first constraint correspond to plot (a) of figure 1 and those satisfying the second constraint correspond to plot (a) of figure 2. The temperature of the event horizon is given by where one may eliminate M using (7). Note that −2M + 2q 2 u eh and 2(3ω + 1)cu eh 3ω are the derivatives, evaluated at the point u eh , of the functions y = 1 − 2M u + q 2 u 2 and y = 2cu 3ω+1 , respectively, which are plotted in figures 1 and 2. From the plots (a) of these two figures one sees that the slopes are such that which yields T eh > 0 for the black holes constrained by (23) or (24). As far as T eh > 0, the evaporation of the black hole proceeds by reducing the mass parameter M . Differentiating (7) with respect to u eh yields or, equivalently, Using the plots (a) of figures 1 and 2 one obtains at u ch and u ah similar inequalities to (26) but with the other order sign: 2(3ω + 1)cu ch 3ω < −2M + 2q 2 u ch and 2(3ω + 1)cu ah 3ω < −2M + 2q 2 u ah . This yields ∂M ∂r r ch < 0 and ∂M ∂r r ah < 0, where we have used the fact that ∂M/∂r ∝ M − q 2 u + (3ω + 1)cu 3ω . Using (28) and (29) we conclude that during the evaporation, the event horizon shrinks and the other two horizons expand. For black holes constrained by (23)  During the evaporation, all the other given parameters (q, c, ω) are held constant but the mass parameter, which decreases from some initial value M , constrained by (23) or (24), to some final value M f satisfying The evaporation ends when the value of c 2 , which depends on M as given in (15), ceases to vary and settles to the given value of c. The final mass of the extreme black hole inside a cosmological horizon is solution to the extremality condition [compare with (15) and (16)] which expresses the equality of the slopes in the plot (c) of figure 2, and u 2 is given by Equations (31) and (32) reproduce the correct result for the extreme Reissner-Nordström black hole. In absence of quintessence, c = 0 and Eq. (31) gives u 2 = M f /q 2 . In order to derive from Eq. (32) this same value for u 2 , which does not depend on ω, the only possibility is to set there q 2 = M f 2 which results in u 2 = M f /q 2 . As we shall see below, the static de Sitter spacetimes have different thermodynamic notions of energy but the above results derived in this section are independent of these different notions.
Rotating or static black hole solutions, of general relativity or extended theories, immersed in flat or antide Sitter spacetimes have well defined physical entities, these are the globally so-called charges (mass, electric and magnetic charges, angular momentum, scalar charges and so on). These solutions may have multihorizons among which one finds only one event horizon but no cosmological horizon.
By the static de Sitter spacetimes we mean all solutions where one of the metric components is of the form 1 − 2M/r + q 2 /r 2 − Λr 2 /3, with Λ > 0 and M and q are constants [in the charged Vaidya-de Sitter black hole, M and q are not constants but the metric is nonstatic [50,51]]. The static de Sitter spacetimes, with a non-extremal event horizon and a cosmological horizon, have the property that the thermodynamic temperatures of these two horizons are not equal [52].
We enlarge the above list of static de Sitter spacetimes by including all black holes with a non-extremal event horizon and a cosmological horizon, as the solutions given by (1), (2), and (5). This is called the set of the de Sitter-like spacetimes 2 . They have the property that, in general, the thermodynamic temperatures of the event and cosmological horizons are not equal (this property is violated by some black holes with conformally coupled scalar field, as the Martínez-Troncoso-Zanelli ones [54], where the two horizons have the same temperature [55]).
In the Euclidean formulation, this means that the imaginary time periods are not equal and, consequently, it is not possible to avoid the conical singularities at both horizons at once. Said otherwise, if one of the two horizons is treated as a thermodynamic system, the other horizon cannot be treated so because of the presence of the conical singularity there; it is, rather, treated as a boundary. Considerations by which both horizons are treated simultaneously as different thermodynamic systems, out of equilibrium, are subtle. An instance of that, one cannot add the entropies of the two horizons to obtain the total entropy of the "thermodynamic system made out of the two horizons", which, in fact, it does not exist.
The other issue with the de Sitter-like spacetimes is the definition and evaluation of the conserved quantities. Because of no spatial infinity accessible to observers, the notions of ADM mass, electric charge and other charges, which are defined in asymptotically flat or anti-de Sitter spacetimes, are no longer valid. Different prescriptions to define conserved charges exist however. These are known as Abbott-Deser (AD) [40], Balasubramanian-Boer-Minic (BBM) [41], and Teitelboim [42,43] prescriptions. A statement made in Ref. [3] asserts that the three prescriptions yield the same conserved charges when applied to static asymptotically de Sitter solutions: "Furthermore, for nonrotating case Teitelboim's charges are in full agreement with BBM/AD charges".
For the static de Sitter-like spacetimes, it is straightforward to generalize the expression of the Teitelboim's energy to these holes provided (1) In classical thermodynamics, the thermodynamic description of a system could be achieved using different forms of energy, which yields the notions of ensembles: micro-canonical, canonical, and grand-canonical. For the de Sitter spacetimes, it seems that the notion of ensembles is larger than what one usually encounters in classical thermodynamics. This is because there is not a universally agreed definition of the energy for this class of black holes with two or more horizons. Even within some elected definition of the energy, say Teitelboim's definition (33), it is possible to have some new emerged ensembles.
Another perception of the notion of energy for the de Sitter spacetimes is due to Padmanabhan [44][45][46]. Using a path integral approach one derives the Padmanabhan's energy E P of the event horizon by where u eh is a solution to (7). Using the latter equation we relate the two energies by The cosmological horizon has its corresponding, and similar, formulas to (33) and (34).
Since the evaporation of the black hole proceeds by reducing the mass parameter M . By (33), this results in a reduction of E T and should yield the same for E P . This is in fact the case since Eq. (28) is just Hence, dE T and dE P have the same signs for the black holes constrained by (23) or (24). The evaporation ends when T eh = 0 and E T reaches its minimal value M f given by Eqs. (31) and (32). E P also reaches its minimal value E P,f at the end of the process where E P,f = 1/(2u 2 ) with u 2 given by (31) and (32).
As claimed earlier in this section, for the de Sitter spacetimes, the notion of ensembles is larger than what one usually encounters in classical thermodynamics. This is why in this paper, we will not employ the classical thermodynamic terminology of micro-canonical, canonical, and grand-canonical; rather, we will describe ensembles by the constancy of the corresponding thermodynamic variable or potential. Instances are provided by the ensembles where either E T or E P is held constant. The ensembles E T = C 1 and E P = C 2 , where (C 1 , C 2 ) are constants, describe different thermodynamic systems since E T held constant is represented in the 3dimensional space (E P , q, c) by the 2-dimensional curved surface (35) where E T = C 1 . This shows that there is, in fact, no ambiguity in the definition of the energy for the de Sitter spacetimes: E T and E P are different energies just because they correspond to different ensembles, exactly in the same way as the Gibbs free energy differs from the Helmholtz free energy.

V. EVENT HORIZON THERMODYNAMICS: ENTHALPY VERSUS INTERNAL ENERGY
Since the Teitelboim's energy E T = M , we will, for time to time, insert E T in a couple of mathematical expressions of this section to show its relation to some thermodynamic potentials.
With the entropy of the event horizon given by we re-write the expression of M , using 1 − 2M u eh + q 2 u eh 2 = 2cu eh 3ω+1 [see Eq. (7)], as Considering (s, q, c) as independent thermodynamic variables and using (38), it is easy to establish 3 the generalized Smarr formula [30] M = 2T eh S + Aq + (3ω + 1)Θc (39) where are the electric potential on the event horizon and the thermodynamic conjugate of c, respectively. It is straightforward to check that T eh as defined in (25) is just The first law of thermodynamics takes then the form The last term in (42) does not have a direct physical meaning; rather, we prefer to introduce a new thermodynamic variable and its conjugate which both have a familiar physical meaning. These variables are the value of the pressure p q evaluated at the event horizon P ≡ p q | r eh and its conjugate, the thermodynamic volume, V . Using (3) along with (37) we obtain [In the cosmological constant case ω = −1 and Λ = 6c, P reduces to constant pressure P Λ = −Λ/(8π)]. In terms of the new independent thermodynamic variables (s, q, P ), Eq. (38) takes the form where it appears as homogeneous in (S, q 2 , P −1 ) of order 1/2. The Euler identity for thermodynamic potentials that are not homogeneous functions of their natural extensive variables yields [56] This shows that E T = M is not a familiar thermodynamic potential.
For ω = −1, Eq. (46) reduces to the following familiar formula: where is interpreted as the enthalpy and the conjugate of P , is the thermodynamic and geometric volume excluded from a spatial slice by the black hole horizon. This interpretation given to the mass parameter M (= E T ) of the static de Sitter spacetimes extends that for the static anti-de Sitter spacetimes [1]- [11]. We aim to extend this interpretation to all the de Sitter-like spacetimes. For the purpose of this paper, we will do that for the de Sitter-like spacetimes given by (1), (2), and (5), that is, we will enlarge the scope of the above-made interpretation to include the cases where ω = −1 by introducing a new thermodynamic potential. The new thermodynamic potential H is defined such that This is achieved upon adding to M , given by (44), the following term: where P is given by (43). Equation (50) reduces to (48) if ω = −1. H is homogeneous in (S, q 2 , P −1 ) of order 1/2. With (∂H/∂S) q,P = T eh and (∂H/∂q) q,P = A, the Euler identity yields where V is the thermodynamic volume, conjugate of P , given by which reduces to the geometric volume (49) if ω = −1.
The differential of H leads to the familiar well known first-law-equivalent formula dH = T eh dS + Adq + V dP, (53) by which H is interpreted as the enthalpy.
We have thus reached the conclusion that the Teitelboim's energy is in general not the internal energy or enthalpy of the event horizon of the de Sitter-like spacetime. The Teitelboim's energy or the mass parameter is related to the enthalpy by and the internal energy U = H − P V is related to M by The first law should read dU = T eh dS + Adq − P dV, but since V depends on S via (52), one of the two variables, S or V , is redundant. This shows that U is not a convenient thermodynamic potential for expressing the first law for the static de Sitter-like spacetimes. With M given by (44), the r.h.s of (55) reduces to a function of (S, q) only as it could be reduced, using (52), to a function of (V, q) only. This has been noticed for the static anti-de Sitter spacetimes where the thermodynamic volume V depends on the entropy S too so that "they cannot be varied independently and so V seems redundant. Indeed this may be the reason why V was never considered in the early literature on black hole thermodynamics. But this is an artifact of the non-rotating approximation, V and S can, and should, be considered to be independent variables for a rotating black hole" [8].
We should be able to do the same upon including rotation; we may pursue that in a subsequent work. It is worth noticing that ω, being dimensionless, cannot be considered as a thermodynamic variable.
It is also worth noticing that the expression of H, as given by (50) and (43), is totally independent of the definition of Teitelboim's energy. Throughout this section we have used the mass parameter M for the derivation of the expression (50) of H. Since E T = M , we have, from time to time, inserted E T to show the relation of E T to the enthalpy, as in (48).
We verify that the conjecture made in Ref. [9] concerning the Reverse Isoperimetric Inequality remains true for Reissner-Nordström black holes surrounded by quintessence. This large inequality reads where D is the dimension of the spacetime, A is the area of the event horizon, and A D−2 is the area of the unit (D − 2)-sphere. In our case, D = 4, A 2 = 4π, A = 4πs, and V is given by (52). This yields which is always true and the equality holds for Reissner-Nordström-de Sitter black hole.

VI. REISSNER-NORDSTRÖM-DE SITTER BLACK HOLE SURROUNDED BY QUINTESSENCE
Up to now, we only considered separately the case of (1) the Reissner-Nordström-de Sitter black hole, that is, the Reissner-Nordström black hole surrounded by a cosmological density, and the case of (2) the Reissner-Nordström black hole surrounded by quintessence. We aim to extend the results of Sec. V to the case of the Reissner-Nordström-de Sitter black hole surrounded by quintessence, that is, to the case where the Reissner-Nordström black hole is surrounded by a cosmological density and quintessence. On doing this we extend the phase space by including two physical constants (Λ, c) the variations of which yields two thermodynamic volumes. This extension should apply to any fundamental theory with many physical constants [9].
The metric f now takes the form The mass parameter M is expressed in terms of s by [compare with (38)] The temperature of the event horizon is no longer given by T eh ; rather, it is given by It is straightforward to generalize the results of Sec. V. For instance, Eq. (39) becomes where Θ ≡ (∂M/∂c) S,q,Λ is as given in (40) and Θ Λ ≡ (∂M/∂Λ) S,q,c = −s 3/2 /6. In a similar way we generalize Eqs. (50), (51) and (53) to where A, P , P Λ , V Λ , and V have the same expressions as in Sec. V. T is either given by (60) or by The internal energy is defined by U = H − V Λ P Λ − V P and retains its expression given by (56) as does the expression of M given by (54).
That the internal energy has the same expression for a Reissner-Nordström black hole and a Reissner-Nordström-de Sitter black hole both surrounded by quintessence seems to be a universal law that applies, not only to all de Sitter-like spacetimes, but to all static charged black holes even in the case where the mass parameter depends r. In Ref [57], we show that the internal energy of any static charged black hole, with possibly a variable mass parameter, is given by This depends only on the entropy of the event horizon and on the electric charge, which are the intrinsic properties of the black hole, and it does not depend on any extrinsic properties, as a cosmological density, quintessence, or any other force that may exert a pressure on the black hole. We have thus the following conclusion: The internal energy of any static charged black hole, with possibly a variable mass parameter, does depend explicitly only on the intrinsic properties of the black hole.
However, U depends implicitly on M and other physical constants through s which is a solution to f (s) = 0. It is worth noticing that the first term √ s/2 = r eh /2 is just Padmanabhan's energy and the second term q 2 /(2 √ s) = q 2 /(2r eh ) is an electric-energy contribution.
For Schwarzschild and Reissner-Nordström black holes, U coincides with the mass M and the pressure P is identically zero.
Finally, the Reverse Isoperimetric Inequality (57) is satisfied separately for V Λ and for V .

VII. CONCLUSION
We have determined the exact general conditions under which extreme solutions exist for the Reissner-Nordström black holes surrounded by quintessence with a negative quintessencial pressure. For q 2 ≤ M 2 , the only existing extreme solutions are generalizations of Nariai black holes. For q 2 > M 2 , but q 2 < 9ω 2 M 2 /(9ω 2 − 1), we may have both extreme solutions: extreme black holes inside cosmological horizons or generalizations of Nariai black holes.
In the limit case q 2 = 9ω 2 M 2 /(9ω 2 − 1) we were led to the extreme cosmo-black-hole solution where all horizons merge. The limit ω → (−1/3) − is characterized by the presence of (1) a huge extreme Reissner-Nordström cosmo-black hole with a huge electric charge but finite mass surrounded by a finite quintessence charge and vanishing event horizon pressure if the mass is held constant as the limit is approached; (2) a massless huge extreme Reissner-Nordström cosmoblack hole with a finite electric charge surrounded by a finite quintessence charge and vanishing event horizon pressure if the electric charge is held constant as the limit is approached; (3) a massless and neutral extreme Reissner-Nordström cosmo-black hole, with a finite radius, surrounded by a finite quintessence charge and nonvanishing event horizon pressure if the mass remains proportional to 3ω + 1 as the limit is approached.
We have shown that during the evaporation, the event horizon shrinks and the other two horizons expand, and that the final mass at the end of the evaporation is independent of the initial order relation between the squares of the electric charge and the mass parameter, provided the three horizons exist at the beginning of the process.
The inclusion of the P -V term has led to a consistent thermodynamic description of the first law of thermodynamics. The results obtained here generalize the results obtained for the anti-de Sitter spacetime, where the pressure exerted on the horizon is positive, as well as the results for the de Sitter one [1]- [11]. This shows that the sign of the pressure is irrelevant. We have commented that the internal energy has a universal expression for any static charged black hole, with possibly a variable mass parameter. We have also shown that the Reverse Isoperimetric Inequality holds.
The results concerning the thermodynamics were easily generalized to the case of the Reissner-Nordström-de Sitter black hole surrounded by quintessence with two physical constants yielding two thermodynamic volumes.
Phase transitions and critical phenomena will be discussed elsewhere.