Quantum Reference Frames at the Boundary of Spacetime

An analysis is given of the local phase space of gravity coupled to matter to second order in perturbation theory. Working in local regions with boundaries at finite distance, we identify matter, Coulomb, and additional boundary modes. The boundary modes take the role of reference frames for both diffeomorphisms and internal Lorentz rotations. Passing to the quantum level, we identify the constraints that link the bulk and boundary modes. The constraints take the form of a multi-fingered Schr\"odinger equation, which determines the relational evolution of the quantum states in the bulk with respect to the quantum reference fields at the boundary.


Introduction
To the best of our knowledge, all matter obeys quantum mechanical laws.The equivalence principle speaks of another such universality of nature.Gravity interacts with all matter in the same universal way.In the linearized regime with small perturbations around Minkowski space, the equivalence principle is a consequence of local gauge invariance and global momentum conservation [1].Besides the particles of the standard model, there is now the graviton propagating over a fixed background.From this perspective, gravity behaves just like another form of matter.Thus, it should have a quantum description as well.The conceptually simplest approach towards such a quantum theory of gravity is to keep the background fixed and treat the perturbations as a quantum field on Minkowski space.There are two major problems with this viewpoint.The first problem is technical-the resulting theory is perturbatively nonrenormalizable, that is, it fails at high frequencies with respect to the fixed classical background.The second problem is conceptual.The perturbative quantization leaves the background fixed, treating it as a classical arena on which the quantum fields evolve.From the perspective of the full nonlinear Einstein equations, this seems highly problematic.The split g ab = η ab + 2f ab of the curved spacetime metric g ab into a classical background η ab , which has no quantum properties at all, and a perturbation f ab is completely arbitrary.If quantum theory is truly universal, it will apply to the gravitational field as a whole.The background would have quantum properties as well.While this possibility has been largely ignored in the perturbative approaches, it is the basic raison d'etre of nonperturbative approaches to quantum gravity.
In this article, we present a novel solution for this second and more conceptual problem within the perturbative regime.The key novelty is that we take into account the quantum properties of the background metric and encode them into a particular form of quantum reference frame at the boundary of spacetime.In this way, our approach sits halfway between perturbative gravity and nonperturbative approaches.Moreover, it connects two converging areas of research: the area of boundary modes, which studies the state space of gauge theories in bounded regions of spacetime, e.g.[2][3][4][5][6][7][8][9][10][11][12][13][14], and the field of quantum reference frames, which has regained much interest in recent years in the quantum foundations community, e.g.[15][16][17][18][19][20][21][22].
The basic idea is motivated by the following observation.For any flat signature (−+++) metric η ab on a four-dimensional spacetime manifold M ≃ R 4 there always exist coordinate fields 1{X µ = (cT, X, Y, Z)} that bring the metric into the standard diagonal Minkowskian form Since the metric is flat, we can refer to the four scalar fields X µ as inertial frames of reference.More concretely, Eq. ( 1) can be read in two ways.If we know η ab , we will also know X µ modulo global symmetries (translations, boosts and rotations).If, on the other hand, we take the fields X µ as given, Eq. ( 1) defines a flat reference metric η ab with respect to which the reference frame {X µ } is inertial.Thus, the metric η ab and the inertial reference fields {X µ } contain the same information.
Combining the universality of quantum theory with the observation that all reference frames are, ultimately, instantiated by material objects, it is then natural to consider the possibility that the reference frames themselves have a quantum description.This is the idea underlying the research on quantum reference frames.By the above argument, the resulting quantum coordinates define a flat Minkowskian quantum geometry η ab = η µν ∂ a X µ ∂ b X ν .Finding a quantum description of the background geometry and understanding quantum inertial frames are thus two sides of the same coin.
From a general relativistic point of view, it seems odd at first to place the focus on coordinates.After all, the observables of the theory are coordinate invariant.General nonlinear coordinate transformations are mere gauge symmetries, mathematical redundancies that do not affect physical predictions.However, this is no longer true when we take into account the boundary conditions [23].When we consider finite subsystems, or spacetimes with asymptotic or inner boundaries such as a black hole horizon, we need to impose boundary conditions to obtain definite predictions.Only those diffeomorphism that preserve the then-fixed boundary conditions remain gauge transformations.The residual diffeomorphisms that change the boundary conditions turn into actual symmetries, transforming one physical state into an inequivalent configuration.Gauge symmetries are thus broken at the boundary that separates the system from its environment.At the boundary, the coordinate fields X µ regain physical meaning and the value of physical observables will depend on them.As a result, we obtain additional boundary degrees of freedom.Since these boundary modes are fields intrinsic to a codimension-two surface bounding a partial Cauchy surface, they are often referred to as edge modes.Classically, these edge modes can be interpreted as reference frames for the symmetry groups of the theory [12][13][14]-an interpretation that extends to the quantum level, where, as we will show, the edge modes manifest themselves as a type of quantum reference frame.The connection between the background metric, edge modes, and quantum reference frames, which is established in our work, is depicted in Fig. 1.

