{'fraction_non_alphanumeric': 0.04341767918849627, 'fraction_numerical': 0.0131813621669825, 'mean_word_length': 4.3511780495079035, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 20, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Cosmological horizons, Noether charges and entropy

It has recently been shown that, in the vicinity of their event horizons, black holes exhibit an infinite-dimensional symmetry. This symmetry captures relevant physical information about the black hole, and in particular about its thermodynamics. Here, we show that the same holds in a cosmological setup. More precisely, we show that around the de Sitter cosmological horizon, an infinite set of diffeomorphisms preserving a sensible set of boundary conditions emerges. These boundary conditions are similar to those considered by Price and Thorne in the context of the membrane paradigm, and permit to accommodate interesting gravity solutions. As for other boundary conditions considered previously, they are preserved by an infinite-dimensional asymptotic isometry algebra that includes supertranslations and two copies of the Virasoro algebra. This symmetry has associated a set of Noether charges that turn out to be conserved, finite, and integrable. We derive these charges explicitly using the covariant formalism and analyze their physical meaning by evaluating them to the case of de Sitter space. In this case, the zero-mode of the charge is found to account for the Gibbons-Hawking entropy of the cosmological horizon. We then consider a much more general set of solutions, including asymptotically de-Sitter black holes and asymptotically Taub-NUT-de Sitter black holes.


Introduction
Observations indicate that we live in a dark energy dominated universe whose late time regime is well approximated by de Sitter space [1]. The accelerating fate of the universe implies the existence of a cosmological event horizon, a surface that surrounds any observer and that defines a point of no return: Distant comoving galaxies scape behind the cosmological horizon and information from there can not be retrieved, as it happens with matter falling into a black hole. In fact, black hole horizons and cosmological horizons share many mathematical and physical properties. In particular, cosmological horizons are also thermodynamical objects; they have temperature and entropy [2]. As for black holes, the entropy of cosmological horizons obeys the area law 1 with A being the area of the cosmological horizon, set by the Hubble scale. This fact has profound implications: In particular, it means that quantum gravity in de Sitter space, if consistent at all, demands the product of the cosmological constant Λ and the Newton constant G to take special, non-continuous values [3]. It also raises the fundamental question about the actual reason for such special functional dependence with the area. In the last two decades, thanks to the holographic principle and its ramifications, we have gained insights about the meaning of the area law and of its generality. However, despite all these advances, a full understanding of it is still missing. The central question about the microscopic derivation of (1) has only been answered in idealized examples. In the case of black hole horizons, the question can be addressed in the context of string theory, but the situation is quite different with de Sitter space, for which an embedding in string theory has shown to be notably elusive [4]. Probably, the most important advances in understanding the thermodynamics of de Sitter space came from the dS/CFT correspondence [5]; but even in that context there are still open questions.
Also in the case of black holes, which certainly offer a much more tractable scenario, there are important conceptual features of their thermodynamics yet to be understood. Such is the case of the information loss paradox, a longstanding problem that has recently been revisited. One of the recent attempts to address this problem is the idea of Hawking, Perry, and Strominger that the study of an intriguing infinite-dimensional symmetries -known as supertranslations-that black holes exhibit in their near-horizon region [6] could shed some light on the information hidden by the horizon (see also [7,8] and references thereof). This study of infinite-dimensional symmetry is inspired by the recent application of similar ideas to the infrared physics in asymptotically flat space; see [9] and references therein. How much information is exactly encoded in this infinite-dimensional symmetry is still an open question. Nevertheless, the existence of such a large symmetry in the vicinity of event horizons is interesting on its own right, and the question remains as to what extent the thermodynamical properties of black holes can be captured by such a mathematical structure. In a more general context, it is interesting to investigate further the connection between the thermodynamical variables of event horizons, their infinite symmetries, and the associate conserved charges. The aim of this paper is precisely to extend to cosmological horizons the study of infinite-dimensional symmetries which have recently been successfully applied to black hole event horizons, and then apply it to investigate their thermodynamics 2 . In order to do this in an efficient way, we will introduce a special set of asymptotic boundary conditions at the horizon, which are conveniently expressed in terms of Gaussian type coordinates. This will facilitate notably the task of finding the explicit change of coordinates needed to write the physically relevant solutions in the suitable near-horizon form. The discussion will be organized as follows: In Section II, we will review the main properties of de Sitter space and of its cosmological horizon. In Section III, we will derive the infinite-dimensional symmetry that emerges in the vicinity of the cosmological horizon and the Noether charges associated to it. In Section IV, we will apply these charges to de Sitter space and show that the zero-mode of the charge reproduces the Gibbons-Hawking entropy of the cosmological horizon. We will then extend the analysis to a much more general set of solutions: We will consider the case of asymptotically de Sitter and asymptotically Taub-NUT-de Sitter spaces, including black holes, for which we also study the infinite symmetries that appear both close to their black hole horizons and to their cosmological horizons. Section V contains our conclusions.

de Sitter space
De Sitter space is the maximally symmetric solution to Einstein equations with positive cosmological constant. It can be defined by a hyperbola embedded in (4+1)-dimensional Minkowski space: On R 4,1 , we define the metric and impose the constraint where ∈ R is the de Sitter radius. This defines a (3+1)-dimensional space with isometry group SO(4, 1), i.e. the remaining Lorentz part of the Poincaré isometry group of (2), after the constraint (3) has broken the affine part of it. The space obtained by these means is non-compact and has positive, constant curvature; namely which in particular implies the Einstein space condition R µν = Λg µν with cosmological constant Λ = 3/ 2 .
In the so-called flat slicing, the metric of de Sitter space takes the form which follows from the embedding This coordinate system only covers the region X 0 + X 4 > 0 of the geometry. In these coordinates, de Sitter space has a clear cosmological interpretation: It describes an exponentially expanding isotropic, homogeneous universe with Hubble constant H ≡ −1 . The space presents a cosmological horizon at the characteristic distance H −1 . The causal structure of the space, and in particular the existence of the horizon, is clearly captured by the static coordinate system, which is defined by with ζ 2 1 + ζ 2 2 + ζ 2 3 = 1. In these coordinates the metric reads where γ AB is the metric of the unit 2-sphere (A, B = 1, 2); ζ i = ζ i (z 1 , z 2 ) with i = 1, 2, 3 give a standard embedding of S 2 . We clearly observe from (8) that the event cosmological horizon appears at the constant radiusr = H −1 . The static patch corresponds to the region r 2 ≤ H −2 ; see Figure 1. In this region, the Killing vector ∂ t is timelike. The time-dependent picture is recovered in the regionr 2 > H −2 ; in the larger limit this is given byr e τ / + ... where the ellipsis stand for terms that are subleading in the large τ limit.

Horizon boundary conditions
Now, let us analyze the region near the cosmological horizon and derive the infinitedimensional symmetries that appear there. We can study this by prescribing a set of boundary conditions that employs Gaussian type coordinates similar to those considered in the context of the membrane paradigm [16]. These coordinates, which are different from those recently considered to study black hole horizon symmetries in [17,18], have the advantage that the explicit change of coordinates from metrics in Boyer-Lindquist form can be easily written down [16]. In particular, this is well suited to study the cosmological case: Consider the following near-horizon boundary conditions where Ω AB and N A are arbitrary functions of the angular coordinates z A , A, B, C = 1, 2, and where O(r n ) stand for z A -dependent functions that near r = 0 damp off as fast as r η with and uv = −e −2τ / . These coordinates permit to see the maximal extension of the geometry, which shows the existence of two pair of causally disconnected static patches separated by a past and a future event horizon.
η ≥ n. The horizon here is located at r = 0. When accommodating de Sitter metric in the boundary conditions (9), one identifies κ with the surface gravity at the horizon. Therefore, if one considers an isolated horizon, one has the condition κ = const.
Boundary conditions (9) are similar to those introduced by Price and Thorne in [16] in the context of the membrane paradigm. These coordinates suffice to express the spacetime metric in a canonical form near the horizon, provided the horizon evolves slowly. In (9), time coordinate t coincides with the universal time coordinate, and for de Sitter it is the time of the static patch close to the future horizon. This is different from the null coordinate v of the boundary conditions considered in [17,18] and which we review in section 4.1.1. Notice that the dependence of the lapse function g tt with the surface gravity κ in (9) is quadratic, in contrast to the linear dependence g vv ∝ κ of the the boundary conditions of [17,18]; see also (26) below. The boundary conditions of [17,18] also differ from (9) in that in the latter we are not fixing any gauge condition, allowing all the components of the metric to acquire subleading falling-off contributions, as in [16]. One may still wonder about the connection between (9) and the boundary conditions of [17,18]. For instance, one might think that both coordinate systems can be related by defining the tortoise coordinate r * = r g rr dr together with the null coordinate v = t + r * . However, although in the case of static spacetimes such change of coordinates may suffice to translate one system into the other, it is not generically the case: For example, in the case of Taub-NUT spacetime, which we will discuss later, the presence of an off-diagonal term in (9) under such a naive change of coordinates translates into a component of the metric that violates 3 the boundary conditions of [17,18]. The change of coordinates relating both systems is therefore more involved. Coordinates (9) certainly 3 In the notation of [17,18], it would result in a component g Aρ O(ρ 0 ), where ρ = 0 represents the horizon location. For details about the coordinate system of [17,18]  facilitate to express several relevant physical solutions in their near-horizon form.

Asymptotic Killing vectors
The set of metrics (9) is preserved by diffeomorphisms generated by vectors ξ µ ∂ µ obeying That is, an asymptotic Killing vector ξ is defined by providing a set of three real functions . This generates the asymptotic isometry group at the horizon which act on the metric functions in the following way: together with δκ = 0. That is, functions Ω AB and N A , which are the leading terms in the small r expansion of components g AB and g At respectively, have variations that depend on the Lie derivatives with respect to the 2-dimensional vectors Y A on the constant-t slices of the horizon. In the case one considers on the horizon the conformal gauge Ω AB ≡ Ωγ AB , with Ω(z A ) being an arbitrary function of z A and γ AB being the constant curvature metric on the unit 2-sphere, the first equation of (11) reduces to with D A being the covariant derivative defined by γ AB . In this case, Y A turn out to be conformal Killing vectors on the 2-sphere, and can be represented by holomorphic functions Y z (z) and anti-holomorphic functions Yz(z).
Function T (z A ), being an arbitrary function defined on the constant-t foliation of the horizon, represents a supertranslation in the time direction; that is, a translation that a priori acts differently in each point (z 1 , z 2 ). Analogously, Y A (z B ) can be thought of as generating superrotations at the horizon.

Noether charges
The asymptotic Killing vectors ξ form an infinite-dimensional Lie algebra, which is defined by the modified Lie bracket introduced in [20]; see [18] for an explicit construction in the case of event horizons. Here, we will directly turn to the algebraic structure formed by the Noether charges associated to vector fields (10). By resorting to the covariant formalism [19] one can derive these charges, which are found to be (see [18] for details on this kind of computation): where γ ≡ det γ AB , and where the integral is over the constant-t section of the horizon. Q 0 is an arbitrary integration constant that sets the value of the charges for the background geometry. Notice that expression (13) is the integrated form of the charge. The covariant formalism rather gives the definition of the functional variation of the charge δQ[ξ], and it is not generally guaranteed that such expression is actually integrable. In fact, if one relaxes the isolated horizon condition and allows, for example, the function Ω AB to depend on time, then δQ[ξ] turns out to be non-integrable. One can thus appreciate the fact that one could arrive to the integrated expression (13), which is moreover conserved and, provided the constant-t foliations of the horizon is compact, finite.
In order to express the charge algebra formed by the charges (13), it is convenient to consider complex coordinates z,z, the conformal gauge Ω zz = Ω(z,z)γ zz , and the mode ex- As explained above, consistency demands ∂zY z = ∂ z Yz = 0. For example, local supertranslations are generated by ξ = T (m,n) z mzn ∂ t , with the standard (rigid) supertranslation in t corresponding to m = n = 0. Superrotations are generated by Y m andȲ m . If one defines the supertranslation charge and the superrotation charges one has an appropriate basis. In this basis, defining the Poisson bracket as , we find the following algebra: and While (14) are the charges associated to the local generalization of translation ∂ t , (15) are the charges associated to the local generalization of rotations ∂ z −∂z. Together, these charges form an algebra that is the semi-direct sum of two copies of the Virasoro algebra (16) and supertranslation (18). In the following section, we will evaluate these charges to the case of de Sitter cosmological horizon. This will allow us to investigate the physical meaning of (13).
In particular, the charge T (0,0) , being related to a rigid translation in t, has to correspond to an energy measure. We will see below what physical quantity such energy corresponds to, and so for different types of solutions exhibiting event horizons.

de Sitter cosmological horizon
A first interesting application of the boundary conditions (9) is the computation of charges associated to the de Sitter cosmological horizon, which is the simplest example of a more complete family of solutions we will consider later. In static coordinates, de Sitter is described by (8) and the cosmological horizon is located atr = H −1 ≡r ++ . Defining the inverse radial coordinate the metric (8) is easily accommodated in the boundary conditions (9) with functions Then, evaluating the charge (13) on de Sitter space yields vanishing charges except for the zero-mode This charge has a clear interpretation: Notice that κ satisfies the equation χ µ ∇ µ χ ν = κχ ν , with χ µ being the components of the timelike Killing vector that generates the rigid translation in the coordinate t and has norm one at the horizon. This confirms that κ = H actually coincides with the surface gravity at the horizon, and therefore it gives the Gibbons-Hawking temperature of the cosmological horizon; namely, On the other hand, the area of the horizon is given by Notice that this is valid for arbitrary horizon topology, regardless Ω AB is the metric of a closed manifold, locally flat or locally of constant curvature.
The zero-mode (21) can be rewritten in terms of these horizon quantities as Therefore, since the boundary conditions (9) are defined for fixed κ, what the zero-mode of the supertranslation charge is actually computing is the Gibbons-Hawking entropy (1) associated to the cosmological horizon (times the temperature). This simple example thus shows the physical meaning of the charge (13): the latter can be considered as a generalization of Wald's definition of entropy as a Noether charge on the horizon [21], but now in the cosmological setup. Notice that, in the case of a horizon with spherical topology, this result can be equivalently expressed as which will turn out to be a convenient form in order to compare with more general solutions we will discuss in the next sections. Notice also that T (0,0) commutes with the charges T (m,n) , Y n , andȲ n ; that is, those charges do not change the horizon entropy.

Comparison with other boundary conditions
For sake of comparison, let us now consider the horizon boundary conditions of [17,18] and apply them to the cosmological case: Close to a horizon located at ρ = 0, consider the following asymptotic conditions together with the gauge fixing conditions g vρ = g ρρ = 0 and g vρ = 1; where Ω AB , λ AB and θ A are arbitrary functions of the angular coordinates z A , A, B, C = 1, 2. Null coordinate v is the advanced time at the horizon and ρ measures the radial separation from the horizon.
In order to compare the computation performed with the boundary conditions (9), take de Sitter space in static coordinates (8) and perform the change of coordinates This yields the form where ∆g ij O(ρ 2 ) with i, j ∈ {v, A}. Then, we observe that (28) satisfies the boundary conditions (26) with These are the boundary conditions considered in [17,18], which are preserved by asymptotic Killing vectors of the form see [18] for details. This yields the conserved charges for which, again, gives the product of the entropy and the temperature of the cosmological horizon, in perfect agreement with (24). This shows that, despite the boundary conditions being different, the physical result is the same. Nevertheless, we emphasize the fact that Gaussian boundary conditions (9) have the advantage of being well suited to analyse more general asymptotically de Sitter solutions, as we explicitly show below.

Schwarzschild-de Sitter event horizons
Let us now consider black holes in de Sitter space. The metric describing such objects in Boyer-Lindquist type coordinates is given by Carter stationary solution where , ρ 2 (r, θ) =r 2 + a 2 cos 2 θ, . This solution describes a Kerr-de Sitter black hole; M represents the mass of the black hole and J = aM its angular momentum. The Carter solution above is a particular case of the Demiański-Plebański solution; the latter including, in particular, the Taub-NUT charge N , which is absent in (32)-(33). We will analyse the effects of the Taub-NUT in section 4.3.

The black hole event horizon
The rotation effects from the horizon point of view were already analyzed 4 in [17,18]. Here, as we are mainly interested in the effects on the cosmological horizon, it will be sufficient to consider the static case a = 0. In this case, one obtains the Schwarzschild-de Sitter metric where F (r) = 1 − 2GM /r − (Λ/3)r 2 , and where we recall that γ AB stands for the metric of the unit sphere. It is worth noticing that, providing M = 0, the cosmological horizon is not longer located at (3/Λ) 1/2 , but at a radius that depends on both Λ and M . In other words, there is a gravitational back-reaction on the cosmological horizon due to the black hole mass. In order to find the new position of the horizon, one needs to solve the polynomial equation This equation has three roots, given by the formulâ with x 2 = 1 − 1/(3ΛGM ) 2 . Notice that (36) actually gives three, in principle different values, each of them corresponding to the different cubic roots of 1 + x. Provided the cosmological constant takes values in the range 0 < Λ < 1/(3GM ) 2 , the three values of (36) are real and obey the hierarchyr ++ >r + > 0 >r −− ; we are also assuming M > 0. The cosmological horizon corresponds to the largest valuer ++ , whiler + is the radius of the black hole horizon. We can easily verify that in the limit Λ → 0 one recovers the Schwarzschild radiusr + = 2GM , while in the limit M → 0 the horizon isr ++ = (3/Λ) 1/2 . Now, we want to investigate whether black hole solutions (34), near their horizonsr + and r ++ , can be accommodated in the boundary conditions (9). In the case of the black hole horizon, this can actually be achieved by considering the new coordinate r defined by where the function k(r + ,r ++ ) is given by The black hole horizon is now located at r = 0, and in terms of this coordinate, the nearhorizon limit of metric (34) takes the form (9) with the functions The value κ = k(r + ,r ++ ) actually coincides with the surface gravity of the black hole event horizon. This means that (24) holds with the Hawking temperature T + = k(r + ,r ++ )/(2π) and Bekenstein-Hawking entropy S = πr 2 + /G. This reproduces the result of [17], now in presence of Λ > 0.

The cosmological horizon
In the case of the cosmological horizon, this can be done in the same way by simply replacingr + ↔r ++ in the definition (37) and in (39). More precisely, one finds where κ now coincides the surface gravity of the cosmological horizon, giving the Gibbons-Hawking temperature T ++ = |k(r ++ ,r + )|/(2π) and entropy S = πr 2 ++ /G; i.e.
In the limit M → 0 one recovers (25). Expression (41) shows that, even in presence of the black hole, where one has M as an additional, independent dimension-1 quantity, the charge T (0,0) computes the product of the temperature and the entropy, and this holds both for the black hole horizon and for the cosmological horizon.

The Nariai limit
The Nariai limit of the Schwarzcshild-de Sitter black hole solution, in whichr + tends tô r ++ , can also be studied in this context. This limit produces a degenerate horizon. As pointed out in [18], the case of degenerate horizons is peculiar in what regards to its near-horizon symmetries and thus it has to be treated separately; this can be done by using the results of reference [22].
The Nariai solution is obtained by taking in metric (34) the particular mass value In this case, we haver which for small ρ yields the metric takes the form where ρ = 0 corresponds tor =r + =r ++ , and where ∆g ij O(ρ 2 ), i, j ∈ {v, A}. This is precisely the form of the metric studied in [22] to work out the symmetries near to extremal horizons. Such asymptotic form for the metric is preserved by Killing vectors of the form [18] Evaluated on (45), this symmetry has associated a conserved charge This result is analogous to what has been obtained in [18] for extremal black holes, where the charge X (0,0) ≡ Q[X = 1, Y A = 0] was shown to reproduce the Bekenstein-Hawking entropy S = A/(4G) multiplied by a simple prefactor 1/(2π) that can be identified as the temperature T L that appears in the Kerr/CFT analysis of the extremal Frolov-Thorne vacuum; see [23,24]. (47) is the cosmological analogue of that.

Taub-NUT-de Sitter horizons
We move now to consider a more general set of solutions, namely the generalization for Λ > 0 of the Taub-Newman-Unti-Tamburino spacetime (hereafter referred to as Taub-NUT-de Sitter). In the larger limit this actually yields a different asymptotic. We will consider here the Schwarzschild-Taub-NUT-de Sitter solution, which can be thought of as a generalization of the Schwarzschild-de Sitter black hole with gravitomagnetic charge. This is also a solution to cosmological Einstein equations in vacuum. Its metric can be written as follows with where, as M , N is an arbitrary integration constant 5 . The solution reduces to Schwarzschildde Sitter solution when N = 0.

The Schwarzschild-Taub-NUT black hole horizon
In order to analyze some features of the Taub-NUT black hole solution (48), let us first consider the case of vanishing cosmological constant. Replacing Λ = 0 in (48) and performing the shift we find with U (r) = (r 2 − 2M Gr − N 2 G 2 )/(N 2 G 2 +r 2 ), where nowr ∈ R, φ ∈ [0, 2π], θ ∈ [0, π], and = ±1. The change of coordinates (50) is needed to cure the wire-type singularity that the solution (48) exhibits at θ = 0 and θ = π. As shown by Misner, this can be done by considering two patches: = +1 in the hemisphere θ < π/2 and = −1 in the hemisphere θ > π/2. Due to the periodicity of φ, (50) demands t to be periodic too, with a period Taub-NUT space can be regarded as the gravitational analog of a magnetic monopole, and the periodicity in t can be thought of as the gravitational analog of the Dirac quantization condition, with N being the gravitomagnetic charge. The space constructed in this way turns out to be spatially homogeneous and of topology R × S 3 . It has isometry SO(3), and thus is spherically symmetric. It contains closed timelike curves, though. Provided N = 0, the surfacer = 0 does not present curvature singularity, and so one can consider the entire rangê r ∈ R. The solution, however, does exhibit coordinate singularities atr =r ± , witĥ being the two solutions of U (r ± ) = 0. The radiir ± can be considered as horizons. For N = 0, the solution reduces to Schwarzschild metric, and in that caser ∈ R ≥0 withr − = 0 andr + = 2GM . For arbitrary values of M and N , the function U (r) behaves at larger like U (r) 1 − 2M G/r + O(N 2 /r 2 ), so that M represents the mass of the solution. On the other hand, N enters in the larger expansion as the charge contribution does in the Reissner-Nordström solution, although it is of purely gravitational origin. The fact that when N = 0 the off-diagonal component g φt does not vanish at infinity implies that, strictly speaking, the solution with Λ = 0 is not asymptotically flat, but it rather exhibits a special kind of asymptotic behavior characterized by the parameter N and known as asymptotically locally flat (ALF). Provided M is positive -as we will consider hereafter-the solution has two horizons, located atr + > 0 andr − < 0. Since M = (r + +r − )/(2G) ≥ 0 and N 2 = −r −r+ /G 2 ≥ 0, one can write the metric function U (r) as follows In terms of the horizon radii, the metric takes the form We are now ready to study the near-horizon region of the Taub-NUT black hole: We define the near-horizon coordinate so that r = 0 corresponds tor =r + . Then, close to the black hole event horizon we find g tt = − r 2 4r 2 + + O(r 4 ), g tr = 0, with A, B ∈ {φ, θ}. This obeys (9), with κ = 1 2r + , From (58), we find the charge which reduces to the result for the Schwarzschild solution found in [17] when N = 0, i.e. when r − = 0. More precisely, when N = 0 (59) gives one half of the ADM mass, T (0,0) = M/2, which coincides with the product T S through the Smarr formula. In the case M = 0, on the other hand, (59) gives one half of the NUT charge, T (0,0) = N/2. Notice that here we are writing (59) in such a way that no special values of the parameters M and N were assumed. This is why one should not expect T (0,0) to give in general the product of the entropy S = A/(4G) times the temperature T = U (r + )/(4π) = 1/(4πr + ) for arbitrarŷ r ± . Having this would demand a special relation between the parameters of the solution for cycles of the Euclidean Taub-NUT geometry to be uniquely defined; i.e. it would demand T = β −1 = 1/(4πN G) that yields the usual condition r + = N G, in which caser + −r − = 2N G and A = 4πr + (r + −r − ) = 8πN 2 G 2 . In fact, for these particular values (59) takes the form which gives the NUT charge and in this case coincides with the inverse of the Euclidean time period multiplied by one quarter of the area in Planck units. That is to say, the charge (59) does reproduce the expected result once the adequate relation between the parameters is considered. We do not find necessary to repeat here the standard discussion 6 on Taub-NUT and Taub-Bolt spaces and the consistency conditions for the real section of the Euclidean solutions withr + ≥ N G to be well defined. We rather prefer to write (59) in such a way that the expression of the horizon charge T (0,0) as a function ofr ± tends to the result for the Schwarzschild black hole when N goes to zero.
Notice that, in the case of the Taub-NUT (-de Sitter) black hole, unlike for the case of the Schwarzschild(-de Sitter) black hole, N A is different from zero. In the case of rotating black holes, the non-vanishing value of N A results in a contribution to the charge Q[T = 0, Y A ] that accounts for the black hole angular momentum [17]; in the case of Schwarzschild-Taub-NUT solution, in contrast, we find That is, the charge associated to the Killing vector ξ = Y A ∂ A vanishes despite N A = 0. This is due to the different values that takes in each hemisphere for the metric to be singularity free, what makes Y 0 −Ȳ 0 to be proportional to the angular integral π/2 0 dθ sin θ(cos θ − 1) + π π/2 dθ sin θ(cos θ + 1) = 0. This is analogous to what happens with the horizon superrotation charge of dyons in Einstein-Maxwell theory [22]. A way of interpreting this result is that Y 0 −Ȳ 0 computes the Misner string contribution to the angular momentum when | | = 1.

The Schwarzschild-Taub-NUT/Bolt-de Sitter cosmological horizon
Now, let us go back to cosmology: consider the case Λ > 0, for which the solution (48) represents a Schwarzschild-Taub-NUT-de Sitter black hole. We can writê wherer ++ ,r + ,r − , andr −− are the four roots of F (r) = 0. In particular, it implies . For a Λ-dependent range of the parameters M , N these four roots are real and obey the hierarchyr ++ >r + > 0 >r − >r −− , withr ++ representing the cosmological horizon andr + the black hole horizon. As in the case of the Schwarzschild-de Sitter black hole, in order to analyze the region near the cosmological horizon we define the coordinate with k(r ++ ,r + ) = Λ 6 such that r = 0 whenr =r ++ . Notice that (64) reduces to the value of k(r ++ ,r + ) for the Schwarzschild-de Sitter black hole when N = 0, i.e. whenr − = 0. This yields the charge which, as expected, agrees with the result for the cosmological horizon of the Schwarzschildde Sitter black hole (41) whenr − = 0. That is, expression (65) comes to generalize (25) and (41) which now appear as the particular casesr ++ >r + =r − = 0 andr ++ >r + =r − = 0, respectively. It also agrees, of course, with the result obtained in [17] for the Schwarzschild black holer ++ = ∞ >r + ≥r − = 0.

Other solutions
It is actually possible to extend the analysis of horizon symmetries to other solutions with regular cosmological horizons. One such example is the McVittie solution of Einstein equations [28], which describes a black hole in an expanding universe. The universe in McVittie solution is described by a fluid with an energy density σ that is constant at constanttime slices and with a pressure p that acquires an inhomogeneous contribution due to the presence of the black hole. This solution has been recently revisited in [29], where it has been noticed that, provided a positive cosmological constant is present, the McVittie solution is regular everywhere on and outside the black hole event horizon and away from the big bang singularity. The metric can be casted in the form where F (r) = 1−2M G/r, with M being the mass of the black hole, and where H(t) obeys the As said, in the case H(t = ∞) > 0 the McVittie metric describes a black hole whose horizon is regular and is embedded in an FLRW spacetime. Its causal structure is quite interesting, as it has a spacelike and inhomogeneous big bang singularity in the remote past.
It also has a spacelike future infinity at large radial distancer and late time t. The surfacê r = 2M G is a curvature singularity at finite t. The solution also contains a null apparent horizon atr =r − and t = ∞, wherer − is the smaller positive root of the cubic, timedependent equation F (r ± ) = H 2 (t)r 2 ± . Provided H(t = ∞) > 0, the solution also exhibits a cosmological event horizon at the pointr =r + (t = ∞) and t = ∞; atr =r ± (t) and t finite, the solution exhibits an apparent horizon. Analyzing the near-horizon symmetries and the associated conserved charges for this type of solutions with the techniques discussed here 7 is also possible. The result is consistent with (65). Other solutions with apparent horizons can also be studied in the same way.

Conclusions
The results of this paper provide an extension to the cosmological setup of previous studies [17,18] of infinite-dimensional symmetries in the vicinity of event horizons: We have shown that, in its proximity, cosmological de Sitter horizons exhibit an infinite-dimensional symmetry. The latter is the asymptotic isometry group that preserves a suitable set of boundary conditions adapted to study cosmological horizons, boundary conditions which happen to be similar to those considered by Price and Thorne in the context of the membrane paradigm. The corresponding isometry group is generated by an infinite-dimensional algebra, which is identified as the semi-direct sum of two copies of the Virasoro algebra (with vanishing central charge) and supertranslations. This symmetry has an associated Noether charge, which we have derived using the Barnich-Brandt covariant formalism. This charge is conserved, integrable, and finite, and we discussed its physical meaning by evaluating it on de Sitter space. The zero-mode of this charge was found to give the Gibbons-Hawking entropy of the cosmological horizon. We then studied a larger class of geometries, considering the case of Schwarzschild-Taub-NUT-de Sitter spaces: we found that their associated Noether charge, in region close to both the black hole horizon and the cosmological horizon, computes the entropy of these horizons.
The fact that symmetries at the cosmological horizon yield conserved charges associated to their thermodynamics is interesting on its own right. A much more interesting, although speculative question is whether these symmetries somehow encode information hidden behind the cosmological horizon. As pointed out in [30], it is an interesting conceptual question whether the matter passing out through the cosmological horizon will leave behind some imprint in the same way that information that fell through a black hole reappears in a scrambled form in the Hawking radiation. It would be interesting to study whether the infinite-dimensional symmetries discussed here have something to say in this regard.
Another interesting application of the near-horizon symmetries of cosmological horizons is the connection of such symmetries with the membrane paradigm [31] and, more specifically, with the Carroll ultra-relativistic symmetries that black hole horizons exhibit [32,33]. To this respect, it is worth emphasizing that the boundary conditions we considered here are well suited to the membrane paradigm formalism [16].