A note on entropy of de Sitter black holes

A de Sitter black hole or a black hole spacetime endowed with a positive cosmological constant has two Killing horizons -- a black hole and a cosmological event horizon surrounding it. It is natural to expect that the total Bekenstein-Hawking entropy of such spacetimes should be the sum of the two horizons' areas. In this work we apply the recently developed formalism using the Gibbons-Hawking-York boundary term and the near horizon symmetries to derive the total entropy of such two horizon spacetimes. We construct a suitable general geometric set up for general stationary axisymmetric spacetimes with two or more than two commuting Killing vector fields in an arbitrary spacetime dimensions. This framework helps us to deal with both the horizons in an equal footing. We show that in order to obtain the total entropy of such spacetimes, the near horizon mode functions for the diffeomorphism generating vector fields have to be restricted in a certain manner, compared to the single horizon spacetimes. We next discuss specific known exact solutions belonging to the Kerr-Newman- or the Plebanski-Demianski-de Sitter families to show that they fall into the category of our general framework. We end with a sketch of further possible extensions of this work.


Introduction and motivation
Since the discovery of the accelerated expansion of our universe, there has been intensified interest in the study of spacetime physics with a positive cosmological constant, Λ. A tiny positive Λ can very satisfactorily explain the data corresponding to the current cosmological evolution of our universe. Thus, black hole spacetimes endowed with a positive cosmological constant are expected to provide us reasonable and physically well motivated models to study the global properties of the black holes living in our current universe.
One of the most exotic features of spacetimes endowed with a positive cosmological constant perhaps is the existence of a cosmological event horizon, when the parameters of that solution obey certain conditions [1]. It can be regarded as a complementary part of the black hole event horizon, arising due to the repulsive effect due to positive Λ, at large length scales. The cosmological event horizon acts as a causal boundary surrounding us, as no communication along a future directed path is possible beyond it [1]. When we have a black hole located inside the cosmological horizon, we call the entire spacetime a de Sitter black hole. Thus for such black holes the natural region of interest is the region between the black hole and the cosmological event horizon.
Just like the black hole event horizon, the cosmological event horizon also creates particles [1], and possesses thermodynamical properties, see. e.g. [2]- [11] and references therein for some recent developments in this direction.
Given that a black hole has thermodynamical properties and it creates particles, it has been an exciting topic to understand this from the symmetry and the microscopic point of view. In particular, for black hole spacetimes, the Bekenstein-Hawking entropy can solely be derived utilizing the near horizon conformal properties [12]- [18]. Precisely, these formalisms use some suitable fall-off conditions near a horizon and investigates under action of which vector fields these near horizon structures are 'preserved'. The algebra of charges corresponding to these symmetry generating vector fields gives a Virasoro algebra with a central extension. Then one uses the Cardy formula [19]- [20] to determine the entropy of the spacetime, which, under suitable choice of the mode functions, coincides with the Bekenstein-Hawking entropy. In other words, the entropy of a black hole spacetime is solely determined by the local symmetry at its boundary, i.e the event horizon.
We further refer our reader to [21] and references therein for a recent review on various approaches to understand black hole thermodynamics, the derivation of the Bekenstein-Hawking entropy and its quantum corrections.
A recent approach for deriving the thermal properties of a Killing horizon can be found in [22]- [25]. This novel formalism uses the the Noether current associated with the variation of the Gibbons-Hawking-York boundary term to obtain the conserved charges associated with the diffeomorphism generating vector fields. It turns out that the variation of the surface term gives a boundary integral of the Noether current located on the event horizon. Then the requirement of the near horizon symmetry gives the Bekenstein-Hawking entropy of the spacetime. This method has also recently been applied to time dependent cosmological black hole spacetimes [26], and to conformal scalar hairy black holes [27].
Let us return to our focus -black holes in the de Sitter spacetime. Since such a spacetime has two Killing horizons, it is natural to expect that its entropy will be the sum of the two horizon areas, e.g. [2]. Our precise goal in this work is to derive this total entropy using the formalism of [22]- [25]. Since such spacetimes has two natural boundaries, i.e. the two horizons, we get two surface integrals located at the two horizons. Then we shall find out the modes of the vector fields that preserve both the horizons' structure to find a Virasoro algebra 'effectively' representing the whole spacetime, giving us the total entropy (Sec. 3). To the best our knowledge, this has not been done before. We refer our reader to [28] for a verification of the second law of thermodynamics with this definition of entropy for some specific cases. See also [29] for some further discussions motivating this definition of the total entropy, from the point of view of area theorem.
Let us now emphasize a rather peculiar feature of this total entropy. We consider the Schwarzschildde Sitter spacetime, where M is the mass parameter. For 3M G √ Λ < 1, we have two Killing horizons located at, where the smaller root r H is the black hole and the larger r C is the cosmological horizon. For 3M G √ Λ = 1, r H and r C merge to 1/ √ Λ, known as the Nariai limit. The entropy and the temperature of the black hole and the cosmological horizons are respectively given by (A H /4G, κ H /2π), and (A C /4G, κ C /2π), where A denotes the horizon area, and κ H and −κ C are respectively the two horizon's surface gravity (with κ C > 0, and κ H ≥ κ C ) [1]. Now, if we consider the variation of the total entropy of the spacetime, S = (A H + A C )/4G, we obtain a Smarr formula with an effective equilibrium temperature, T eff = κH κC 2π(κH +κC ) , e.g. [7,11]. This implies that even though the two horizons have different characteristic temperatures, there can be an effective thermal equilibrium state when their entropies are combined. This has been demonstrated earlier in [30] via semiclassical tunneling method, but to the best of our knowledge, any field theoretic derivation of this is yet unknown.
Thus, it is highly motivating to understand this effective thermal equilibrium state, and quite naturally, a first step would be to actually derive the total entropy of such spacetimes.
The derivation of the entropy of the cosmological event horizon in any dimension can be found in [31], where suitable fall-off condition on an asymptotic de Sitter spacetime and the near horizon symmetry has been used. A discussion on the relation between the Friedman equation and the Cardy formula can be seen in [32]. We also refer our reader to [33] for a discussion on phase transition of de Sitter black holes using the aforementioned effective thermal equilibrium state.
The paper is organized as follows. In the next section we outline the general geometric scheme in which we work in. This will help us to deal with general stationary axisymmetric spacetimes with arbitrary number of commuting Killing vector fields in arbitrary spacetime dimensions. Precisely, this will provide us a timelike (non-Killing) vector field that foliates the spacetime between the two horizons. This vector field becomes null and Killing on both the horizons. This will enable us to identify Rindler geometries in the vicinity of both the horizons, and hence to treat them in an equal footing. Apart from the existence of the cosmological event horizon, we shall not assume any further explicit form or fall-off for the metric there. In other words, our method works well for the Nariai class de Sitter black holes as well, where the two horizons have comparable length scales.
In Sec. 3, we use this general framework to derive the total entropy of such spacetimes extending the formalism of [22]- [25]. We point out there that the mode functions corresponding to the symmetry generating vector fields near the horizons have to be restricted in a certain manner, compared to the single horizon spacetimes, in order to obtain the total entropy. Sec. 4 is devoted to address known non-trivial explicit examples from four and higher dimensions, to demonstrate they all fall under the general framework we built. We also address two cases of non-minimal couplings here and show in particular, some qualitative difference for the Brans-Dicke field for our case, when compared to the asymptotically flat spacetimes. Finally, we discuss our results in Sec. 5.
We shall work with mostly positive signature of the metric and set c = k B =h = 1 throughout, but will retain Newton's constant, G. In different spacetime dimensions, different values of it will be understood.

