Large Gauge Symmetries of an Asymptotically de Sitter Horizon: An Extended First Law of Thermodynamics
Fatemeh Mahdieh Hossein Shojaie








f_mahdieh@sbu.ac.ir
 Department of Physics, Shahid Beheshti University, G.C., Evin, Tehran 1983969411, Iran



h-shojaie@sbu.ac.ir
Department of Physics, Shahid Beheshti University, G.C., Evin, Tehran 1983969411, Iran




In this paper, we show that a universe with a dynamical cosmological constant approaching pure de Sitter at timelike infinity, enjoys an infinite dimensional symmetry group at its horizon. This group is larger than the usual SO(4,1) of pure de Sitter. The charges associated with the asymptotic symmetry generators are non-integrable, and we demonstrate that they promote an extended version of the first law of thermodynamics. This contains four pairs of conjugate variables. The pair (Θ,Λ) corresponds to the change in the cosmological constant and its conjugate volume Θ. The contribution of the surface tension of the horizon and its conjugate parameter surface area make a pair (σ, A). The usual conjugate variables (T, S), (Ω, J) and a term ∂_v δ S corresponding to entropy production, are included. In addition, this extended first law describes the non-conservative behaviour of the asymptotic charges in non-equilibrium.



Large Gauge Symmetries of an Asymptotically de Sitter Horizon:
 An Extended First Law of Thermodynamics 
    Hossein Shojaie
    December 30, 2023
========================================================================================================







§ INTRODUCTION

The extra symmetries of asymptotically flat spacetimes were first discussed by Bondi, van der Burg, Metzner and Sachs in the early 60's <cit.>. In fact they were discovered in a “fortunate" failure of trying to obtain the Poincaré group as the exact symmetry group of asymptotically flat spacetimes.

The mathematical machinery for asymptotic symmetries has been studied further and developed extensively by different research groups <cit.> and references therein. But the physical concept of what these extra symmetries describe has received additional attention after Strominger et.al. <cit.> related three notions from completely different areas together. They introduced a triangle with asymptotic symmetries being one of the corners, and the other two are soft theorems and the memory effect. This area of research has thereafter developed in work done by various groups at null infinity, spatial infinity <cit.> and the black hole horizon <cit.>. The asymptotic limit is usually a boundary in spacetime, and these symmetries are generated by vector fields that are asymptotically Killing on the boundary. In asymptotically flat spacetimes, null infinity or spacelike infinity is taken to be the boundary, and asymptotically Killing vector fields on this boundary are derived. These generators can be written as an expansion of a parameter that approaches zero in the asymptotic region. Accordingly, ℒ_ξ g is also of order of this parameter, and not exactly zero, but asymptotically vanishes. 

This is the case in asymptotically flat spacetimes at their boundaries with the parameter being 1/r. This gives ℒ_ξg_μν = 𝒪(1/r). Likewise, if the horizon of a black hole is taken to be the boundary, the same expansion can be written in terms of a radial parameter ρ, that describes the horizon through the equation ρ = 0.

In asymptotically flat spacetimes, the BMS group is the semi-direct product of the Poincaré group and super translations, i.e. BMS = Poincare× ST. This group is manifest when studying the symmetries of an asymptotically flat spacetime at null infinity. Here, the supertranslations have been derived by Bondi, Metzner and Sachs <cit.>. This group can additionally be extended to give Virasoro-like generators referred to as superrotations <cit.>. In the case of spacelike infinity in an asymptotically flat spacetime, it has been shown that in the canonical formalism, by considering suitable matching conditions, supertranslations can be recovered at spacelike  infinity of a dynamical asymptotically flat spacetime<cit.>. In this context, electromagnetism has also been included <cit.>. The study of asymptotic symmetries at spatial infinity was worked out covariantly in <cit.>, where the relation between the symmetries at past and future null infinity was discussed and again electromagnetism was included  <cit.>. 

The physical concept of the asymptotic symmetries can be understood as the transition between different asymptotically flat states of an evaporating black hole. In other words, these symmetries characterize the vacuum transitions in the evaporation of a black hole. A boosted radiating black hole <cit.> and as shown recently, a rotating black hole <cit.>, are other examples of spacetimes possessing these extra symmetries.  