Minkowskian quantum geometry quantum reference frames local quantum geometry boundary modes & asymptotic symmetries
Figure 1: Basic themes of this paper.We consider linearized quantum gravity in finite regions (local quantum geometry).At the boundary, we obtain additional edge modes that carry a representation of the boundary symmetry group (left corner).Going to an asymptotic region, we recover asymptotic Bondi-Metzner-Sachs symmetries.In the absence of radiation, the boundary modes take the form of inertial coordinate frames, which give rise to a notion of Minkowskian quantum background geometry (upper corner).Finally, the boundary modes serve as quantum reference frames for the relational evolution in the bulk (right corner).
To access all possible gravitational boundary modes, we consider a compact region of space on an initial hypersurface.While there is a large and growing literature characterizing the phase space of such regions in general relativity and the emergence of additional reference degrees of freedom at the boundary [2][3][4][5][6][7][8][9][10][11][12][13][14]24] most of the analysis is performed within the full nonperturbative framework.For the present work, on the other hand, we utilize perturbative methods.This viewpoint has two main advantages.First of all, we obtain a clear physical interpretation of the bulk and boundary modes.Second, we can take the analysis to the quantum level without going into the mathematical difficulties of the nonperturbative approaches.Starting at the classical level, we consider the perturbative expansion of the presymplectic two-form with respect to the gravitational coupling constant, which is √ 8πG.To understand the effects of matter in a physically transparent model, we consider N point particles contained within the region.Upon expanding the metric to second order in the perturbation, the presymplectic two-form splits into contributions from particles, gravitational waves, Coulomb fields sourced by both the particles and the gravitational waves themselves, as well as additional reference fields at the boundary.The reference frames take the form of coordinate fields and local SO(1, 3) transformations at the boundary, which provide a reference frame for internal boosts and rotations of the tetrads.The conjugate variables find their interpretation as boundary currents for energy, momentum and spin.At the quantum level, our analysis has three important implications.First of all, the partition of the phase space into matter, radiative, and boundary modes, together with the resulting Poisson commutation relations, informs us about the basic constituents of the algebra of observables at the quantum level.We will argue, in particular, that in any quantum theory of gravity the boundary modes must admit a quantum description as well, which leads immediately to the notion of quantum reference frames.
Second, as we will see below, the field equations place constraints on the relation between the bulk and boundary modes.The then-necessary constraints are a remnant of the Wheeler-DeWitt equation of the full nonperturbative theory and describe the relational evolution of the wavefunction with respect to the reference fields at the boundary.Finally, we take the limit to an asymptotic boundary and show how unitary representations between different asymptotic reference frames will carry a quantum representation of the group of asymptotic Bondi-Metzner-Sachs (BMS) transformations [25][26][27][28] at spatial infinity.Our results go beyond the existing framework for quantum reference frames, because they provide an extension to pointwise transformations at the boundary and thereby establish the connection to asymptotic symmetries, which are infinite dimensional.At every point of the boundary, there sits a physical frame of reference for Poincaré translations, boosts and rotations.Since any full theory of quantum gravity should agree at least at low energies with the linearized theory, the insights gained from our analysis are expected to hold across different approaches.
As for the remainder of this article, we will start by explaining in Sec. 2 the emergence of boundary modes in gauge theory and gravity using the example of electromagnetism.Sec. 3 contains the main technical part of the paper.We review the perturbative expansion of the Einstein equations in the tetrad formulation before applying our analysis to linearized gravity in a bounded region.Readers that are not interested in the technical derivation may skip the first part of Sec. 3 and start with the summary at the end of Sec.3.2.After discussing the physical characterization of the bulk degrees of freedom-matter and radiation-we turn to the boundary modes in Sec. 4. We close our analysis by summarizing and comparing our results to previous works and provide a short outlook in Sec. 5

