## Symmetries, charges and conservation laws at causal diamonds in general relativity

Venkatesa Chandrasekaran 1, ∗ and Kartik Prabhu 2, †

1 Center for Theoretical Physics and Department of Physics, University of California, Berkeley, CA 94720

2 Cornell Laboratory for Accelerator-based Sciences and Education (CLASSE), Cornell University, Ithaca, NY 14853

## Abstract

We study the covariant phase space of vacuum general relativity at the null boundary of causal diamonds. The past and future components of such a null boundary each have an infinite-dimensional symmetry algebra consisting of diffeomorphisms of the 2-sphere and boost supertranslations corresponding to angle-dependent rescalings of affine parameter along the null generators. Associated to these symmetries are charges and fluxes obtained from the covariant phase space formalism using the prescription of Wald and Zoupas. By analyzing the behavior of the spacetime metric near the corners of the causal diamond, we show that the fluxes are also Hamiltonian generators of the symmetries on phase space. In particular, the supertranslation fluxes yield an infinite family of boost Hamiltonians acting on the gravitational data of causal diamonds. We show that the smoothness of the vector fields representing such symmetries at the bifurcation edge of the causal diamond implies suitable matching conditions between the symmetries on the past and future components of the null boundary. Similarly, the smoothness of the spacetime metric implies that the fluxes of all such symmetries are conserved between the past and future components of the null boundary. This establishes an infinite set of conservation laws for finite subregions in gravity analogous to those at null infinity. We also show that the symmetry algebra at the causal diamond has a non-trivial center corresponding to constant boosts. The central charges associated to these constant boosts are proportional to the area of the bifurcation edge, for any causal diamond, in analogy with the Wald entropy formula.

∗ ven\_chandrasekaran@berkeley.edu

† kartikprabhu@cornell.edu

## CONTENTS

| 1. Introduction                                                |   2 |
|----------------------------------------------------------------|-----|
| Notation and conventions                                       |   4 |
| 2. Null boundary symmetries and charges                        |   5 |
| 1. Universal structure and symmetries on null boundaries       |   5 |
| 2. Gaussian null coordinates                                   |   6 |
| 3. Null boundary symmetries, charges and fluxes in GNC         |   8 |
| 3. Causal Diamonds                                             |  12 |
| 1. Reduced symmetry algebra g CD at causal diamonds            |  14 |
| 2. Behavior of the fields and charges near the corners         |  15 |
| 4. Conservation laws at causal diamonds                        |  18 |
| 5. Central charges and area of the bifurcation edge            |  20 |
| 6. Discussion                                                  |  22 |
| Acknowledgements                                               |  23 |
| A. Covariant phase space formalism and the Wald-Zoupas charges |  23 |
| B. Structure of g CD as a central extension                    |  26 |
| References                                                     |  27 |

## 1. INTRODUCTION

Understanding the description of finite subsystems in diffeomorphism invariant theories is an important problem in both classical and quantum gravity. Over the years there has been a lot of progress in uncovering various aspects of gravitational subsystems by studying the covariant phase space formalism in the presence of boundaries [1-10]. The presence of a boundary promotes a subset of the boundary preserving diffeomorphisms to symmetries of the covariant phase space. These boundary symmetries then result in non-trivial boundary charges which can be thought of as capturing aspects of the degrees of freedom contained within the subregion.

Null boundaries are particularly important as they play a fundamental role in gravitational thermodynamics [11-14], as well as in holography and quantum gravity [15-17]. Moreover, it

was recently conjectured that the symmetries and charges at stationary event horizons are relevant to the black hole information problem [18-20]. A particularly important class of null surfaces are the boundaries of causal diamonds, which are fundamental to the description of gravitational subregions. Black hole thermodynamics [2, 11, 21, 22] and entanglement entropy in AdS/CFT [23] have demonstrated that geometric properties of causal horizons are deeply related to the thermodynamic and statistical properties of spacetime. This strongly suggests that these deep connections generalize to any causal diamond in any spacetime. While there have been important insights in this direction [24-26], a complete understanding for arbitrary gravitational subregions remains elusive. A potential avenue of progress lies in the covariant phase space formalism applied to gravity at the boundaries of causal diamonds.

In [27] the covariant phase space of general relativity at null boundaries was studied in detail. There it was shown that there exists an infinite-dimensional symmetry algebra for general relativity at all null boundaries, including non-stationary event horizons, in any spacetime. For null boundaries of the form N = Z × R , where Z is the space of null generators, the algebra takes the form diff( Z ) glyph[multicloseleft] s with the 'supertranslation' subalgebra s consisting of angle-dependent translations and rescalings of affine parameter along the null surface. The charges (at cross-sections of the null surface) and fluxes associated to the symmetries were computed from the covariant phase space formalism using the prescription given by Wald and Zoupas [3].

In this paper we use the results of [27] at the boundaries of causal diamonds. We consider causal diamonds obtained from the intersection of the chronological past and future of timelike separated points in a convex normal neighborhood. The boundary of such a causal diamond is a null surface N with a 2-sphere bifurcation edge B . We use Gaussian null coordinates adapted to the causal diamond boundary to show that both N and B for any causal diamond in any spacetime can be identified. The resulting reduced symmetry algebra which preserves B takes the form diff( S 2 ) glyph[multicloseleft] b with b consisting of the angle-dependent rescalings along the null generators. The angle-dependent translations along the null generators are eliminated by requiring the bifurcation edge B to be preserved.

By considering the behavior of geometric fields on N near its corners, we show that the boundary charges and the boundary presymplectic potential vanish in the limit to the corners of the causal diamond. From this we show that the Wald-Zoupas fluxes are Hamiltonian generators on the covariant phase space, which in particular provides an infinite family of boost generators for any smooth causal diamond in general relativity. This is similar in spirit to the boost generator at Killing horizons, which is also the vacuum modular Hamiltonian, where in the present context the boost generators act on the gravitational data associated to the causal diamond. Furthermore, we show that the reduced symmetry algebra at N has a non-trivial center. The charges associated to the elements of the center are precisely the boost generators, whose values are proportional to the area of B . Thus there exists a Wald

entropy [2], and a quasi-local first law, for any smooth causal diamond in general relativity.

Using the smoothness of the spacetime metric and the vector fields representing the null boundary symmetries we then show that the Wald-Zoupas fluxes associated to the symmetries are conserved between the past and future components of N . This gives an infinite set of conservation laws for finite subregions in general relativity on any spacetime. This is analogous to the conservation laws between past and future null infinity [28-31] except, in this case, the smoothness of fields at B are much simpler to analyze. Just as the asymptotic conservation laws between past and future null infinity place an infinite number of constraints on gravitational scattering (conjectured to hold even in the quantum theory [28]), the conservation laws we derive for finite causal diamonds likely place important constraints on the properties of scattering in local gravitational subsystems.

The rest of the paper is organized as follows. In Sec. 2 we review the formulation of symmetries and the associated charges and fluxes in general relativity from [27], and reformulate them in terms of Gaussian null coordinates. In Sec. 3 we adapt this formalism to the boundaries of causal diamonds and detail the reduction of the symmetry algebra of a general null surface to a subalgebra which preserves the structure on a causal diamond. We also investigate the behavior of the fields and charges near the corners of the causal diamond and show that the fluxes associated to the symmetries at a causal diamond are also Hamiltonian generators on the corresponding phase space. In Sec. 4 we show that the smoothness of the spacetime metric at the causal diamond implies an infinite number of conservation laws for the fluxes through the null boundaries. In Sec. 5 we compute the charges associated to the central elements of the symmetry algebra and show that these take the form of a 'first law' . We end with Sec. 6 summarizing and discussing the potential applications of our results. In Appendix A we collect the essential ingredients of the covariant phase space formalism and the Wald-Zoupas prescription for calculating charges and fluxes. In Appendix B we analyze the structure of the symmetry algebra at a causal diamond and show how it arises as a non-trivial central extension.

## Notation and conventions

We follow the conventions of [32]. We use abstract indices a, b, . . . to denote tensor fields, e.g. g ab is the spacetime metric, and indices A, B, . . . to denote components of tensor fields in some coordinate system on S 2 , e.g. q AB is a metric on S 2 . Boldface quantities like ω will denote differential forms.

We also use the following terminology for the charges associated to the symmetry algebra at a null boundary N . Quantities associated to null boundary symmetries evaluated as integrals over cross-sections of the null boundaries will be called 'charges', while the difference

of these charges on two cros-sections evaluated as an integral over a portion of the null boundary will be called 'fluxes'. When certain conditions are satisfied the fluxes can also be considered as Hamiltonian generators on the null boundary phase space (see Eqs. A.10 and A.11).

## 2. NULL BOUNDARY SYMMETRIES AND CHARGES

In this section we briefly review the basic formalism and results of [27], namely the symmetries and charges at a null boundary in general relativity. We then recast the null boundary phase space, and the resulting symmetries and charges, in terms of Gaussian null coordinates. This will prove to be useful when considering causal diamonds.

The relationship between the covariant approach of [27] which is intrinsic to the null boundary and the coordinate-based approach in Sec. 2.2 and 2.3 is the same as that between the intrinsic universal structure approach [33-35] and one based on Bondi coordinates [36-38] or the conformal Gaussian null coordinates [39, 40] at null infinity in asymptotically-flat spacetimes.

## 1. Universal structure and symmetries on null boundaries

Consider a spacetime ( M,g ab ) with null boundary N . For now we will assume the null generators of N are complete, i.e. N ∼ Z × = R , where Z is the space of null generators. Later we will consider null surfaces with boundary, which is the relevant setting for causal diamonds. The null boundary N is naturally equipped with the equivalence class [ /lscript a , κ ] where /lscript a is the null generator of N κ , is the non-affinity defined by 1

$$\ell ^ { b } \nabla _ { b } \ell ^ { a } & \cong \kappa \ell ^ { a },$$

and the equivalence class [ /lscript a , κ ] is defined by the rescaling freedom

↦

$$\ell ^ { a } \mapsto e ^ { \beta } \ell ^ { a } \,, \quad \kappa \mapsto e ^ { \beta } ( \kappa + \mathcal { L } _ { \ell } \beta )$$

↦

where β is a smooth function on N .

In [27] it was shown that the structure [ /lscript a , κ ] is universal in the sense that different such structures on N , as induced by different background metrics, are all related by diffeomorphisms (we shall show this explicitly in Sec. 2.2 in a Gaussian null coordinate system). We can then define the field configuration space F to be the set of smooth metrics g ab on a

1 We use the notation = to mean 'equality on ̂ N ' throughout the paper.

manifold M with null boundary N which is equipped with the universal structure [ /lscript a , κ ]. The covariant phase space is the subset F ⊂ F consisting of on-shell metrics satisfying the vacuum Einstein equation.

The group of symmetries on a null boundary N is the subgroup of diffeomorphisms on M which preserves the null boundary N and the universal structure on it. It will be easier to work with the symmetry algebra instead of the group. The symmetry algebra on the null boundary consists of vector fields ξ a on M which are tangent to N and preserve the linearized version of the equivalence relation Eq. 2.2. This results in the conditions