The general near horizon geometry
We shall derive below the general geometric framework we will be working in. The first part of which essentially deals with the construction of Killing horizons in stationary axisymmetric spacetimes of general dimensions with two or more commuting Killing vector fields. This will help us to deal with both the black hole and the cosmological horizon in an equal footing and in a much convenient manner than dealing with exact solutions case by case. The essential details of this can be found in e.g. [34,35] and references therein. For the sake of self consistency and convenience of the reader, we shall briefly outline them here.
The next part consists of identifying an (1 + 1)-dimensional geometry in a general way, such that it becomes the Rindler on any of the two Killing horizons.
We assume that the spacetime is an n-dimensional torsion-free manifold and satisfies Einstein's equations. We assume that the spacetime is stationary, axisymmetric and is endowed with two or more than two commuting Killing vector fields, where ξ a and φ (i) respectively generates stationarity and axisymmetries and the vanishing commutators are represented by the vanishing Lie derivatives. Since the spacetime is stationary, not static, we take ξ a φ (i) a = 0 for all i = 1, 2, . . . m. Also, to allow sufficient generality in our method, we further assume that the axisymmetric Killing vector fields φ (i) 's are not mutually orthogonal as well. There can be additional spatial isometries orthogonal to the stationary Killing vector field ξ a , but for our present purpose we need to worry only about isometries non-orthogonal to ξ a .
For convenience, we shall first discuss the case of three commuting Killing vector fields (ξ, φ (1) , φ (2) ). Generalization to higher numbers or specialization to two such commuting fields will be clear from this, as we shall see below.
We assume that the (n−3)-dimensional spacelike surfaces orthogonal to these three commuting Killing vector fields form integral submanifolds, which essentially means the vector fields spanning the subspace form a Lie algebra between themselves, which in turn implies Frobenius-like conditions [35], Clearly, the chief difference between the static and stationary axisymmetric spacetime is that, for the later the timelike Killing vector field is not hypersurface orthogonal. For our convenience, we shall now construct a foliation of the spacetime by constructing a hypersurface orthogonal (non-Killing) vector field. To do this, we define a 1-form χ a as so that χ a φ (1)a = 0 = χ a φ (2)a identically everywhere, giving where we have written for the norms, φ (1) · φ (1) = +f 2 1 and φ (2) · φ (2) = +f 2 2 . Let ξ a ξ a = −λ 2 . Then the norm of χ a is given by Since for any two spacelike vector fields A and B, we always have A · B ≤ |A||B|, the denominator of the second term on the right hand side is negative. This shows that χ a is timelike when β 2 > 0. The price we have paid doing this orthogonalization is that, χ a is not a Killing field in general, In terms of χ a , the first of Eq. (4) can be written as φ where µ ia (i = 1, 2, 3) are 1-forms orthogonal to χ a , φ (1)a and φ (2)a , and ν i (x) (i = 1, 2, 3) are functions. These functions and 1-forms can be determined exactly, chiefly using the commutativity of the Killing vector fields [34], which shows χ a satisfies the Frobenius condition, χ [a ∇ b χ c] = 0, and hence is orthogonal to the family of (n − 1)-dimensional spacelike hypersurfaces containing φ (i) a 's. This is a crucial result for proceeding further.
Having obtained the foliation of the spacetime, we now proceed to define the Killing horizons. Eq. (10) shows by the torsion-free condition that, for any β 2 = 0 hypersurface (say H), so that on any such hypersurface H, we may write where κ(x) is a smooth function defined on H. The next step is to prove that any such compact surface H is a Killing horizon, in the sense that the functions α i (x) (Eq.s (5), (6)) becomes constant on H. This involves constructing a null geodesic congruence for k a = e −κ(x)τ χ a (with χ a ∇ a τ := 1), and then solving for the Raychaudhuri equation on H. We shall not go into the details of this proof here referring our reader to [34] and references therein for this.
Then following similar steps as in four spacetime dimensions [35], we can show that κ is a constant on H, and is given by known as the surface gravity of the Killing horizon. We shall assume in the following κ = 0, always. Thus we have seen that the foliation timelike vector field χ a , smoothly becomes the horizon Killing vector field (say, χ a H ). For the de Sitter black hole spacetimes we wish to deal with, we have two such compact β 2 = 0 surfaces. The smaller one is the black hole event horizon and the larger one is the cosmological event horizon, and the vector field χ a smoothly becomes null and Killing on both of them.
The next step is to show using this general framework that, we can select a part of the near horizon geometry, which is Rindler-like. In order to see this, we define a 1-form which is orthogonal to χ a and φ (i)a , as can be seen by using (3). Also, Eq. (13) shows, X a X a = +1. Let X be the parameter along X a , such that X a ∇ a X := 1. We have from the action of a vector field on functions [35], which gives, X = β κ . If we choose X a to be one of the basis vectors orthogonal to χ a , it is clear that the metric infinitesimally close to H takes the form where the spacelike compact (n − 2)-section γ ab is orthogonal to both χ a and X a . It is clear that while χ H is tangent to H, the vector field X a defines orthogonality or 'away from' H. The 'χ − X' part of the near horizon coincides with the Rindler metric. We assume that the basis vectors spanning γ ab have neither vanishing nor diverging norms. This is just because otherwise we will have either vanishing or diverging horizon 'area' (:= (det γ ab )).
We note that the function κ 2 (x) = (∇aβ 2 )(∇ a β 2 ) 4β 2 , where β 2 is not necessarily vanishing this time, smoothly coincides with κ 2 (Eq. (13)) on H. Then it is clear that with the vector field X a = ∇aβ κ(x) , we may write the general spacetime metric as which smoothly coincides with (16) on any of the two Killing horizons. This helps us to deal with the black hole and the cosmological horizon in an equal footing.
To summarize, we have found a foliation of an n-dimensional stationary axisymmetric spacetime with three commuting Killing vector fields along a timelike vector field. Whenever that vector field becomes null on a compact surface, it becomes Killing as well, making the null surface a Killing horizon. For our concern, we have two such Killing horizons. We have also shown that the general spacetime metric (17) coincides with (16), in an infinitesimal neighborhood of any of the horizons.
For two commuting Killing vector fields, we set any one of φ (1) and φ (2) (and hence f 12 , Eq.s (6)) to zero. Also, there is only two conditions analogous to (4) now. They include ξ and any one of the axisymmetric Killing vector fields.
For more than three commuting Killing vector fields, say four, we assume that the (n − 4)subspace orthogonal to those Killing vector fields form integral submanifolds. Accordingly, we have four Frobenius-like conditions (6). We proceed then as earlier to define χ a as χ a = ξ a + α i φ (i) a with i = 1, 2, 3. We find α i by imposing the orthogonality between χ a and the axisymmetric Killing vector fields and proceed as earlier for the rest of the construction. This process may go on and clearly can accommodate arbitrary number of commuting Killing vector fields.
For our convenience we further define a new coordinate X = κ X/2 in (16) to get where X a is the tangent vector field associated with the new coordinate X, and X a X a = (2κX) −1 .
Finally, we note that since the Killing vector fields commute, we may specify coordinates along them, at least locally. On the other hand, since the horizon Killing vector field χ a H is a linear combination of those Killing fields with constant coefficients, it is clear that we can treat χ a H as a coordinate Killing vector field. We shall denote the surface gravities (κ) of the black hole and the cosmological event horizon by κ H and −κ C (with κ C > 0) respectively and will always work with the absolute value of the cosmological event horizon's surface gravity, in order to maintain the correct signature of the metric (18).
With all these necessary geometric ingredients, we are now ready to go into the derivation of entropy.