Boundaries Break Gauge Symmetries
In the following, we choose a finite domain D ⊂ Σ on a Cauchy surface Σ and characterize its phase space.The boundary C = ∂D separates the inside of the domain from the rest of the Cauchy surface.The presence of the boundary breaks the gauge symmetries.This is true both for the usual Yang-Mills type of gauge theories with internal symmetry groups, but it is also true for generally covariant theories, such as general relativity.In either case, new edge modes are excited at the boundary [2,29,30].In quantum theory, edge modes provide a neat tool to characterize the quantum entanglement of physical states across a codimension-two surface.Therefore, such a boundary C = ∂D is also often referred to as an entangling surface.An intuitive and often cited way to understand the emergence of such boundary modes is to consider Wilson loops in electromagnetism.This will be the focus of this section.In Sec. 3, we turn to general relativity.
Consider thus a closed loop α : S 1 → Σ on the initial hypersurface.The Wilson loop observable is invariant under U (1) gauge transformations that shift A into A − dλ.This is no surprise because the gauge parameter λ is not among the physical degrees of freedom in the bulk Σ. If, however, we restrict ourselves to a finite domain D and assume (see Fig. 2) that the entangling surface C cuts the loop α into two separate paths γ ⊂ D and γ ⊂ Σ − D, the restriction of the Wilson loop observable to the compact domain D is no longer invariant under the U (1) transformations.Instead, we now find where ε p is the relative orientation between the loop α and the entangling surface C at the puncture p ∈ α∩C.In this sense, the gauge parameter λ becomes physical at the boundary once we introduce an entangling surface that breaks the Wilson lines apart.The gauge parameter becomes a boundary mode dual to the charge excited at the puncture.Another way to make the same observation is to consider the classical phase space of a bounded region in the covariant phase-space formalism [31].While we point the reader to [24,32,33] for pedagogical introductions, we recall here the most important concepts and definitions that we need for the remainder of this work.If we start from the Lagrangian of the theory and introduce the phase space P via a Legendre transformation, the entire construction depends on a choice of time.This seems to be at odds with general covariance.The covariant phase-space formalism provides a way to reconcile the Hamiltonian description with covariance.Phase space P is built from the infinite-dimensional manifold of solutions to the field equations equipped with a symplectic form Ω, which is a two-form on the infinite-dimensional space.The two-form Ω provides a mapping between vectors and one-forms on field space (i.e., it can be used to lower indices).It does not, however, define a positive-definite inner product because it is skew-symmetric and may have degenerate null (i.e., gauge) directions.Given a Lagrangian L, the symplectic form can be derived from the first variation of L. Consider a theory of dynamical fields Φ I on a d-dimensional manifold.The variation of the Lagrangian is given by where E I = 0 are the Euler-Lagrange equations of motion and ϑ denotes the symplectic potential.Notice that ϑ is a one-form on field space and a (d − 1)-form on spacetime.It will be convenient to distinguish vector fields δ from exterior derivatives d on field space.When they act on the fundamental fields, the two are related via dΦ I (δ) = δΦ I .From the symplectic potential, one then constructs the symplectic current by taking the exterior derivative dϑ.The symplectic two-form Ω derives from the latter through Ω = Σ dϑ, where the integral is defined over a Cauchy hypersurface Σ ⊂ M. Note that although the integral is defined over a Cauchy hypersurface, it is independent of which Cauchy hypersurface we choose, because ddϑ = d 2 L = 0 on shell.Thus, the choice of Σ does not break the manifest covariance of the formalism.Finally, let us point out one important caveat in the derivation of the symplectic form given the Lagrangian: in the presence of symmetries, Ω is not invertible.It is therefore more correctly referred to as a pre-symplectic form.To obtain an invertible symplectic form, one has to quotient out the symmetries of the phase space [34].We will, however, simply work with the presymplectic form for the remainder of our discussion.
Returning to the example above, we are interested in the symplectic structure of a bounded region in electrodynamics.Considering the variation of the Lagrangian for Maxwell theory, ), where F = dA is the field strength and ⋆ is the Hodge dual, one finds that To describe the phase-space structure of a bounded region, we further restrict the integral of the presymplectic current to a compact domain D. The resulting presymplectic two-form is Consider now an arbitrary gauge-fixing condition such that we can write the vector potential A as the sum of a gauge-fixed solution A gf , which captures all the physical degrees of freedom in the bulk, and some arbitrary gauge transformation λ.In other words, Inserting this decomposition back into (6), and taking into account the field equations d ⋆ F = 0, we obtain The first term describes the geometry of the gauge-invariant phase space for the electromagnetic field in the region D. The boundary term describes the geometry of the edge modes.Its simple functional form tells us that the gauge parameter λ is conjugate, in the phase-space sense, to the electric field through the boundary.
The goal of the next two sections is to provide an in-depth analysis of the gravitational case for linearized fields, where the analog of the boundary modes λ will provide us with a local Poincaré frame at the boundary.
There is, however, a fundamental difference between the gravitational and the electromagnetic case.In electromagnetism, we can consider situations where there are no charged particles and the radiation field is confined to a finite pulse of compact support, localized within D. For any such configuration, the boundary term in Eq. (8) clearly vanishes.This is so because photons are charge neutral.They do not act as a source for radiation themselves.The gravitational case is fundamentally different.Gravitational radiation carries gravitational charge (energy, momentum, spin) and thus sources gravity as well.As soon as there is gravitational radiation in D, it will also excite a surface charge at the boundary C.As we will see below, this surface charge survives as a finite contribution to the presymplectic structure at the perturbative level.
3. Modes of Gravity: Radiative Fields, Coulombic Fields, Boundary Modes 3.1.Review: Perturbative expansion of the field equations in the tetrad formalism Let us now turn to the general relativistic case.Our starting point is gravity in the first-order tetrad formalism.The analysis could be done equally well in the metric formalism.One important difference is that working with tetrads, we obtain an additional gauge redundancy in the bulk, namely internal SO(1, 3) frame rotations of the frame field.This leads to additional reference fields for the Lorentz group at the boundary.These would be absent in metric gravity.In the context of this paper, the main advantage of working with tetrads rather than the metric is that tetrads are tailor-made to deal with how boundaries couple to the fields inside the manifold. 2In tetrad gravity, the geometry is encoded into an orthonormal frame e α that determines the metric tensor: In addition, there is the spin connection A α β with corresponding covariant exterior derivative Here, internal Lorentz indices α, β, . . .are raised and lowered with the flat (−+++) Minkowski metric η αβ , a, b, c, . . .are abstract spacetime indices and [•, •] is the Lie bracket for the Lie algebra so (1,3).In this language, solving the Einstein equations amounts to imposing the equations where T α is the energy-momentum three-form of matter, which is the Hodge dual of the usual stress energy tensor T ab , i.e., In addition, * F α β denotes the internal Hodge dual of the so(1, 3)-valued curvature two-form, i.e., β denoting the curvature of the connection and ϵ αβγδ is the Levi-Civita tensor with respect to internal indices with conventions ϵ 0123 = 1.Furthermore, ε abcd = ϵ αβγδ e α a e β b e γ c e δ d is the usual Levi-Civita tensor in spacetime indices.In flat space, with vanishing energy momentum tensor, i.e., T ab = 0, the general solution to the field equations is simply given by where the coordinate functions X µ : M → R 4 determine an inertial frame of reference on spacetime M and Λ α µ : M → SO(1, 3) is a local Lorentz transformation.In the following, we consider a perturbative expansion of the field equations around such Minkowskian solutions.We are interested in a decoupling limit with weak interaction and self-interaction among gravity and matter.
The perturbative expansion is written as a formal power series with respect to the coupling constant ℓ = √ 8πG, which is the Planck length.We thus introduce a perturbation and expand the remaining dynamical fields, i.e., the tetrad perturbation f µ , the difference tensor ∆ µ ν , and the energy momentum three-form T µ into a formal power series with respect to the coupling constant ℓ.We thus write Notice that the expansion of the energy-momentum three-form starts out one order lower than the expansion for f µ and ∆ µ ν .This is so because at lowest order in the expansion, the energymomentum tensor is independent of the metric perturbation.
It is then also useful to decompose the tensorial coefficients that appear in this expansion, ( 17) and (18), in terms of the zeroth order of the frame field, i.e., we write Here, all indices µ, ν, ρ, . . .are treated as internal, i.e., they are raised and lowered with the signature (−+++) flat Minkowski metric η µν , e.g.f Note that, at each order in the perturbation, any antisymmetric part of f (n) µν can be reabsorbed back into a redefinition of the local Lorentz transformation Λ α µ .Hence, we can always assume without loss of generality We are now ready to consider the perturbative expansion of the field equations.At lowest order n = 1, the field equations are * d ∆ where the operator * only acts on internal indices, i.e., * d ∆ where we introduced the effective energy-momentum three-form and torsion two-form of the order Here, [•, •] denotes the Lie bracket between so(1, 3)-valued p-forms, e.g.
The source terms on the right-hand side of ( 25) and ( 26) are a consequence of the self-interaction of gravity-perturbations of the gravitational field at first order are themselves a source of gravity at second order.
The linearized field equations of the coupled system can now be solved in the usual way, see, e.g., [35] for a detailed analysis with pointlike sources.Taking advantage of the gauge freedom of the field equations, 3 the general solution of (23,24) is a superposition of plane waves with amplitudes a ± ( ⃗ k), polarization tensors m µ m ν , and four-momentum k µ , where, in terms of, e.g., stereographic coordinates (z, z) on the two-sphere in momentum space, To solve the second-order equations ( 25) and ( 26) by the method of Green's functions, it is useful to split the perturbation of the spin connection into two terms that take values in so(1, 3), Here, the first terms is defined as the solution of the equation which can be solved algebraically by In the same way, the components of C (2) µ ν = C (2) µ νρ dX ρ can be obtained algebraically by inverting (28).Writing θ We thus have Given the standard Lorentz gauge conditions Eq. ( 35) amounts to the Einstein equations at second order in the perturbation, i.e.,

□ f
(2) 3 A gauge transformation of fµν can be always reabsorbed (order by order in the coupling constant) into a shift of the coordinate scalars X µ and a redefinition of the Lorentz gauge transformation Λ α µ .The perturbation (29) satisfies the usual transverse traceless gauge conditions ∂µ f where are the components of the effective order n = 2 energy-momentum three-form and □ = ∂ µ ∂ µ is the D'Alembertian.Assuming no incoming radiation from past null infinity at order n = 2, the solution to (35) for the gauge conditions (36) is thus given by where G R (X) = δ(X 0 − | ⃗ X|)/| ⃗ X| is the retarded Green function.