$$\mathcal { L } _ { \xi } \ell ^ { a } & \hat { = } \beta \ell ^ { a } \\ \mathcal { L } _ { \xi } \kappa & \hat { = } \beta \kappa + \mathcal { L } _ { \ell } \beta,$$

where β is some smooth function on N which depends on the vector field ξ a . The detailed structure of the resulting symmetry algebra g was derived in [27] and can be summarized as follows. The vector fields of the form ξ a ̂ = f/lscript a with £ £ /lscript ( /lscript + κ f ) ̂ = 0 form an infinitedimensional abelian Lie ideal s ⊂ g of supertranslations . The quotient algebra g s / is isomorphic to diff( Z ), the algebra of smooth diffeomorphisms of the space of null generators Z . There is an additional Lie ideal s 0 ⊂ s of affine supertranslations given by ξ a ̂ = f/lscript a with ( £ /lscript + κ f ) ̂ = 0. Hence the symmetry algebra g can be written as

$$\mathfrak { g } \cong \text{diff} ( \mathcal { Z } ) \ltimes ( \mathfrak { b } \ltimes \mathfrak { s } _ { 0 } )$$

where b ∼ = s / s 0 .

The charges and fluxes associated to this symmetry algebra were also derived in [27] using the Wald-Zoupas procedure. Writing the covariant expression for these charges would require introducing a significant amount of formalism and notation. Instead we will derive the symmetry algebra, and express the charges and fluxes, using Gaussian null coordinates in Sec. 2.3.

## 2. Gaussian null coordinates

It will be convenient to introduce coordinates adapted to the null surface N called Gaussian null coordinates (GNC) [41], which have been used in a variety of contexts [40, 42-45]. We briefly review the construction of GNC below. Since our main interest is in causal diamonds, we will now restrict (for convenience) to the case of 4-dimensional spacetimes where the space of null generators is a 2-sphere Z = S 2 , but our results can be readily generalized.

Let /lscript a be an affinely-parameterised ( κ = 0) null normal to N and let v be an affine

parameter along these null generators i.e., v is some function on N such that /lscript a ∇ a v ̂ = 1. Now let S ∼ = S 2 be a cross-section of N such that v | S = 0, and let x A be a coordinate system on S . We extend the coordinate functions x A to all of N by parallel-transport along the null generators, /lscript a ∇ a x A ̂ = 0. This defines a coordinate system ( v, x A ) on N .

To define a coordinate system in a neighborhood of N , let u be a function in such a neighborhood so that u | N = 0 on N . Then, /lscript a ≡ -du is the normal to N and the vector field n a ≡ ∂/∂u is transverse (i.e. not tangent) to N . To fix coordinates away from N we choose u such that n a is an affinely-parameterised null vector field i.e. n n a a = 0 and n b ∇ b n a = 0. Then we extend the coordinates ( v, x A ) away from N by parallel transport along n a . The coordinate functions ( u, v, x A ) define a GNC in a neighborhood of the null surface N . It follows from the above definition of the GNC that in these coordinates the spacetime metric satisfies [45]

$$g _ { u u } & = g _ { u A } = 0 \,, \quad g _ { u v } = - 1 \\ g _ { v v } & \cong g _ { v A } \, \hat { = } \, \partial _ { u } g _ { v v } \, \hat { = } \, 0$$

and thus we can write the line element in the form (this is equivalent to the form in [41])

$$& d s ^ { 2 } = - W d v ^ { 2 } - 2 d u d v + q _ { A B } ( d x ^ { A } - W ^ { A } d v ) ( d x ^ { B } - W ^ { B } d v ) \\ \text{where} \quad & W | _ { u = 0 } = \partial _ { u } W | _ { u = 0 } = W ^ { A } | _ { u = 0 } = 0$$

and W W , A , q AB are functions of ( u, v, x A ), and can be considered as tensors on S 2 which depend on ( u, v ). The tensor q AB defines a Riemannian metric on the 2-spheres of constant u and v . The extensions of the null generator /lscript a and the auxilliary null vector n a in the neighborhood of N are given by

$$\ell ^ { a } \equiv \partial _ { v } - \frac { 1 } { 2 } W \partial _ { u } + W ^ { A } \partial _ { A } \,, \quad n ^ { a } \equiv \partial _ { u }$$

The shear σ AB and expansion θ of N are given by the relation

$$\frac { 1 } { 2 } \partial _ { v } q _ { A B } \, \widehat { = } \, \sigma _ { A B } + \frac { 1 } { 2 } q _ { A B } \theta$$

where σ AB q AB = 0, while the Há´ iček rotation 1-form of the u = constant cross-sections is given by

$$\omega _ { A } \, \widehat { \equiv } \, - q _ { A } \, ^ { a } n _ { b } \nabla _ { a } \ell ^ { b } = - \frac { 1 } { 2 } \partial _ { u } ( q _ { A B } W ^ { B } )$$

We emphasize that the above construction of the GNC can be carried out in any spacetime in a neighborhood of any null surface. Now let ( M ,g 1 1 ) and ( M ,g 2 2 ) be two spacetimes with null surfaces N 1 and N 2 along with the GNCs ( u , v 1 1 , x A 1 ) and ( u , v 2 2 , x A 2 ), as constructed above, respectively. Without any loss of generality we can identify a neighborhood of N 1 in ( M ,g 1 1 ) with that of N 2 in ( M ,g 2 2 ) by identifying the corresponding GNCs ( u , v 1 1 , x A 1 ) =

( u , v 2 2 , x A 2 ). Thus, we can identify all the spacetimes under consideration, and work on a single manifold M with boundary N such that the configuration space F consists of all the metrics in M for which N is a null surface and the metrics take the form Eq. 2.6 in a GNC in a neighborhood of N . From the above construction we see that while the induced metric q AB depends on the particular choice of the spacetime metric in F , the null generator /lscript a is common to all metrics in F . Thus, the null surface N along with the affine null generator /lscript a are universal . Note that we can construct the GNC even with a non-affinely parametrized /lscript a , which leads to the universal structure used in [27] as described in Sec. 2.1.

## 3. Null boundary symmetries, charges and fluxes in GNC

The symmetry algebra at the null boundary N consists of the vector fields ξ a generating infinitesimal coordinate transformations which preserve the GNC form of the metric in Eq. 2.6. We expand the vector field ξ a in the GNC to first order in u as,

$$\xi ^ { a } \equiv ( f _ { 0 } + u f _ { 1 } ) \partial _ { v } + ( \beta _ { 0 } + u \beta _ { 1 } ) \partial _ { u } + ( X _ { 0 } ^ { A } + u X _ { 1 } ^ { A } ) \partial _ { A } + O ( u ^ { 2 } )$$

To preserve the location u = 0 of the null surface N ξ , a must be tangent to N and hence β 0 = 0. To preserve the form of the metric Eq. 2.6 we have

$$\mathcal { L } _ { \xi } g _ { u u } = \mathcal { L } _ { \xi } g _ { u v } = \mathcal { L } _ { \xi } g _ { u A } = 0$$

$$\mathcal { L } _ { \xi } g _ { v A } & = O ( u )$$

$$\mathcal { L } _ { \xi } g _ { v v } = O ( u ^ { 2 } )$$

Evaluating Eq. 2.11a at u = 0 we have

$$\mathcal { L } _ { \xi } g _ { u u } = 0 \implies f _ { 1 } = 0$$

$$\mathcal { L } _ { \xi } g _ { u v } = 0 \implies \beta _ { 1 } = - \partial _ { v } f _ { 0 }$$

$$\mathcal { L } _ { \xi } g _ { u A } = 0 \implies q _ { A B } X _ { 1 } ^ { A } = \partial _ { A } f _ { 0 }$$

The conditions Eq. 2.11b imply

$$\partial _ { v } X _ { 0 } ^ { A } = 0$$

while Eq. 2.11c, evaluated to O u ( ) using Eq. 2.13, gives

$$\partial _ { v } \beta _ { 1 } = 0 \implies \partial _ { v } ^ { 2 } f _ { 0 } = 0$$

where the second condition follows from Eq. 2.12b. Similar conditions were derived independently in [45].

From Eqs. 2.13 and 2.14 we conclude that the symmetries on the null boundary N are characterized by ( α, β, X A ) where α and β are functions and X A is a vector field on S 2 . In a neighborhood of N this symmetry is represented by a vector field ξ a , which in GNC takes the form

$$\xi ^ { a } \equiv ( \alpha - v \beta ) \partial _ { v } + u \beta \partial _ { u } + \left [ X ^ { A } + u q ^ { A B } \partial _ { B } ( \alpha - v \beta ) \right ] \partial _ { A } + O ( u ^ { 2 } ) \, \quad \ \ ( 2. 1 5 )$$

Note that the vector field ξ a at N , i.e. u = 0, is parametrized entirely by ( α, β, X A ) and is independent of the choice of metric in the configuration space F .

We now analyze the structure of the symmetry algebra g generated by such vector fields. Consider two symmetries ξ 1 = ( α , β , X 1 1 A 1 ) and ξ 2 = ( α , β , X 2 2 A 2 ) in g . Their Lie bracket can be computed using their representations in terms of vector fields as in Eq. 2.15. A straightforward computation gives

$$\left [ ( \alpha _ { 1 }, \beta _ { 1 }, X _ { 1 } ^ { A } ), ( \alpha _ { 2 }, \beta _ { 2 }, X _ { 2 } ^ { A } ) \right ] = ( \alpha, \beta, X ^ { A } )$$

$$\text{where} \quad \alpha = - \alpha _ { 1 } \beta _ { 2 } + \alpha _ { 2 } \beta _ { 1 } + X _ { 1 } ^ { A } \partial _ { A } \alpha _ { 2 } - X _ { 2 } ^ { A } \partial _ { A } \alpha _ { 1 } & & ( 2. 1 6 b )$$

$$\beta = X _ { 1 } ^ { A } \partial _ { A } \beta _ { 2 } - X _ { 2 } ^ { A } \partial _ { A } \beta _ { 1 }$$

$$X ^ { A } = \left [ X _ { 1 }, X _ { 2 } \right ] ^ { A } = X _ { 1 } ^ { B } \partial _ { B } X _ { 2 } ^ { A } - X _ { 2 } ^ { B } \partial _ { B } X _ { 1 } ^ { A }$$

where the last line is the Lie bracket of vector fields on S 2 . Note that the sign of β in Eq. 4.10b [27] is incorrect and has been corrected in Eq. 2.16c above.

From Eq. 2.16 it is easy to deduce the following structure of g . If X A 1 = 0 then X A = 0, i.e. symmetries of the form ( α, β, 0) form an abelian Lie ideal s ⊂ g of supertranslations . The quotient g s / is then isomorphic to the Lie algebra diff( S 2 ) represented by symmetries of the form (0 0 , , X A ). There is an additional Lie ideal in g which is given as follows. In Eq. 2.16, taking β 1 = X A 1 = 0 we get β = X A = 0, that is, symmetries of the form ( α, 0 0) , are also an abelian Lie ideal s 0 ⊂ g called affine supertranslations . The quotient b ∼ = s / s 0 of all the supertranslations by s 0 is represented by symmetries of the form (0 , β, 0). Thus, the symmetry algebra g on any null boundary has the structure (same as in Eq. 2.4)