One should mention that, due to the lack of a timelike Killing vector field in the geometry of a dynamical spacetime, there are difficulties in defining the entropy and surface gravity, which result in issues with assigning a temperature to the black hole horizon. There are conserved charges associated with the asymptotic symmetries of dynamical spacetimes, that are given by Noether's theorem. Using the conserved charges, an analogue of the first law of thermodynamics can be obtained <cit.>.

To highlight our choice of setup for a dynamical system, we note that maximally symmetric spacetimes are of interest to theoretical physicists as well as cosmologists. Regarding the n-dimensional pseudo-Reimannian manifolds, these are Minkowski, de Sitter and anti-de Sitter spacetimes, with zero, positive and negative curvatures respectively. These spacetimes can be seen as the “ground states" of general relativity <cit.>. On the other hand, observations provide that the fate of our Universe is similar to de Sitter spacetime. According to the current values of the density of matter and cosmological constant, it seems that the locally gravitationally bound group of galaxies, i.e. cluster of galaxies, are the typical structures that will remain bound as the Universe expands. Therefore, one can take the Universe to be undergoing a process that eventually tends to pure de Sitter. 

In this work we aim to find the extra symmetries of a dynamical de Sitter-like spacetime at its horizon in order to study thermodynamic aspects of this dynamical system.[The asymptotic symmetries of de Sitter spacetime have been studied in <cit.> and references therein. However, in this paper we consider the de Sitter horizon as the asymptotic boundary, whereas the previous work mentioned has mainly studied the asymptotic symmetries at null infinity.] It is worth noting that extra symmetries on a black hole horizon have been studied in  <cit.>. Here we consider a different geometry and constraints in order to retrieve pure de Sitter at timelike infinity. Furthermore, we do not dismiss the time dependency of the metric components to ensure we have a dynamical system settling in pure de Sitter.  

In a dynamical spacetime such as Friedmann-Lemaître-Robertson-Walker (FLRW), there is no timelike Killing vector field in general and the surface gravity cannot be defined in the usual manner. However, there are several ways to look at the surface gravity, in which not all methods reproduce the expected temperature in the static limit <cit.>. Therefore there is not yet a conclusion to a preferred or correct approach to defining the surface gravity in the dynamical case, and therefore we lack a description for the thermodynamics of such a spacetime.
    
The static patch of pure de-Sitter has a temperature of 2π/κ <cit.>, which is the analogue of the black hole temperature. Since all dynamical black holes ultimately settle in a stationary state, in the language of equilibrium and non-equilibrium thermodynamics, a stationary black hole with temperature T can be considered as the analogue of an equilibrium state. Similarly, the static patch of pure de Sitter spacetime with temperature T is also the analogue of an equilibrium state in an inflating universe. To present an analogue of an extended first law of thermodynamics in a de Sitter-like dynamical spacetime, we turn to non-equilibrium thermodynamics. Here, by a de Sitter-like spacetime, we mean a geometry that tends to pure de Sitter at timelike infinity. 

We show that the near-horizon geometry of a de Sitter-like spacetime admits asymptotic symmetries at its horizon. In other words, this spacetime possesses a larger symmetry group than SO(4,1) of pure de Sitter. The conserved charges related to these symmetries can be obtained by means of Noether's theorem and constructing the symplectic form <cit.>. In a similar manner, for a non-stationary perturbation to the near-horizon geometry of a de Sitter-like spacetime generated by its asymptotic symmetries, the variation of the conserved charges can be calculated. This variation  can be interpreted as an extended first law of thermodynamics that provides an insight into the non-conservation and apparent time-reversal symmetry breaking of this dynamical system. 

This manuscript is organized as follows. In Sect. <ref>, we give a brief description of the general definition of asymptotic symmetries, then proceed by portraying how this definition manifests in general relativity. The remainder of this section is dedicated to specify the near-horizon geometry that we wish to study, and then to derive the generators of the asymptotic symmetries of this geometry. 

In Sect. <ref>, the equations of motion for this geometry are given where the timelike boundary conditions have been taken into account in the solutions. Sect. <ref> consists of the derivation of the variation in surface charge. This calculation is presented as an analogue of an extended first law of thermodynamics in Sect. <ref>. A conclusion and some remarks for possible future work are provided in <ref>.





§ ASYMPTOTIC SYMMETRIES