Perturbative expansion of the symplectic structure in a finite region
We now turn to the perturbative expansion of the presymplectic two-form.The gravitational sector is governed by the Hilbert-Palatini action: To add explicit matter sources, we add the action for N point particles minimally coupled to the gravitational field.The resulting field equations are at best formal.The reason is intuitively clear [36].If we try to concentrate a given mass into a point, we will eventually create a black hole before reaching the point particle limit.A possibility to circumvent this problem is to shrink down the body and take at the same time the ultrarelativistic limit, where the mass goes to zero, but the energy, with respect to some preferred observer at infinity, is kept fixed.In this way one can find exact distributional solutions to Einstein's equations; see for instance [37,38].In the following, we are dealing with perturbative gravity, where we can safely ignore this problem and work with distributional point sources.We can then formally solve the Einstein equations order by order in the coupling constant.
The action for the coupled system is the sum of the gravitational action and the action for the N point particles.The total action is thus given by where p i α is the four-momentum of the ith particle, γ i denotes its worldline in M, and N i denotes a Lagrange multiplier that imposes the mass shell condition with rest mass m i .The resulting energy momentum three-form for the matter Lagrangian is distributional and given by the variation of (40) with respect to the tetrad.That is, where δ(4) p is the four-dimensional Dirac delta distribution peaked at p ∈ M (a scalar density of weight one) and εabcd is the inverse Levi-Civita tensor density.Next, we consider a perturbative expansion around the flat Minkowskian solution given in Eqs. ( 13) and (14).We write p α = Λ α µ p µ , letting indices α, β, γ, . . .refer to the tetrad e α , wile indices µ, ν, . . .from the second half of the Greek alphabet refer to the flat background coordinates {X µ }.We expand the particle action in powers of the coupling constant ℓ = √ 8πG and obtain where q µ (τ ) denotes the coordinates of the trajectory γ : τ ∈ R → M ∋ γ(τ ) with respect to the coordinate frame {X µ }.If we want to solve the field equations order by order in ℓ, we would now also need to consider a perturbative expansion of the particle trajectory q µ (τ ) in powers of the coupling constant, i.e., q µ (τ ) = q (0) µ (τ ) + ℓ q (1) µ (τ ) + . . ., see for example [36], where such an approach is developed in great detail to first order in ℓ.In this paper, we are only interested in the zeroth leading order of the symplectic structure.The first contribution to the symplectic structure from the particle trajectory is O(ℓ 0 ).Hence, it is not required for the present purpose to consider higher-order perturbations of the trajectories.
The first variation of the action gives the Einstein equations with a distributional energy momentum tensor (41) in the bulk plus a boundary term.The boundary term determines the presymplectic current ϑ.Schematically, we have where the term EOM(δ), which is linear in the field variation δ, vanishes provided the equations of motion (10,11) are satisfied.Given the total action (40) of the coupled system, the presymplectic current ϑ is given by4 where p i a = p i α (τ )e α a γ i (τ ) is the four-momentum of the ith particle and denotes the variation of its position, which defines a T γ i M-valued one-form on configuration space.
Our goal in this section is to establish a full characterization of the gravitational phase space in a compact domain D in perturbative gravity.We choose the region D large enough such that it contains all particles i = 1, . . ., N and such that the n = 1 radiative modes (1) f have compact support within D. This simplifying assumption is illustrated in Fig. 3. Next, we integrate the presymplectic potential ϑ along D and identify the contributions from matter and radiation and the additional boundary modes.Taking into account the torsionless condition (11) and the expansion (15,16) of the configuration variables around the Minkowski solution, we can collect all variations of Λ α µ into a boundary term where S µν is the so(1, 3) spin current at the boundary: (a) (2+1)-dimensional representation.
Figure 3: A subregion of spacetime, delimited by a timelike cylinder and two subregions D and D ′ on Cauchy hypersurfaces Σ and Σ ′ respectively.The worldline of a single point particle within this region is depicted by γ.In this section, we determine the bulk and boundary phase space on D in the linearized regime.In the derivation of the presymplectic form, the first-order perturbation of the tetrad field (1) f µ , which accounts for gravitational waves, is assumed to have compact support within D. An example for such a subregion of spacetime is represented by the shaded area in the Penrose diagram.In the white areas, f The fact that the SO(1, 3) gauge element Λ α µ appears only through a boundary term is a consequence of the SO(1, 3) gauge invariance of the action.
Let us now simplify the expression (46) further.Our initial assumption was that the free radiative data f (1) µ has compact support within D. All boundary terms containing f (1) µ will thus vanish.It is for this reason, in fact, that no O(ℓ −1 ) term appears in the definition of the SO(1, 3) spin current S µν ; see (47).For the same reason, the second term in (46) also vanishes.Indeed, taking into account [d, d] = 0, and repeatedly using Stokes's theorem, we find where we set all terms to zero that vanish thanks to the Einstein equations ( 23) at order n = 1.
We are now in the position to turn our attention from the presymplectic potential to the perturbative expansion of the presymplectic two-form on D. The presymplectic two-form Ω D is the exterior derivative of the presymplectic potential, i.e., Going back to the perturbative expansion of the presymplectic potential (46), we obtain where S µ ν is the so(1, 3)-valued spin current (47) and we introduced the Maurer-Cartan form in field space, which is dual to it, i.e., The bulk terms can be neatly reorganized.First of all, we introduce a field-space covariant derivative D given by the ordinary exterior derivative d on field space shifted by a diffeomorphism,5 Here, X a is a vector-valued one-form on field space, explicitly, and L is the Lie derivative satisfying L ξ (•) = ξ⌟(d•) + d(ξ⌟•) for all ξ a ∈ Γ(T M).Next, we reorganize the various terms.Consider first the second term in (50).Taking into account the linearized Einstein equations at order n = 2, we can rewrite it as We insert this expression back into the presymplectic structure and identify three separate terms at order O(ℓ 0 ) or lower.Let us first list them and then comment on their physical significance in the next subsection.We have The first term is the matter contribution: The second term is the radiative part of the presymplectic potential: with t (2) µ denoting the vector-valued three-form (27).A straightforward algebraic manipulation (see Appendix A) simplifies this expression, which can be cast into Finally, there is the boundary term where S µ ν is the so(1, 3) spin current (47) and P µ is the local momentum -Summary.Before we turn to the interpretation of these terms, let us briefly summarize what we have done so far.Starting from general relativity in terms of tetrads, we performed an expansion in terms of the coupling constant ℓ = √ 8πG up to second order.We considered gravity coupled to N massive point particles in a bounded region D and derived the perturbative expansion of the presymplectic two-form Ω D on phase space.Assuming that the region is big enough such that gravitational radiation (1) f does not reach the boundary and using the linearized equations of motion, we simplified the expression to obtain three terms, corresponding to matter modes (56), radiation modes (58), and additional boundary modes (59).

Modes of gravity
Let us now explain the physical interpretation of the various contributions to the presymplectic two-form Ω D .First of all, we have the contribution from the matter degrees of freedom (56).This contribution describes the phase space of the point particles.The particles' states (points on phase space) are labeled by the particles' position and their four-momentum p i µ on the t = const initial surface Σ.Notice that the variation of the particles' position is dressed by a coordinate transformation replacing the differential dq µ i by Dq µ i , see (56).This dressing ensures that the symplectic structure of the point particles is invariant under simultaneous diffeomorphisms φ ∈ Diff(M : M) transforming both tetrads and trajectories into e α → φ * e α and γ i → φ −1 • γ i .As the coordinate functions transform as scalar fields, i.e., X µ → X µ • φ, the dressed differential Dq is invariant under such combined transformations of the trajectories and the fields in spacetime.
Next, there are the radiative modes (58).Variations of the first-order radiative perturbations f (1) µ and ∆ (1) µ ν appear in the presymplectic two-form only through the covariant field differential D = d − L X .Going back to the mode expansion (29), we find and similarly for ∆ (1) µ ν .The symplectic structure of the radiative modes thus depends only on the variations of the true physical degrees of freedom, i.e., the modes a ± ( ⃗ k) that characterize the strength and polarization of gravitational waves in the given gauge.From an operational point of view, the gauge condition that we implicitly used to parametrize the radiation field (29) states that the X i = const coordinate lines are trajectories of test particles at rest.These trajectories are geodesics up to terms of order n = 1 in the perturbative expansion.Gravitational waves will change the relative physical distance between these geodesics, but leave their coordinate distance fixed.
Finally, we have the boundary contribution (59).This is where the otherwise irrelevant gauge modes become physical.There are additional boundary modes X µ dual to diffeomorphisms and SO(1, 3)-valued boundary modes that are dual to the Lorentz charges S µν .These will be the focus of the next section.