$$\mathfrak { g } \cong \text{diff} ( \mathbb { S } ^ { 2 } ) \ltimes ( \mathfrak { b } \ltimes \mathfrak { s } _ { 0 } )$$

It was shown in [27] that this symmetry algebra coincides with the definition given by Wald and Zoupas [3], reviewed below Eq. A.4.

↦

Remark 2.1 (Symmetry group at N ) . The symmetry group can also be obtained by considering finite coordinate transformations of the GNC which preserve the metric form Eq. 2.6. In particular at N , i.e. u = 0, we have the coordinate transformations ( v, x A ) → ( v, x A ) with

$$\bar { v } = \alpha ( x ^ { A } ) + e ^ { - \beta ( x ^ { A } ) } v \,, \quad \bar { x } ^ { A } = \bar { x } ^ { A } ( x ^ { B } )$$

Thus, the symmetry group consists of arbitrary diffeomorphisms of S 2 along with angledependent translations (given by α x ( A )) and angle-dependent rescalings (given by β x ( A )) along the null generators.

When the null surface N has additional structure which is also universal - i.e., common to all the spacetimes under consideration - the symmetry algebra can be reduced further. For instance, when all the spacetimes have a Killing vector field in a neighborhood of N which becomes tangent to N , the symmetries proportional to this Killing field provide a preferred 1-dimensional subalgebra of g (see Sec. 4.4 [27]). Similarly, if the null surface is stationary for all spacetimes, so that the shear and expansion of N vanish, then the symmetry algebra can be reduced so that β = constant and X A is a conformal Killing field on S 2 i.e., an element of the Lorentz algebra (see Sec. IV.B [46]). 2 We show in Sec. 3.1 that when N is the null boundary of causal diamonds a similar reduction of the symmetry algebra occurs due to the presence of a preferred cross-section corresponding to the bifurcation edge. Specifically, since the bifurcation edge is a preferred cross-section of the null boundary of a causal diamond, only those symmetry vector fields which preserve its location i.e. are tangent to the bifurcation surface are permitted in the symmetry algebra. This has the effect of eliminating the affine supertranslations ξ a ≡ α∂ v ∈ s 0 from the symmetry algebra (see also Sec. 4.5 [27]).

The charges and fluxes associated to the null boundary algebra g were computed in [27] using the covariant phase space formalism along with the Wald-Zoupas prescription. It was also shown that the ambiguities in the symplectic current and the Wald-Zoupas prescription do not affect the resulting charges and fluxes. We do not repeat the full analysis of [27], but below we write down the relevant expressions for the boundary presymplectic potential Θ ( g ; δg ), the Wald-Zoupas (WZ) charges Q ξ and fluxes F ξ for vacuum general relativity, derived in [27], in terms of GNC.

The boundary presymplectic potential on N is given by

$$\Theta ( g ; \delta g ) = \frac { 1 } { 1 6 \pi } \varepsilon _ { 3 } \left ( \sigma ^ { A B } - \frac { 1 } { 2 } q ^ { A B } \theta \right ) q _ { A } ^ { a } q _ { B } ^ { b } \delta g _ { a b }$$

where ε 3 ≡ ε abc is the 3-volume element on N .

Let S be any cross-section of N with area-element ε 2 ≡ ε ab and ∆ N be a region of N bounded by two cross-sections. The charges (on S ) and fluxes (through ∆ N ) associated to a

2 Note that the reduction of diff( S 2 ) to the Lorentz algebra in the stationary case relies crucially on the S 2 topology of the cross-sections.

supertranslation ξ a ̂ = f/lscript a = ( α -vβ /lscript ) a are

$$\mathcal { Q } _ { f } [ S ] & = \frac { 1 } { 8 \pi } \int _ { S } \varepsilon _ { 2 } \left [ ( \alpha - v \beta ) \theta + \beta \right ] \\ \mathcal { F } _ { f } [ \Delta N ] & = \frac { 1 } { 8 \pi } \int _ { \Delta N } \varepsilon _ { 3 } \left ( \alpha - v \beta \right ) \left ( \sigma _ { A B } \sigma ^ { A B } - \frac { 1 } { 2 } \theta ^ { 2 } \right )$$

while those associated to a diff( S 2 ) generator X A (taken to be tangent everywhere to the v = constant cross-sections of N ) are given by

$$\mathcal { Q } _ { X } [ S ] & = \frac { 1 } { 8 \pi } \int _ { S } \varepsilon _ { 2 } \left ( - \omega _ { A } X ^ { A } \right ) \\ \mathcal { F } _ { X } [ \Delta N ] & = \frac { 1 } { 8 \pi } \int _ { \Delta N } \varepsilon _ { 3 } \left ( \sigma _ { A B } - \frac { 1 } { 2 } \theta q _ { A B } \right ) D ^ { A } X ^ { B }$$

where the Há´ iček rotation 1-form ω A (in the GNC foliation) is as defined in Eq. 2.9, and D A is the derivative operator with respect to q AB in the foliation given by the GNC.

Note that the supertranslation charge expression Eq. 2.20 can be evaluated on any choice of cross-section, while the charge expression Eq. 2.21 only holds on the cross-sections of the v = constant foliation. 3 If the diff( S 2 ) charge is to be evaluated on some arbitrary foliation of N with normal ˆ n a such that /lscript a ˆ n a ̂ = -1, then we have instead

$$\mathcal { Q } _ { \hat { X } } [ \hat { S } ] = \frac { 1 } { 8 \pi } \underset { \hat { S } } { \int } \hat { \varepsilon } _ { 2 } \left ( \beta _ { \hat { X } } - \hat { \omega } _ { A } \hat { X } ^ { A } \right )$$

where now ˆ X A is taken to be tangent to the cross-sections of the chosen foliation, ˆ ω A is the corresponding Há´ iček rotation 1-form, while β ˆ X is given by

$$\beta _ { \hat { X } } = - \hat { X } ^ { A } \hat { q } _ { A } ^ { a } \mathcal { L } _ { \ell } \hat { n } _ { a }$$

Remark 2.2 (Fluxes from the charges) . Given the charge expressions Eqs. 2.20 and 2.21, the corresponding fluxes can be obtained using the vacuum Einstein equations R ab = 0 on N . Specifically we have

$$R _ { a b } \ell ^ { a } \ell ^ { b } & = 0 \implies \partial _ { v } \theta \quad = - \frac { 1 } { 2 } \theta ^ { 2 } - \sigma _ { A B } \sigma ^ { A B } \\ R _ { a b } q _ { A } ^ { a } \ell ^ { b } & = 0 \implies \partial _ { v } \omega _ { A } \quad = - \theta \omega _ { A } - D ^ { B } \sigma _ { A B } + \frac { 1 } { 2 } D _ { A } \theta$$

where D A is the derivative operator with respect to q AB in the foliation given by the GNC. These are the Raychaudhuri equation and the Damour-Navier-Stokes equation respectively

3 The dependence of the diff( S 2 ) charge expression on the foliation is a result of the semidirect structure of g in Eq. 2.6, that is, there does not exist any unique choice of a diff( S 2 ) subalgebra of g . This is similar to the status of the Lorentz algebra within the BMS algebra at null infinity.

## [47, 48].

## 3. CAUSAL DIAMONDS

In this section we apply the above results to the boundaries of causal diamonds, with the appropriate modifications necessary for null surfaces with boundary. We begin by recalling the definition of causal diamonds and the structure on their null boundaries.

In a given spacetime ( M,g ), consider two points p + , p -∈ M such that p + is in a convex normal neighborhood of p -and is in its chronological future, i.e., p + is 'inside the future light cone' of p -. The intersection of the chronological past of p + with the chronological future of p -defines a causal diamond or a double cone (see for instance [25, 26, 49, 50]). We assume that the causal diamond is 'small enough' so that the null generators emanating from p ± form smooth null surfaces N ± respectively, which intersect at a smooth 2-surface B , the bifurcation edge , which is topologically S 2 . We denote the null boundary of the causal diamond by N = N + ∪ N -(see Fig. 1).

FIG. 1: Diagram of a causal diamond in a spacetime ( M,g ). The points p ± denote the corners of the causal diamond and B is the bifurcation edge while N ± denote the future/past null surfaces joining B to p ± , respectively. The functions v and u are affine coordinates with affine null normals /lscript a and n a on N ± .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00017v2-figure-d4709ab0-92cf-4ab4-b951-567d530acf78.jpg)

We now investigate the null boundary symmetries and charges defined in [27] for the null boundary N of a causal diamond. It will be convenient to use the formulation in terms of GNC as detailed in Sec. 2.2, which we adapt to a causal diamond as follows.

Unlike a general null surface, the boundary of a causal diamond has a preferred crosssection determined by the bifurcation edge B . At B there exist unique null vector fields /lscript a and n a both future-directed such that /lscript a is tangent to N + , n a is tangent to N -, and /lscript a n a B | = -1. We can extend these vector fields to N ± so that /lscript a is the affine null generator of N + , n a is the affine null generator of N -, and /lscript a n a N | = -1. Let x A be some coordinates on B ; we pick the affine parameter v of /lscript a ≡ ∂/∂v on N + , and similarly u of n a ≡ ∂/∂u on N -, such that B ≡ ( u = 0 , v = 0). Since, /lscript a is future-directed, the coordinate v increases moving towards p + from v = 0 at B . Note that in a general spacetime N + will have both a shear and an expansion, which depend on the space of generators, hence the value of the affine-parameter v at the corner p + will depend on the null generator along which we

approach p + i.e. v | p + = V ( x A ). Similarly, on N -the affine parameter u along n a decreases moving towards p -from u = 0 at B and p -lies at u | p -= U x ( A ). These are depicted in Fig. 1.

As described in Sec. 2.2, we can extend ( u, v, x A ) to form a GNC in a neighborhood of the causal diamond. 4 Since we have two null surfaces we obtain two different GNCs, one based on N + which we denote by ( u + , v + , x A + ) and another based on N -denoted by ( u -, v -, x A -). In general, these two coordinate systems will not agree in a neighborhood of B and will be related by a coordinate transformation that preserves neither GNC. We will not need the explicit form of the transformations between these coordinates but we note that (by construction)

$$( u _ { \pm }, v _ { \pm } ) | _ { B } = ( u = 0, v = 0 ) \,, \quad x _ { + } ^ { A } \Big | _ { B } = x _ { - } ^ { A } \Big | _ { B }$$

and

$$\ell ^ { a } \equiv \partial _ { v } = \partial _ { v _ { + } } \Big | _ { B } = \partial _ { v _ { - } } \Big | _ { B } \, \ \ n ^ { a } \equiv \partial _ { u } = \partial _ { u _ { + } } \Big | _ { B } = \partial _ { u _ { - } } \Big | _ { B }$$