Considering a gauge theory on a principal fiber bundle (P, ℳ,π), a gauge transformation is defined as an automorphism ϕ : P→ P such that the following diagram commutes


    
        P [r]^ϕ [rd]_π     P [d]^π 

                   M 

A gauge is therefore fixed by choosing a section S: M→ P of this fiber bundle, and a gauge transformation maps sections to sections. 

One can define a large gauge transformation as a gauge transformation that is not homotopic to the identity map. In terms of physical systems that we deal with, trivial gauge transformations do not change the physical states of a system, whereas large gauge transformations in fact map physical states to different ones. In other words, while a trivial gauge transformation can not affect the physical parameters of the system, a large gauge transformation can influence the physical state.

The gauge transformations in general relativity are diffeomorphisms, and the physical quantities are defined through geometry. A trivial gauge transformation on a manifold (ℳ,g) can be constructed by small diffeomorphisms, generated by exact Killing vector fields of the geometry. On the other hand, large gauge transformations in a diffeomorphism invariant theory such as general relativity, are defined by the quotient space of all possible gauge transformations modulo the trivial gauge transformations of the theory. If the group of trivial gauge transformations is a normal subgroup of all possible transformations, then its quotient space is the group of large gauge transformations. 

To proceed in finding large gauge transformations in a spacetime with geometry g_μν, we first fix a coordinate system to discard the gauge freedom related to trivial gauge transformations in GR. Thereafter, to find the quotient space, we need to find all remaining diffeomorphisms that leave the fixed gauge untouched.



 §.§ Gauge fixing

In the case of the de Sitter horizon defined by the hypersurface ρ = 0, we can construct a near-horizon geometry that takes into account non-stationary perturbations of order ρ in the metric field. To define the perturbation parameter, it should be noted that the geometry will not change at the horizon by adding terms of order ρ and higher. Hence, we fix the gauge by defining a near-horizon geometry using the definition in <cit.> defined for any null hypersurface. For the sake of maintaining some level of self-consistency, in the following, we give a brief description of the construction of this coordinate system[See 
<cit.> for a rigorous description of this construction and further examples]. 



 §.§ The near-horizon geometry as the gauge


Let (ℳ,g) be an 4-dimensional pseudo-Riemannian manifold. By constructing a suitable coordinate system near the horizon, we wish to make the asymptotic symmetries on the horizon manifest.

According to the theorem stated in <cit.>, the geometry near a smooth null hypersurface ℋ (⊂ℳ) can always be described by the metric g as follows

    g =   ρϕ(v,ρ,x)dv^2 + 2dv dρ + 2ρ h_A(v,ρ,x)dvdx^A 
       + g_AB(v,ρ,x)dx^Adx^B

where v is the time coordinate, ρ is the radial coordinate, x=x^A (A=θ,ϕ) are the angular coordinates on the two-sphere and ℋ is given by the equation ρ = 0.

Specifically, according to the theorem, this geometry is formed by considering a coordinate system x^A on a two-sphere 𝒮(⊂ℋ). Then a coordinate v can be constructed by extending the coordinates x^A to a neighbourhood ℋ'(⊂ℋ) of 𝒮, by demanding ∇_χχ = 0 for some vector field χ = χ^μ∂_μ with initial value χ|_𝒮, and taking
χ(x^A) = χ^μ∂_μ x^A = 0. Thus, χ |_ℋ'≡∂_v is a null generator on ℋ'(⊂ℋ), which means g(χ,χ)|_ℋ' = g_vv|_ℋ' = 0. Henceforth, it takes the form g_vv≡ρϕ(v,ρ,x), where ϕ is a function of coordinates and must be finite on the horizon. So far, this gives a coordinate system (v,x) on the horizon. Furthermore, one can deduce g(χ,∂_A)|_ℋ'=g_vA|_ℋ' = 0 from the fact that χ |_ℋ' is null and thus normal to ℋ'(⊂ℋ). On this account, g_vA≡ρ h_A(v,ρ,x), where h_A is also finite at the horizon. 

To extend this coordinate system, we can consider a neighborhood 𝒰(⊂ℳ) of ℋ'  and demand ∇_ΥΥ=0 with initial value Υ|_ℋ'. Hence, the coordinate system (v,x) on ℋ' is extended to (v,ρ,x) on 𝒰, by solving Υ(v) = Υ(x^A) = 0, and defining ρ to be the solution of the equation Υ(ρ)=1 with initial value ρ=0 on ℋ', which is the equation defining the null surface ℋ. As a result of this construction, we have Υ = ∂_ρ, and the vρ component of the metric is fixed to g_vρ = 1. 