General derivation of the entropy
We shall use below the formalism developed in [22]- [25] using the Gibbons-Hawking-York surfacecounterterm in order to calculate the total entropy of a stationary axisymmetric de Sitter black hole spacetime.
Let us briefly review the formalism first. The Gibbons-Hawking-York surface term subject to the variation at the boundary is given by, where G is the Newton constant in dimension n, and [d n−1 x], [d n x] stand for the invariant volume measures in respective dimensions. K is the trace of the extrinsic curvature of an (n−1)-dimensional boundary hypersurface, ∂M. In the second integral, which is valid on the entire spacetime manifold M, K(x) and N a (x) are respectively a function and vector field that smoothly coincide with K and N a on ∂M, where N a is the unit normal to ∂M. The conserved Noether charge corresponding to the variation of the integrand of the second integral of (19) under infinitesimal diffeomorphism generated by a vector field ζ a is given by where J a is the conserved Noether current, J ab is an antisymmetric tensor field given by, . Σ is a suitable hypersurface, and the second integral is the (n − 2)-dimensional boundary of Σ. The choice of Σ is made in such a way that its boundary coincides with the (n − 2)-section of ∂M. √ γ is the determinant of the induced metric on that boundary, and is the area element. The vector fields N a and M a are chosen to be unit spacelike and timelike, respectively.
For a black hole spacetime, the natural choice of the (n − 1)-dimensional hypersurface in (19) is clearly the event horizon. The hypersurface Σ in (20) is a spatial hypersurface, and the (n − 2)dimensional subspace in (20) is the compact spatial section of the event horizon, spanned by angular coordinates.
The bracket algebra of charges (20) generated by different vector fields is given by where δ ζm Q[ζ n ] := Σ dΣ a £ ζm ( √ det gJ a [ζ n ]). The next step is to identify an infinite discrete set of diffeomorphism generating vector fields {ζ a m } which leave the near horizon geometry invariant. It can then be shown that for such vector fields, (21) can be identified with the Virasoro algebra with a central extension. One then uses the Cardy formula [19,20] in order to compute the entropy of the spacetime.
For the de Sitter black holes, we have two natural boundaries -the black hole horizon along with the cosmological event horizon. In other words, the so called 'bulk' of the de Sitter black hole spacetimes is the region between these two horizons. For such two natural boundaries, the first integral in Eq. (19) splits into two pieces -on the two hypersurfaces located at the two horizons.
Accordingly, the Noether charge in (20) will consist of two integrals at the two Killing horizons, similarly for the algebra satisfied by the charges, Eq. (21). The hypersurface Σ is orthogonal to the foliation vector field χ a , derived in the previous section. Since χ a smoothly becomes null and Killing (χ a H ) on both the horizons, we have obtained the two boundary integrals in a quite natural manner.
Let us first evaluate the charge corresponding to the horizon Killing vector field χ a H in Eq. (20). We have seen in the previous section that the spacetime metric formally looks the same in the neighborhood of both the Killing horizons, (18). We choose for the black hole horizon, N a = √ 2κ H XX a , and M a = χ a H √ 2κH X . The trace of the extrinsic curvature at the black hole event horizon is given by K H = − κ H /2X. Likewise, we get the unit vectors and the trace of the extrinsic curvature on the cosmological horizon, by replacing κ H with κ C . We find where we have defined the (n − 2)-dimensional 'area' as A = √ γd n−2 x, which correspond to the compact spatial sections at the two horizons.
We shall next obtain a Virasoro algebra for the charges generated by vector fields which preserves the near horizon structures. Eq. (21) for our case becomes, We shall look for the set of vector fields {ζ a }, which has only non-vanishing 'time' and X components, with respect the spacetime metric (18). Let us collectively denote the spatial coordinates and basis vectors tangent to the horizon by {Θ i } and {Θ a } respectively, spanning γ ab in Eq. (18). Let Φ(x) be the norm of any of the basis vectors {Θ a }. Clearly, in order to have the area of the horizon to be finite and non-vanishing, Φ has to finite and non-vanishing, too. Thus, in the neighborhood of a Killing horizon, we may expand Φ = Φ 1 (Θ i ) + O(β) + higher order terms in β. Also, since Θ a 's are tangent to the horizon, where β = 0, we must have Θ a ∇ a β → 0 in the infinitesimal neighborhood of any Killing horizon. Putting these all in together, using Eq.s (15), (16), (18) and the chain rule for the partial derivatives, we get to be vanishing in the neighborhood of any Killing horizon. This is analogous to the static and spherically symmetric case : the metric functions spanning the 2-sphere (at r = r H ) do not depend upon the spacelike Rindler coordinate, X, on or in the infinitesimal neighborhood of the horizon. In order to preserve the near horizon geometry, we must impose £ ζ g τ τ = 0 = £ ζ g XX , where τ represents the parameter or coordinate along the horizon Killing vector field χ a H . Solving these two equations involves only the Rindler part (τ − X) of the metric, and one obtains [22,23], where T is some smooth but otherwise arbitrary function. In terms of this vector field, Eq. (22) reads We now expand the function T in terms of infinite number of discrete eigenmodes as, T = m A m T m with m integer, so that for each m, we call the corresponding vector field as ζ a m . It is usual to choose T m = 1 l0 e im(l0τ +liφ i +g(x)) , where φ i 's are the parameters along the axisymmetric Killing vector fields, l 0 and l i 's are constants, g(x) = −l 0 dX 2κX . With this choice of the modes, ζ m 's satisfy an infinite dimensional discrete Lie algebra over a circle, where the subscript 'LB' denotes the Lie bracket. Thus the modes formally look the same on both the horizons, but the eigenvalue l 0 may be different. We shall call them as l 0H and l 0C respectively, for the black hole and the cosmological event horizon. Moreover, we have to fix them uniquely as well in order to derive the entropy of the whole system, as we shall see below. This is qualitatively different from the single horizon system discussed in [22,23], where one can leave l 0 completely arbitrary. Now, due to the axisymmetric geometry, the mode functions must be periodic in the Killing parameters φ i of the axisymmetric Killing vector fields. Then, since we have assumed the horizons to be compact, φ i 's are tangent to them and (26) becomes Likewise, the algebra of the charges, Eq. (21) gives for our two horizon spacetimes, which is a Virasoro algebra, effectively encompassing both the boundaries, and hence the bulk. We can identify the zero mode energy or the Hamiltonian and the central charge, C from Eq.s (28), , According to the Cardy formula [19,20], the entropy of the system is given by S = 2π CQ0 6 . Then it is clear that in order to get the entropy, we must set l 0H = κ H and l 0C = κ C in (30). This choice gives the same modes as one obtains via the method of 'asymptotic' fall-off near the horizon, e.g. [17].
With this choice, we get via the Cardy formula i.e. the total Bekenstein-Hawking entropy of the de Sitter black hole spacetimes. Before proceeding further, let us summarize what we have done so far. Since the de Sitter black holes have two Killing horizons, the Gibbons-Hawking-York surface counter-term splits into two pieces, corresponding to the two horizons. The hypersurface Σ appearing in Eq. (20) is the one which is orthogonal to the timelike vector field χ a , derived in Sec. 2. In terms of the two boundary integrals and choice of appropriate mode functions on the horizons, we have actually derived the total entropy of such spacetimes.
We have also seen that as long as the derivation of the total entropy is concerned, the choice of the mode functions are much more restricted than the single horizon spacetimes. Clearly, doing so is absolutely justified. We start with considering vector fields generating diffeomorphism in the entire spacetime, and due to the existence of boundaries, we only consider their explicit forms on the boundaries themselves. Since our spacetime has two horizons, the diffeomorphism generating vector fields ζ m 's assume different forms on them. In other words, this analysis can be regarded as 'doubly local' instead of 'local' [22,23], as the single horizon systems. Most importantly, both the surface gravities can be formally expressed by (13), so that κ H and κ C are nothing but the β 2 → 0 limits of the smooth function κ(x) = (∇ a β 2 )(∇ a β 2 )/4β 2 . Then it is clear that even though l 0H and l 0C have different values numerically, they are formally exactly the same.
We shall consider below some explicit and non-trivial exact solutions in order to demonstrate that they indeed fall under the scope of the general analysis we have done so far.