The spacetime metric g ab , which we assume is smooth written in either coordinate system, coincides at B .

We define the 3-volume elements ε ± abc on N ± and the 2-area elements ε ± ab on the crosssections of N ± as follows:

$$\text{on } N ^ { + } & \colon \quad \varepsilon _ { a b c } ^ { + } = n ^ { d } \varepsilon _ { d a b c } \,, \quad \varepsilon _ { a b } ^ { + } = - \ell ^ { c } \varepsilon _ { c a b } ^ { + } = - \ell ^ { c } n ^ { d } \varepsilon _ { d c a b } \\ \text{on } N ^ { - } & \colon \quad \varepsilon _ { a b c } ^ { - } = - \ell ^ { d } \varepsilon _ { d a b c } \,, \quad \varepsilon _ { a b } ^ { - } = - n ^ { c } \varepsilon _ { c a b } ^ { - } = n ^ { c } \ell ^ { d } \varepsilon _ { d c a b }$$

Note that on N + these conventions are the same as those of [27] while on N -the sign of ε -abc is the opposite of that in [27]. We have chosen these conventions so that the area elements on the bifurcation edge B induced from N ± coincide, that is,

$$\varepsilon _ { a b } ^ { + } \Big | _ { B } = \varepsilon _ { a b } ^ { - } \Big | _ { B }$$

Similar to the case of a general null surface, we can now identify the boundaries of any two causal diamonds in any two spacetimes by identifying the GNCs ( u ± , v ± , x A ± ). Note that with this identification the bifurcation edge B ≡ ( u = 0 , v = 0) is common to all causal diamonds and is universal , but the corners u = U x ( A ) and v = V ( x A ) depend on the specific choice of causal diamond and spacetime metric, and are thus not universal. Henceforth we will work with the covariant phase space F of general relativity at the boundary N of a causal diamond where the bifurcation edge B is a common universal surface for all spacetimes in F .

4 To define the GNC in a neighborhood of B , we need to extend the null surface N + smoothly 'a little' to the past of B , and similarly extend N -to the future of B . We assume, henceforth, that this has been done.

## 1. Reduced symmetry algebra g CD at causal diamonds

Since the bifurcation edge B is universal the symmetry algebra for a causal diamond must preserve B . Consider the null boundary symmetry algebra on the future null surface N + . From the form of the vector fields ξ a in Eq. 2.15 we see that the symmetries on N + which preserve the surface B ≡ ( v = 0) are the ones which satisfy α x ( A ) | N + = 0. In other words, the bifurcation edge B breaks the affine supertranslation symmetry. Similarly, the affine supertranslations of the symmetry algebra at the past surface N -are also broken.

A priori it seems we have two independent symmetries for the causal diamond: one induced from N + and the other from N -, with the respective affine supertranslations set to vanish. However, there is a natural isomorphism between the future and past symmetries which follows from the smoothness of the vector field ξ a in spacetime. To see this let ξ a be a smooth vector field in the spacetime M which is a representative of a symmetry on N ± respectively, preserving the bifurcation edge B . In the GNCs ( u ± , v ± , x A ± ) based on the null surfaces N ± we have (see Eqs. 2.15 and 3.2)

$$\xi ^ { a } & \equiv \beta _ { + } ( - v _ { + } \partial _ { v _ { + } } + u _ { + } \partial _ { u _ { + } } ) + X _ { + } ^ { A } \partial _ { A _ { + } } + \dots \\ & \equiv \beta _ { - } ( - u _ { - } \partial _ { u _ { - } } + v _ { - } \partial _ { v _ { - } } ) + X _ { - } ^ { A } \partial _ { A _ { - } } + \dots$$

where as before . . . denotes the subleading terms in the respective GNCs. Note that while the GNCs do not match in a neighborhood of B , from Eqs. 3.1 and 3.2 and the smoothness of ξ a at B we can conclude that

$$\beta _ { + } | _ { B } = - \, \beta _ { - } | _ { B } \, \quad X _ { + } ^ { A } \Big | _ { B } = X _ { - } ^ { A } \Big | _ { B }$$

This implies a natural isomorphism between the symmetry algebras on N + and N -given by the matching conditions Eq. 3.6 at B . Thus the elements of the symmetry algebra g CD on the boundary of a causal diamond are given by ( β, X A ); for definiteness we choose ( β, X A ) = ( β + , X A + ) = ( -β -, X A -) to represent an element in g CD .

The Lie brackets of the algebra g CD can be derived from Eq. 2.16 by setting α 1 = α 2 = 0. We have

$$& \left [ ( \beta _ { 1 }, X _ { 1 } ^ { A } ), ( \beta _ { 2 }, X _ { 2 } ^ { A } ) \right ] = ( \beta, X ^ { A } ) \\ & \text{where} \quad \beta = X _ { 1 } ^ { A } \partial _ { A } \beta _ { 2 } - X _ { 2 } ^ { A } \partial _ { A } \beta _ { 1 } \\ & X ^ { A } = [ X _ { 1 }, X _ { 2 } ] ^ { A } = X _ { 1 } ^ { B } \partial _ { B } X _ { 2 } ^ { A } - X _ { 2 } ^ { B } \partial _ { B } X _ { 1 } ^ { A } \\.$$

If X A 1 = 0 then X A = 0 hence symmetries of the form ( β, 0) form an infinite-dimensional abelian Lie ideal b of boost supertranslations . Thus,

$$\mathfrak { g } _ { \text{CD} } & \cong \text{diff} ( \mathbb { S } ^ { 2 } ) \ltimes \mathfrak { b }$$

Further, if β 1 = constant and X A 1 = 0 then β = X A = 0, that is, the symmetries of the form ( β = constant 0) commute with any element of , g CD and thus form a 1-dimensional Lie subalgebra b 0 of central elements which we call boosts . 5 Consider the quotient g CD / b 0 ∼ = diff( S 2 ) glyph[multicloseleft] ( b b / 0 ). Then g CD has the structure of a central extension of g CD / b 0 by the abelian Lie algebra b 0 ; the fact that this is a non-trivial central extension is shown in Appendix B. 6

We show in Sec. 5 that the charges associated to the central elements in b 0 can be interpreted as providing a 'first law' for causal diamonds.

## 2. Behavior of the fields and charges near the corners

In this section we analyze the behavior of the relevant fields and charges on N + near the corner p + ; similar results hold also for N -at p -. Essentially, for causal diamonds in smooth spacetimes the behavior of the fields of interest near p + is the same as that of a light cone in Minkowski spacetime with some subleading corrections away from p + . We invoke the results of [51, 52] below. The main consequence of interest for our purposes is that in the limit to the corner p + along N + we have

$$\mathcal { Q } _ { \xi } \rightarrow 0 \,, \quad i _ { \xi } \Theta ( \delta g ) \rightarrow 0$$

for all symmetries ξ ∈ g CD , and all metric perturbations δg ab which are smooth at p + in any spacetime. The limit Q → ξ 0 near the corner ensures that the total flux F ξ associated to the symmetries through all of N + is finite and is, in fact, equal to the charge at B . As we will show in Sec. 4, along with smoothness of the spacetime at B , this gives infinitely many conservation laws between the incoming and outgoing fluxes through any causal diamond. The implication of the limit i ξ Θ ( δg ) → 0 is as follows: a perturbation of the total flux on N + can be written as (see Eq. A.10)

$$\delta \mathcal { F } _ { \xi } = \underset { N ^ { + } } { \int } \, \omega ( g ; \delta g, \mathcal { L } _ { \xi } g ) + \underset { B } { \int } i _ { \xi } \Theta ( \delta g ) - \underset { p ^ { + } } { \int } i _ { \xi } \Theta ( \delta g )$$

where the integral over the corner p + should be interpreted as a limit of integrals on crosssections of N + which suitably limit to p + as described below. Since i ξ Θ ( δg ) → 0 in the limit to p + and i ξ Θ ( δg ) | B = 0 (as ξ a is tangent to B for any symmetry in g CD ), we have

$$\delta \mathcal { F } _ { \xi } = \underset { N ^ { + } } { \int } \omega ( g ; \delta g, \mathcal { L } _ { \xi } g )$$

5 The terminology 'boosts' for elements of b 0 is motivated by the fact that if one considers a bifurcate Rindler horizon in Minkowski spacetime, instead of a causal diamond, then Lorentz boosts which preserve the Rindler horizon are precisely elements in b 0 .

6 Note that if one eliminates the non-constant boost supertranslations, for instance by imposing a weakly isolated horizon structure when N is stationary, then the central extension also becomes trivial (see for instance Sec. IV.B [46]). 15

for all symmetries ξ ∈ g CD and all perturbations δg ab which are smooth at p + and B . Thus the flux F ξ , viewed as a function on the covariant phase space on N + , generates a Hamiltonian flow associated to the symmetry ξ a (see Eq. A.11).

We now describe the arguments leading to the above described result, referring to [51, 52] for the details. Since the value of the GNC coordinate v at the corner is direction-dependent, v not a regular coordinate at p + . Similarly, the cross-sections v = constant do not limit to p + . Thus, to analyze the behavior of the symmetries and charges we need a coordinate system on N + which is more suited to the structure near p + . In a neighborhood of p + such a coordinate system can be constructed as follows. Consider the tangent space Tp + at p + , and let y i = ( y , y 0 1 , y 2 , y 3 ) be coordinates in Tp + so that y i | + p = 0 and the coordinate vector fields ∂ i are orthonormal, with ∂ 0 being timelike and future-directed. Defining

$$r ^ { 2 } = ( y ^ { 1 } ) ^ { 2 } + ( y ^ { 2 } ) ^ { 2 } + ( y ^ { 3 } ) ^ { 2 } \,, \quad u = y ^ { 0 } - r$$

the past-directed light cone in Tp + from p + is then given by u = 0, and coordinatized by ( r, x A ) where x A are coordinates on the space of past-directed null directions at p + isomorphic to S 2 .

There exists an exponential map from Tp + to a sufficiently small neighbourhood of p + so that y i are coordinates in this neighborhood, called Riemann normal coordinates . In such a neighborhood, using ( u, r, x A ) as coordinates, the metric g ab takes the form (see [51, 52], note we have changed some signs to conform to our orientation conventions)

$$d s ^ { 2 } = \mu d u ^ { 2 } - 2 \nu d u d r - 2 \nu _ { A } d u d x ^ { A } + q _ { A B } d x ^ { A } d x ^ { B }$$

The analysis of [51, 52] then shows that near the corner p + the metric components in Eq. 3.13 behave as

$$\mu = 1 + O ( r ^ { 2 } ) \,, \quad \nu = 1 + O ( r ^ { 4 } ) \,, \quad \nu _ { A } = O ( r ^ { 3 } ) \,, \quad q _ { A B } = r ^ { 2 } q _ { A B } ^ { 0 } + O ( r ^ { 4 } ) \, \quad \quad ( 3. 1 4 )$$