For our purpose we work with the following form of the near-horizon geometry to describe the geometry near the de Sitter horizon to first order in ρ,

    ds^2 =    2ρκ(v,x) dv^2 + 2dvdρ + 2ρ h_A(v,x)dvdx^A 
       + g_AB(v,ρ,x)dx^Adx^B,

where 

    g_AB(v,ρ,x) = (ρ + l(v))^2 q_AB(x),

and q_AB(x) is the metric on the unit two-sphere. We later constrain the variable fields κ, h_A and l through the equations of motion, to recover pure de Sitter in the limit v→∞. 

This construction provides four equations that fix the diffeomorphism gauge freedom appropriate for finding the form of the asymptotic symmetry generators. 





 §.§ Asymptotic Symmetry generators

According to this geometry, the gauge equations are by definition

    δ g_ρρ = 0, δ g_ρ A = 0, δ g_ρ v = 0

with boundary conditions

    δ g_vv = 𝒪(ρ^2),δ g_vA = 𝒪(ρ^2),δ g_AB = 𝒪(ρ^3).


Solving these equations, one can find the form of the generators of the transformations that leave the gauge fixed. They are consequently the generators of the asymptotic symmetry group. This yields, 

    0    = δ g_ρρ = ℒ_ξg_ρρ
       ⟹ξ^v = f(v,x) 
    
            0    = δ g_ρ A = ℒ_ξg_ρ A
       ⟹ξ^A = Y^A(v,x)- ∂_B f(v,x)∫_0^ρg^ABdρ' 
    
            0    = δ g_ρ v = ℒ_ξg_ρ v
       ⟹ξ^ρ = F(v,x) - ρ∂_vf + ∂_Bf∫_0^ρg^ABg_vAdρ',

where f, Y^A and F are arbitrary functions of the coordinates (v,x). 




 §.§ Dynamics of Generators


The variation of the metric field now gives the non-stationary perturbations in the geometry. For the non-zero components of the metric, we thus have 

    δ g_vv = 2ρδκ
       δ g_vA = ρδ h_A
       δ g_AB = 2(ρ + l)δ l q_AB.

Using the boundary conditions and the field variations through δ g_μν = ℒ_ξg_μν, leads to 

    𝒪(ρ^0)_δ g_vv = 2∂_v F + 2κ F = 0.

Noting that the asymptotic symmetry generators are taken to be fixed to first order, they do not depend on the dynamical fields [In <cit.>, they mention this assumption as the state independence of the boundary conditions.]. On this account, we take F = 0, and as a result, 

    𝒪(ρ^0)_δ g_vA = Fh_A + ∂_A F + l^2 q_AB∂_vY^B = 0

gives ∂_v Y^A = 0. To next order, one has 

    𝒪(ρ)_δ g_vv   = 2ρ[ ∂_v(fκ) + h_A∂_vY^A -∂_v^2f ] 
       ⟹δκ =  ∂_v (fκ) -∂_v^2f.

From this equation and the boundary conditions  we can write 

    f(v,x) = T(x) + Z(v,x),

where T(x) is an arbitrary function on the two-sphere at the horizon. This generates a supertranslation-like transformation along the horizon, that causes the transition between different states of the dynamical process. Also, 

    𝒪(ρ)_δ g_vA   = ρ[ f∂_v h_A + Y^B∂_B h_A + h_B∂_AY^B - ∂_A∂_v f ] 
       = ρδ h_A
    ⟹   δ h_A = f∂_v h_A + Y^B∂_B h_A + h_B∂_AY^B - ∂_A∂_v f.

The last components yield, 

    𝒪(ρ^0)_δ g_AB   = 2fl∂_v l q_AB + l^2 (D_AY_B + D_BY_A) = 2lδ l q_AB
       ⟹δ l = f∂_v l + 1/2lψ,

where ψ≡ D_AY^A. 



§ EQUATIONS OF MOTION


