Local subsystems in gauge theory and gravity

We consider the problem of defining localized subsystems in gauge theory and gravity. Such systems are associated to spacelike hypersurfaces with boundaries and provide the natural setting for studying entanglement entropy of localized subsystems. We present a general formalism to associate a gauge-invariant classical phase space to a spatial slice with boundary by introducing new degrees of freedom on the boundary. In Yang-Mills theory the new degrees of freedom are a choice of gauge on the boundary, transformations of which are generated by the normal component of the nonabelian electric field. In general relativity the new degrees of freedom are the location of a codimension-2 surface and a choice of conformal normal frame. These degrees of freedom transform under a group of surface symmetries, consisting of diffeomorphisms of the codimension-2 boundary, and position-dependent linear deformations of its normal plane. We find the observables which generate these symmetries, consisting of the conformal normal metric and curvature of the normal connection. We discuss the implications for the problem of defining entanglement entropy in quantum gravity. Our work suggests that the Bekenstein-Hawking entropy may arise from the different ways of gluing together two partial Cauchy surfaces at a cross-section of the horizon.


Introduction
Entanglement plays an increasingly prominent role in fundamental physics. It has long been recognized that the entanglement naturally present in the vacuum could account for the Bekenstein-Hawking entropy of black holes [1][2][3] (see [4] for a review). In fact it is conjectured that the scaling of entanglement entropy with area is a universal feature of near-vacuum states in quantum gravity [5]. More ambitiously, it has been proposed that the geometry of spacetime [6] and its topology [7] may emerge from entanglement. Moreover, if the Bekenstein-Hawking entropy really is entanglement entropy, then all regions of space should possess an entropy and therefore be subject to the laws of thermodynamics. Applying the laws of thermodynamics to local Rindler horizons [8], or to small spherical regions of spacetime [9], one can derive the Einstein equation as an equation of JHEP09(2016)102 state. Together these results point toward idea that entanglement, an essentially quantum phenomenon, is deeply connected to classical spacetime geometry and its dynamics. Despite the central role that entanglement is conjectured to play in quantum gravity, there is so far no sharp definition of entanglement between regions of space that would apply in quantum gravity. In order to define entanglement entropy, we require that for any Cauchy surface divided into two disjoint regions Σ and Σ, there is a tensor product factorization of the total Hilbert space (1.1) We can then define for any density matrix ρ, a reduced density matrix ρ Σ := Tr Σ ρ and an entanglement entropy We do not expect the Hilbert space of quantum gravity to be naturally equipped with a tensor product structure such as (1.1). An essential property of the tensor product is that operators within the different factors commute: [A ⊗ I, I ⊗ B] = 0 for any operators A and B. There are strong arguments that such locally supported commuting operators simply cannot exist in quantum gravity [10,11]. The basic obstruction is diffeomorphism invariance: physical operators in quantum gravity should be diffeomorphism invariant, but the action of a diffeomorphism is to permute the points of spacetime. Thus an invariant operator with support at any spacetime point must also have support on every other spacetime point. Whatever notion of local subsystems survives in quantum gravity must clearly be a generalization of (1.1). The goal of this paper is to unravel the mechanism behind such a generalisation and the basis for its construction. The absence of tensor factorisation is expected for any theory with gauge symmetry. In a gauge theory the initial data must obey constraints, which take the form of differential equations on each time slice. In Yang-Mills theory the constraint is Gauss' law, and in gravity it is the diffeomorphism constraints (often further divided into a spatial diffeomorphism constraint and a timelike Hamiltonian constraint). Because of these constraint equations, initial data in the two regions Σ and Σ cannot be specified independently. This fact is reflected at the quantum level by the failure of the factorisation (1.1).
Before proceeding to describe our construction, we discuss some alternative perspectives on the problem of entanglement entropy and localized subsystems in field theory. The most straightforward way to define entanglement entropy for a field theory is to regularize on a lattice. There each site carries a Hilbert space, which could be for example that of a harmonic oscillator, or a finite-dimensional system. Then the Hilbert space of each region is simply the tensor product of all the sites inside the region, and the factorisation structure (1.1) is manifest. One can therefore define entanglement entropy for any quantum field theory that arises as a continuum limit of such a lattice system. Although the entanglement entropy defined in this way is divergent and depends on the precise lattice realization of the theory, certain universal quantities can be extracted from it. One might hope to define entanglement entropy in quantum gravity via such a bottom-up approach. Of course, this requires gravity to emerge from a theory with a local tensor product struc-JHEP09(2016)102 ture, as in string-net models [12]. It has recently been suggested that the bulk gravity theory in AdS/CFT must have an ultraviolet completion that does factor in this way [13].
In such a lattice model, gauge invariance is not fundamental. Instead, it emerges only within a low-energy sector of the Hilbert space. This low-energy sector does not admit a factorization into local subsystems; only the full ultraviolet completion does. Entanglement entropy of the low-energy gauge-invariant states is then defined by embedding them into a larger Hilbert space which does factorize. We will see that a similar embedding can be defined directly in the gauge theory, without specifying a lattice regularization.
An alternative approach is to identify a region of spacetime by its algebra of local observables, as is standard in the algebraic approach to QFT [14]. One can define the entropy of a region to be the entropy of its observable algebra, as advocated in ref. [15]. This approach to subsystems is problematic for several reasons [11]. First, in quantum field theory the algebra of observables associated to a spacetime region is a type III von Neumann algebra [16]. Unlike a type I von Neumann algebra, a type III von Neumann algebra is not associated with a factorization of the Hilbert space. Since the type III property depends on the ultraviolet behavior of quantum field theory, it is not clear whether this property will persist in quantum gravity. A much more serious issue is the absence of compactly supported gauge-invariant observables in gravity [10]. Thus it would seem that the algebraic approach, which is based on algebras of compactly supported gauge-invariant observables, does not know how to assign an entropy to systems without local gauge-invariant operators.
Instead we will follow an approach that has been called the extended Hilbert space construction [17][18][19][20][21]. We can divide the constraints of the theory into three classes: those with support in Σ, those with support in Σ, and those with nontrivial support on the boundary. To the region Σ we assign an extended Hilbert space H Σ , consisting of states that satisfy the gauge condition away from the entangling surface S = ∂Σ, i.e. which are annihilated by the constraints with support only inside Σ. These states of the extended Hilbert spaces H Σ and H Σ generally transform nontrivially under transformations with nontrivial support on the boundary, so the Hilbert space carries a nontrivial representation of the boundary gauge group. We can then identify the physical Hilbert space as a subspace of the tensor product H ⊂ H Σ ⊗ H Σ . (1. 3) The space H consists of those states that are invariant under boundary gauge transformations. Equivalently, it is the space of states annihilated by the constraints with nontrivial support on the boundary. The extended Hilbert space construction can be defined precisely in Yang-Mills theory on a lattice [18] and in two dimensions [22]. Here the extended Hilbert space H Σ consists of Wilson loops (or more generally, generalized spin network states [18]) which are allowed to have open ends on the boundary S. This is in a sense the minimal extension of the Hilbert space: in order for H to contain gauge-invariant operators such as a Wilson loop in a representation R, the space H Σ has to contain a state that looks like a Wilson loop in the interior of Σ, but which terminates on S. The endpoint of the Wilson line acts as a surface charge transforming in the representation R. The generators of this gauge transformation JHEP09(2016)102 are the nonabelian electric fields along links that cross the boundary. The subspace H consists of the gauge-invariant states for which the surface charges associated on one side of the boundary match up with the surface charges on the other side to form singlets. Equivalently, the normal component of the electric field must match across the boundary.
The result of this construction is that the gauge symmetry of the boundary is promoted to a physical symmetry on H Σ . The entanglement entropy that results from the embedding (1.3) includes a contribution from the surface charge degrees of freedom. While it may seem like overcounting to count gauge-related states in the entanglement entropy, multiple lines of evidence point to this being the physically relevant definition of entanglement entropy. For example, in ref. [22] it is shown that such states are necessary for consistency of the thermal description of the de Sitter vacuum. Additional evidence for this is provided by the calculation in ref. [23], where the entanglement entropy of a quark-antiquark pair was found to contain a contribution from surface charges where the nonabelian electric string pierces the entangling surface. In abelian gauge theory, the surface charges are needed to provide agreement with Euclidean methods, and to obtain agreement between the conformal anomaly and the universal logarithmic term in the entanglement entropy [24,25]. Thus the surface charge degrees of freedom are an essential component of the horizon entropy.
Our goal is twofold. First, we want to revisit and improve the previous construction of the extended Hilbert space in gauge theory by analysing its classical analog, the extended phase space. Second, we would like to define an analogous extended phase space for gravity.
We first show how to assign to each region of space Σ an extended phase space P Σ , which contains additional degrees of freedom at the boundary. This extended phase space is the classical analog of the extended Hilbert space. Our analysis shows that gauge invariance dictates what type of additional degrees of freedom are needed and their commutation relations. These additional degrees of freedom are designed to make the phase space invariant under all gauge transformations of the original system, whether or not they vanish at the boundary. In addition to the gauge symmetries extended to the boundary, we discover that the boundary degrees of freedom transform covariantly under a new group of physical symmetries, which we call surface symmetries, that are implemented as Hamiltonian transformations on the extended phase space. The appearance of a new type of boundary symmetry group G S where S = ∂Σ as a result of gauge invariance is one of the key results of this paper. We find, quite remarkably, that the boundary symmetry group commutes with the gauge transformations; this implies that the charges generating G S are gauge-invariant boundary observables. The second goal of this paper is to identify these charges for Yang-Mills and gravity. At the quantum level this implies that the extended Hilbert space H Σ carries a representation of G S .
We can use this Hamiltonian group action to define the the phase space associated to a union of regions Σ ∪ Σ. It is given by a fusion product where we take the Cartesian product of the phase spaces (analogous to the tensor product of Hilbert spaces) and mod out by the action of the surface symmetry group. At the