where q 0 AB is the unit-metric on S 2 . Here, for any tensor T AB... we use O r ( k ) to denote that T AB... = r k t AB... for some t AB... which, in general, has a non-vanishing limit as a tensor field on S 2 as r → 0. Roughly speaking, to leading order the metric g ab near p + behaves as the Minkowski metric at the corner of a light cone.

The expansion and shear of N + have the behavior

$$\theta = - \frac { 2 } { r } + O ( r ^ { 3 } ) \,, \quad \sigma _ { A B } = O ( r ^ { 3 } )$$

The normal to the foliation by r = constant surfaces is ˆ n a = dr . The Há´ iček rotation 1-form on N + relative to the foliation by r is essentially the quantity denoted by ξ A in [51, 52],

which satisfies

$$\hat { \omega } _ { A } = O ( r ^ { 2 } )$$

We have put a 'hat' on the rotation 1-form to emphasize its dependence on the foliation.

To consider the limit of the charges associated to the symmetries on N + , we now relate these coordinates to the GNC used in the main arguments above. The non-affinity ˆ of the κ null generator ˆ /lscript a ≡ -∂ r is given by

$$\hat { \kappa } = - \partial _ { r } \ln \nu = O ( r ^ { 3 } )$$

and thus ˆ /lscript a is an affine null generator of N + up to O r ( 3 ). Thus near p + , we can identify ˆ /lscript a with the GNC null generator /lscript a ≡ ∂ v and the coordinate r with GNC coordinate v as

$$\ell ^ { a } = \hat { \ell } ^ { a } + O ( r ^ { 3 } ) \,, \quad v - V ( x ^ { A } ) = - r + O ( r ^ { 4 } )$$

Note that, as is to be expected, the cross-sections of N + given by r = constant and those given by v = constant do not coincide. In particular, their normals ˆ n a ≡ dr and n a ≡ -dv are related by

$$n _ { a } = \hat { n } _ { a } + \partial _ { A } V d x ^ { A } + O ( r ^ { 3 } )$$

Now consider a symmetry ξ a = -vβ∂ v + X ∂ A A in GNC where, as before, X A is tangent to the v = constant cross-sections. We rewrite this vector field as ξ a = ˆˆ f/lscript a + ˆ X ∂ A A so that ˆ X A is tangent to the r = constant cross-sections. From Eqs. 3.18 and 3.19 we have

$$\hat { f } = - ( V - r ) \beta - X ^ { A } \partial _ { A } V + O ( r ^ { 3 } ) \,, \quad \hat { X } ^ { A } = X ^ { A } + O ( r ^ { 3 } )$$

and also

$$\beta _ { \hat { X } } = - \hat { X } ^ { A } \hat { q } _ { A } ^ { b } \mathcal { L } _ { \hat { \ell } } \hat { n } _ { b } = O ( r ^ { 2 } )$$

The WZ charge Q ξ (see Eqs. 2.20 and 2.22) evaluated on some cross-section S r with r = constant is then

$$\mathcal { Q } _ { \xi } [ S _ { r } ] = \frac { 1 } { 8 \pi } \int _ { S _ { r } } \hat { \varepsilon } _ { 2 } \, \left [ \hat { f } \theta + \beta + \hat { \beta } _ { \hat { X } } - \hat { \omega } _ { A } \hat { X } ^ { A } \right ]$$

From Eqs. 3.20 and 3.21 we see that

$$\mathcal { Q } _ { \xi } [ p ^ { + } ] = \lim _ { r \rightarrow 0 } \mathcal { Q } _ { \xi } [ S _ { r } ] = 0$$

where we have used that ˆ ε 2 = r 2 ε 0 2 + O r ( 4 ) with ε 0 2 the area-element of the unit-sphere.

Next consider the integral of i ξ Θ ( δg ) on the cross-sections S r

$$\int _ { \tilde { S } _ { r } } i _ { \xi } \Theta ( \delta g ) = - \frac { 1 } { 1 6 \pi } \int _ { \tilde { S } _ { r } } \hat { f } \left ( \sigma ^ { A B } - \frac { 1 } { 2 } q ^ { A B } \theta \right ) \hat { q } _ { A } ^ { \, a } \hat { q } _ { B } ^ { \, b } \delta g _ { a b }$$

Any metric perturbation δg ab which is smooth at p + has smooth components in the Riemann normal coordinates y i described above, and its spherical components behave as ˆ q A a ˆ q B b δg ab = O r ( 2 ). Thus, we have

$$\lim _ { r \to 0 } \int i _ { \xi } \Theta ( \delta g ) = 0$$

## 4. CONSERVATION LAWS AT CAUSAL DIAMONDS

We now show that there exist an infinite-number of conservation laws associated to the symmetry algebra g CD between fluxes through N -and N + for any causal diamond. These conservation laws follow directly from the smoothness of the relevant fields at the bifurcation edge B .

First we show that the smoothness of the spacetime at B implies that the charges corresponding to the symmetries in g CD evaluated at B are equal. From Eqs. 2.20, 2.21 and 3.6, the charges at B induced from N ± are

$$\mathcal { Q } _ { \xi } [ B ] & = + \frac { 1 } { 8 \pi } \int \varepsilon _ { 2 } ^ { + } \ \left ( + \beta - \omega _ { A } ^ { + } X ^ { A } \right ) \\ \mathcal { Q } _ { \xi } [ B ] & = - \frac { 1 } { 8 \pi } \int \varepsilon _ { 2 } ^ { - } \ \left ( - \beta - \omega _ { A } ^ { - } X ^ { A } \right )$$

where the difference in the sign of these expressions is due to our conventions for the area elements on N ± given in Eq. 3.3 and the matching conditions on the symmetries Eq. 3.6. To show that these charges are equal we need to consider the relation between the Há´ iček rotation 1-forms ω + A and ω -A which can be obtained as follows. Let /lscript a ± and n a ± be the extensions in the respective GNCs of the null vector fields /lscript a and n a on B . Then we can compute

$$\omega _ { A } ^ { + } | _ { B } & = - ( q ^ { + } ) _ { A } ^ { c } n _ { b } ^ { + } \nabla _ { c } ^ { + } \ell _ { + } ^ { b } = ( q ^ { + } ) _ { A } ^ { c } \ell _ { b } ^ { + } \nabla _ { c } ^ { + } n _ { + } ^ { b } \\ & = ( q ^ { + } ) _ { A } ^ { c } \ell _ { b } ^ { - } \nabla _ { c } ^ { + } n _ { - } ^ { b } \\ & = ( q ^ { - } ) _ { A } ^ { c } \ell _ { b } ^ { - } \nabla _ { c } ^ { - } n _ { - } ^ { b } \\ & = - \omega _ { A } ^ { - } | _ { B } \\. &. &. &. &. &. &. &. &. &.$$

where in the first line we have used n + a /lscript a + | B = n /lscript a a = -1, in the second line we have used the fact that /lscript a and n a are continuous at B (see Eq. 3.2), in the third line the continuity

of the induced metric q ab and the spacetime derivative operator ∇ (which follows from the smoothness of the metric g ab ) and in the last line the definition of ω -A at B . Thus, from the smoothness of the spacetime metric and the continuity of the GNCs at B we have 7

$$\omega _ { A } ^ { + } \Big | _ { B } = - \, \omega _ { A } ^ { - } \Big | _ { B }$$

Combining Eq. 4.3 with Eqs. 3.4, 3.6 and 4.1 we have

$$\mathcal { Q } _ { \xi } [ B ] \, \text{from } N ^ { + } = \mathcal { Q } _ { \xi } [ B ] \, \text{from } N ^ { - }$$

Next, we consider the fluxes through N ± given by

$$\mathcal { F } _ { \xi } [ N ^ { + } ] & = \mathcal { Q } _ { \xi } [ B ] - \mathcal { Q } _ { \xi } [ p ^ { + } ] \\ \mathcal { F } _ { \xi } [ N ^ { - } ] & = \mathcal { Q } _ { \xi } [ B ] - \mathcal { Q } _ { \xi } [ p ^ { - } ]$$

Note that the flux on N + is outgoing while that on N -is incoming relative to the causal diamond (in accordance with our conventions Eq. 3.3). As shown in Sec. 3.2 the charges at the corners p ± vanish and thus from Eq. 4.4 we have

$$\mathcal { F } _ { \xi } [ N ^ { + } ] = \mathcal { F } _ { \xi } [ N ^ { - } ]$$

That is, the incoming flux through N -is equal to the outgoing flux through N + for any symmetry in g CD . Thus, there are infinitely-many conservation laws associated to the symmetry algebra on any causal diamond in any spacetime in general relativity.

/negationslash

/negationslash

Remark 4.1 (Affine supertranslations) . Note that in Sec. 3.1 we eliminated the affine supertranslations α x ( A ) = 0 from the symmetry algebra of the causal diamond by demanding that the bifurcation surface B be preserved under the symmetries. If we had kept α then i ξ Θ ( δg ) | B = 0 - since such vector fields are not tangent to B - and thus, the flux of the affine supertranslations is not a Hamiltonian generator on the phase space of the null boundary. Furthermore, the affine supertranslations α /lscript + a defined on N + and α n -a defined on N -cannot be matched at B , as the corresponding vector fields are not continuous. Even if one imposes the condition α + ( x A ) = α -( x A ) by hand, the charges corresponding to the affine supertranslations at B (see Eq. 2.20) do not match since the expansions θ ± along N ± need not be equal at B in general. Thus, there do not exist any conservation laws at a causal diamond in general spacetimes analogous to Eq. 4.6 for the affine supertranslations.

Remark 4.2 (Non-affine parametrization of the null generators) . For convenience we chose the null generators of N ± to be affinely-parametrized, but our result is invariant under this

7 In the Newman-Penrose notation [53] Eq. 4.3 is simply the identity β + α = - -( β -α ), while in the Geroch-Held-Penrose notation [54] it is β -β ′ = - -( β + β ′ ), which follow from the continuity of the spin-coefficients of the spacetime derivative operator ∇ at B .

/negationslash choice. One can construct a GNC on a null surface relative to an arbitrarily parametrized null generator (with κ = 0). The resulting symmetry algebra is then as described in [27] and Sec. 2.1. The affine supertranslations are eliminated by the condition f | B = 0, in which case the boost supertranslations in b are parametrized by the function -( £ /lscript + κ f ) which is invariant under arbitrary rescalings of the null generators. For the boosts in b 0 we have -( £ /lscript + κ f ) = constant. The remainder of our analysis can also be generalized in a similar fashion; we only note that since f | B = 0, the non-affinities κ ± of the generators of N ± do not enter into the matching of the symmetries and charges at B and the resulting conservation laws.

## 5. CENTRAL CHARGES AND AREA OF THE BIFURCATION EDGE

As discussed in Sec. 3.1, the symmetry algebra g CD at the boundary of the causal diamond can be viewed as a non-trivial central extension of g CD / b 0 by the 1-dimensional abelian subalgebra b 0 of boosts. One expects the charges of such central elements to be of special significance. We show below that these charges are directly related to the area of B , in analogy with the Wald entropy formula for black holes [2]. A similar result was found in [6] through different considerations.