Explicit examples 4.1 The Kerr-Newman-and the Plebanski-Demianski-de Sitter families
Let us begin with the Kerr-Newman-de Sitter spacetime in four spacetime dimensions, whose metric in the Boyer-Lindquist coordinate reads where ∆ r = r 2 + a 2 1 − Λr 2 /3 − 2M Gr + q 2 , ∆ θ = 1 + Λa 2 cos 2 θ/3 with M , a and q are respectively the mass, rotation parameter and charge. For a Dyonic black hole, q 2 is understood as the sum of the square of the electric and magnetic charges. Firstly, since the timelike and axisymmetric Killing vector fields are coordinate fields (ξ a = (∂ t ) a , φ a = (∂ φ ) a ), they trivially commute 1 . The 2-planes orthogonal to these Killing vectors are spanned by the coordinate vector fields, (∂ r ) a and (∂ θ ) a . Hence they also commute to give a trivial Lie algebra and thus form integral submanifolds [35]. This was a crucial assumption made in Sec. 2. Also, as we have seen in Sec. 2, this guarantees the existence of the hypersurface orthogonal timelike vector field χ a = ξ a − (ξ · φ/φ · φ)φ a = ∂ a t − (g tφ /g φφ )∂ a φ . The black hole and the cosmological event horizon of (32) correspond to the largest and the next to largest positive roots of ∆ r = 0. The inner or the Cauchy horizon will not concern us for our present purpose.
The norm −β 2 of χ a , close to any of the horizons (∆ r → 0) is given by which is null on both the horizons. Also, it is eassy to check from the metric functions (32) that the function α = ξ·φ φ·φ becomes a constant whenever ∆ r = 0. Thus the Kerr-Newman-de Sitter spacetime falls into the general geometric category we discussed in Sec. 2.
Let us now discuss the derivation of the Rindler coordinate following Sec. 2. According to Eq. (14), infinitesimally close to any horizon, we choose where the subscripts or superscripts 'H, C' denote black hole an the cosmological horizon respectively. Using Eq.s (32), (34), we have (∂ θ ) a X a ∼ O( √ ∆ r ) in the infinitesimal vicinity of any of the horizons. Likewise, for the Rindler coordinate X a (Eq. (18)), it is easy to check that (∂ θ ) a X a ∼ O(∆ 1/4 r ). Thus infinitesimally close to the horizon, the vector fields ∂ a θ becomes orthogonal to X a or X a .
So, the Kerr-Newman-de Sitter spacetime takes the form given in Eq. (18) infinitesimally close to the horizons with γ ab spaanned by (∂ θ ) a and (∂ φ ) a . Thus, it falls into the general scheme discussed in Sec.s 2 and 3.
An exact asymptotically anti-de Sitter black hole solution with two independent rotation parameters in five dimensional minimal supergravity can be seen in [36]. From this solution we can obatain a de Sitter black hole spacetime via analytic continuation [10], where ρ 2 = r 2 + a 2 cos 2 θ + b 2 sin 2 θ , ∆ θ = 1 + a 2 g 2 cos 2 θ + b 2 g 2 sin 2 θ , Ξ a = (1 + a 2 g 2 ), The parameters M, a, b, q specify respectively the mass, two independent rotations and the charge of the black hole and g 2 is the positive cosmological constant. Clearly, being coordinate fields, all the three Killing vector fields (∂ t ) a , (∂ φ ) a , (∂ ψ ) a commute. The 2-planes orthogonal to them are spanned by coordinate fields along r and θ, which commute and trivially form a Lie algebra. Thus those two planes are integral submanifolds. This gurantees the existence of the timelike vector field χ a orthogonal to the (r, θ, φ, ψ) family of hypersurfaces, which from the discussions of Sec. 2 is written as The horizons of this spacetime are the positive roots of ∆ r = 0. A discussion on thermodynamics of this spacetime can be seen in [10]. Infinitesimaly close to any of the horizons, the norm of χ a takes the form, Thus χ a becomes null on the horizons. It is easy to check from the metric functions that infinitesimally close to any of the horizons, coefficient functions in (38) become constants where r H and r C denote respectively, the balck hole and the cosmological horizon radii. As earlier, it can also be checked easily that the near horizon Rindler coordinates can be defined and γ ab in Eq. (18) for this case is spanned by (θ, φ, ψ). Thus the general analysis of Sec.s 2 and 3 holds perfectly for (36).
We shall next consider the Plebansky-Demianski-de Sitter black hole spacetimes. This class is a geneailazation of the Kerr-Newman family, in the sense that apart from mass, charge, rotation and the cosmological constant, it contains additional parameters. The complete family of the Plebanski-Demianski-de Sitter class spacetimes which might represent de Sitter black holes has the metric [37], where Ω = 1 − α ω (l + a cos θ) r, ρ 2 = r 2 + (l + a cos θ) 2 , P = sin 2 θ 1 − a 3 cos θ − a 4 cos 2 θ The parameters α, ω, q, q m , ǫ and k are arbitrary, and a 3 and a 4 are determined via them. These parameters have their physical meaning in certain special sub-classes only. For example, for α = 0, the above metric becomes the Kerr-Newman-NUT-de Sitter solution [37], where ρ 2 = r 2 + (l + a cos θ) 2 , P = sin 2 θ 1 + 4Λal 3 cos θ + Λa 2 cos 2 θ 3 where q and q m are electric and magnetic charges and l is the NUT parameter. We shall consider the most general class given by (41), implicitly assuming it represents de Sitter black holes. The black hole and the cosmological horizons are the two largest roots of ∆ r = 0. It is easy to argue the existence of the hypersurface orthogonal timelike vector field χ a as earlier. The norm of χ a behaves as ∆ r → 0 as, which is null. It can also be verified from the metric functions that the function (ξ·φ)/(φ·φ) becomes a constant on ∆ r = 0. Thus the most general class of Plebanski-Demianski metrics (41) falls into our geometrical assumptions, and clearly, when it represents a de Sitter black hole spacetime (such as (43)), we can derive its total entropy. Finally, we shall briefly discuss the Kerr-de Sitter spacetime in generic spacetime dimensions n [38,39]. The metric of which reads in the Boyer-Lindquist like coordinates [10,38,39], where where N is the integer part of (n − 1)/2. The constant ǫ is +1 (0) for even (odd) spacetime dimensions. φ i 's are the coordinates along the axisymmetric Killing vector fields, and a i are the corresponding independent rotation parameters and g 2 is the cosmological constant. The coordinates µ i and ν are not independent, but are related via the constraint, ensuring the correct dimensionality of the spacetime.
Discussion of thermodynamic properties of this spacetime including that of the variation of the cosmological constant can be found in [10].
The spacetime is endowed with total (N +1) Killing vector fields. Since all of them are coordinate Killing vector fields, they commute. The (n−N −1)-dimensional subspace is spanned by coordinate vector fields {(∂ µi ) a , (∂ ν ) a }. Since the spatial coordinates along these vector fields satisfy the constraint (48), we may replace by using the chain rule of the partial derivatives, any specific ∂ µi as a linear combination of partial derivatives of the set {∂ µj , ∂ ν }, j = i. Then it is clear that the Lie brackets like [∂ µi , ∂ µ k ] a LB or [∂ µi , ∂ ν ] a LB , will always be a linear combination of the coordinate vector fields spanning this subspace. This is a Lie algebra, but not trivial as the earlier examples. Nevertheless, this ensures that the subspace orthogonal to the Killing vector fields form an integral submanifolds [35].
Then from this, the existence of the hypersurface orthogonal vector field χ a , and its Killing property when it becomes null follows as earlier.