In a universe dominated by a cosmological constant, we consider a mass locally centered within its cosmological horizon. The equation of state near the horizon, can be approximated by: 

    w=p_tot/ϵ_tot   =(-1+1/3ϵ_rad/ϵ_Λ)(1+ϵ_rad/ϵ_Λ)^-1
       ≈-1+4/3ϵ_rad/ϵ_Λ>-1,

where ϵ and p are the energy and pressure densities respectively. In equation (<ref>) and by the indices “Λ" and “rad", we label the cosmological constant and the radiation emitted out from the mass inside the bulk. Here, the radiation density is not comparable in value with the vacuum energy density, i.e. ϵ_rad/ϵ_Λ≪ 1. Since this ratio is also decreasing as the constituents of the mass in the bulk gradually settle down in their ground states, one can approximate T_μν at the horizon defined by ρ = 0 in (<ref>), as 

    T_μν|_ρ=0≈ -ϵ_Λ(1-1/3ϵ_rad/ϵ_Λ)g_μν.

As a consequence, by introducing the function λ(v,x) such that

    lim_v→∞λ(v,x) = Λ_0 = const,

the Einstein equations on the horizon can be written as

    G_μν + λ(v,x)g_μν≈ 0.

The function λ(v,x) accounts for the effects caused by legitimately small amounts of matter passing through the horizon, radiation, or any other type of perturbation. The constraint on λ(v,x) is to ensure pure de Sitter is recovered in timelike infinity[Although we have considered the fate of the Universe to settle in pure de Sitter at timelike infinity, other scenarios also seem to be possible if one loosens the timelike boundary conditions.].

Continuing with the Einstein equations, the vv component at the horizon is, 

    G_vv|_ρ = 0 = -2κ∂_v l/l - 2∂_v^2 l/l = -λ(v,x)g_vv|_ρ = 0 = 0.

Since we have taken l to be only a function of v, this component of the equations results in ∂_Aκ = 0. For the ρ v component we solve G_ρ v = -λ(v,x)g_ρ v at the horizon to find

    G_ρ v|_ρ = 0 =    -1/l^2(1 - 2 ∂_v l + 2l κ - 1/4h_Ah^A + 1/2q^AB∂_Ah_B 
       + 1/4 q^ABh^C∂_Cq_AB) = -λ(v,x).

As mentioned above, in order to retrieve pure de Sitter as v→∞, we wish lim_v→∞λ(v,x)= 3/l_0^2, and lim_v→∞ l(v) = l_0, where l_0 is the cosmological horizon in pure de Sitter. This gives lim_v→∞κ = 1/l_0. We also have,

    lim_v→∞[ -h_Ah^A + 2q^AB∂_Ah_B + q^ABh^C∂_Cq_AB] = 0.

The vA components of the EOM  at the horizon reduce to

    G_vA = ∂_A κ - ∂_vl/lh_A - 1/2∂_vh_A = -λ(v,x)g_vA|_ρ = 0 = 0.

Since ∂_Aκ = 0 and l=l(v), this equation gives the v functionality of h_A(v,x), 

    ∂_vl/lh_A = -1/2∂_v h_A
       ⟹ h_A(v,x) = a_A(x)b(v) = a_A(x)/l(v)^2,

where separation of variables is applied to h_A(v,x)=a_A(x)b(v). In turn, the AB components at the horizon are

    G_AB|_ρ = 0 =    -1/2h_Ah_B+3/4q^CDh_Ch_Dq_AB
       + 1/2(∂_Ah_B + ∂_Bh_A) - q^CD∂_Ch_Dq_AB
       - 2lκ q_AB - 1/2q^CDh_B∂_Aq_CD
    
            =    -λ(v,x) g_AB|_ρ=0
    
            =    - λ(v,x) l^2 q_AB.

From relation (<ref>) we have, 

    -h_Ah^A + 1/2q^AB(∂_Ah_B+∂_Bh_A + q^CDh_B∂_Aq_CD) + 4lκ
        = 2λ(v,x) l^2.

As a side note, the reason that the coefficient of lκ in (<ref>) does not match to give lim_v→∞κ = 1/l_0, is that we have not considered the higher order in the g_vv. If we take g_vv = 2ρκ_1 + ρ^2κ_2, the G_AB equations will give 

    G_AB =    -1/2h_Ah_B+3/4q^CDh_Ch_Dq_AB + 1/2(∂_Ah_B+∂_Bh_A) 
       - q^CD∂_Ch_Dq_AB - 2lκ_1 q_AB - l^2κ_2 q_AB
       - 1/2q^CDh_B∂_Aq_CD,