From the results of Sec. 4 we have that 8

$$\mathcal { F } _ { \beta } [ N ^ { + } ] = \mathcal { Q } _ { \beta } [ B ] = \frac { \beta } { 8 \pi } \text{Area} ( B )$$

for any ( β = constant , X A = 0) ∈ b 0 . This can be written in a more illuminating form as follows: consider the vector field ξ a | N + = -vβ/lscript a corresponding to an element in b 0 , so that

$$\xi ^ { b } \nabla _ { b } \xi ^ { a } = \kappa _ { ( \beta ) } \xi ^ { a } \,, \quad \kappa _ { ( \beta ) } = - \beta = \text{constant}$$

and thus

$$\mathcal { F } _ { \beta } [ N ^ { + } ] = \mathcal { Q } _ { \beta } [ B ] = - \frac { \kappa _ { ( \beta ) } } { 2 \pi } \times \frac { 1 } { 4 } \text{Area} ( B )$$

If we interpret the charge Q β [ B ] as an 'energy', the 1 4 Area( B ) as an 'entropy' and -κ ( β ) 2 π as a 'temperature', relative to the vector field ξ a , then Eq. 5.3 takes the form of a 'first law' [2, 32, 47]. Note that if ξ a is future-directed on N + we have β &lt; 0 and so κ ( β ) &gt; 0 and the temperature is negative. This difference in sign compared to the temperature of bifurcate Killing horizons essentially arises due to the fact that the future-directed null generator /lscript a points 'inwards' on N + . Such a negative temperature was also found for causal diamonds in maximally symmetric spacetimes in [25]. Also note that for asymptotically-flat stationary

8 It can be verified that for a causal diamond in Minkowski spacetime where B is a sphere of radius R , our conventions Eq. 3.3 give ∫ B ε + 2 = Area( B ) = +4 πR 2 .

black holes the scaling of the horizon Killing vector field K a is fixed by the requirement that at spatial infinity K a asymptotes to a future-directed unit-normalized timelike Killing vector field (plus a rotational vector field). This completely fixes the scaling of the surface gravity, and hence the temperature, of the black hole. In contrast, there is no natural normalization for the boost vector fields at a causal diamond so we get an entire 1-dimensional family of surface gravities κ ( β ) and temperatures corresponding to the symmetries b 0 .

The charges Q β [ B ] associated to elements of b 0 are central even in the sense of a Poisson bracket on the phase space of N + , i.e. the boost charges Poisson-commute with all the other charges. This can be seen as follows: since the fluxes F ξ are Hamiltonian on the phase space of N + (see Eq. 3.11), we can define their Poisson bracket as (see also Sec. 8 [27])

$$\{ \mathcal { F } _ { \xi _ { 1 } }, \mathcal { F } _ { \xi _ { 2 } } \} = - \int _ { N ^ { + } } \omega ( g ; \mathcal { L } _ { \xi _ { 1 } } g, \mathcal { L } _ { \xi _ { 2 } } g )$$

for any two symmetries ξ , ξ 1 2 ∈ g CD . Since i ξ Θ | B = 0 and F ξ [ N + ] = Q ξ [ B ] for any such symmetry it follows from the analysis of Sec. 8 [27] that

$$\{ \mathcal { Q } _ { \xi _ { 1 } } [ B ], \mathcal { Q } _ { \xi _ { 2 } } [ B ] \} = \mathcal { Q } _ { [ \xi _ { 1 }, \xi _ { 2 } ] } [ B ] - \underset { B } { \int } \mathcal { L } _ { \xi _ { 1 } } \mathcal { Q } _ { \xi _ { 2 } }$$

where in the final term Q ξ is the 2-form whose integral on B gives the charges Eq. 4.1. If this term is vanishing then the Poisson algebra of charges is isomorphic to the Lie algebra of spacetime symmetries in g CD (see Sec. 8 [27] for details). In the present case this term does indeed vanish; we have

$$\int _ { B } \mathcal { L } _ { \xi _ { 1 } } \mathcal { Q } _ { \xi _ { 2 } } = \int [ i _ { \xi _ { 1 } } d \mathcal { Q } _ { \xi _ { 2 } } + d \left ( i _ { \xi _ { 1 } } \mathcal { Q } _ { \xi _ { 2 } } \right ) ] = 0$$

where the first term vanishes since any ξ 1 ∈ g CD is tangent to B , while the second term vanishes upon integration over B . Thus the Poisson algebra of the charges on B is isomorphic to the Lie algebra g CD . The boost charges Eq. 5.3 associated to the central elements in b 0 are then also central charges on the phase space on N + in the sense of the Poisson algebra.

We emphasize that the appearance of central charges in the above analysis is quite different from that of previous approaches. In particular, we work in full nonlinear general relativity at any causal diamond in any spacetime satisfying the vacuum Einstein equations without any restriction to 'near horizon' geometries of stationary black holes (as done in [4, 55, 56]). Since the Einstein equations are not conformally invariant there is no conformal symmetry or Virasoro algebra at the causal diamond in the general case we have considered. In fact, since we have allowed for non-vanishing shear, the induced 2-metrics on the cross-sections of the causal diamond are not conformally related (see Eq. 2.8). We always work with smooth

vector fields as representatives of the symmetries (as opposed to the singular vector fields considered in [20]). Furthermore, as discussed above, the Poisson algebra of the charges in our case is isomorphic to the Lie algebra of symmetries with no additional central extension (in contrast with [4, 7, 20]). The central extension we obtain already exists in the structure of the spacetime symmetry algebra g CD .

## 6. DISCUSSION

We studied the covariant phase space formalism at the boundaries of causal diamonds in vacuum general relativity. In suitable Gaussian null coordinates, we showed that one can identify all causal diamonds and their bifurcation edges across all spacetimes, and that the symmetry algebra at the null boundaries of the casual diamond takes the form g CD ∼ = diff( S 2 ) glyph[multicloseleft] b where diff( S 2 ) maps different null generators of the causal diamond boundary into each other and b consists of angle-dependent rescalings of the affine parameter along the null generators. Suitable smoothness conditions at the corners of the causal diamond imply that the Wald-Zoupas charges vanish at the corners - so that the total flux across the null boundary is equal to the charge at the bifurcation edge - and that the Wald-Zoupas fluxes define Hamiltonian generators of the symmetries on the null boundary phase space. The smoothness of the symmetry vector fields and the fields at the bifurcation edge then give rise to an infinite-number of conservation laws for the Wald-Zoupas fluxes between the past and future components of the causal diamond boundary. We also showed that the charge associated to the central elements of the symmetry algebra - i.e. the elements of the subalgebra b 0 consisting of the angle-independent supertranslations - is related to the area of the bifurcation edge through a 'first law' similar to the Wald entropy formula for stationary black holes.

While our analysis focused on causal diamonds in classical vacuum general relativity we expect that it can be generalized to include matter fields described by a suitable QFT on curved spacetimes. For instance, it was shown in [50] that a comparison of suitable states defined on causal diamonds in different spacetimes can be used to extract properties of the local curvature of the spacetimes. Similarly, in [57] it was shown that the relative entropy of quantum states in linearized general relativity in an asymptotically-flat black hole spacetime is related to the area of the black hole and Bondi flux at null infinity. The infinite-dimensional symmetry algebra at the causal diamond could be useful to analyze other properties of a QFT on curved spacetimes.

It has been conjectured that the conservation laws at null infinity strongly constrain the scattering matrix of quantum gravity in asymptotically-flat spacetimes [28]. Similarly, we expect the conservation laws derived in Sec. 4 can be used to constrain the transition

amplitudes in quantum gravity on local causal diamonds. To do this one needs to suitably quantize the gravitational degrees of freedom on the null boundary (see for instance [15, 17, 58]) and promote the charges and fluxes to operators with the bracket structure Eq. 5.5 in the corresponding quantum theory. We leave further investigation of this problem to future work.

We also expect that our analysis can be extended to causal diamonds at an asymptotic boundary in a spacetime, an interesting example of which arises in the AdS/CFT duality. In this context, for an asymptotically-AdS spacetime, the entanglement entropy of a CFT state defined on a causal diamond lying on the asymptotic boundary (conformal to Minkowski spacetime with one fewer dimensions) is dual to the area of the Ryu-Takayanagi surface in the bulk spacetime [23, 59]. The Ryu-Takayanagi surface can itself be considered as the bifurcation edge of an 'entanglement wedge'. Presumably, our analysis can be suitably generalized to this case, taking into account the asymptotic AdS boundary conditions. We expect that the resulting symmetries are related to the boundary modular Hamiltonian, and that the associated charges and fluxes could provide further insight into the bulk dual of boundary modular flow following [60-62].

## ACKNOWLEDGEMENTS

V.C. is supported in part by the Berkeley Center for Theoretical Physics; by the Department of Energy, Office of Science, Office of High Energy Physics under QuantISED Award de-sc0019380 and contract DE-AC02-05CH11231; and by the National Science Foundation under grant PHY1820912. K.P. is supported in part by NSF grants PHY-1404105 and PHY-1707800 to Cornell University.

## Appendix A: Covariant phase space formalism and the Wald-Zoupas charges

The computation of boundary charges and their fluxes makes use of the covariant phase space formalism [1] and the Wald-Zoupas prescription [3]. We quickly review the main ingredients needed for null boundaries, and refer the reader to [27] for details.

↦

Consider a diffeomorphism covariant theory of the metric g ab as a dynamical field which is described by a Lagrangian 4-form L ( g ) that depends locally and covariantly on g ab . Under perturbations g → g + δg the Lagrangian changes as

$$\delta L = E ^ { a b } \delta g _ { a b } + d \Theta ( g ; \delta g )$$

where E ab is a 4-form presenting the equations of motion of the theory and the 3-form

θ ( g, δg ) is the presymplectic potential. The 3-form presymplectic current is defined by

$$\omega ( g ; \delta _ { 1 } g, \delta _ { 2 } g ) = \delta _ { 1 } \theta ( g ; \delta _ { 2 } g ) - \delta _ { 2 } \theta ( g ; \delta _ { 1 } g )$$

where δ g 1 and δ g 2 are any two independent perturbations.

Given a vector field ξ a , one can then show that

$$\omega ( g ; \delta g, \mathcal { L } _ { \xi } g ) = d [ \delta Q _ { \xi } - i _ { \xi } \theta ( \delta g ) ]$$

where the 2-form Q ξ is the Noether charge [1-3].

Consider now a null boundary N in the spacetime and a spacelike hypersurface Σ that intersects N at some cross-section S . Integrating Eq. A.3 we get

$$\underset { \Sigma } { \int } \omega ( g ; \delta g, \mathcal { L } _ { \xi } g ) = \underset { S } { \int } \delta Q _ { \xi } - i _ { \xi } \theta ( \delta g )$$