JHEP09(2016)102
quantum level this construction is implemented by an entangling product where the entangling product ⊗ G S is a tensor product of modules: H Σ acts as a right module for G S and H Σ acts as left module for G S and we make the identification (g |ψ ; g |ψ ) ≃ (|ψ ; |ψ ), for g ∈ G S . In other words the entangling product restricts to states that are singlets under the diagonal action of surface symmetries. We begin by constructing an extended phase space for Yang-Mills theory with a general gauge group G in section 2. We follow the covariant canonical formalism, which has been developed by many authors, and presented e.g. in [26]. The key step is to include the choice of gauge as a dynamical variable in the phase space. 1 In the mathematical language this is a choice of bundle trivialization, but locally it is simply a G-valued function ϕ : S → G. Following the standard algorithm we define a presymplectic potential form θ, which can be integrated over a Cauchy surface Σ to define the presymplectic potential of the theory. If the Cauchy surface Σ has a boundary, the presymplectic potential obtained from this procedure is not gauge-invariant. We show how it can be made invariant by adding a boundary term that depends on ϕ and its conjugate momentum: the non-abelian surface charge E, which is a Lie algebra valued density. When this is done we obtain a presymplectic potential that is fully invariant under gauge transformations, including those which act at the boundary.
The introduction of the additional variable ϕ allows us to introduce a new symmetry of the system, which rotates ϕ without changing the other fields. Thus we distinguish between two different types of transformations: Gauge symmetry: Surface symmetry: The former is pure gauge, acts on the right of the trivialisation, and has vanishing generator H[α]. This generator can be expressed as a bulk term proportional to the Gauss constraint, and a boundary term that enforces the identification of the normal electric field with the surface charge (see section 2.5). The latter is a symmetry, acts on the left of the trivialisation, and its generator E[α] is nonvanishing and is given by the the non-abelian surface charge. The vanishing of the boundary component of the gauge generator H[α] can be understood as a boundary condition that identifies the non-abelian surface charge density E with the normal component of the electric field in the frame ϕ. Explicitly this reads n a is the unit time normal to the slice Σ and s b is the unit normal of S within Σ, F is the Yang-Mills curvature. These normal electric fields therefore generate an algebra of

JHEP09(2016)102
non-abelian surface charges. The Poisson brackets of this extended phase space coincide with the commutators of the extended Hilbert space of lattice gauge theory; we have thus succeeded in constructing a continuum version of the lattice construction [18]. The most significant advantage of our covariant formalism is that it extends naturally to general relativity, as we show in section 3. Here in addition to a metric g ab , we introduce a variable X a : R D → M that specifies the coordinate location of the entangling surface S = X(s), where s is a fixed (D−2)-dimensional closed surface in R D . X a plays a role analogous to ϕ which specifies the choice of boundary gauge in Yang-Mills. The non-temporal variable X a drops out of the presymplectic potential when integrated over a manifold without boundary. As shown in [28] the temporal component δT does not in general, unless we restrict to hypersurface normal deformations. On a surface Σ with boundary, we show that the addition of an X a -dependent boundary term makes the presymplectic potential diffeomorphism-invariant, even for diffeomorphisms supported on the boundary.
The additional phase space variable X a again leads to physical symmetries of the extended phase space. Thus again we distinguish between two types of transformations: Surface symmetry: V is a vector field on M that labels infinitesimal diffeomorphism of spacetime: these are gauge. w is a vector field on R D that generates an infinitesimal change of the surface S: these are the surface symmetries. The gravitational case is more subtle than Yang-Mills because of the different ways that a transformation of X a can act in a neighbourhood of the codimension-2 surface S. Thus we distinguish between three classes of surface transformations: • Surface boosts linearly transform the normal plane of S in a position-dependent way.
• Surface diffeomorphisms are reparametrizations of S that move S tangent to itself.
• Surface translations move the surface S in the normal direction.
Previous works [29,30] and especially [31][32][33] (see [34] for a more recent overview) have focused on the surface translations. In order to implement the surface translations as canonical transformations, boundary conditions are usually imposed at the surface S, and the correct way to choose those boundary conditions remains somewhat mysterious. Here we focus instead on the first two classes of transformations, which we call the surfacepreserving symmetries. These are the boundary symmetries that preserve the location of the surface as a whole. We find canonical generators for these symmetries, without the need to specify the location of the surface or any associated boundary conditions. However it is clear that the full story must include both classes of transformations, especially if we are to consider dynamically evolving entangling surfaces. While we hope that our formalism will lend insight to the quantization of the surface translations, we leave this task for future work.

JHEP09(2016)102
The central result of this paper is twofold. First, we identify the surface-preserving symmetry group which is given by the semi-direct product of diffeomorphisms of the entangling surface, times the algebra of surface boosts represented as group valued maps S → SL(2, R). We also show that the action of this symmetry group is Hamiltonian on the extended phase space of gravity and identify the Hamiltonian generators in terms of intrinsic and extrinsic geometry of S. These Hamiltonian generators play the role of gravitational surface charges for the entangling surface S. They are the analog of the dressed normal electric field (1.8) for Yang-Mills. Remarkably they can be described explicitly as follows.
To evaluate the Hamiltonian in terms of components of the metric we introduce a 2 + (D − 2) decomposition and adapted coordinates (x i , σ µ ). σ µ are coordinates on S, and (x i ) i=0,1 are coordinates in the two normal directions. In these coordinates we can parametrize the metric as Here q µν is the induced metric on S, h ij a generalized lapse which defines the normal geometry and A µ i a generalized shift. A µ i can be viewed as a connection on the normal bundle of S. The generators of diffeomorphisms and SL(2, R) symmetries are integrals of surface charge densities, given by the normal curvature and the conformal normal geometry, respectively: Here √ q is the volume form on S, while √ h is the square root of the determinant of the normal metric, and ǫ ij is the two-dimensional Levi-Civita symbol.
We see that the algebra of surface boosts is generated by conformal components of the metric normal to S and form a local SL(2, R) algebra whose Casimir is the area element of the surface S. The diffeomorphisms are generated by the curvature of the normal connection, which vanishes if the normal planes are integrable. This is our main result. This result, along with the identification of the fusion product of phase space and the entangling product, form the basis of understanding gauge systems with boundaries.
While quantization is beyond the scope of this work, our approach shows that the Hilbert space associated with Σ must carry a unitary representation of the boundary symmetry group G S . In section 4 we mention some general results in representation theory which may be useful. A key point is the identification of the area element of the surface with the Casimir of SL(2, R). The representation theory of this algebra therefore singles out the area of the entangling surface S as the physical quantity that labels the dimension of the irreducible representations, and area preserving diffeomorphisms as a key subgroup to represent unitarily. We conclude in section 5 with many areas for future work. Some technical details of our formalism have been exiled to two appendices.