Non-minimal couplings
Before we end, we shall address two cases of non-minimal couplings, the hairy black hole with conformal scalar field and the Brans-Dicke theory in the Jordan frame. Since for any non-minimal coupling the Ricci scalar term in the action gets modified as, f (ϕ) R, where ϕ is the scalar field, the presence of such coupling modifies the surface counterterm in Eq. (19), as well.
For static and spherically symmetric de Sitter black hole spacetimes, a derivation of Bekenstein-Hawking-Wald entropy using the surface counterterm has recently been done in [27], for the black hole event horizon. The surface term reads in this case, We shall briefly discuss this case for stationary axisymmetric spacetimes. It was shown in [40] that a black hole with conformal scalar hair and a positive cosmological constant cannot have a slow rotation in four spacetime dimensions. As was also argued there, it may be possible though, that a solution with generic rotation exists -at least one cannot rule out the possibility. In the following, we shall assume that such a solution indeed exists, and it falls into the geometric category described in Sec. 2. Following [27], the conserved Noether charge Q corresponding to the diffeomorphism, instead of Eq. (26), in this case becomes We note that for axisymmetric spacetimes, the field ϕ may depend on the non-Killing coordinates tangent to the horizons (such as the polar angle θ). Thus unlike the spherically symmetric spacetimes [27], we cannot pull out the scalar field out of the integration in Eq. (50).
Using the suitable mode decomposition as described in Sec. 3, we find where on any of the horizons we have defined, Likewise, we find the algebra of charges Setting l 0H = κ H and l 0C = κ C , we obtain the entropy. Finally, we shall discuss thermodynamics of black holes in the Brans-Dicke theory (see, e.g. [41] and references therein, for discussions on asymptotically flat spacetimes). The action of the Einstein-Brans-Dicke theory with a cosmological constant in the Jordan frame reads (see e.g. [42] and references therein) where ϕ is the Brans-Dicke scalar field, and ω is called the Brans-Dicke parameter. The inverse of the Brans-Dicke field, ϕ −1 acts as a spacetime dependent or dynamical Newton's 'constant'. Also, for ω = ∞, the field becomes a constant and the theory coincides with the Einstein gravity, ϕ = ϕ (0) = 1 16πG . The derivation of the surface couterterm for the Brans-Dicke theory is similar to that of (49), giving Now, for stationary asymptotically flat black hole spacetimes in four spacetime dimensions, the Brans-Dicke theory obeys a no hair theorem [43], which states that ϕ is necessarily a constant in the exterior of such black hole spacetimes. However, it does not constrain the parameter ω anyway. Thus, it was argued in e.g. [41] that the entropy of such black holes should be 4πϕ 0 A, where A is the horizon area and ϕ 0 is constant (corresponding to some finite value of ω), i.e not necessarily it equals ϕ (0) = 1 16πG . However, when we consider the de Sitter black hole spacetimes, it turns out that ϕ is not only a constant between the black hole and the cosmological event horizon, but also we must have ω = ∞ [42]. In other words, for such spacetimes, we must have ϕ 0 = ϕ (0) = 1 16πG . It is then clear that the total entropy of the de Sitter black holes in Brans-Dicke theory equals (31). Clearly, this is qualitatively different from the asymptotically flat black hole spacetimes.