Boundary phase space
As we have seen in Sec. 2, by the example of electromagnetism, restricting the phase space of a gauge theory to a bounded region leads to additional boundary modes.The U (1) gauge parameter λ, which is redundant in the description of the physics of the bulk, becomes physical at the boundary.Mathematically, this is realized by the emergence of boundary contributions to the presymplectic two-form.
In the case of general relativity in the tetrad formulation, we have two such gauge groups.The field equations ( 10) and ( 11) are invariant under four-dimensional diffeomorphisms and internal SO(1, 3) frame rotations.Inside a compact region D on a Cauchy surface Σ, it is impossible to distinguish configurations that are connected by small gauge transformations.Any such small gauge transformation can be written as a finite product of diffeomorphisms generated by the exponential map of vector fields ξ a 1 , ξ a 2 , . . .times the exponential of a local so(1, 3) gauge element λ α µ .For such transformations to be small, both λ α µ and the vector fields ξ a have compact support.If we extend them to the boundary, they become physical, generating a symmetry that maps a given state into an inequivalent one.The characterization of the physical states of the theory will thus include a description of additional boundary modes.In our setup, where we consider the perturbative expansion of the gravitational field in a neighborhood of Minkowski space, the boundary modes show up as the coordinate fields X µ and the SO(1, 3) gauge elements Λ α µ , arising from the internal Lorentz symmetry.The corresponding contribution to the presymplectic form is given by Eq. ( 59).
The coordinate fields X µ , which are inertial with respect to the background metric, act as physical reference frames with respect to which we can construct gauge-invariant (dressed) observables (see Fig. 4).This becomes particularly evident in our setup, where we consider point particles coupled to the gravitational field.In a generally covariant theory, the embedding of a trajectory γ i into the abstract spacetime manifold is meaningless without giving a physical prescription for how to measure the location of the trajectory.The physical reference fields X µ allow us to resolve this issue.The relative position of the trajectory with respect to the coordinate system X µ is gauge invariant.To leading order in the perturbation, the particles move along straight lines with respect to the reference frame X µ , i.e., (X µ • γ i )(τ ) = q µ o + m −1 i p µ i τ .Along the particles' trajectories, the reference frames become physical and serve as coordinate fields that give physical meaning to the otherwise gauge-dependent position γ i (τ ) of the particles.This relational definition is important, because it allows us to access the gauge-invariant part of the presymplectic two-form of the matter modes through the dressed variation of the particles' positions, The differential Dq µ i is gauge invariant, because a diffeomorphism will shift both the trajectories γ i of the particles and the coordinate fields X µ , in such a way that the dressed quantity X µ • γ i does not change.In the same way, the fields Λ α µ provide an internal reference frame for the Lorentz group at each point of the boundary.Take, for example, the SO(1, 3) curvature tensor F α β ; its Λ-dressed version Λ µ α F α β Λ β ν is clearly SO(1, 3) gauge invariant.The coordinate fields X µ have an additional role.Besides serving as dressing fields that give physical meaning to radiative modes and particle trajectories, they also define the location of the boundary itself.To the zeroth order of the expansion, the coordinate fields X µ define an inertial frame of reference, which embeds the abstract boundary of D into flat Minkowski space R 4 .In this sense, the boundary modes X µ provide a concrete realization of the embedding fields of [5].Given the presymplectic two-form at the boundary, see (59), we define an extended phase space that contains the conjugate boundary observables that are dual to the boundary reference frame X µ and the Lorentz fields Λ α µ .First of all, we have the coordinate fields X µ at the boundary, which are dual to the local momentum If we send the boundary to infinity, the integral of the local momentum (62) along the asymptotic two-sphere returns the Arnowitt-Deser-Misner (ADM) linear momentum [28,[40][41][42][43][44].The probably simplest way to see this is to use the standard representation of the ADM momentum in terms of the electric part of the Weyl tensor [28].Take hyperbolic coordinates in the vicinity of spacelike infinity, Next, we integrate (62) along a t = const and ρ = const surface.Performing a partial integration and taking into account the falloff conditions of the Coulombic fields, we obtain where t µ = ∂ t X µ and n µ = ρ −1 X µ are the timelike and spacelike normal vectors to the asymptotic two-sphere, d 2 Ω is the round area element d 2 Ω = sin ϑ dϑ dφ, and C (2) µνρσ n ν n ρ denotes the n = 2 leading order of the electric part of the Weyl tensor, which falls off like O(ρ −3 ) in the outside region Besides the local momentum, there is also a spin current (47) at the boundary, which is dual to the local Lorentz elements Λ α µ .In the linearized regime considered here, the local momentum and spin current P µ and S µν depend only linearly on the Coulombic modes f (2) µ , which are themselves determined by the effective energy momentum tensor via (38) and (26).
A few further observations.Using the equations of motion, it is straightforward to see that both the local momentum and angular momentum are conserved.Throughout this work, we assumed that the n = 1 radiative modes f (1) µ have compact support in D. Together with the second-order Einstein equations (25), this implies immediately that In the same way, we have conservation of the local angular momentum, which is the sum of the local orbital angular momentum and the intrinsic spin density (47), To demonstrate that the local angular momentum J µν is conserved in the outside region Σ = Σ−D, see Appendix B, where we show that It is straightforward to see that 0 returns the Beig-Ó Murchaddha relativistic center of mass (the global boost generator) [44].
Finally, we invert the boundary symplectic structure and obtain the commutation relations, where ⃗ ζ and ⃗ ζ ′ label points on the boundary surface C. In here, the local four-momentum Pµ and the spin current Sµν are now treated as surface densities intrinsic to C. The connection to the fields in the interior of the manifold is imposed by the momentum and spin-momentum constraints where the boundary Hamiltonian density Hµ and the spin surface density Sµν are given by the pull back of the two-forms ( 62) and ( 47) to C, i.e., We have thus found an extended phase space of bulk and boundary modes and a set of additional constraints ( 71) and ( 72) that establish the connection between the theory at the boundary and the theory in the interior.

Implications for quantum gravity
The structure of the phase space that has emerged from the discussions in the previous sections has important implications for the quantum theory.As any full nonlinear and nonperturbative theory of quantum gravity must agree at low energies with the linearized theory, there is no need to select a particular model or approach to quantum gravity.At low energies, foundational questions are shared across different approaches.
Through our preceding analysis, we found a neat separation of the classical state space into matter modes (p µ i , q µ i ), radiative modes a ± ( ⃗ k) and boundary modes X µ and Λ α µ with corresponding momentum currents Pµ and Sµν .A similar decomposition must also happen at the quantum level.At the linearized level, 6 we expect a kinematical state space that admits a tensor product structure The first term is the Hilbert space for the matter degrees of freedom, i.e., the tensor product of the one-particle Hilbert spaces for each constituent system.The second term describes the radiation modes.Inverting the presymplectic two-form at leading order in the coupling constant, we obtain the Poisson bracket between the positive and negative frequency modes a ± and ā± : Notice that the G → 0 decoupling limit removes the nonlinear higher-order terms on the righthand side.In this limit, the Hilbert space of the radiation modes is generated from the usual Fock vacuum |0⟩, Gravitons are excited from the vacuum by creation operators a To obtain an element of H rad D , the test functions f s ( ⃗ k) must be such that their Fourier transform has compact support in D.
Besides matter and radiation, there is the boundary Hilbert space H ∂D , which carries a representation of the canonical commutation relations (68), ( 69), (70).At a formal level, we may choose a functional Schrödinger representation, for which the momentum operators act as functional derivatives, i.e., whereas X µ ( ⃗ ζ) and Λ α µ ( ⃗ ζ) act by multiplication and ⃗ ζ = (ϑ, φ) are fiducial coordinates at ∂D.So far, the entire construction is formal.The obvious difficulty is that it is not enough to have a representation of the canonical boundary Heisenberg commutation relations.We also need an inner product.The formal infinite-dimensional Lebesgue measure gives rise to a vast Hilbert space, which is nonseparable.This is hardly surprising.The algebra of boundary charges is ultralocal.There is no boundary-intrinsic Hamiltonian H ∼ X 2 + P 2 with respect to which we could built a (separable) boundary Fock space.At first sight, this seems to be a major problem ahead, but we will argue below that the presence of constraints will render this problem rather fictitious.At the level of the kinematical state space, no Fock space structure is required.Only at the level of the physical phase space, do we expect to obtain a separable Hilbert space.This is reminiscent of a similar construction in the context of Loop Quantum Gravity (see [45,46]), where the constraints are imposed on an auxiliary kinematical Hilbert space, which is nonseparable.
Let us now turn our attention to the constraints ( 73) and ( 74) that establish the correlation between the bulk and boundary fields.In quantum theory, there are two ways to look at such constrained systems: solving the constraints before quantizing (reduced quantization) or quantizing first and then imposing the constraints at the level of the wavefunction [47][48][49].In the latter case, one starts with a kinematical state space with some auxiliary inner product.Physical states are identified by the requirement that the matrix elements of the constraints vanish between any two such physical states.If the constraints are first class, we impose instead the even stronger condition that the constraints annihilate all physical states [50].Given the specific form of the constraints ( 71) and ( 72) as a difference of a canonical momentum and a complicated function of the bulk configuration variables, the constraints can be given a particularly neat interpretation.They provide a relational evolution for the wavefunction of the modes in the bulk with respect to the boundary modes X µ ( ⃗ ζ) and Λ α µ ( ⃗ ζ).More concretely, consider the quantization of the constraint (71) on the physical states in the configuration space basis.Assuming, for a moment, that all constraints are first class, we obtain a local Schrödinger equation at the boundary, i.e., Here, the dependence of the wavefunction Ψ phys [X, Λ] ∈ H bulk ⊗H ∂D on the bulk degrees of freedom is left implicit.Notice that there are four such equations for every point ⃗ ζ at the codimension-two boundary C = ∂D generating local translations in X µ ( ⃗ ζ).Such multi-fingered evolution equations frequently appear whenever we consider the relational evolution of gravitational observables in terms of material frames of Refs.[49,51,52].However, there are two crucial differences to the usual formalism.First, the reference fields are now built from gravity itself, 7 without introducing any matter.Second, the constraints live in one dimension less, at the codimension-two boundary of D, where they describe the relational evolution of the wavefunction in the bulk relative to the coordinate fields X µ ( ⃗ ζ) at the boundary.There is, however, a subtle difficulty with imposing (80) strongly.It is not at all obvious whether the constraints ( 71) and (72) are first class, i.e., whether the (classical or quantum) algebra of the constraints closes under the bracket. 8In fact, our expectation is that they do not close at finite distance and that there is an anomaly due to cross terms such as At spacelike infinity, when we send ρ = X µ X µ to infinity, this problem disappears due to the falloff conditions-the cross terms go to zero.However, this limit also creates new difficulties.The boundary mode X µ diverges as ρ → ∞.To remove this divergence, we proceed as follows.Split a generic reference frame {X µ } into a divergent part and a finite shift where ⃗ ζ = (ϑ, φ) are coordinates at the asymptotic two-sphere.The divergent part consists of a global (yet arbitrary) Lorentz transformation Ω µ ν in SO(1, 3), i.e., dΩ = 0, dΩ ̸ = 0, and a fiducial reference frame, which is kept fixed once and for all, i.e., dX µ o = 0.The global SO(1, 3) transformation is required to have access to the global ADM angular momentum and boost charges.The finite translation Q µ ( ⃗ ζ), on the other hand, is an arbitrary Lorentz vector-valued function on the asymptotic two-sphere.The resulting IR-regularized variation of the boundary modes is then given by For any resulting reference frame X µ , there is a corresponding metric tensor In addition, we also have the fiducial reference metric Clearly, both metric tensors are flat.Taking into account the falloff condition for the coordinate transformation (81), the two line elements can only differ by a tensor which is O(ρ −1 ).The map from X µ o to X µ defines therefore an asymptotic symmetry, namely an asymptotic BMS transformation [25][26][27][28] in standard Minkowskian coordinates.
The analogous regularization for the SO(1, 3) spin frame of reference is where λ µ ν ( ⃗ ζ) is an so(1, 3) gauge element intrinsic to the two-surface boundary C = ∂D.It is possible to show that, given this specific falloff, the boundary presymplectic two-form (59) has a finite limit to spacelike infinity. 9The resulting boundary phase space is the cotangent bundle of the BMS group times an ultralocal contribution [T * so(1, 3)] C .A detailed analysis of this phase space is beyond the scope of this paper and will be presented elsewhere.For the following discussion, the important point is that this construction leaves us with an algebra of constraints that is first class.
For concreteness, let us restrict our final remarks to the BMS part of the asymptotic symmetries.The case of asymptotic and internal frame rotations (83) follows in complete analogy, but more care is needed because of the IR divergence induced by Ẽµν , which diverges quadratically as ρ → ∞, see (74).A generic state in the kinematical Hilbert space K will now be entangled between the bulk and boundary degrees of freedom.It will be constructed from a quantum superposition of bulk states Ψ bulk [Q µ , Ω µ ν ) ∈ H bulk that depend parametrically on the asymptotic BMS frame of reference,10 i.e., where dµ Ω is the left-and right-invariant measure dµ Ω = Tr( 6 Ω −1 dΩ) on SO(1, 3) and D[Q µ ] is a formal Lebsgue measure for the translational subgroup, which is infinite dimensional.At this point, the construction is formal, as we have not specified the Hilbert space inner product between the states |Q µ , Ω µ ν ⟩.Still, we can proceed to discuss the basic structure of the resulting physical Hilbert space.Physical states are those that satisfy a multi-fingered Schrödinger equation with respect to asymptotic BMS translations, where Hµ is the momentum density (73).In the same way, there will also be a constraint for the global Lorentz charge, where Λ▷Q µ = Λ µ ν Q ν is a left action and Ω µ ν ◁Λ = Ω µ ρ Λ ρ ν is a right action for any Λ µ ν ∈ SO(1, 3).The corresponding Hamiltonian is J o µν , whose space-space and spacetime components are the ADM angular momentum and Beig-Ó Murchaddha center of mass.In our representation, where Jo µν is the two-surface density At the full nonperturbative and background-independent level, the dynamics will be encoded into a generalized projector P that maps kinematical states Ψ ∈ K into physical states; see [45,[58][59][60].As we see now, a remnant of this generalized projector survives the decoupling limit that we have studied in this paper.The projector onto physical states is the formal delta distribution on the constraints Formally, the physical inner product is defined by the matrix elements of the generalized projector on the kinematical state space Consider now two such kinematical states.Take, for example, a simple product state between the bulk degrees of freedom and a quantum reference frame peaked perfectly on the configurations Q 0 and Ω 0 , and an entangled state where the reference frame is in a quantum superposition of several such configurations.Assume a delta function normalization for the kinematical inner product, i.e., ⟨Q, Ω|Q ′ , Under this assumption, the formal matrix elements of the projector (89) become where defines a representation of the BMS transformation on the bulk Hilbert space that maps the BMS frame {Q i , Ω i } onto {Q 0 , Ω 0 }.This viewpoint resonates with approaches to quantum gravity such as loop quantum gravity [58,59,61], where there is a kinematical boundary Hilbert and an amplitude map that is defined via a generalized projector P .This projector can act on widely different boundary states.Following this logic, Eq. ( 94) can be used to construct the corresponding matrix elements of BMS transformations.The entire structure is, of course, formal, but, even at the formal level, the emerging geometry is clear.The kinematical state space K is a principal fiber bundle (K, H phys , BMS, P ), the base space is the physical Hilbert space (H phys , ⟨•|•⟩ phys ), which is isomorphic-but not in a canonical way-to the bulk Hilbert space, i.e., the tensor product between the relativistic Hilbert space of matter and the Fock space of the radiation modes.The standard fibers are the BMS frames (Q µ , Ω µ ν ) on which the BMS group acts freely; the projector onto the base manifold is P , which defines the physical inner product.Consider then a specific model of quantum gravity, defined by a projector P , see, e.g., [58,59,61] and calculate its matrix elements to infer, via Eq.( 94).If we take the matrix elements of P between asymptotic boundary states, it would be possible, in principle, to define via (94) the matrix elements of asymptotic symmetries at the full nonperturbative level.Our long-term goal is to understand in this way the IR structure of quantum gravity from within a full nonperturbative approach.
Furthermore, the sum implements what is called a quantumcontrolled reference frame transformation, transforming the quantum state of the bulk degrees of freedom with a different unitary map U BMS i→0 that depends on the configuration of the reference frame (Q i , Ω i ).It thus constitutes a quantum reference frame transformation in the sense of e.g., [16,18,20].To make this connection more explicit and to get a better intuition of this transformation in general, let us restrict ourselves to a change between two reference frames that differ only by translations.That is, consider states for which Ω i = Ω 0 .In this case, the expression simplifies to implementing a pointwise translation by Q µ i − Q µ 0 at each point of the boundary S 2 , whose action on the bulk degrees of freedom is defined by Hµ , which we assume is self-adjoint.The operator given in Eq. (95) generalizes the existing quantum reference frame transformations for the translation group [16,18] and Lorentz group [62] in two ways: by going beyond global translations to pointwise translations at the boundary and by identifying their action on the gravitational degrees of freedom within the linearized regime.Moreover, these transformations provide a further step towards constructing general quantum-controlled coordinate transformations and their action on the geometry, as considered abstractly in [63].In this respect, they are complementary to the transformations considered in [64,65], which allow one to change into the local inertial frame of a quantum particle, and [66], which implement general coordinate transformations but rely on an auxiliary system to identify the different branches of the superposition and do not include the action on the gravitational field.

Outlook and Discussion
Physics describes correlations among observable facts in terms of fundamental laws of nature.These laws frequently make use of auxiliary elements that are not themselves observable.In electrodynamics, the charge current j µ couples to the electromagnetic field via the vector potential A µ .In linearized gravity, the energy momentum tensor T µν couples to the gravitational field via the metric perturbation f µν .Both A µ and f µν are not directly observable.They are the channels through which component systems interact.The most elementary way to access these otherwise invisible channels appears when we cut the gauge system along a boundary.At the surface that separates the two component parts, gauge invariance is broken.The then necessary boundary modes carry a representation of the fundamental gauge group.
The purpose of this study is to provide a complete characterization of gravitational boundary modes in the simplified setting of linearized gravity in a finite domain.Our discussion started out in Sec. 2 with a brief overview of the situation in electrodynamics.In Sec. 3, we considered the perturbative expansion of the Einstein-Cartan equations coupled to a system of N point particles.At zeroth order in the coupling constant ℓ = √ 8πG, all propagating degrees of freedom, gravitational waves and matter, are washed away.At this order, we are left with empty Minkowski space.We saw that there are infinitely many ways to parametrize Minkowski space in terms of inertial coordinates {X µ : D → R 4 } and local Lorentz frames {Λ α µ : D → SO(1, 3)}.Are these different representations of Minkowski space all gauge equivalent?They are not.They are not, because there is a boundary ∂D = C, which breaks gauge invariance.At the boundary, the coordinate fields X µ are dual to the local momentum (60) while the Lorentz frame Λ α µ is conjugate to the so(1, 3)-valued boundary area two-form (47).For any finite domain D, these boundary currents have a non-vanishing limit as G → 0. In this sense, even empty space has infinitely many physical boundary modes, which provide reference frames for the gauge symmetries of gravity.
Next, we studied the perturbative expansion.At order n = 1, we have the free radiation field.At the next to leading order, we have the Coulombic fields, which are sourced by the effective energy-momentum tensor.Assuming no incoming radiation from past infinity at order n = 2, the second-order perturbation f (2) µ is fully determined by matter and radiation.Upon inserting the perturbative expansion back into the presymplectic two-form, we identified two terms.The symplectic structure is the sum of a bulk contribution and a boundary contribution.At lowest order in the expansion, the bulk symplectic structure depends only on the radiation field and the states of matter.The boundary symplectic structure is where the Coulombic part of the gravitational field enters.It enters the boundary term through the local momentum P µ and the boundary spin current S µν ; see Eqs. (62) and (47).The local momentum P µ ( ⃗ ζ) and the spin current S µν ( ⃗ ζ) are conjugate to the local reference frame X µ ( ⃗ ζ) and the internal spin frame of reference Λ α µ ( ⃗ ζ) at the boundary.Let us also stress that in our analysis the boundary modes emerge already from the zeroth-order flat metric.In other approaches, the boundary modes are introduced by hand.For example, in [4,5], they arise as embedding fields that give physical meaning to the boundary of the subregion.Another approach was discussed in [12,13], where they emerge from within the full nonperturbative theory through post-selection.In [2], edge modes are added to render the total presymplectic two-form gauge invariant.In a similar way, it is possible to extend the phase space such that the Komar charges form a closed algebra under the Poisson bracket [4][5][6].In our approach, nothing is added.The boundary modes (X µ , Λ α µ ) are the fields that parametrize the flat metric, which is treated as a q number as well.
Finally, we considered the quantization of the phase space.In a full nonperturbative theory of quantum gravity, all geometry will be subject to quantum fluctuations.In Sec. 4, we saw how to realize the low-energy regime of such a quantum geometry within perturbative gravity.Besides matter and radiation, there are the boundary contributions to phase space, which need to be taken into account in the quantum theory, too.Quantization then proceeds in a rather conventional way.Ignoring higher orders in the coupling constant, the positive and negative frequency modes of the radiation field turn into creation and annihilation operators of a free field on Minkowski space.Since all elements of phase space ought to be quantum, we also need a quantum representation of the reference frames (X µ , Λ α µ ).At the kinematical level, the boundary charges (P µ ( ⃗ ζ), S µν ( ⃗ ζ)) are fields intrinsic to the boundary, dual to the boundary reference frames (X µ ( ⃗ ζ), Λ α µ ( ⃗ ζ)).Physical states lie on the hypersurface defined by the momentum and spin-momentum constraints, i.e., (71) and (72).These boundary constraints are what remains of the Wheeler-DeWitt scalar, vector and Gauss constraint of the Einstein-Cartan Hamiltonian system in the linearized regime.The constraints lead to a multi-fingered Schrödinger equation, which describes a relational evolution of the quantum state with respect to the reference fields at the boundary.Since the reference fields are part of the phase space, different subregions are now realized as different configurations of the bulk plus boundary phase space.Depending on the boundary and falloff conditions at ∂D, the system may be open or closed.If it is closed, there exists a Hamiltonian that drives the time evolution with respect to the reference frames at the boundary. 11If the system is open, we can work with a modified bracket that takes into account dissipation [57].
At finite distance, a further subtlety arises.We expect that the charges are not conserved under the Poisson bracket.To resolve this issue of anomalies at finite distance, we considered the ρ → ∞ limit.On phase space, this limit is subtle, because X µ ( ⃗ ζ) diverges as ρ = X µ X µ goes to infinity.To remove this infrared divergence, we made the following observation.Although X µ ( ⃗ ζ) may diverge, a variation δX µ ( ⃗ ζ) ≡ δQ µ ( ⃗ ζ) may still be finite.In this way, we found a rather immediate derivation of the asymptotic BMS translations and their dual charges P µ ( ⃗ ζ).Given the conservation law (65) for the local momentum and taking into account the asymptotic expansion of the Green's function near spacelike infinity, it is easy to check that P µ is, in fact, conserved along the ∂ t X µ generators of the asymptotic hyperboloid ρ = X µ X µ → ∞.Thus, we obtain a familiar result: there are an infinite number of angle-dependent conserved charges at spacelike infinity [67].At the quantum level, these conservation laws turn into a multi-fingered Schrödinger equation.The quantum states in the bulk evolve with respect to a BMS quantum reference frame at the boundary.Thus, our analysis extends the existing work on quantum reference frames (QRF) to asymptotic symmetry groups, which are infinite dimensional.Moreover, it contributes to a growing body of recent results that demonstrate the close connection of QRFs to gravitational physics, and boundary modes, e.g., [13,64,66,[68][69][70][71].In particular, QRF transformations that realize global translations and rotations [66], global Lorentz transformations [62], and local coordinate transformations [64,65] have been developed and applied to gravitational physics.Our work strengthens the connection between these different fields.First of all, we established the identity between boundary modes and QRFs explicitly in the linearized regime.Moreover, we extended the existing framework by providing the QRF transformations for the BMS group, which amount to pointwise translations and Lorentz rotations at the boundary.
Let us close with a few final remarks on how our results relate to nonperturbative approaches to quantum gravity.We mentioned above that even an empty region of space has infinitely many boundary modes.This observation is crucial.In nonperturbative approaches such as Spinfoams, Loop Quantum Gravity, Group Field Theory and Quantum Regge Calculus [45,58,59,[72][73][74], a quantum spacetime consists of a discrete number of flat building blocks that are glued together with deficit angles around codimension-two corners [59,[75][76][77][78][79][80].A smooth geometry emerges at large scales from the coarse-grained behavior of many such building blocks [81][82][83].By imposing flatness inside, all physical degrees of freedom are the boundary modes at the now internal two-surface defects.There is no radiation inside the fundamental building blocks; only the boundary modes are quantized.In perturbative gravity, one takes the opposite extreme-the radiation modes are quantum but the boundary modes are classical.Here, we steer a middle course.There is only a single boundary and we consider the quantum theory of both the bulk and boundary modes in the linearized regime.

Figure 2 :
Figure 2: The entangling surface C cuts a Wilson line α into two parts γ and γ, breaking gauge invariance, introducing boundary modes at the puncture.

Figure 4 :
Figure 4: The additional boundary modes for the subregion D: coordinates X µ and Lorentz frames Λ α µ at each point of the boundary.The coordinate fields serve two purposes: in the bulk, they allow us to define the diffeomorphisminvariant, relative position of the point particle moving along γ by expressing the path in terms of the coordinate fields as X µ •γ (indicated by the dotted line).In addition, they give physical meaning to the position of the boundary itself by providing an embedding into the flat background spacetime.