JHEP09(2016)102
2 Yang-Mills symplectic potential We first construct an extended phase space for bounded regions in Yang-Mills theory. In addition to being of interest in its own right, Yang-Mills theory provides a useful warmup exercise for gravity and allows us to review some aspects of the covariant canonical formalism. Moreover, it will allow us to see that the method we have developed is a continuum generalization of the extended Hilbert space approach to lattice Yang-Mills theory [17][18][19][20][21].
Consider Yang-Mills theory with gauge group G. The field variable is a connection A, which can be represented in a trivialization as a g-valued 1-form. The Lagrangian density is given by the D-form where ⋆ denotes the D-dimensional Hodge duality operation 2 and F = dA + 1 2 [A, A] is the curvature two-form.
The first step of deriving the canonical formalism is to find the symplectic potential of the theory, which we do following the method of refs. [26,[35][36][37][38]. The space of classical solutions of a Lagrangian system is naturally equipped with a presymplectic structure. The symplectic potential itself is obtained from an on-shell variation of the action. We introduce a differential δ on field space which satisfies Leibniz rule and squares to zero: δ 2 = 0. It acts on the space of solutions of the Lagrangian system mapping infinitesimally solutions to solutions and it is such that the product of field differentials anti commute δAδB = −δBδA for A, B two scalar fields. Following the conventions of ref. [36] we do not include an explicit wedge product between the 1-forms δ g ; we reserve the notation ∧ for spacetime forms and keep the antisymmetrization of δ implicit. This field space differential δ should not be confused with the spacetime differential denoted d. A general field is a (p, q) form where p denotes the degree of the form in field space and q the degree as a spacetime form.
We begin by computing the symplectic potential current density θ by varying the Lagrangian, θ depends on the field A and is linear in the variation δA; it is a (1, D − 1)-form, a oneform in field space in the sense of ref. [36] and a (D − 1)-form in spacetime. To obtain the presymplectic potential of the theory, we integrate θ over a Cauchy surface. The expression δA is analogous to the finite-dimensional expression dq where q is a canonical coordinate, and the integral of the current θ is analogous to the finite-dimensional symplectic potential i p i dq i .

JHEP09(2016)102
It is important to note that θ depends explicitly on A, and is not gauge-invariant. This is potentially a problem, because we would like to define θ as a one-form on the gauge orbits, and this is only possible if θ is gauge-invariant. As we will see this poses no problem when θ is integrated over a closed surface Σ without boundary, but when Σ has a boundary the symplectic structure does change nontrivially under gauge transformations. We will then show how to correct for this non-invariance by addition of a boundary term to the symplectic structure.
Let g * denote the gauge transformation by g, which acts on A and F as (2.5) Now consider how θ changes under this gauge transformation. After transformation the new symplectic potential current density is given by We now introduce some notation that will clarify this calculation and generalize more readily to the gravitational case. The quantity δg is a somewhat unnatural object: it is an element of the tangent space to G at the point g. A more natural object is obtained by translating it back to a neighbourhood of the identity: δ g is an element of the Lie algebra, and the transformation of the presymplectic potential form (2.6) can be expressed in terms of the gauge covariant derivative of δ g : The presence of the second term in eq. (2.9) shows precisely the degree to which the symplectic potential current density is not gauge invariant. The reason that θ is not gauge invariant is that the variation δ does not commute with the gauge transformation g * . Instead the difference can be expressed as an infinitesimal gauge transformation by δ g : Note that δA being a difference of two connections transforms as g * (δA) = g δA g −1 .
Similarly, for the field strength F we have δg * (F ) = g * (δF + [δ g , F ]), (2.11) which is again the transformation rule for F under an infinitesimal gauge transformation by δ g .

JHEP09(2016)102
Later we will need the variation of δ g itself, which is given by Recall that the two δ's in this formula correspond to two different (anticommuting) variations.
To obtain the presymplectic potential Θ of the theory, we integrate the symplectic potential current over a Cauchy surface Σ: (2.13) This is not manifestly gauge invariant because of the nontrivial transformation law for θ in eq. (2.9). But the latter term can be integrated by parts, using the equation of motion d A ⋆F= 0, with= denoting equality on shell. We find that the presymplectic potential changes by a boundary term, (2.14) When θ is integrated over a Cauchy surface without boundary, or when the boundary conditions are such that the boundary term in (2.14) vanishes, the presymplectic potential is gauge invariant. However we are interested in associated a phase space with a region Σ whose boundary S is the entangling surface. In this case we do not wish to impose any particular boundary conditions at S, since the entangling surface is not a physically distinguished location as far as the gauge field is concerned.
In order to define a gauge-invariant symplectic potential we need to add to Θ Σ a term whose gauge transformation compensates that in eq. (2.14). To construct such a potential we recall that a connection on a bundle is not specified just by the 1-form A but also by a local trivialization ϕ, which we may think of, at least locally, as a G-valued function ϕ : S → G. Under a gauge transformation ϕ transforms as g * (ϕ) = ϕg −1 . We can now augment the symplectic structure with the term where analogous to δ g we define δ ϕ := ϕ −1 δϕ. Under a gauge transformation δ ϕ transforms as which leads to the gauge transformation law for Θ S The second term in (2.17) is exactly what is needed to cancel the second term in (2.14).
Combining these two terms we obtain the presymplectic potential

JHEP09(2016)102
This presymplectic potential is gauge invariant for any spacelike surface Σ, even for gauge transformations with support on the boundary.

Gauge invariance
The above presymplectic potential satisfies, by construction, the gauge invariance property Note that this depends not only on the gauge field A everywhere on Σ (including the boundary S), but also on the trivialization ϕ on S. Because Θ is gauge-invariant, we can pull it back to define a symplectic potential on the space gauge orbits, which is the true physical phase space of the theory. This means that for any α valued in the Lie algebra g, the infinitesimal gauge transformation is a null direction of the total symplectic form, for any α including those with support on the boundary. We will now show explicitly that the canonical generator for this transformation vanishes on-shell. To find the generators of these transformations, we take the variation of the symplectic potential to obtain the presymplectic variational 2-form Ω = δΘ. This presymplectic form is the sum of a bulk term Ω Σ and a boundary term Ω S , where the latter depends on the boundary degree of freedom ϕ: We are now interested in constructing the generators of gauge transformations associated with the bulk and boundary symplectic structure in order to determine the Poisson brackets between boundary observables. Recall that a Hamiltonian H determines a variational vector field δ H via Since we have already introduced a field differential δ it is natural to also introduce a field "interior product": given any infinitesimal field transformation δ α A, we denote I δα the inner product along this variation that maps n-form in field space to (n − 1)-forms. I δα acts trivially on functionals of A (field space 0-forms), satisfies the graded Leibniz rule (see appendix A) and is defined on field one-forms by The Poisson bracket (2.24) can then be expressed in terms of the interior product and the symplectic structure Ω via

JHEP09(2016)102
where I is the interior product. Thus the Hamiltonian H and the variational vector field it generates are related by We can use the relation (2.27) to find the generator of gauge transformations. Let α be a Lie-algebra-valued function on Σ and δ α denote the gauge transformation (2.21). We introduce the bulk and boundary Hamiltonian generators H Σ [α] and H S [α] respectively such that I δα Ω Σ = δH Σ [α], and I δα Ω S = δH S [α]. One finds that they are given by For instance we can explicitly check the boundary variation We have used that the second term vanishes since (ϕ −1 δ α ϕ) = α. We can now check that vanishes on-shell: This is the expression of the gauge invariance of our formulation.

Formal derivation of gauge invariance
Since we would like to generalize our results to the case of gravity, and since direct computations in gravity such as the one done in the previous section are very cumbersome, it is convenient to develop a formal but simpler proof of the gauge invariance of the presymplectic potential. Such a proof will then generalize easily to the gravity case and also be conceptually clearer. Moreover, this proof elucidates the key mechanisms behind the restoration of gauge invariance by the addition of boundary degrees of freedom.
Using the field interior product we can define a field Lie derivative 3 L δα := δI δα + I δα δ. (2.30) We can now evaluate The symplectic potential we have constructed is invariant under gauge transformations, which was the main motivation to introduce the trivialisation degree of freedom ϕ. Indeed the infinitesimal version of (2.20) when This shows that the Hamiltonian is simply given by H[α] = −I δα Θ. This can be evaluated directly Here we recognize the contraction of the Yang-Mills equation of motion with the gauge parameter α pulled back to Σ. As we will see the same derivation applies in the gravity case.

Boundary observables
The addition of ϕ to the phase space allows us to construct new gauge-invariant operators. These are operators that depend on the choice of trivialization on the boundary. One is the electric field normal to the boundary, ϕ⋆F ϕ −1 . Another is a Wilson line that begins and ends on the boundary, where γ is a curve with endpoints x 1 and x 2 on the boundary. Let us now focus on the gauge invariant normal electric field observable: Its variation is given by where we define the variational vector field ∆ α such that I ∆α δ ϕ = −ϕ −1 αϕ. ∆ α acts only on the boundary trivialisation, so we have therefore generates a local rotation of the boundary trivialisation. This rotation acts on the left of the trivialisation, unlike the gauge transformation which acts on the right. This equation is equivalent to the statement that the Poisson bracket of {E[α], ·} = ∆ α with a boundary observable generates the transformation ∆ α . Accordingly the previous equation determines the Poisson brackets between E α 's via Thus we find that the different components of E := ⋆(ϕF ϕ −1 ) pulled back on S do not commute, but satisfy the current algebra of the group G. Expressed in terms of a set of This is the current algebra of the group G, i.e. the algebra satisfied by nonabelian gauge charges on the surface S. Thus the dressed normal electric field E acts as a surface charge.

Gauge versus symmetry
As a summary, we see that the introduction of the boundary degree of freedom ϕ has two effects. First, it restores the gauge invariance of the theory. Second, it introduces a new boundary symmetry which rotates the boundary observables (ϕ, E). The gauge transformations acts on the bulk variables and on the right of the boundary trivialisation vanishes on-shell. The boundary symmetry leaves the bulk variables invariant and acts on the left of the boundary trivialisation

Extended phase space and gluing
We are now able to define precisely the extended phase space P Σ associated with an open region Σ and the gluing property of phase space associated with two regions Σ L , Σ R . A key notion in order to describe these two procedures is the concept of symplectic reduction. Consider a phase space (P, Ω) that admits a group action: G × P → P. This action is Hamiltonian if it preserves the symplectic structure. For an infinitesimal transformation δ V this translates to the condition From the Cartan definition of L δ V and the closure of Ω we conclude that δI δ V Ω = 0, hence locally in field space we have I δ V Ω = δH V . The action is said to be globally Hamiltonian when the identity I δ V Ω = δH V is valid globally on field space. 4 H is the Hamiltonian of the group action. Given such a Hamiltonian group action we can define the symplectic quotient which is denoted by a double quotient labelled either by the group or by the set of Hamiltonians: The procedure is well-known from the theory of constrained Hamiltonian systems [39]. First one considers the constraint subspace C H ⊂ P on which the constraint H = 0 is imposed. C H is naturally equipped with a presymplectic structure Ω H = i * H Ω obtained by the pullback of the embedding map i H : C H ֒→ P. To obtain a symplectic manifold we need to quotient out the action of the group G, so that the phase space is then given by C H /G. In the more general case where we are given a family of Hamiltonians {H} we JHEP09(2016)102 only need to divide out C H by the kernel Ker(Ω H ), so that the phase space is P//{H} := C H /Ker(Ω H ). The advantage of this description is that it treats first and second class constraints simultaneously. The vector fields in Ker(Ω H ) are the infinitesimal generators of a group of transformations we denote G H . The double quotient notation refers to the fact that the phase space is obtained by a double operation: first a restriction onto the kernel C H of constraints, and then then a quotient of this set by the group G H generated by the Hamiltonian vectors preserving C H . The group G H and the Hamiltonians H therefore play a dual role in this double quotient.
We can now use the symplectic reduction to define a classical analog of the entangling product. Suppose that we have two symplectic spaces (P 1 , Ω 1 ) and (P 2 , Ω 2 ) and consider the constraints H 12 := H 1 − H 2 = 0, which generate the group of transformations G 12 . We can then define the fusion product This fusion product or classical entangling product identifies the two Hamiltonians H 1 = H 2 and divides out by their flow. Thus in the fusion product the phase space variables are those acting on the space in which H 1 = H 2 that are also invariant under its action: observables We can now describe precisely the extended phase space for Yang-Mills. This description formalises the construction that was first given in [40], in the context of first order gravity. One starts with a bulk phase space (P Σ ,Ω Σ ), whose phase space variables are the pull-back of (⋆F, A) on Σ and a boundary phase space (P S ,Ω S ) whose phase space variables are given by the pairs (E, ϕ) on S. Here E is a (D − 2)-form on S valued in the Lie algebra and ϕ is a G-valued function on S. The symplectic structures are given bŷ On the entangling surface S we recognise the natural symplectic structure of the group cotangent bundle T * (G S ) ≃ (G⋉g) S . The fusion of these symplectic structures is obtained by imposing the boundary condition for α ∈ g S . The extended phase space associated with Σ with boundary ∂Σ = S is simply the fusion product What is remarkable is the fact that the boundary condition (2.48) which determine the fusion product is in fact dictated by gauge invariance. Indeed as we have seen the generator of gauge invariance H[α] contains two contributions, one from the bulk and one from the boundary. The boundary part of the Hamiltonian generator is given by (2.48)

JHEP09(2016)102
The gluing constraints (2.48) can therefore be understood as the boundary component of the Gauss law and results from demanding that H[α] = 0 which imposes gauge invariance. Note that in this language the generator of boundary symmetries G S is simply We can now describe precisely the gluing procedure of two regions Σ L and Σ R sharing a common boundary S = Σ L ∩ Σ R . It is given by the fusion product associated with the constraints This constraint generates a diagonal action of the boundary symmetry group on P Σ L and P Σ R . The extended phase space associated with the region Σ L ∪ Σ R is simply given by the fusion product We have described the procedure of extension and gluing at the classical level in terms of phase space. This procedure has a direct, and in some ways simpler, analog at the quantum level. At the quantum level the phase space P Σ is replaced with a Hilbert space H Σ . The quantum analog of the symplectic quotient P//G is simply the quotient H/G of the Hilbert space by the group action. A simple quotient is enough since demanding invariance under the action of G is equivalent to imposing the Hamiltonian constraint H |ψ = 0, since the Hamiltonian H generates the action of G. That is at the quantum level we have the equality Ker(H) = H/G, which makes the symplectic reduction a conceptually simpler operation. 5 Accordingly, the extended Hilbert space is simply given by Similarly the gluing of two extended Hilbert spaces is obtained by considering the space of singlets under the boundary symmetry, that is by considering the entangling product

Gravity symplectic potential
We now turn to general relativity. The formalism developed in section 2 generalizes with some additional features. We will see that again the symplectic potential can be made gauge-invariant by addition of a boundary symplectic potential. Just as the boundary symplectic potential in Yang-Mills depends on the choice of gauge at the boundary, the symplectic potential of general relativity depends on a choice of coordinates describing the location of the boundary.

JHEP09(2016)102
In gravity we have the Lagrangian density (in units where 8πG = 1) where R is the Ricci scalar of the metric g ab , ǫ := the volume form, and s g is a sign factor that is −1 when g ab is Lorentzian, and +1 when g ab is Euclidean. The symplectic potential current density is given by -form defining the Hodge duality. Here we used the abbreviated notation δg = g ab δg ab , and δg ab = g ac g bd δg cd (i.e. δg ab is the variation of the metric with the indices raised, not the variation of the inverse metric, which is given by δ(g ab ) = −δg ab ). This follows from the direct calculation of the variation of the Ricci tensor given by In the following we will use three notions of equality: = when equality holds off-shell,= when the vacuum equation of motion for the background field is satisfied G ab= 0 but the variation δg ab is arbitrary and . = when the variation is also restricted to be on shell δG ab . = 0. We now consider how θ transforms under the gauge symmetries of general relativity, the diffeomorphism group. Let Y : M → M be a diffeomorphism of spacetime and denote by Y * : T * M → T * M the pullback under this diffeomorphism. We also denote by Y * : T M → T M the push forward. Since Y is invertible, the pullback map can be applied to vectors by defining Y * = (Y −1 ) * on vectors and then to tensors with both upper and lower indices (see appendix C of [42]). In particular, the pullback of the metric is given by Just as in Yang-Mills, the variation δ does not commute with the pullback Y * due to terms that involve the variation of Y . Instead we have the following relation: where L V denotes the spacetime Lie derivative along the vector field V and where we have introduced the vector field Eq. (3.6) is the generalization of the equations (2.10) and (2.11) of Yang-Mills to the diffeomorphism group. The relation (3.5), which we will use repeatedly, is proved in appendix B for general tensor-valued variational forms. Now using (3.5) we see that under a diffeomorphism of the underlying fields θ transforms as: (3.7)

JHEP09(2016)102
This is not simply the pullback of the symplectic potential, because of the second term in (3.7). However we can evaluate it on shell: In the last line we evaluate the symplectic potential on-shell, and in the hatted equality= we have made use of the vacuum Einstein equation R ab = 0. For the last term we use the identity where we have used the notation g for the operator of lowering an index on a vector to obtain a 1-form, gV = V a dx a , and ⋆ is the Hodge star acting on 2-forms 6 which is given by Here ǫ ab = i ∂ b i ∂a ǫ is the (D−2)-form obtained by double contraction with the volume form.
Note that the previous equality can be summarised by introducing a variational vector field I V δg = L V g. Then the equality becomes . (3.12) This is an identity that we will repeatedly use. Here π g [V ] is a (D − 2)-form which depends linearly on V and defined to be π g [V ] := 1 2 ⋆ dgV. (3.13) This shows that under diffeomorphism the symplectic potential current changes by shift under Y and a total derivative: (3.14) We now consider the symplectic potential Θ Σ obtained by integrating the current θ over the slice Σ with boundary ∂Σ = S: (3.15) 6 We recall that the convention for Hodge duality on a p-form α = 1 p! αa 1 ···ap dx a 1 ∧ · · · ∧ dx ap is where ǫ is the volume form.

JHEP09(2016)102
To arrive at this expression we have used the integral identity for a differential form ω. There are two sources of diffeomorphism non-invariance in this expression. The first is that we integrate over the surface Y (Σ) with boundary Y (S): this simply reflects that the diffeomorphism moves the coordinate location of the surface over which we integrate. The second source of diffeomorphism non-invariance is the presence of the second term in (3.16), which can be cancelled by addition of a boundary symplectic potential as we will now show. As in the case of Yang-Mills, we introduce a new variable into the phase space whose transformation law will cancel the boundary term in the symplectic potential. In gauge theory, this extra degree of freedom was a trivialization of the bundle; in a diffeomorphisminvariant theory, the natural analog is a coordinate system X which we view as a mapping X : U → M where U ⊂ R D is an open set. X is assumed to be invertible on its image. We assume for simplicity that Σ can be covered with one open set and we denote by σ ⊂ U the preimage of Σ, with boundary ∂σ = s. Under this map we have Σ = X(σ) and S = X(s). Under a diffeomorphism Y : M → M , the coordinate system X changes to Y * (X) = Y −1 • X. We can use this coordinate system to define the symplectic potential Θ Σ [g, δg] as the integral over the slice Σ = X(σ). Under a diffeomorphism its transformation is given by (3.18) Thus when integrating over the surface Σ the only change in the symplectic structure under diffeomorphism is via the boundary term.
In order to cancel this boundary term in the transformation law of Θ Σ , we introduce an additional boundary symplectic potential that depends on δX. The pullback of X by Y does not commute with the variation δ. Instead we find, defining δ X := δX • X −1 as in (3.6) and using the chain rule, that under diffeomorphism δ X transforms as Again this equation is directly analogous to the Yang-Mills case (2.16). Using δ X we can construct a boundary symplectic potential Θ S to cancel the action of the diffeomorphism on Θ Σ . Let us define where we recall that S = X(s). Under a diffeomorphism it will transform as

JHEP09(2016)102
Addition of this boundary term to the symplectic structure will therefore cancel the unwanted term in the transformation of the bulk symplectic potential (3.14). We therefore have the diffeomorphism-invariant symplectic potential given by The X dependence of this expression is implicit in Σ = X(σ), S = X(s), and in the definition of δ X . This symplectic potential satisfies by construction the invariance We can express the boundary symplectic potential in a more familiar tensorial form using the binormal to the codimension-2 surface S. Let n a be a timelike unit normal to S, and s a a spacelike unit normal. We can then form the unit binormal n ab = (n a s b − n b s a ), which is independent of the choice of normals n and s. The densitised unit binormal is ǫ ab = √ q n ab d D−2 σ where σ α are coordinates on S, and q is the determinant of the induced metric q αβ on S. The boundary symplectic potential Θ S can be expressed as

Gauge invariance and Hamiltonian generators
We now want to explore the gauge invariance of the extended symplectic structure. This follows closely the formal derivation given in section 2.2 for Yang-Mills. Let V be a vector field and consider the following variation This is the infinitesimal version of the transformation (g, . The infinitesimal version of diffeomorphism invariance of the symplectic potential shown in (3.23) translates into the identity where the equality is valid off-shell. This shows that the Hamiltonian is given by −I V (Θ) since We are therefore left with the evaluation of where we have used the identity (3.12). This shows that the generator of gauge transformation vanishes on shell even for gauge transformations with nontrivial support at the boundary.

Bulk and boundary symplectic forms
Having found Θ we can now evaluate the symplectic form Ω = δΘ, which determines the Poisson brackets between physical observables. There will be two terms in this expression, one coming from Θ Σ , and the other from Θ S . Here there is a new subtlety in the derivation that was not present in the Yang-Mills case. Let us first recall that σ is a surface in R D with boundary s and that we define X(σ) = Σ and X(s) = S as their images in spacetime.
The new subtlety is that the bulk symplectic potential is an integral over the surface X(σ) and this leads to nontrivial variation of the bulk symplectic potential under variations of X. This is in contrast with Yang-Mills where the bulk symplectic potential (2.13) does not depend on the choice of gauge. This dependence on X leads to δ X δg terms in the variation of Θ Σ . However we will see that these terms are all total derivatives, and hence can be expressed as boundary integrals. Once this is done the symplectic form consists of a bulk integral containing δgδg terms, and a boundary integral containing δgδ X and δ X δ X terms. This structure is similar to that of Yang-Mills theory. In order to vary an expression such as Θ Σ , which is an integral over the surface Σ = X(σ), we first use the pullback by X to express it as an integral over σ: (3.29) Then we need only vary the integrand. To vary the pullback of a tensor-valued variational form such as θ we make use of the identity (B.1): Recall that δ a X = δX a • X −1 is a vector field associated with the variation δX. Using (3.30) to vary the bulk symplectic potential Θ Σ , we find (3.31) The first integral is the usual bulk symplectic form for gravity [26]. It is obtained as the metric variation of the symplectic potential Θ Σ and therefore contains only δgδg terms: Explicitly, it is given by (see [43]): Here we use again the abbreviated notation used in (3.2). The expression can then be evaluated in terms of the metric tangent to the slices and the normal accelerations [28,44,45]. If we choose a slice Σ defined by the condition T = constant and equipped with a unit normal one-form n a dx a = −N dT , where N is the lapse we can, following [28], write the bulk component of the symplectic structure explicitly

JHEP09(2016)102
in terms of the induced metric h and extrinsic curvature K ab n := h ac h bd (∇ c n d ). If Σ has no boundary it is given by Ω Σ = Σ δΠ ab δh ab + δΠ T δT, (3.34) where the momenta conjugate to (h ab , T ) are given by where ∆ = D a h ab D b is the spacelike Laplacian. The second term in (3.34) says that the Laplacian of the lapse is canonically conjugate to the time variable. This seems unusual to promote the lapse to a kinematical variable. The reason it is usually not considered is because in the ADM approach, one fixes the foliation from the outset and focuses solely on hypersurface-normal deformations. These are only a subset of the deformationsthose that preserve the foliation -hence they correspond to deformations for which δT = 0. Once we do this and ignore the canonical pair (∆N, T ) the lapse becomes simply a Lagrange multiplier for the Hamiltonian constraint. Note however that the restriction to hypersurface-normal deformations is the reason why the ADM algebra is not equivalent to the algebra of spacetime diffeomorphisms, since in general a diffeomorphism changes the foliation. The second term in (3.31) is a mixed term that depends on δ X and δg. To analyse it we use Cartan's formula, The first term in (3.36) is proportional to dθ = ⋆ 1 2 (∇ a ∇ b δg ab − ∇ 2 δg) which vanishes for on-shell variations by the trace of the linearized vacuum Einstein equation. It is therefore a pure boundary term when δg is an on-shell variation. In general, when we do not restrict to on-shell variations, we note that the term i δ X dθ in the symplectic structure is an unfamiliar feature of the canonical formalism for gravity. Because δ X appears only in the interior product, this term depends on the normal component of δ X and hence arises from changes of the foliation. As such, it is similar to the δT term discussed in the preceding paragraph and disappears if one restricts to hypersurface-normal deformations only. Clearly, this issue of general deformations deserves further investigation, but for now, and because it doesn't play any role in our boundary study, we are content to consider variations for which this term vanishes.
We can now focus on the boundary terms in the symplectic form. The second term in (3.36) is manifestly a total derivative, and so can be written as an integral over S: (3.37) We also have to vary the boundary symplectic form, which can be accomplished using (3.30) again:

JHEP09(2016)102
We can simplify this expression by expanding the middle term using Cartan's formula and discard the total derivative, since the surface S has no boundary. The result is Combining this variation of Θ S with the boundary integral from the variation of Θ Σ we arrive at the boundary symplectic form: The first and third terms in this equation are the mixed δgδ X terms, while the second and fourth are the δ X δ X terms. The bulk piece contains only terms of the form δgδg if we restrict either to on-shell or to hypersurface-normal deformations. Thus we find a structure similar to that of Yang-Mills theory; the boundary terms in the symplectic structure lead to nontrivial Poisson brackets between X and certain components of the metric in the neighbourhood of S. Note that we can write this symplectic structure using the map π g [V ] = 1 2 ⋆ dg(V ) as A significant difference between Yang-Mills and gravity is that the gravity symplectic potential depends on derivatives of δX on the boundary: this shows that there will be a larger group of physical transformations that transform the metric in a first-order neighbourhood of the surface S. Moreover we expect the Poisson brackets between observables to contain derivatives of the generators, unlike in Yang-Mills.

Surface symmetries
As in gauge theory, we can distinguish between two types of transformations. The transformations considered in section 3.1 are diffeomorphisms, they simply relabel the points of M and hence are pure gauge. This is reflected in the fact that these transformations are null directions in the symplectic form. However, there is a different class of transformations we can consider that transforms the reference frame X a . These transformations are genuine physical symmetries of the phase space and have nontrivial generators that we will derive in this section. We will refer to them as surface symmetries, because they act at the entangling surface. Before describing the surface symmetries, we first recall the action of gauge transformations. A diffeomorphism acts on the metric by the pullback, and acts on X on the left: . This means that infinitesimally we have I V δX = −V • X, and so I V δ X = −V . Thus we have the infinitesimal transformation

JHEP09(2016)102
Recall that X is a map from an open set U ⊂ R D into M , and a diffeomorphism Y is a map M → M that relabels the points. Because of the transformation rule for X, the point X(x) still refers to the same physical point after a diffeomorphism. This transformation is just a relabelling of points in M , and this explains why it is pure gauge. Instead we will consider a transformation of the reference system Z : U → U . This corresponds to a change of the reference surface s which is mapped onto S, and of the surroundings of that surface. It acts by changing the labelling of the points, but keeps the dynamical fields unchanged: Note that like in Yang-Mills the symmetry acts on the opposite side of X from the gauge transformation.
To find the infinitesimal version of this transformation, we let Z = I +w+O(w 2 ), where w is a vector field on U . Then we have in components (X This defines a vector field ∆ w on phase space where Here we use upper case W to denote the pushforward of the vector field w to a vector field on M , i.e. W = X * (w), where X * : T U → T M is the pushforward. The infinitesimal transformation (3.45) is the analog of the transformation generated by the electric field in Yang-Mills (2.37). Note that because W depends implicitly on X, it has a nontrivial variation δW , which we will need below. This follows from the identity X * (W ) = w, and our formula for the variation of the pullback. Varying this identity we find (3.47)

Classification of surface symmetries
We now have a boundary symplectic structure for general relativity analogous to that of Yang-Mills theory. We have seen that as in Yang-Mills, the extension of the phase space via the introduction of boundary coordinates allows us to construct two types of canonical transformations: the gauge transformations whose Hamiltonian generators vanish on shell and the surface symmetries that act purely on the boundary variables. The latter are associated with a nonvanishing Hamiltonian generator that we will construct in the following section. The key difference between Yang-Mills and general relativity is that the Yang-Mills symmetry generators are spacetime scalars, while the symmetry generators in general relativity are vectors. This leads to a classification of different types of symmetries, depending on how the vector field W meets the surface S. To fix the idea it is convenient to introduce a local coordinate system (x i , σ µ ) where (x i ) i=0,1 denote coordinates normal

JHEP09(2016)102
to the surface and (σ µ ) µ=1,··· ,D−2 denote coordinates tangential to the surface which is located at x i = constant. Given an infinitesimal symmetry generated by a vector field W , we can divide its restriction to the surface S into its normal part W ⊥ = W i ⊥ ∂ x i and tangential part W = W µ ∂ σ µ as W = W + W ⊥ . The important point is that the contraction of the symplectic structure with the generator ∆ w depends not only on the value of W on S but also on the normal derivatives of W ⊥ , so that even when the normal component vanishes on S, it will have a nontrivial generator. Vector fields whose parallel component vanish W ⊜ 0 and whose normal component and its first derivative in the normal direction also vanishes W i ⊥ ⊜ ∂ j W i ⊥ ⊜ 0, where ⊜ denotes equality on the surface S, are in the kernel of the symplectic structure and therefore pure gauge. Thus we only implement as Hamiltonian transformations those vector fields that are nonvanishing in a first-order neighbourhood of S.
We can classify these surface symmetries into three classes: • Surface boosts These are generated by vector fields with W ⊜ 0 and W ⊥ ⊜ 0 but ∂ i W j ⊥ ⊜ 0. They generate position-dependent linear deformations of the normal plane of S. We call them surface boosts because they are a close analog of the boost transformations in the normal plane (x 0 , x 1 ). Like the usual boost, which leaves invariant a codimension-2 plane, the surface boosts leave S invariant. Unlike the usual boost, the surface boosts are not isometries of any particular background. Moreover we will see that these surface boosts generate a larger SL(2, R) subalgebra of the two-dimensional linear group at each point of S.
• Surface diffeomorphisms Vector fields with W ⊥ ⊜ 0 and W ⊜ 0 generate infinitesimal diffeomorphisms that map the surface S onto itself.
• Surface translations The transformations with W ⊥ ⊜ 0 transform the surface S normal to itself.
Here we will focus on the class of surface-preserving transformations, i.e. the transformations such that W ⊥ ⊜ 0. These consist of the surface boosts and surface diffeomorphisms, which do not move the surface S. This is appropriate to the study of an entangling surface which is kept at a fixed location in spacetime, while we allow the angle at which the foliation Σ meets S to change. Previous works have focused on the surface translations (which have been called "would-be pure gauge" degrees of freedom [46]), but have largely overlooked the surface-preserving generators. We leave a full treatment (treating all surface symmetries together) as an interesting area for future work, and comment on some of the technical challenges in the discussion.

Surface symmetry generators
We now consider the generators H w of the surface symmetries ∆ w . These are obtained by solving the equation δH w = I ∆w Ω. We first find I ∆w Ω, which has both a bulk component and a boundary component. The bulk component I ∆w Ω Σ vanishes, because the

JHEP09(2016)102
transformation ∆ w vanishes when acting on the metric. We are left with the boundary term, which is: Here we have made use of the identity I ∆w δ(δ X ) = − 1 2 I ∆w [δ X , δ X ] = [δ X , W ]. Here (δπ g ) denotes the variation of the metric terms that appear in π g without variation of the argument.
We see that this interior product contains two types of terms, one associated with surface-preserving symmetries and one associated with translations. The first one is the translation contribution (3.49) This term vanishes for surface-preserving symmetries. Indeed for surface-preserving symmetries, W is tangent to S, and hence the interior product i W α vanishes when pulled back to S, for any form α.
From now on we specialize our analysis to the case of surface-preserving symmetries. In this case the previous term vanishes and we are left only with the surface-preserving contribution where we have used the definition of π g and discarded a total derivative term in order to write the second term as a Lie derivative. Using the identity (3.47) for δW , and the variational identity (3.30), we clearly see that this contraction is Hamiltonian that is δH w = I ∆w Ω, where H w is therefore the Hamiltonian generator of the surface-preserving transformations.
We can now easily find the algebra of the surface-preserving transformations. If v is another surface-preserving transformation, then by direct substitution we have (3.52) Note that the surface-preserving condition is preserved by the Lie bracket: if v and w are surface-preserving then so is [v, w]. We also observe that the surface-preserving symmetries do not acquire a central charge at the classical level. We will give a more explicit characterization of this algebra and its Poisson brackets in the following subsection. We also note that our Hamiltonian (3.51) coincides with the Noether charge of refs. [47,48]. However the context here is different: W is not assumed to be an isometry of any fixed background, and S need not be a section of a Killing horizon. Instead W = X * (w) parametrises an infinitesimal change of coordinate frame of the boundary surface.

Algebra of surface symmetries
In order to describe explicitly the algebra of surface symmetries and its generators, we make a "2+2" decomposition 7 of the metric in a neighbourhood of S, and evaluate the generator H w of a surface-preserving symmetry w. In terms of the timelike unit normal n and spacelike unit normal s to S, the generator takes the form To evaluate the Hamiltonian explicitly in terms of components of the metric we introduce coordinates (x i , σ µ ) adapted to the 2 + 2 decomposition. Let σ µ be coordinates on S, and let (x i ) i=0,1 be two coordinates in the normal directions. In these coordinates we can parametrize the metric as Here q µν is the induced metric on S, h ij is a generalized lapse which defines the normal geometry and A µ i is a generalized shift. We can view A µ i as a normal connection, since under relabelling δσ µ = ϕ µ it transforms as a connection δA µ i = ∂ i ϕ µ + [A i , ϕ] µ , where the bracket is the Lie bracket for vectors tangent to S. The relationship between the normals n A and the normal coordinates x i is given by a 2 × 2 matrix with coefficients n a i : where η AB = diag(−1, +1) is a flat 2D normal metric. This expresses the fact that n A is normal to the the level sets of constant x i . It is also necessary to express the normal vector fields n A = n a A ∂ a in terms of the metric coefficients: Indeed we can readily check that i n A (n B ) = δ B A and that i n A (dσ µ − A µ i dx i ) = 0.

JHEP09(2016)102
We can now further evaluate the Hamiltonian (3.53) in terms of the 2+2 decomposition of the metric where we have used that W ⊥ ⊜ 0. Contracting with ǫ AB we find It will be convenient to introduce the densitised normal metric It is traceless by construction. Moreover, since det(h) = −(det n) 2 , and det(ǫ) = 1, we find that Using the form (3.56), we can also evaluate the twist density vector It is proportional to the curvature of the normal connection. The proportionality coefficient √ q/det(n) is the same we used to rescale the normal metric. From (3.59), (3.62) we see that the Hamiltonian can be written in terms of the densitised normal metric and curvature as We can now see that the normal component of the Lie bracket for surface-preserving transformations becomes the sl(2, R) algebra while the tangential component reduces to the (D − 2)-dimensional diffeomorphism group acting on S: from which we find the Poisson brackets between generators: We see that the components of the densitised curvature F act as generators of tangential diffeomorphisms. They have a nontrivial action on H i j , which transforms as a scalar under JHEP09(2016)102 diffeomorphism, while H i j itself generates a local sl(2, R) algebra. Thus we see that the group of surface-preserving symmetries is Diff(S) ⋉ SL(2, R) S . Let us introduce sl(2, R) generators: K 0 is an elliptic generator while K 1 and K 2 are hyperbolic. They satisfy the sl(2, R) commutation relations The determinant of H is which is precisely the Casimir of sl(2, R). The unitary representations of SL(2, R) come in a discrete series for which the Casimir is positive and continuous series for which it is negative. Since we require det H to be negative, we are interested in the continuous series. Note that the previous relation means that the Casimir of SL(2, R) is the area element Since the local SL(2, R) symmetry is a new feature, we briefly comment on the relation between this symmetry and the more familiar boost symmetry of Carlip and Teitelboim [49]. In that work it was argued that the near horizon symmetry consists of a single boost, and that the quantity canonically conjugate to the boost is the horizon area density √ q. The difference between the two symmetry groups comes because our SL(2, R) transformations are not isometries of the normal geometry. The relation between the two symmetries can be understood as a particular gauge-fixing. The normal metric h ij transforms in the adjoint representation of SL(2, R), and since the metric must be nondegenerate it will spontaneously break this symmetry. It is natural to fix the metric to the flat metric h ij = η ij , in which case one should consider the little group of transformations that preserves h ij . This gauge corresponds to the conditions H 0 0 = 0 and H 0 1 = H 1 0 = −1. The little group has a single unbroken generator K 2 which generates boosts. In this gauge, the generator becomes the area element of the entangling surface. Thus we recover Carlip and Teitelboim's result as a gauge fixing of ours. The situation can be summarized with a simple analogy to the rotation group SO(3). The Casimir J 2 is of course different from the rotation generator J z . If we fix the axis of rotation then classically J z and √ J 2 coincide, but from the point of view of the representation theory they are of course distinct and we should not conflate the two. In particular, they have different spectra at the quantum level and the representation theory is controlled by the Casimir J 2 .

JHEP09(2016)102 4 Quantization and entanglement
So far we have constructed the algebra of boundary symmetries in Yang-Mills theory and in general relativity. At the quantum level, we expect this algebra to be realized in the usual way as an algebra of operators on a Hilbert space. In Yang-Mills theory, this is simply the algebra of surface charges on the entangling surface. In a lattice realization it is the 'electric center' algebra of ref. [15], and it is represented on a Hilbert space spanned by generalized spin network states used in ref. [18].
In gravity, we do not have an explicit construction of the Hilbert space as we do in Yang-Mills. In constructing such a Hilbert space it is extremely useful to have a group of symmetries as a guide. Here we have constructed a large group of surface symmetries whose representation theory can guide the construction of the Hilbert space of a local region. We therefore give a brief sketch of how to construct representations of the algebra of surface symmetries. This appears to be a rather challenging problem whose full solution we are forced to leave to future work.

Quantization
Having found the algebra satisfied by the surface symmetries in gravity we can look briefly at its quantization, which amounts to studying representations of this algebra. Our surface preserving symmetry group Diff(S) ⋉ SL(2, R) S has the structure of a semidirect product, whose projective representations have been extensively studied by Mackey [50]. This approach is a generalization of the familiar classification of projective representations of the Poincaré group SO(1, 3) ⋉ R 4 . First one starts with irreducible representations of the translation group R 4 , which are labelled by momenta, and then defines the little group as the subgroup of the Lorentz group that preserves the given momentum vector.
In our case we have diffeomorphisms playing the role of the Lorentz group and SL(2, R) S playing the role of translations. We first fix an irreducible representation of SL(2, R) S , which amounts to choosing the value of the Casimir operator at each point of S. We will assume that this choice is made smoothly which means that by the relation (3.72) we specify an area form on S. This area form is analogous to the choice of momentum, and the total area (which is diffeomorphism-invariant) is analogous to the mass. Next we consider irreducible representations of the little group, which is the subgroup of diffeomorphisms that preserve the given area form, which we denote APD(S), where APD stands for Area Preserving Diffeomorphism. The choice of APD representation is analogous to the spin in the classification of Poincaré group representations. Thus the irreducible representations of our algebra are labelled by a total area, and an irreducible representation of APD(S). Using this information, one can then construct the irreducible representations of the surface symmetry algebra by the method of induced representations.
The algebra of area-preserving diffeomorphisms APD(S) of a surface S has been studied particularly in the case where S = S 2 (which fortunately is also the most physically relevant case, as it would apply to both black hole and cosmological horizons in a four-dimensional spacetime). In particular its structure constants are known [51], and it is known that for JHEP09(2016)102 a surface of genus g there are 2g central charges. It would therefore be of great interest to understand more completely the representations of APD(S 2 ).

Entangling product
Provided we can construct Hilbert spaces associated to local regions, the next question is how to glue together two such Hilbert spaces along a common boundary to make the Hilbert space of a larger region. The diffeomorphism-invariance of general relativity would appear to be an obstacle, since diffeomorphisms naively can move the location at which the two Hilbert spaces are glued. Here we point out that such a gluing construction has effectively been carried out for a large class of rigorously-constructed diffeomorphisminvariant quantum field theories: topological quantum field theories (TQFTs). Although these theories have no local degrees of freedom, they provide useful models of how we can glue subsystems together in a diffeomorphism-invariant theory.
In the axiomatic approach to TQFT, every closed codimension-1 surface Σ is assigned a Hilbert space Z(Σ), and each codimension-0 surface M with boundary ∂M is associated with a state in the Hilbert space Z(∂M ) [52]. These satisfy a set of identities that allow one to calculate the partition function of a closed manifold M by cutting along codimension-1 surfaces and gluing them back together. This structure is sufficient to classify two-dimensional TQFTs [53], but in higher dimensions this is not sufficient. This is essentially because codimension-1 surfaces are more complicated objects as the spacetime dimension increases.
One of the axioms of TQFT states that a disjoint union of two surfaces is assigned a tensor product of Hilbert spaces, Z(Σ 1 ⊔ Σ 2 ) = Z(Σ 1 ) ⊗ Z(Σ 2 ). One can augment this structure by allowing codimension-1 surfaces to be cut along codimension-2 surfaces, the resulting structure is called an extended TQFT [54]. This additional structure is used to classify higher-dimensional extended quantum field theories in terms of n-categories [55].
In general, the structure associated to a codimension-2 surface is a category, and a codimension-1 surface with boundary is a morphism in its boundary category. However one can understand this better in a simple example, such as two-dimensional Yang-Mills theory [56]. Two-dimensional Yang-Mills is not quite a TQFT according to the usual axioms. First, its Hilbert space is infinite-dimensional, and second it is not invariant under the full diffeomorphism group but only under the subgroup of area-preserving diffeomorphisms. However it is close enough to a topological quantum field theory that much of the axiomatic framework carries over to this case [57].
In 2D Yang-Mills, the Hilbert space associated to a closed codimension-1 manifold (a circle) is the space of square-integrable class functions of the gauge group G. These are functions ψ(u) of the holonomy u = exp(i A) ∈ G that are invariant under the gauge transformation u → g −1 ug, i.e. ψ(u) = ψ(g −1 ug). One can also construct the Hilbert space associated to an interval, which is the space of square-integrable functions of the connection integrated along the interval u = exp(i A) ∈ G. In the case of the interval, the states are not restricted to be class functions. The Hilbert space of an interval then carries extra structure: it carries two representations of G corresponding to left and right multiplication JHEP09(2016)102 of the group element u acting as u → gu and u → ug −1 respectively. This endows the Hilbert space of an interval with the structure of a module over the group algebra of G.
The Hilbert spaces of intervals are not glued together with the usual tensor product of Hilbert spaces. Recall that one defines the tensor product over C as the cartesian product H 1 × H 2 modulo the equivalence (αv, w) ≃ (v, αw) where α ∈ C. In the case of a tensor product over an algebra A, we define H 1 ⊗ A H 2 such that (va)⊗w = v⊗(aw) for all elements a ∈ A. When the algebra is the group algebra A = C(G) we denote this entangling product H 1 ⊗ G H 2 . As we have seen this entangling group and the corresponding product naturally appear in Yang-Mills and gravity as the surface symmetry group. Interestingly it has been conjectured in [58] that the slice of a spacetime outside a black hole horizon should have precisely this structure. We see here that it appears in fact generically for an entangling surface in a gauge theory and in gravity.

Discussion
We have defined a formalism for associating a classical phase space to a bounded region that applies to both Yang-Mills theory and general relativity. We have restricted our analysis to entangling surfaces: codimension-2 surfaces which are fixed in spacetime. By including the choice of gauge as a canonical coordinate, we obtain an extended phase space that is invariant under gauge transformations and which carries a representation of an induced boundary symmetry group G S . In Yang-Mills theory, the boundary degrees of freedom consist of a section of a principal G-bundle and the normal component of the electric field, which act like surface charges on the boundary. In gravity, the boundary degrees of freedom consist of the location of an embedding surface s → X a (s) and its normal and tangential geometry, both intrinsic and extrinsic. The remarkable feature of such a construction is that it selects a subset of the boundary observables to be the charges for this emergent boundary symmetry. This is the normal electric field in Yang-Mills, and the normal conformal metric, normal curvature and area density in gravity. While these components of the fields are commuting variables from the point of view of the bulk theory they become noncommutative in the presence of a boundary, forming a Lie subalgebra. We have seen that the appearance of this boundary symmetry group allow us to define at the quantum level the total Hilbert space associated with the gluing of two region as the entangling product. The entangling product is a generalisation of the tensor product where states are restricted to be singlets under the boundary symmetry. We have described a classical version of this entangling product, which takes the form of a symplectic reduction on the product of two phase spaces.
It is well-known in the case of lattice Yang-Mills that the identification of boundary charges as the normal electric field can be understood as a consequence of the discretisation. In a discretisation two facts are clear. First, the discrete electric field generates the action of a nonabelian group and hence must be noncommutative. Second, one must demand continuity of the electric field, but not of the magnetic field across the discrete surface. This is because the electric and magnetic fields are non-commuting variables and the identification E L = E R generates a canonical transformation acting on the magnetic fields JHEP09(2016)102 B L and B R under which admissible magnetic field observables must be invariant. This means that the discretisation requires a splitting of the boundary data into two halves. One half becomes the nonabelian charges and the other half corresponds to the bulk conjugate variables. We can now appreciate why putting gravity on a lattice in a way that respects gauge symmetry is so challenging: it requires an understanding of how to split in half the metric variables along the discrete lattice, and it also requires us to understand what commutation relations these discrete variables should satisfy. What is remarkable in our analysis is that we recover this splitting in two halves of the boundary variables and their commutation relations from an analysis of gauge invariance. This suggest that our analysis should provide new clues into the problem of discretising gravity in a covariant manner. If one views the process of discretisation as a necessary step in order to define a quantum field theory nonperturbatively, this could be relevant for the construction of a theory of quantum gravity.
In fact this reflexion leads to an interesting conjecture: in the continuum we have seen that preserving gauge symmetry in the presence of boundaries implies the existence of a boundary symmetry on entangling surfaces. One can view the vertices of a discrete spacetime as a set of regions bounded by entangling surfaces, which are glued together along links of the lattice. If we take this point of view we can ask what are the necessary conditions on the discrete variables to ensure that the resulting theory manifests diffeomorphism symmetry in the continuum limit? A natural conjecture in view of our analysis is that each vertex of the discretisation needs to carry a representation of the entangling group G S . In other words, not only does the boundary symmetry appear as a result of gauge invariance, but gauge invariance could be seen to emerge as the result of preserving the local surface symmetry. This is a new and very tantalising possibility that deserves a proper analysis.
Our analysis so far has been restricted to entangling surfaces at a fixed position in the bulk of the spacetime. We have therefore restricted our attention to surface symmetries that preserve the location of the entangling surface. But the local symmetry group also includes surface translations that move the entangling surface. These transformations introduce new complications; in particular we have not yet found a Hamiltonian generator of surface translations analogous to that for surface-preserving transformations (3.51). We have laid the groundwork for such an analysis, but it needs to be completed. Previous work suggests that one will have to impose boundary conditions that fix the location of the surface in order to have a nondegenerate symplectic structure, as found in [33]. It remains to be seen whether such an invariant specification of the boundary region is possible in general and whether it is necessary and sufficient to complete the symplectic structure.
It would be very interesting to generalise our results to the case of more general surfaces and to understand how our results connects to the results obtained when the boundary is at infinity. Much recent work has focused on asymptotic symmetries of Yang-Mills theories [59,60] and gravity [61,62]. In particular, [63] showed how the asymptotic symmetries of electrodynamics may arise by augmenting the phase space with edge modes living at the "edge of infinity". However, so far there is no construction of an asymptotically flat phase space for gravity that contains a representation of the BMS group, except in 2+1 dimen-

JHEP09(2016)102
sions [64]. The methods developed in this paper may be useful to construct a phase space that includes the soft photon and graviton modes on which the action of the asymptotic symmetries is Hamiltonian. We have seen so far that our method constructs a phase space containing a representation of a large nontrivial subgroup of the diffeomorphism group. It would therefore be interesting to impose asymptotically flat boundary conditions within our formalism, and to see what asymptotic symmetry group is obtained.
Another related question concerns the relationship of our analysis with the usual Israel junction conditions. Our analysis shows that we need to identify the normal geometry across the codimension-2 entangling surface but not the tangential geometry (except for the area density), while the Israel junction condition (in the absence of a surface stress tensor) demands the continuity of the tangential metric and extrinsic curvature across a codimension-1 boundary. While there is no immediate contradiction between these results since they apply to radically different boundary conditions (codimension-2 versus codimension-1 boundaries) they create an interesting tension that has to resolve itself when we consider dynamical boundaries. In order to understand this we first need to complete the canonical analysis of surface translations.
In particular, this matching raises the following issue. In electrodynamics, the division of the Hilbert space requires splitting electric field lines into positive and negative surface charges, which must cancel in pairs to create gauge-invariant states. Does this mean that the corresponding matching condition in general relativity requires the introduction of both positive and negative energy? 8 Here we can provide a resolution of this tension. Consider for example a black hole spacetime with a bifurcation surface divided into two regions along this surface. Then the energy is related to the generator of the Killing flow, which we have found is simply a component of the normal metric. These metric components are matched up to a sign coming from the convention of choosing n to be an outward-facing normal vector. Thus there is no contradiction because the generators of surface boosts do not have any fixed sign. On the other hand, the area element does need to be positive in order for the surface S to be spacelike. However we have identified this as the Casimir of SL(2, R), and so the matching condition imposes that it must be equal on both sides. Thus we see that our matching condition gives the expected positivity properties, and there does not seem to be any problem of negative energy. We expect to find no further contradiction when considering the generators of surface translation, but this remains to be explored.
Here we have focused on Yang-Mills theory and general relativity, but our formalism appears to be flexible enough to handle a very general class of field theories. It would be interesting to apply our formalism to more general theories of gravity, for example with higher curvature couplings. As a simplest example, inclusion of a nonzero cosmological constant would allow us to make contact with gravity in asymptotically AdS space. The algebra we derive is reminiscent of the commutation relations of a surface stress tensor, which is an important part of the AdS/CFT correspondence, though we have so far not imposed any boundary conditions. Another interesting case for our formalism is Chern-Simons theory, where the edge states and their relation to entanglement entropy have been understood JHEP09(2016)102 for some time. There we anticipate nontrivial Poisson brackets between ϕ, which do not appear in pure Yang-Mills theory. This may help to explain the appearance of quantum group symmetry, and may shed new light on the open problem of discretizing Chern-Simons theory, which has so far only been solved for particular lattices and gauge groups [65,66].
More generally, our work seems to be pointing toward a certain extension of the formalism of quantum field theory. In topological quantum field theory (TQFT) there is a useful notion of an extended TQFT [54]. In extended TQFT, one associates Hilbert spaces not only to closed codimension-1 surfaces, but also to codimension-1 surfaces with boundary. The codimension-2 boundaries of these surfaces are then given some additional structure that encodes how the Hilbert spaces of two regions with boundary may be glued together to make a larger region. This additional structure seems to be precisely what is needed in order to define entanglement entropy, which suggests the utility of considering extended field theory beyond the TQFT context. In addition to formalizing the notion of extended QFT, it also raises the interesting question of whether a quantum field theory can always be extended in this way and whether such an extension is unique. Theories with electric-magnetic duality may provide an interesting example, as they can be viewed as gauge theories in two different ways, suggesting multiple extensions and hence multiple ways of defining the entropy for such theories. The existence of different ways to define entanglement entropy in the algebraic definition of entropy were noted in [15], and there may be an analogous freedom in the extended Hilbert space definition. Another interesting case are certain quantum Hall states admitting distinct edge phases [67], further suggesting different ways of extending the same TQFT.
Ultimately, we would like to apply these classical considerations to the problem of entanglement entropy in quantum gravity. In order to do this, we need a Hilbert space which is the quantization of our classical phase space. This is a hard problem, but it may be easier than the usual quantization problem because we have a large physical symmetry group. Thus an important step in the quantization is to find the irreducible representations of Diff(S) ⋉ SL(2, R) S . These irreducible representations play an important role in the entanglement entropy, as we will review below. In particular, the irreducible representations label superselection sectors of the Hilbert space H Σ . Note that the entanglement entropy for compact groups contains a term that depends on the dimension of the irreducible representations of the gauge group. Since SL(2, R) is noncompact, all irreducible representations are infinite-dimensional; presumably this leads to a divergence that must be regulated. However the regularization of this divergence via heat kernel methods has been understood in other models with SL(2, R) symmetry [68], so we do not expect this to pose any great difficulty.
Another important question that affects the quantization is the existence of central charges. These central charges appear at the classical level in AdS 3 gravity [69], and underlie the derivation of BTZ black hole entropy [70] and its generalization to higher dimensions [32]. The derivation of (3.52) shows that there is no classical central term in the algebra of surface-preserving transformations. However the central charge in [32] comes from diffeomorphisms that move the surface, so we may expect central terms to appear once we include the surface translations.

JHEP09(2016)102
Finally, our work can be seen as a first step towards the problem of defining entanglement entropy in nonperturbative quantum gravity. Once we have identified the Hilbert spaces H Σ carrying an action of the group G S of surface symmetries, we can define the entanglement entropy by the embedding We can now define a density matrix ρ Σ by embedding states of H into H Σ ⊗ H Σ , and tracing over H Σ . Any density matrix ρ Σ defined in this way commutes with the action of G S . This leads to a large symmetry of the problem, and a significant simplification of the entanglement entropy. To make use of this symmetry, we first decompose the Hilbert space into irreducible representations of G S . Letting R denote the irreducible representations, we have H Σ = R H Σ,R . Since the reduced density matrix commutes with the action of the surface symmetries, we can also decompose it as ρ Σ = R p(R)ρ Σ,R . Here p(R) is the probability that the state lies in the superselection sector described by the representation R; this factor is needed to keep the normalization condition Tr(ρ Σ ) = 1 and Tr(ρ Σ,R ) = 1. The entanglement entropy in gauge theory then takes the form [18]: S ent = R p(R) (− log p(R) + log dim R + S(ρ Σ,R )) .

(5.2)
Thus the representation theory of the surface symmetry group G S and the identification of its generators already tells us a lot about the entanglement entropy. It is natural to compare formula (5.2) for the entanglement entropy with the generalized entropy In both of these equations the entropy is a sum of an entropy associated with degrees of freedom localized away from the entangling surface (the term S out of the generalized entropy and the term S(ρ Σ,R ) of the entanglement entropy) and the expectation value of an operator (the Bekenstein-Hawking entropy term in the generalized entropy, and the term log dim R in the entanglement entropy). The quantity log dim R is a gauge-invariant observable that depends on the fields at the entangling surface -for example, in Yang-Mills we can determine it explicitly as a functional of the normal electric field. This suggests the possibility that the Bekenstein-Hawking entropy arises as the part of the entanglement entropy that is determined by the representation theory of the diffeomorphism group. In other words, the Bekenstein-Hawking entropy of a black hole could arise from counting the different ways in which the interior can be glued to its exterior; our boundary symplectic form gives the appropriate measure for this counting. There is of course much to be done to establish or refute this conjecture. In particular, we have to study the representation theory of the surface symmetry group to establish the relationship between dim R and the geometry of the surface Σ. While the appearance of the area density as the Casimir of local SL(2, R) symmetry is a suggestive, the real test of quantum gravity is to obtain the area law with the correct finite prefactor of 1/4G.

JHEP09(2016)102
For the Lie derivative along δ X , we must define the Lie derivative with the opposite sign, This leads to the graded Liebniz rule L δ X (αβ) = (L δ X α)β + (−1) p α(L δ X β) (A.8) which reflects the fact that L δ X increases the form degree by one. We now show that the field-space Lie derivative is related to the spacetime Lie derivative as: We first show that it holds for scalars and the coordinate one-forms δφ, and then extend this to arbitrary products. For a field-space 0-form f we have and a coordinate one-form, we have Now suppose that forms α and β both obey the identity (A.9). Then their product does also, which follows from the fact that both L V and L V + I δV satisfy the Liebniz rule: This completes the proof of (A.9). For a transformation such as L δ X , the identity (A.9) is slightly modified. In this case L δ X is an infinitesimal transformation of the fields that depends linearly on δ X , so it is a vector-valued differential form. The Lie derivative can be naturally extended to vectorvalued forms (see [71]) and we can carry out the same generalization to phase space. Since δ X is a field-space 1-form, L δX increases the degree by one and hence unlike L V obeys an anti -Liebniz rule. This leads to the identity: (A. 13) To see this note that on a scalar f we have (as before) However since L δ X and δ are both antiderivations, we have L δ X (δφ) = I δ X (δδφ) − δ(I δ X δφ) = −δ(L δ X φ) = L δ X δφ − L δ(δ X ) φ = (L δ X − I δ(δ X ) )δφ. (A.15) Thus (A.13) holds for field-space 0-forms and coordinate 1-forms. We can now extend this to all forms, using the fact that L δ X , L δ X and I δδ X are all antiderivations, and hence satisfy a graded Leibniz rule. Thus the relation (A.13) is proved for all (p, q) forms.

JHEP09(2016)102
B Variation of the pullback Let X : M → M be a diffeomorphism of spacetime, T be a spacetime tensor and fieldspace p-form, and let X * T denote the pullback under diffeomorphism. The purpose of this appendix is to prove the following variational formula: δX * (T ) = X * (δT + L δ X T ) (B.1) To prove this formula we first establish it on spacetime tensors, and then on a basis of variational 1-forms. We can then extend to arbitrary tensor-valued variational forms by the product rule, using the fact that both δ and L δ X act as derivations on spacetime tensors, and as antiderivations on variational forms. To prove (B.1) for spacetime tensors, it suffices to prove it for scalars, vectors and covectors and then to extend to tensors of arbitrary rank by taking products. For a scalar function f : where we have used that δ A X = δX A •X −1 . Here we temporarily introduce upper-case letter indices on X A to distinguish them from the lower-case indices on x a . For a vector v a : For a covector w a : δX * w a = δ(∂ a X A w A (X)) = ∂ a X A (δw A )(X) + ∂ a δX A w A (X) + ∂ a X A δX B ∂ B w A (X) = X * (δw a + ∂ a δ b X w b + δ b X ∂ b w a ) = X * (δw a + L δ X w a ). (

B.4)
This establishes the identity for all spacetime tensors, using the fact that both δ and L δ X are derivations.
To extend the identity to variational forms, it suffices to prove it on a basis of variational 1-forms. Our basis of variational 1-forms will consist of δg ab and δ a X . For the metric variation δg ab : δX * δg ab = δ(∂ a X A ∂ b X B δg AB (X)) ∂ a δX A ∂ b X B δg AB (X) + ∂ a X A ∂ b δX B δg AB (X) + ∂ a X A ∂ b X B δX C ∂ C δg AB (X) = X * ((∂ a δ c X )δg cb + (∂ b δ c X )δg ac + δ c X ∂ c δg ab ) = X * (L δ X δg ab ). (B.5)