Discussions
In this work we have considered the thermodynamics of stationary axisymmetric de Sitter black hole spacetimes. We have utilized the formalism developed using the Gibbons-Hawking-York surface couterterm and near horizon symmetries [22]- [25] to derive the total entropy of such two horizon spacetimes. To the best of our knowledge, this has not been done before. Since the spacetimes we have considered are endowed with two Killing horizons, the surface counterterm method provides us a natural convenience to deal with them. We have used a very general geometric framework to perform the derivation of the entropy in Sec.s 2 and 3. We have also seen in Sec. 3 that in order to do that, we have to choose uniquely, constants in the mode functions preserving the near horizons' geometries. Such restriction is not present for single horizon spacetimes.
This entire analysis can be thought of as doubly-local, instead of local as the single horizon spacetimes. Since nowhere in our analysis we assumed any precise asymptotic behaviour of the metric, our analysis is absolutely valid for the Nariai class de Sitter black holes, where the black hole and the cosmological event horizons have comparable sizes.
We also emphasize here that considering two boundary integrals in such two horizon spacetimes is the most natural choice.
After performing this general analysis, we have considered known non-trivial stationary axisymmetric solutions belonging to the Kerr-Newmann-or more general Plebanski-Demianski-de Sitter classes, and have demonstrated that for all of them the general derivation of the total entropy holds. We have also considered two cases of non-minimal couplings. In particular, we have pointed out a qualitative difference of the de Sitter black hole entropy in the presence of a Brans-Dicke scalar field, compared to the asymptotically flat spacetimes.
The analysis we have done above holds well as long as the spacetime is not exactly extremal, as the very use of Eq. (18) requires, even though κ could be 'small', but could never be vanishing. This is just because of the fact that the extremal black holes are qualitatively different objects from the usual non-extremal or even near-extremal ones. The vanishing of the surface gravity indicates vanishing temperature and this leads to a debate about the entropy of such black holes : should the entropy of such black holes vanish? We refer our reader to [44] for a recent discussion and a novel proposal about the calculation of entropy of extremal Schwarzschild-de Sitter spacetime, where the black hole and cosmological event horizon coincide. In that spacetime, there exists no timelike Killing vector field. It was shown in [44] that there exists a vector field, which is not Killing on the bulk, but becomes so and also null on the coinciding horizons, thereby enhancing the symmetry of the spacetime near the horizons. This is in accordance with the usual requirement that the near horizon structure must contain the conformal group. Rather curiously, this particular vector field gives bifurcate Killing horizons, spatially coincident with the original horizons we started with, which might give rise to a nonvanishing entropy. It should be interesting to attempt to relate the current formalism we have used, with that of [44]. Precisely, for non-extremal spacetimes, where the black hole and the cosmological horizon are separated, one might look for such a vector field and check its near horizon properties. The enhanced algebra containing this vector field might yield meaningful result for all values of κ, and hence might give a meaningful notion of the yet not-well understood κ → 0 limit. This procedure seems to have some qualitative similarity with what we did for the stationary axisymmetric spacetimes -we found a specific timelike non-Killing vector field χ a in the bulk, becoming Killing and null on both the horizons.
As we also have emphasized in the beginning, the current work is a step towards understanding aspects of de Sitter black hole thermodynamics when we treat the entire two horizon spacetime as a thermodynamic system as a whole. The variation of the total entropy of such spacetimes gives a Smarr formula predicting an effective equilibrium temperature. Thus the next step is to understand how or under what circumstances or with what choices of the field mode functions one can actually derive such thermal states using field theory. We hope to return to this issue in future works.