Two vector fields ξ a and ˜ ξ a are equivalent representatives of symmetries on N if ξ a | N = ˜ ξ a | N and the right-hand-side of Eq. A.4 evaluated with ξ a and ˜ ξ a are equal for all backgrounds g ∈ F , all perturbations δg within F and all cross-sections S of N . The boundary symmetries on N are then given by vector fields ξ a factored out by the above equivalence relation.

From the above identity it would be 'natural' to define a charge at S associated to a symmetry ξ a as a function Q S ξ [ ] on phase space so that

$$\delta Q _ { \xi } [ S ] = \underset { S } { \int } \delta Q _ { \xi } - i _ { \xi } \theta ( \delta g )$$

for all perturbations δg within F and all cross-sections S . Unfortunately, in general, the right-hand-side of Eq. A.5 is not integrable in phase space and no function Q ξ satisfying Eq. A.5 exists on the phase space. As shown in [3] the integrability condition for the existence of a charge Q S ξ [ ] for some symmetry ξ a is

$$0 = ( \delta _ { 1 } \delta _ { 2 } - \delta _ { 2 } \delta _ { 1 } ) Q _ { \xi } = - \underset { S } { \int } i _ { \xi } \omega ( g, \delta _ { 1 } g, \delta _ { 2 } g )$$

for all perturbations δ g, δ 1 2 g within F and all cross-sections S . This above criteria is not satisfied, even at null infinity in general relativity, except in very special cases [3].

Nevertheless, Wald and Zoupas [3] developed a prescription for defining a modified charge which is always integrable. Define a boundary presymplectic potential Θ ( g ; δg ) for the

pullback to N of the presymplectic current,

$$\overline { \omega } ( g, \delta _ { 1 } g, \delta _ { 2 } g ) = \delta _ { 1 } \Theta ( g, \delta _ { 2 } g ) - \delta _ { 2 } \Theta ( g, \delta _ { 1 } g )$$

where ω denotes the pullback to N . Then define the Wald-Zoupas charge (WZ charge) by

$$\delta \mathcal { Q } _ { \xi } [ S ] = \underset { S } { \int } \delta Q _ { \xi } - i _ { \xi } \theta ( \delta g ) + i _ { \xi } \Theta ( \delta g )$$

It can be shown using Eqs. A.4 and A.7 that δ Q ξ [ S ] is integrable in phase space. Thus Eq. A.8 determines a function Q ξ [ S ] up to a constant of integration on F , which can be fixed by choosing a reference solution g 0 such that Q ξ [ S ] ∣ ∣ ∣ g 0 = 0 for all symmetries ξ a and all cross-sections S .

The prescription given by Wald and Zoupas is to choose the 3-form Θ such that Θ ( g ; δg ) vanishes for all perturbation δg ab for any background g ab which is stationary, and to choose the reference solution g 0 to also be stationary. The consistency conditions for such choices and the ambiguities in them are detailed in [3].

If the above choices can be made then the flux F ξ [∆ N ] of the WZ charge Q ξ through a part of the null boundary N is given by [3]

$$\mathcal { F } _ { \xi } [ \Delta N ] = \underset { \Delta N } { \int } \, \Theta ( g ; \mathcal { L } _ { \xi } g )$$

It can also be shown that the perturbed flux δ F ξ for any symmetry ξ a and any perturbation δg ab satisfies (see Eq. 29 [3])

$$\delta \mathcal { F } _ { \xi } = \underset { N } { \int } \omega ( g ; \delta g, \mathcal { L } _ { \xi } g ) + \underset { \partial N } { \int } i _ { \xi } \Theta ( \delta g )$$

If i ξ Θ ( δg ) → 0 on ∂N for all perturbations δg ab then F ξ is a function on the covariant phase space F satisfying

$$\delta \mathcal { F } _ { \xi } = \underset { N } { \int } \omega ( g ; \delta g, \mathcal { L } _ { \xi } g )$$

for all perturbations δg ab , that is, F ξ defines a Hamiltonian which generates the flow on the covariant phase space F associated to the symmetry ξ a .

In [3] this procedure was applied to the asymptotic symmetries at null infinity in general relativity to derive the charges and symmetries for the BMS algebra. The case of finite null boundaries in vacuum general relativity was handled in [27], where it as shown that the notion of symmetries defined below Eq. A.4 coincides with those defined in Sec. 2 and the Wald-Zoupas prescription gives the charges and fluxes described in Eqs. 2.20 and 2.21.

In [27] the reference solution (used in the Wald-Zoupas prescription) was chosen to be the horizon of a Schwarzschild black hole in the limit that the mass tends to zero. It was shown that this reference solution satisfies all the criteria given by Wald and Zoupas [3]. In this paper, we simply adopt the formulae for the charges and fluxes from [27] and do not analyze the choice of reference solution in detail.

## Appendix B: Structure of g CD as a central extension

In this section we explore the structure of the summetry algebra g CD of the causal diamond as a non-trivial extension of g CD / b 0 by the boosts b 0 .

Recall from Sec. 3.1 that the elements of g CD ∼ = diff( S 2 ) glyph[multicloseleft] b are of the form ( β, X A ) where β is a smooth function and X A a vector feld on S 2 . The central elements (i.e. those which commute with all other elements in g CD ) of boosts in b 0 are the ones given by ( β = constant 0). Consider the quotient , g CD / b 0 which consists of equivalence classes given by the relation ( β, X A ) ∼ ( β +constant , X A ). Thus, g CD is a central extension of g CD / b 0 by the abelian algebra b 0 . We now show that this central extension is, in fact, a non-trivial central extension. What this means is the following:

Does the bracket of two representative elements belonging to g CD / b 0 , computed in g CD , have a non-vanishing b 0 -part?

If the answer is 'no' then the central extension is trivial and g CD will be a direct product of g CD / b 0 and b 0 . If the answer is 'yes' then g CD has the structure of a non-trivial central extension of g CD / b 0 by b 0 . 9

Since the symmetry algebra g CD is independent of the metric q AB on N the null boundary of the causal diamond, we can deduce its structure in any choice of metric, in particular it is convenient to choose q AB to be the metric of a unit-sphere in the standard ( θ, φ ) coordinates on S 2 . We now compute the bracket of any two elements in g CD / b 0 , with the only relevant case being the bracket between an element of b b / 0 with an element of diff( S 2 ) which gives (see Eq. 3.7)

$$[ ( \beta _ { 1 }, 0 ), ( 0, X _ { 2 } ^ { A } ) ] = ( \beta, 0 ) \quad \text{with} \quad \beta = - X _ { 2 } ^ { A } \partial _ { A } \beta _ { 1 }$$

To answer the above question we expand in terms of spherical harmonics Y l,m ( θ, φ ). Note that elements of b 0 are purely l = 0 spherical harmonics. Now let β 1 be a l 1 -harmonic with l 1 ≥ 1 and X A 2 be a l 2 -vector harmonic. We can write X A 2 = ∂ A X + /epsilon1 AB ∂ B ˜ X where X,X ˜ are l 2 -harmonics with l 2 ≥ 1. We want determine whether β can have a non-trivial l = 0 part, that is a non-vanishing constant piece.

9 In a more mathematical language, every central extension of g CD / b 0 by the abelian algebra b 0 corresponds to a 2-cocycle in the cohomology group H 2 ( g CD / b 0 , b 0 ) (see Sec. IV.2 [63]). Our computation in this section amounts to showing that the cocycle which gives the Lie algebra structure of g CD is non-trivial.

Before considering the general case, we note the following example

$$\beta _ { 1 } = \cos \theta \,, \ \ X _ { 2 } ^ { A } \equiv - \sin \theta \partial _ { \theta } \implies \beta = \sin ^ { 2 } \theta$$

Thus, β has non-vanishing l = 0 2 parts in terms of spherical harmonics. , This already shows that g CD is a non-trivial central extension of g CD / b 0 by b 0 .

For the general situation, first consider the case X = 0, X A 2 = /epsilon1 AB ∂ B ˜ . X Then we have, by integrating-by-parts on the unit-sphere (we leave the area element implicit for notational convenience)

$$\int \beta \, \overline { Y } _ { l = 0, m } \, \infty \, \int \beta = - \int \epsilon ^ { A B } \partial _ { B } \tilde { X } \partial _ { A } \beta _ { 1 } = - \int \partial _ { B } ( \epsilon ^ { A B } \tilde { X } \partial _ { A } \beta _ { 1 } ) = 0$$

Thus, β always has a vanishing l = 0 component for X A 2 = /epsilon1 AB ∂ B ˜ . X

Next consider the case ˜ = 0, X X A = ∂ A X we have

$$\int \beta \, \overline { Y } _ { l = 0, m } \, \infty \int \beta & = - \int \partial ^ { A } X \partial _ { A } \beta _ { 1 } = \int X \partial ^ { 2 } \beta _ { 1 } \\ & = - l _ { 1 } ( l _ { 1 } + 1 ) \int X \beta _ { 1 }$$

Expanding the functions X and β 1 in terms of the corresponding spherical harmonics, and using the orthonormality and completeness of the spherical harmonic basis we conclude that, the right-hand-side is non-vanishing if and only if

$$l _ { 1 } = l _ { 2 } \geq 1 \,, \quad m _ { 2 } = - m _ { 1 }$$

Thus, β has a non-vanishing constant ( l = 0) part whenever Eq. B.5 is satisfied, which has many solutions; taking l 1 = l 2 = 1 and m 1 = m 2 = 0 gives the above example Eq. B.2. Thus, g CD is a non-trivial central extension of g CD / b 0 by b 0 as we wished to show.

[1] J. Lee and R. M. Wald, Local symmetries and constraints , J. Math. Phys. 31 (1990) 725.

[2] V. Iyer and R. M. Wald, Some properties of Noether charge and a proposal for dynamical black hole entropy , Phys. Rev. D50 (1994) 846 [ gr-qc/9403028 ].

[3] R. M. Wald and A. Zoupas, A general definition of 'conserved quantities' in general relativity and other theories of gravity , Phys. Rev. D61 (2000) 084027 [ gr-qc/9911095 ].

[4] S. Carlip, Black hole entropy from conformal field theory in any dimension , Phys. Rev. Lett. 82 (1999) 2828 [ hep-th/9812013 ].

- [5] A. P. Balachandran, L. Chandar and A. Momen, Edge states in gravity and black hole physics , Nucl. Phys. B461 (1996) 581 [ gr-qc/9412019 ].