which results in lim_v→∞κ_1 = 1/l_0 and lim_v→∞κ_2 = 1/l_0^2 and we retrieve pure de Sitter as desired. 



§ SURFACE CHARGES

Having found the generators of the asymptotic symmetries, we can calculate the surface charge using <cit.>

    δ Q_ξ[h,g] =    - 1/16π∫_S(d^2x)_μν√(-g) [D^ν(ξ^μh) + D_σ(h^μσξ^ν) 
       + D^μ(h^νσξ_σ) + 3/2hD^μξ^ν
       + 3/2h^σμD^νξ_σ + 3/2h^νσD_σξ^μ - (μ↔ν)],

where h_μν≡δ g_μν, D_μ is the covariant derivative of the metric, and

    (d^2x)_μν = 1/4ϵ_μνσδdx^σ∧ dx^δ, ϵ_vρθϕ=1.

Note that δ indicates that the integrability of δ Q is not determined. 

In our case, the symmetry generators are given by (<ref>), and the corresponding change in the asymptotic charges is given accordingly by

    δQ_ξ[h,g] =
               +1/8π∂_v∫_ℋ dΩ^2 √(q) fδ l^2 + 1/8π∫_ℋ dΩ^2 √(q)f∂_vlδ l 
       + 1/16π∫_ℋ dΩ^2 √(q) Y^Aδ (l^2h_A) 
        - 1/4π∫_ℋ dΩ^2 √(q)δ l^2∂_v f + 1/8π∫_ℋ dΩ^2 √(q)fκδ l^2.


A related point to consider is that in the limit v→∞, κ_0 = 1/l_0, ∂_vl=0, and for ξ_0 = ∂_v (f=1,Y^A=0), the expression (<ref>) leads to

    δ Q_0 = 1/8πδ∫_ℋ dΩ^2 √(q)(l_0^2κ_0)
       ⟹ Q_0 = l_0/2.

Considering T_0 ≡ (2π l_0)^-1 and S_0=A_ℋ/4 = π l_0^2, Q_0 = S_0T_0 can be interpreted as the heat observed by the static patch observer in pure de Sitter  <cit.>. In other words, the static patch observer assigns an apparent temperature to its horizon as a result of its inability to survey the entire spacetime. 





§ EXTENDED FIRST LAW OF THERMODYNAMICS FOR A DYNAMICAL INFLATING UNIVERSE


It is well known that there is an analogy between equilibrium thermodynamics and gravity in black holes <cit.>. In non-equilibrium thermodynamic cases, the end state of dynamical gravitational systems such as mergers or radiating black holes, can be either a stationary or Minkowski spacetime. In other words, these dynamical systems decay to settle into an analogue of an equilibrium configuration. Undeniably, thermodynamic properties of these systems cannot be described completely through equilibrium thermodynamics. 

For instance, one can refer to <cit.> for non-equilibrium thermodynamic properties of an asymptotically flat spacetime containing gravitational radiation. There, gravitational waves have been shown to be an entropy producing mechanism of a viscous dynamical screen.   

In an irreversible thermodynamic process, we have 

    dS >dQ/T,

which is due to entropy production in the system. In addition, dS/dt≠ 0, which makes any process with a gravitational wave production mechanism, irreversible and non-equilibrium. Here, we show that the relation (<ref>) can be taken as an analogue of an extended first law of thermodynamics, hence indicating the dynamics of our setup to include entropy production. 

Proceeding with the first term in (<ref>), we find that 

    1/2π∂_v δ[1/4∫_ℋdΩ^2√(q)fl^2 ] = 1/2π∂_v δ S

can be understood as a consequence of the non-conservation of entropy in a dynamical system in an irreversible process. It has a close resemblance to dS/dt which is not zero in non-equilibrium thermodynamics. The existence of f again is due to the infinite dimensional symmetry in the system. 

Taking Λ(v)=3/l^2(v) as a state parameter, the second term

    1/8π∫_ℋdΩ^2 √(q)f∂_vlδ l    = -1/8π∫_ℋdΩ^2 √(q)f l^3/6∂_vlδΛ
       =-∫_ℋdΩ^2 √(q)fl^3/6∂_vl δρ_Λ
       ≡Θδρ_Λ

resembles Θ dΛ in <cit.>[A similar analysis has been done for anti-de Sitter in <cit.>, where the negative cosmological constant is considered as a thermodynamic variable]. Here Θ is the generalized force conjugate to Λ, and corresponds dimensionally to a generalized volume. Specifically, Θ is the generalized time-dependent volume inside the dynamical cosmological horizon up to a factor. Since ρ_Λ≡Λ / 8π, the term Θδρ_Λ is the contribution of vacuum energy production which is consistent with the increase of vacuum energy in a dynamically inflating universe. As long as ∂_v l>0, this term is the contribution of the increase in the volume enclosed by the horizon, due to the time dependent cosmological constant. In the limit v→∞, this term vanishes as expected for pure de Sitter spacetime with a constant cosmological constant. 

The term

    1/16π∫_ℋdΩ^2 √(q)Y^Aδ(l^2h_A) ≡Ωδ J

is the contribution of angular momentum to the thermodynamic potential.

The fourth term can be rewritten as

    -∫_ℋdΩ^2 √(q)∂_vf/4π f fδ l^2 ≡σδ A

and can be thought of as an analogue of the contribution of surface tension when considering a screen or bubble. It has been shown that this contribution is σδ A, where σ is the surface tension and A is the area of the horizon <cit.>. In equation (<ref>), the surface tension is defined as σ_ℋ≡ -∂_v f / 4π f, and can be understood as a negative pressure. This term acts as the work done by a repulsive force, expanding the horizon. Our result is different from that of <cit.> since we have assumed T_μν = 0 near the horizon, while their definition of surface tension is proportional to T^r_r|_r_c, where r_c is the radius of the cosmological horizon in their coordinates. Accordingly, that contribution is not included in our case. The reason we have not combined this term with (<ref>) is based on the assumption of taking the surface gravity to be only a function of v. 

Finally, in the term 

    κ/2π1/4∫_ℋdΩ^2√(q)f δ l^2 ≡ Tδ S,

we interpret κ/2π as the temperature of the dynamical inflating spacetime and the expression ∫ dΩ^2√(q)fδ l^2 can be taken as the change in the area. Consequently, this term can be interpreted as the equivalent of Tδ S. The coefficient f(v,x) is defined by (<ref>), and this can be a proposal for the dynamical entropy and dynamical surface gravity. Note that κ is restricted by the equations of motion and is not an arbitrary function. It is also worth pointing out, this term contributes to non-integrability due to the fact that δκ≠ 0. This issue will become crucial in the study of charge algebras, which will be studied in future work and is not a subject in this work.

In a gravitational system such as a merger, or a ringing black hole, entropy production is implemented by gravitational waves. Here we have found, the stretching of spacetime itself contributes to entropy production, when dealing with an otherwise empty spacetime with positive curvature due to a positive cosmological constant. 

The extended first law of thermodynamics for such a system is, therefore,

    δ Q = 1/2π∂_v δ S + Θδρ_Λ + Ωδ J + σδ A + Tδ S .


As a result, what is found from (<ref>) and consequently (<ref>), is that the dynamical spacetime of our setup is an analogue of a non-equilibrium process.  




§ CONCLUSION

In this paper, we propose a process described by the equations of motion for the near-horizon geometry of a universe with a dynamical cosmological constant. We consider this system to approach pure de Sitter at timelike infinity, and find that such a system enjoys a larger symmetry than the usual SO(4,1) symmetry of pure de Sitter given by (<ref>). This larger symmetry consists of supertranslation-like generators on the horizon. The supertranslations generate the dynamical transition between different states of an irreversible process. By calculating the change in the charges associated to these symmetry generators, we find an extended version of the first law of thermodynamics, expressed in (<ref>). This displays the non-conservation of the charges in a non-equilibrium process. The apparently non-conservation of the charges is due to the dynamical aspect we have assumed and the incapability of the observer to examine the degrees of freedom in the entire spacetime beyond its horizon. 




§ ACKNOWLEDGEMENT

Part of this project was done when F.M. was a fellow at Center for the Fundamental Laws of Nature, Harvard University in 2018. F.M. would like to thank Andrew Strominger for providing this incredible opportunity.


Ferreira16