- [6] W. Donnelly and L. Freidel, Local subsystems in gauge theory and gravity , JHEP 09 (2016) 102 [ 1601.04744 ].
- [7] A. J. Speranza, Local phase space and edge modes for diffeomorphism-invariant theories , JHEP 02 (2018) 021 [ 1706.05061 ].
- [8] J. Kirklin, Unambiguous Phase Spaces for Subregions , JHEP 03 (2019) 116 [ 1901.09857 ].
- [9] A. Strominger, Lectures on the Infrared Structure of Gravity and Gauge Theory , 1703.05448 .
- [10] D. Harlow and J.-Q. Wu, Covariant phase space with boundaries , 1906.08616 .
- [11] A. C. Wall, A proof of the generalized second law for rapidly changing fields and arbitrary horizon slices , Phys. Rev. D85 (2012) 104049 [ 1105.3445 ], errata: Phys. Rev. D87 (2013) 069904.
- [12] T. Padmanabhan, General Relativity from a Thermodynamic Perspective , Gen. Rel. Grav. 46 (2014) 1673 [ 1312.3253 ].
- [13] F. Hopfmüller and L. Freidel, Null Conservation Laws for Gravity , Phys. Rev. D97 (2018) 124029 [ 1802.06135 ].
- [14] T. Jacobson, Thermodynamics of space-time: The Einstein equation of state , Phys. Rev. Lett. 75 (1995) 1260 [ gr-qc/9504004 ].
- [15] W. Wieland, New boundary variables for classical and quantum gravity on a null surface , Class. Quant. Grav. 34 (2017) 215008 [ 1704.07391 ].
- [16] R. Bousso, A Covariant entropy conjecture , JHEP 07 (1999) 004 [ hep-th/9905177 ].
- [17] W. Wieland, Generating functional for gravitational null initial data , 1905.06357 .
- [18] S. W. Hawking, M. J. Perry and A. Strominger, Soft Hair on Black Holes , Phys. Rev. Lett. 116 (2016) 231301 [ 1601.00921 ].
- [19] S. W. Hawking, M. J. Perry and A. Strominger, Superrotation Charge and Supertranslation Hair on Black Holes , JHEP 05 (2017) 161 [ 1611.09175 ].
- [20] S. Haco, S. W. Hawking, M. J. Perry and A. Strominger, Black Hole Entropy and Soft Hair , JHEP 12 (2018) 098 [ 1810.01847 ].
- [21] J. M. Bardeen, B. Carter and S. W. Hawking, The Four laws of black hole mechanics , Commun. Math. Phys. 31 (1973) 161.
- [22] J. D. Bekenstein, Black holes and entropy , Phys. Rev. D7 (1973) 2333.
- [23] S. Ryu and T. Takayanagi, Holographic derivation of entanglement entropy from AdS/CFT , Phys. Rev. Lett. 96 (2006) 181602 [ hep-th/0603001 ].
- [24] T. Jacobson, Entanglement Equilibrium and the Einstein Equation , Phys. Rev. Lett. 116 (2016) 201101 [ 1505.04753 ].
- [25] T. Jacobson and M. Visser, Gravitational Thermodynamics of Causal Diamonds in (A)dS , 1812.01596 .
- [26] J. de Boer, F. M. Haehl, M. P. Heller and R. C. Myers, Entanglement, holography and causal diamonds , JHEP 08 (2016) 162 [ 1606.03307 ].
- [27] V. Chandrasekaran, E. E. Flanagan and K. Prabhu, Symmetries and charges of general

relativity at null boundaries , JHEP 11 (2018) 125 [ 1807.11499 ].

- [28] A. Strominger, On BMS Invariance of Gravitational Scattering , JHEP 07 (2014) 152 [ 1312.2229 ].
- [29] C. Troessaert, The BMS4 algebra at spatial infinity , Class. Quant. Grav. 35 (2018) 074003 [ 1704.06223 ].
- [30] K. Prabhu, Conservation of asymptotic charges from past to future null infinity: Maxwell fields , JHEP 10 (2018) 113 [ 1808.07863 ].
- [31] K. Prabhu, Conservation of asymptotic charges from past to future null infinity: Supermomentum in general relativity , JHEP 03 (2019) 148 [ 1902.08200 ].
- [32] R. M. Wald, General Relativity . The University of Chicago Press, 1984.
- [33] R. Geroch, Asymptotic structure of space-time , in Asymptotic structure of space-time (F. P. Esposito and L. Witten, eds.). Plenum Press, New York, 1977.
- [34] A. Ashtekar and M. Streubel, Symplectic Geometry of Radiative Modes and Conserved Quantities at Null Infinity , Proc. Roy. Soc. Lond. A376 (1981) 585.
- [35] A. Ashtekar, Geometry and Physics of Null Infinity , 1409.1800 .
- [36] H. Bondi, M. G. J. van der Burg and A. W. K. Metzner, Gravitational waves in general relativity. 7. Waves from axisymmetric isolated systems , Proc. Roy. Soc. Lond. A269 (1962) 21.
- [37] R. K. Sachs, Gravitational waves in general relativity. 8. Waves in asymptotically flat space-times , Proc. Roy. Soc. Lond. A270 (1962) 103.
- [38] R. Sachs, Asymptotic symmetries in gravitational theory , Phys. Rev. 128 (1962) 2851.
- [39] S. Hollands and A. Thorne, Bondi mass cannot become negative in higher dimensions , Commun. Math. Phys. 333 (2015) 1037 [ 1307.1603 ].
- [40] S. Hollands, A. Ishibashi and R. M. Wald, BMS Supertranslations and Memory in Four and Higher Dimensions , Class. Quant. Grav. 34 (2017) 155005 [ 1612.03290 ].
- [41] R. Penrose, Techniques of differential topology in relativity , Regional Conference Series in Applied Math. 7. SIAM, Philadelphia, 1972.
- [42] L. Donnay, G. Giribet, H. A. Gonzalez and M. Pino, Supertranslations and Superrotations at the Black Hole Horizon , Phys. Rev. Lett. 116 (2016) 091101 [ 1511.08687 ].
- [43] E. Morales, On a Second Law of Black Hole Mechanics in a Higher Derivative Theory of Gravity , Ph.D. thesis, University of Gottingen, 2008.
- [44] K. Parattu, S. Chakraborty, B. R. Majhi and T. Padmanabhan, A Boundary Term for the Gravitational Action with Null Boundaries , Gen. Rel. Grav. 48 (2016) 94 [ 1501.01053 ].
- [45] T. Mädler, Affine-null metric formulation of general relativity at two intersecting null hypersurfaces , Phys. Rev. D99 (2019) 104048 [ 1810.04743 ].
- [46] A. Ashtekar and S. Bahrami, Asymptotics with a positive cosmological constant. IV. The no-incoming radiation condition , Phys. Rev. D100 (2019) 024042 [ 1904.02822 ].
- [47] E. Gourgoulhon and J. L. Jaramillo, A 3+1 perspective on null hypersurfaces and isolated

horizons , Phys. Rept. 423 (2006) 159 [ gr-qc/0503113 ].

- [48] I. Booth, Spacetime near isolated and dynamical trapping horizons , Phys. Rev. D87 (2013) 024008 [ 1207.6955 ].
- [49] E. Witten, Light Rays, Singularities, and All That , 1901.03928 .
- [50] C. Dappiaggi, N. Pinamonti and M. Porrmann, Local causal structures, Hadamard states and the principle of local covariance in quantum field theory , Commun. Math. Phys. 304 (2011) 459 [ 1001.0858 ].
- [51] Y. Choquet-Bruhat, P. T. Chrusciel and J. M. Martin-Garcia, The Cauchy problem on a characteristic cone for the Einstein equations in arbitrary dimensions , Ann. Henri Poincaré 12 (2011) 419 [ 1006.4467 ].
- [52] Y. Choquet-Bruhat, P. T. Chrusciel and J. M. Martin-Garcia, An Existence theorem for the Cauchy problem on a characteristic cone for the Einstein equations , in 4th International Conference on Complex Analysis and Dynamical Systems, Nahariya, Israel, May 18-22, 2009 , 2010, 1006.5558 .
- [53] E. Newman and R. Penrose, An Approach to Gravitational Radiation by a Method of Spin Coefficients , J. Math. Phys. 3 (1962) 566 errata: J. Math. Phys. 4 (1963) 998.
- [54] R. P. Geroch, A. Held and R. Penrose, A space-time calculus based on pairs of null directions , J. Math. Phys. 14 (1973) 874.
- [55] S. Carlip, Black Hole Entropy from Bondi-Metzner-Sachs Symmetry at the Horizon , Phys. Rev. Lett. 120 (2018) 101301 [ 1702.04439 ].
- [56] M. Guica, T. Hartman, W. Song and A. Strominger, The Kerr/CFT Correspondence , Phys. Rev. D80 (2009) 124008 [ 0809.4266 ].
- [57] S. Hollands and A. Ishibashi, News vs Information , 1904.00007 .
- [58] W. Wieland, Fock representation of gravitational boundary modes and the discreteness of the area spectrum , Ann. Henri Poincaré 18 (2017) 3695 [ 1706.00479 ].
- [59] X. Dong, D. Harlow and A. C. Wall, Reconstruction of Bulk Operators within the Entanglement Wedge in Gauge-Gravity Duality , Phys. Rev. Lett. 117 (2016) 021601 [ 1601.05416 ].
- [60] D. L. Jafferis, A. Lewkowycz, J. Maldacena and S. J. Suh, Relative entropy equals bulk relative entropy , JHEP 06 (2016) 004 [ 1512.06431 ].
- [61] N. Lashkari, J. Lin, H. Ooguri, B. Stoica and M. Van Raamsdonk, Gravitational positive energy theorems from information inequalities , PTEP 2016 (2016) 12C109 [ 1605.01075 ].
- [62] R. Bousso, V. Chandrasekaran and A. Shahbazi-Moghaddam, Ignorance is Cheap: From Black Hole Entropy To Energy-Minimizing States In QFT , 1906.05299 .
- [63] A. W. Knapp, Lie Groups, Lie Algebras, and Cohomology . Princeton University Press, 1988.
Venkatesa Chandrasekaran Kartik Prabhu
Symmetries, charges and conservation laws at causal diamonds in general relativity
We study the covariant phase space of vacuum general relativity at the null
boundary of causal diamonds. The past and future components of such a null
boundary each have an infinite-dimensional symmetry algebra consisting of
diffeomorphisms of the $2$-sphere and boost supertranslations corresponding to
angle-dependent rescalings of affine parameter along the null generators.
Associated to these symmetries are charges and fluxes obtained from the
covariant phase space formalism using the prescription of Wald and Zoupas. By
analyzing the behavior of the spacetime metric near the corners of the causal
diamond, we show that the fluxes are also Hamiltonian generators of the
symmetries on the phase space. In particular, the supertranslation fluxes yield
an infinite family of boost Hamiltonians acting on the gravitational data of
causal diamonds. We show that the smoothness of the vector fields representing
such symmetries at the bifurcation edge of the causal diamond implies suitable
matching conditions between the symmetries on the past and future components of
the null boundary. Similarly, the smoothness of the spacetime metric implies
that the fluxes of all such symmetries is conserved between the past and future
components of the null boundary. This establishes an infinite set of
conservation laws for finite subregions in gravity analogous to those at null
infinity. We also show that the symmetry algebra at the causal diamond has a
non-trivial center corresponding to constant boosts. The central charges
associated to these constant boosts are proportional to the area of the
bifurcation edge, for any causal diamond, in analogy with the Wald entropy
formula.