From Shockwaves to the Gravitational Memory Effect

We study the relationship between shockwave geometries and the gravitational memory effect in four-dimensional asymptotically flat spacetime. In particular, we show the 't Hooft commutation relations of shockwave operators are equivalent to the commutation relation between soft and Goldstone modes parametrizing a sector of the gravitational phase space. We demonstrate this equivalence via a diffeomorphism that takes a shockwave metric to a metric whose transverse traceless component is the gravitational memory. The shockwave momentum in 't Hooft's analysis is related to the soft graviton mode, which is responsible for the memory effect, while the shift in the shockwave position is related to the Goldstone mode. This equivalence opens new directions to utilize the gravitational memory effect to explore the observational implications of shockwave geometries in flat space.


Introduction
The gravitational memory effect measures the net relative displacement between asymptotic detectors induced by the passage of gravitational waves [1][2][3][4][5][6][7].More recently, in the context of high energy physics, it was shown that the memory effect can equivalently be understood as the leading soft graviton theorem in particle scattering amplitudes, as well as a vacuum transition between two inequivalent vacua in asymptotically flat spacetimes [8] (see [9] for a review).
In a complementary series of developments, 't Hooft proposed that incoming and outgoing modes near black hole horizons obey a set of commutation relations [10] and used these to evaluate a universal contribution to the black hole S-matrix.This takes the simple form of a phase shift and is directly related to the time delay acquired by a particle crossing a shockwave.The same contribution was later shown to be reproduced by resumming an infinity of graviton exchanges in the eikonal/high-energy limit of any 2-to-2 gravitational scattering process [11][12][13][14].These results triggered a fruitful, ongoing program of understanding black holes from shockwave dynamics [15][16][17][18][19][20][21][22][23][24].More recently, the 't Hooft commutation relations were utilized in a proposal for how quantum gravitational fluctuations could be observed in table-top interferometers [25][26][27].
A priori, the gravitational memory effect and the time delay experienced by a particle propagating in a shockwave background appear to be very different physical phenomena.Indeed, through its relation with the soft theorems, the gravitational memory qualifies as an infrared effect, while the time delay is revealed in the high-energy limit of scattering processes and therefore more naturally associated with the ultraviolet.The goal of this paper is to provide evidence that the two effects are nevertheless closely related.In particular, we find a diffeomorphism that relates a family of shockwave metrics to a family of Bondi metrics exhibiting a leading order memory effect.Restricting to the linearized theory is sufficient to establish this relation.This allows us to identify the shockwave variables appearing in the 't Hooft commutation relations with the soft and Goldstone1 variables [28] parametrizing the gravitational phase space of asymptotically flat spacetimes.We show that the 't Hooft commutation relations then follow directly from the commutation relations of the soft modes.
More specifically we consider a shockwave background in four spacetime dimensions.
Such shocks have been previously studied in a variety of contexts [10,20,[29][30][31][32], and are solutions to the Einstein field equations with a delta function source.We will focus on the radial shock sourced by a radially expanding shell of massless particles. 2The spherical shock reduces to the more familiar planar shock for two points with small transverse separation on the celestial sphere.
In the context of Schwarzschild horizons, 't Hooft argued that quantum mechanically there is an inherent uncertainty between the locations of ingoing and outgoing modes of Hawking particles [10,33].Here the ingoing probe particle experiences a shift X + due to the shockwave generated by an outgoing particle with momentum P + , while the outgoing shockwave similarly experiences a shift X − due to the ingoing probe particle's momentum P − .This uncertainty is expressed by a nontrivial commutation relation between ingoing and outgoing modes, as shown in Figure 1, which takes the form [P − (z, z), X − (z ′ , z′ )] = −iγ z z δ (2) (z − z ′ ) . (1.1) Note that, due to index raising and lowering conventions of lightcone coordinates, these commutation relations imply that P + is conjugate to X − ; the ingoing and outgoing conjugate and x + = t + r = u + 2r.The particle (or rather, shell of particles) generating the shockwave, depicted in red, is localized at x − = x − 0 and z = z 0 (the transverse directions are suppressed).The probe particle crossing the shockwave experiences a null shift X + that is related to the shockwave momentum P + .The shockwave also experiences a X − shift due to the momentum P − of the probe.
pairs in Figure 1 are accordingly shown in colored pairs.It has been previously suggested that the commutation relation (1.1) applies near generic horizons [13,25,26].In this paper, we will show that (1.1) is equivalent to the canonical commutation relations of soft modes in four-dimensional asymptotically flat spacetimes.To establish this equivalence, we show that a linearized Bondi metric having a leading memory effect can be obtained by a diffeomorphism from a shockwave solution to the linearized Einstein's equations sourced by T M uu ∝ δ(u − u 0 ).The equations of motion then imply a relation between the shockwave momentum P + and the soft graviton mode N, namely where is the transverse Laplacian (to be defined explicitly later).Physically, N measures the gravitational memory effect [8].We therefore establish a relation between the time delay acquired by a probe propagating in the shockwave background and the gravitational memory effect encapsulated concretely in the memory mode N(z, z). 3 We also find that the shockwave metric enjoys a supertranslation symmetry, which allows for shifts in X − to be traded for shifts of the Goldstone mode C parametrizing the memory metric, thereby allowing us to identify Given the identifications (1.2) and (1.3), we can immediately demonstrate that the commutation relation between soft and Goldstone operators, given by [28] z This equivalence has far-reaching consequences.The 't Hooft commutation relations in flat space were the foundation of a derivation in [26], which demonstrated that metric perturbations due to quantum effects are related to fluctuations of the modular Hamiltonian K in the vacuum state of an empty causal diamond, recovering the well-known area law , where A is the area of the bifurcate horizon of the causal diamond.
We view our results as evidence that such modular relations can indeed be applied to lightsheet horizons in flat space, and that modular fluctuations are related to soft effects in celestial holography (related directions were explored in [35]).These modular fluctuations, in turn, are at the center of a proposal of detectably large spacetime fluctuations in quantum gravity [25,36,37].Thus, our work opens new directions in utilizing gravitational memory to study observational signatures of quantum gravity.
The organization of this paper is as follows.In Section 2.1, we give a quick overview of a spherical shockwave solution, and show how by taking an appropriate limit we recover the more familiar Aichelburg-Sexl planar shockwave solution.In Section (2.2), we review the derivation of the memory effect in asymptotically flat spacetimes and its relation to asymptotic symmetries.In Section (2.3), we discuss a class of Bondi metrics with relaxed falloff conditions in the large r limit.In Section 3, we demonstrate that under a diffeomorphism, the shockwave metric considered in Section 2.1 is transformed into a metric that agrees with the Bondi metric up to overleading contact terms.This allows us to determine in Section 3.2 that the standard asymptotically flat memory metric is obtained via the same diffeomorphism from a shockwave solution to the linearized Einstein equations with a delta function perturbation turned on.The symplectic form of this new shockwave metric agrees background, and is therefore not in contradiction with the results of [34], where the "memory effect" for a planar shock is shown to vanish.
with that of the standard memory metric subject to an identification between the shockwave and memory variables.This allows us to determine in Section 3.3 how the soft graviton mode N and the Goldstone mode C are respectively related to the momentum P + and position X − of the outgoing shockwave and establish the 't Hooft commutation relations from the commutators of N(z, z) and C(z, z).Finally, we make some concluding remarks and speculate on future directions in Section 4.

Preliminaries
In this section, we give a brief overview of the background material necessary for our analysis.We begin by discussing some basics of shockwave metrics and reviewing the 't Hooft commutation relations in Section 2.1.In Section 2.2, we review the relation between supertranslations and the memory effect in the Bondi gauge.Finally, in Section 2.3, we introduce a class of metrics obeying slightly weaker boundary conditions compared to asymptotically flat metrics and hence containing the shockwave backgrounds.

Shockwaves and the 't Hooft Commutation Relation
An energetic massless particle in Minkowski spacetime sources a gravitational field referred to as an impulsive wave or a shockwave. 4Such solutions were first found by Aichelburg and Sexl [29] and later studied by Dray and 't Hooft in the context of an ultra-boosted particle, e.g. one moving along the horizon of a black hole [30].
Here, we consider a shockwave metric in retarded spherical coordinates X µ = (u, r, x A ), where u = t − r and x A = (z, z) are the stereographic coordinates on the sphere.This describes an outward propagating spherical shock, and is given by where 2 is the metric on the unit sphere, and D A is the covariant derivative with respect to the transverse space metric γ AB .Note that α(z, z) has dimensions of length.We will throughout this paper adopt the convention where lowercase Greek letters µ, ν, . . .are spacetime indices, and uppercase Latin letters A, B, . . .are transverse spatial indices.This metric describes a gravitational shockwave exiting future null infinity I + at retarded time u = u 0 .
To ensure that the metric (2.1) obeys the Einstein equations, the function α(z, z) is necessarily constrained.The only non-vanishing components of the Ricci tensor are where ≡ D A D A is the transverse Laplacian, and the Ricci scalar is Therefore, the Einstein equations, given by have only two nonzero components, which are Writing and substituting h uu from (2.1) into (2.5),we conclude that the metric (2.1) is a solution to the Einstein equations with a matter source where O(δ ′ , δ 2 ) indicates terms that are proportional to δ(u − u 0 ) 2 and δ ′ (u − u 0 ).We will be mainly interested in linearized solutions, in which case quadratic terms in the delta functions can be ignored. 5As for δ ′ (u − u 0 ), we argue that such terms integrate to zero on any definite u interval,6 and hence will not contribute to the total energy given by the u integral of T M uu .It is straightforward to show that ∇ µ T M µν = 0, so the stress tensor is conserved.In Section 3.2 we will introduce a closely related shockwave solution to the linearized Einstein equations with vanishing T M (2) ur .
We now briefly discuss how to obtain the familiar Aichelburg-Sexl planar shock studied in [29] from the metric (2.1).To this end, we want to zoom into a particular patch of the celestial sphere by taking z, z ≪ 1 (in this region γ z z → 2).If we then introduce the lightcone coordinates (x − , x + ) and transverse coordinate (x 1 , x 2 ), related to the retarded coordinates in (2.1) via the metric becomes where O(z, z) indicates terms that are higher order in z, z and r = 1 2 (x + − x − ).We recognize (2.9) as the Aichelburg-Sexl planar shock, which satisfies the Einstein equation where ⊥ is the transverse Laplacian in the x 1 , x 2 coordinates.
We conclude this subsection with a brief review of the 't Hooft commutation relations proposed to capture quantum dynamics of particles scattering near black hole horizons.
We will largely follow the presentation given in Section 11 of [10].Recall the lightcone coordinates (x − , x + ), which are related to the retarded Bondi coordinates (u, r) in (2.1) by (2.8).In these coordinates, the shockwave metric (2.1) takes the form which falls into the class of metrics analyzed in [30].However, whereas the analysis in [30] was done in a Schwarzschild black hole background, we emphasize that no black hole is present in our setup.We expect our analysis to apply to any sufficiently large causal diamonds with metrics that admit a large r expansion of the form (2.26) near the boundary [25,26,36,37,[40][41][42][43][44].Note that (2.12) Substituting (2.11) into the Einstein equations, we find the −− component to be7 To solve for α, let us define the Green's function for the transverse Laplacian8 It is clear then from (2.13) that a generic matter distribution in the transverse directions with stress tensor leads to a shockwave profile where we used (2.12) to switch p − to p + .For a single source particle of constant outgoing In the presence of the shock (2.16), a probe particle propagating along the future horizon (or equivalently along x − ) will acquire a time delay δx + for x − > x − 0 : This setup is illustrated in Figure 1.The magnitude of the time delay is determined by the wave equation in the shockwave background (2.11), and is given by [11] δx + = α(z, z) . (2.18) The propagation of this particle in the shockwave background is therefore characterized by the phase shift where in the second line we substituted δx + in terms of p + (z, z) according to (2.16) and (2.18).
Thus far, all the equations we derived above follow from the classical equations of motion.
By (2.18), we have a relation between δx + and p + , and there is no classical shift of the shockwave position x − .However, in the special case where the shockwave is sourced by a point particle, (2.19) may be regarded as the high-energy limit of a 2-to-2 scattering amplitude [11,[45][46][47].As such, it is clear that one would arrive at the same result (2.19) by changing reference frame, or equivalently, by exchanging the roles of the source and probe particles.The metric in this case takes the same form as in (2.11) with x + ↔ x − , namely a shock located at x + = x + 0 , and the propagating particle in this new background will acquire a shift paralleling (2.16).
This motivated 't Hooft to promote p ± , δx ± to operators P ± , X ± and postulate the commutation relations where we have suppressed the dependence of all operators on z for clarity.These commutation relations can then be used to recover the phase shift (2.19), which we review in Appendix A. Remarkably, we will later show in Section 3.3 that P − = − 1 2 P + is directly related to the leading gravitational memory operator [8,48,49], while X − = −2X + will be identified with the Goldstone operator involved in discussions of infrared divergences [50][51][52][53][54][55].

Asymptotically Flat Metrics and Gravitational Memory
In this subsection, we review the standard framework used to describe gravitational radiation from isolated sources in spacetime [8].We will assume that observers are located at a relative distance much larger than the extent of the sources, where spacetime is approximated by Minkowski space with metric g (0) µν .Gravitational radiation is characterized by linear perturbations h µν around this flat background, so that As appropriate for the propagation of outgoing gravitational waves, we will as in the shockwave case adopt retarded spherical coordinates X µ = (u, r, x A ).A convenient gauge that the metric can be expressed in is the Bondi gauge [56,57], which is defined to obey the condition as well as the determinant condition [58] ∂ r det Physically these correspond to requiring the normal vectors ∂ µ u to be null and the angular coordinates to be constant along null rays, or equivalently that gravitational waves propagate along the radial direction and have spherical wavefronts.
At large r, the perturbations h µν admit expansions of the form9 with the fall-offs of each component determined by the choice of boundary conditions.Solving the asymptotic radial Einstein equations subject to the standard fall-offs proposed in [56,57] (for a review see [9] and references therein) leads to the metric for an asymptotically flat spacetime near future null infinity I + Here m B is the Bondi mass aspect (note that, despite its label as a mass, it has dimensions of length, effectively having absorbed a factor of G N ), and • • • denote subleading corrections in the large r expansion.The gravitational radiation is characterized by the news tensor, related to the shear C zz by Therefore, spacetimes that can be written in the form of (2.26) are often referred to as radiative spacetimes, and if N zz = 0 then the metric describes a vacuum spacetime.Note from (2.27) that N zz = 0 does not imply the vanishing of C zz .In particular, imposing N zz = m B = 0, we obtain from the uz component of the vacuum Einstein equations that As such, vacuum spacetimes are parametrized by an arbitrary function C(z, z) on the sphere.
This variable is known as the Goldstone mode canonically conjugate to the zero mode of the news tensor or the soft graviton mode [9], as we will review below.Moreover, these different vacua are related by large diffeomorphisms called supertranslations.
Further imposing the constraint equations, broadly defined to be all the remaining Einstein equations, order by order in a large r expansion provides relations between m B and C zz , and also among further subleading metric components.For example, at O(r −2 ) the uu component of the Einstein equations takes the form where the stress tensor T uu (by convention dimensionless here) is defined to be with uu being the coefficient of the r −2 component of the matter stress tensor defined in (2.6).Consider now the ℓ = 0, 1 spherical harmonic modes of N zz .Recalling Y ℓ=0 ∝ 1 and 1+z z , it is easy to see that Thus, the linear terms in N zz and N z z drop out for the ℓ = 0, 1 modes, which means ∂ u m B cannot generically vanish unless T M uu fails to satisfy positive energy conditions in the ℓ = 0, 1 modes.Nevertheless, we will henceforth assume this, and leave the inclusion of a nontrivial Bondi mass for future work.
The asymptotic symmetries of asymptotically flat spacetimes are the transformations that preserve the form of the metric (2.26) and are known as the BMS symmetries [56,57].
A particular subset of BMS symmetries are the supertranslations, which are parametrized by a function f (z, z) on the celestial sphere and are generated by the vector fields Under a diffeomorphism parametrized by vector ξ, the coordinates and metric transform according to where we are using ∇ µ as opposed to D A for spacetime covariant derivatives, and symmetriz- Explicitly, under a supertranslation the coordinates shift via and the shear C zz and news N zz transform as [9,62] The first equation in (2.35) implies that supertranslations act non-trivially on vacua parametrized by (2.28), namely Gravitational flux characterized by N zz = 0 for u ∈ [u i , u f ] generally induces transitions between the early and late vacua at u < u i and u > u f .In particular, because this is no longer a vacuum spacetime, we do not expect (2.28) to hold, and the shear C zz is now a function of u.A sharp probe of this is the gravitational memory effect [1][2][3][4][5][6][7][8], where gravitational radiation induces a net relative transverse displacement in the trajectories of nearby detectors or probe particles.This displacement is directly related to the change in the shear between times u i and u f before and after the gravitational pulse, given by where we can view the last equality as the definition of the field N(z, z).Note that N(z, z) is directly related to the zero mode of the news and can be shown to correspond to the leading soft graviton mode upon quantization [28].We can determine N(z, z) by integrating the constraint (2.29), so that For the special case of a delta function localized stress tensor T uu ∝ δ(u−u 0 ) and m B = 0, we have N zz ∝ δ(u − u 0 ) by (2.29), which in turn implies by applying (2.37) over an integration from −∞ to u that where Θ(x) is the unit step function.Here N(z, z) is determined in terms of the flux by (2.38), and C(z, z) can be viewed as an initial condition, related to the vacuum C vac zz by (2.28).We emphasize that, according to (2.36), supertranslations induce shifts in C(z, z) while leaving N(z, z) invariant.On the other hand, the vacuum transition between early and late times induced by flux is itself parametrized by a supertranslation in the sense that Consequently, one can regard supertranslations as relating physically inequivalent configurations [8].
Solutions analogous to (2.39) in QED were discussed in [63].Choosing C(z, z) and N(z, z) such that they diagonalize boosts towards a point (z 0 , z0 ) on the celestial sphere yields solutions related to the conformally soft wavefunctions in celestial holography [51,64,65].
As argued in [28], the phase space of gravity in asymptotically flat spacetimes has to be augmented by N(z, z) and C(z, z).Upon quantization, the former correspond to insertions of soft gravitons while the latter are key constituents of the gravitational dressings [52][53][54][55].
As we review below, N(z, z) and C(z, z) are also canonically conjugate to one another [28,51].
Consequently, the role of the gravitational dressings is to supply the asymptotic states with the amount of soft charge (related to soft gravitons) necessary to ensure that the net large gauge charge is conserved in any scattering process [50,63].
The commutator between the fields C(z, z) and N(z, z) was first worked out in [28] by starting from the canonical commutation relations in gravity given by (also see Appendix C) [66, 67] and then imposing suitable boundary conditions.Note that from (2.28) and (2.39), it is clear that where S = (z−w)(z− w) (1+z z)(1+w w) and we used the identities (2.44) In the planar limit, the sphere is flattened to a plane (γ z z → 2) and S → |z − w| 2 .Consequently, (2.43) simply becomes (2.45)

Relaxing the Boundary Conditions
In this short subsection, we introduce a class of metrics obeying (2.23) and (2.24), but with h uu allowed to be finite in the large r limit.Notice that this class of metrics contains the shockwave metric (2.1), and it is more general than metrics obeying the fall-off conditions (2.26), where h uu ∼ O(r −1 ).These generalized boundary conditions, as well as a partial gauge fixing where (2.24) is not imposed and (2.25) is allowed to include logarithmic corrections, were extensively studied in [61].
Metrics obeying (2.23) take the form where in particular V, U A and β are allowed to obey weaker fall-offs at large r than those specified by (2.26) in a way consistent with the Einstein equations.Before imposing the determinant condition (2.24), such metrics enjoy an enhanced diffeomorphism invariance parametrized by where F = F (u, z, z) and ξ r is an arbitrary function of (u, r, z, z), and • • • indicate subleading terms in the large r expansion.
Imposing the determinant condition (2.24) fixes ξ r in terms of F , namely Substituting this back into (2.47),we see that a generic diffeomorphism preserving (2.46) is parametrized by a single function F , given by Notice that we recover (2.34), but with f promoted to a function F that may depend on u, due to not having demanded the vanishing of O(r 0 ) component of h uu .In this case, the vector fields (2.49) have non-vanishing divergence ∇ • ξ = ∂ u F = 0, which implies that the associated diffeomorphisms will in general change the trace h = h µ µ .The diffeomorphism considered in the next section falls into this class. 10 Relating Shockwaves to Gravitational Memory In this section we establish an equivalence on the future null horizon away from u = u 0 between a shockwave metric closely related to (2. and N implies the 't Hooft commutation relation (2.21).This extends the analysis of 't Hooft near a Schwarzschild horizon to null horizons at large r.

The Diffeomorphism
In this subsection, we show that (2.1) is diffeomorphic to a metric that satisfies the gauge conditions (2.23) and (2.24) and where the shear takes the form (2.39).Recall that under a diffeomorphism parametrized by vector field ξ, the metric transforms linearly via (2.33).
Requiring that this transformation eliminates the h uu component of the metric (2.1), we find where g(z, z) is an arbitrary function on the sphere.It is easy to see that Together with the requirement (2.49) that Bondi gauge is preserved under the diffeomorphism, F (u, z, z) is determined in terms of the shock profile α(z, z) and g(z, z) to be We show in Appendix B that diffeomorphisms parametrized by (2.49) with F (u, z, z) defined by (3.3), transform the shockwave metric (2.1) to a form closely related to the Bondi metric introduced in (2.26), namely Here C zz takes the form (2.39), and from (B.9) we see that N(z, z) is determined by the shockwave background to be and C(z, z) is related to the integration constant g(z, z) of the vector field in (3.1) via Notice that (3.5) differs from the usual Bondi metric (2.26) by terms proportional to delta functions, given in the last line of (3.5).As described in Section 2.3, (3.4) transforms the shockwave metric (2.1), which is traceless, to one with a nonzero trace, specifically h ur = 0.
The O(r) component of h uz is then necessarily nonzero to ensure that the Einstein equations are obeyed.
We conclude this subsection with a comment on the physical interpretation of diffeomorphisms parametrized by (3.4).Both the time-dependent and the time-independent compo-

Modifying the Shockwave Metric
The analysis in the previous section allows us to easily connect the standard asymptotically flat memory metric (2.26) with m B = 0 to a shockwave metric closely related to (2.1).
This can be achieved by subtracting the delta function contact terms from (3.5), which leads us to consider a shockwave metric of the form where The first line is precisely the spherical shockwave metric (2.1).Note that (3.8) differs from (2.1) both by the choice of Minkowski vacuum and by having a non-vanishing trace h ur = 0. Allowing the trace to be non-vanishing implies the Einstein's equations are obeyed without a ur-component of the matter stress tensor (see (3.10) below).The trace will be removed by the diffeomorphism generated by (3.4) with F given by (3.3).Indeed, because the diffeomorphism generated by the vector field ξ (3.4) is linear, the metric we obtain after applying the (linearized) diffeomorphism parametrized by ξ generating a small gauge transformation is precisely where C zz is defined in (2.39), and we used the identifications (3.6) and (3.7) with g(z, z) = 0 to ensure that ξ is generating a small gauge transformation.We refer to this metric as a memory metric, which are Bondi metrics with nonzero memory mode N(z, z).It is clear this metric is exactly of the form (2.26), which means the commutation relation (2.43) holds.
Evaluating Einstein's equations for (3.8) and using (3.6), the only nontrivial equations are Thus we see that as in 't Hooft's original planar shockwave analysis, the ur-component of Einstein's equation vanishes.Furthermore, integrating the uu-component of Einstein's equation yields for u > u 0 .This is precisely the momentum operator defined in 't Hooft's original analysis, given in (A.3).

From Soft to 't Hooft Commutation Relations
In this subsection, we demonstrate that (2.21) is in fact implied by the canonical commutation relations (2.43) obeyed by C(z, z) and N(z, z).Since the metric (3.8) is diffeomorphic to the Bondi metric (2.26), the canonical commutation relations (2.43) remain unchanged.
First, note that by substituting (3.6) into (3.11),we are able to immediately express P − in terms of N(z, z) for x − > x − 0 , namely This result is not as surprising as it may seem.Under our diffeomorphism, the stress tensor T M −− remains the same to linear order, while the Coulombic degrees of freedom of the shockwave metric are effectively shifted into radiative ones, as one can see by comparing the −− components of the Einstein equations (2.13) and (2.38).As such, the relation between P − and N in (3.12) could have been directly deduced from the equations of motion before and after the diffeomorphism.A completely analogous relation between the Coulombic and memory degrees of freedom in the high-energy limit of QCD was pointed out in [71].
To identify the shockwave variable associated to the Goldstone C(z, z), note that diffeomorphisms with non-vanishing g(z, z) in (3.7), and hence C(z, z), change the boundary vacuum configuration.The metric (3.8) generalizes the shockwave metrics considered by 't Hooft by allowing for C vac zz = 0.However, one can revert to the configuration originally considered by 't Hooft with C vac zz = 0 by a supertranslation (2.47) with F = C(z, z) acting on coordinates while keeping the metric fixed. 11In particular, recalling from (2.8) that x − = u, by (2.34) x − transforms by the simple shift Promoting this field-dependent coordinate shift to an operator, we are led to write δx − = X − , suggesting the identification of the 't Hooft variable X − with the Goldstone mode C via up to a field-independent variable.Note that the relative sign in this identification is consistent with the action of supertranslations (2.34) and (2.36) on X − and C, respectively.
Finally, since N is proportional to the transverse momentum mode P − via (3.12), we see that the commutator between C and N (2.43) implies a nontrivial commutation relation between P − and X − , which for x − > x − 0 is given by This exactly reproduces the 't Hooft commutation relation (2.21).We could have obtained the 't Hooft commutation relation involving P + and X + had we begun with an ingoing shockwave localized at x + = x + 0 rather than an outgoing shockwave.We can recast the commutation relation (3.15) into a more symmetric form.Recalling the classical equations of motion (2.16) and (2.18) and promoting them to operator equations, we have where G(z − z ′ ) is the Green's function for the transverse Laplacian given in (2.14).Substituting this relation between X + and P − in (3.15), we obtain (3.17) This reflects an uncertainty in the X − , X + coordinates of the shockwave and probe respectively in the quantum theory.

Discussion
In this paper, we established a relation between the time delay acquired by a particle propagating in a shockwave background and the gravitational memory effect.This implies that the memory mode, or equivalently the leading soft graviton mode, is related to the shockwave momentum introduced by 't Hooft in the context of scattering near a black hole horizon.Moreover, the canonically conjugate operator he postulated to correspond to the location of the black hole horizon may be identified with the Goldstone mode.Therefore, we are able to reinterpret the 't Hooft commutation relations as a simple consequence of the algebra of boundary operators derived from the covariant phase space formalism applied to asymptotically flat spacetimes.In particular, no black hole horizon needed to be present in our analysis.This was foreseen in [25,26] through the application of the 't Hooft commutation relations to a generic horizon of a causal diamond in flat space.
It is interesting to interpret this result from the perspective of soft gravitons.It is wellknown that the memory effect can be thought of as a vacuum transition [8], and this is associated with the production of soft gravitons.Here we see that this effect is directly related to the classical time delay of a particle crossing a shock, which is constrained by causality to be positive [72,73].By (2.16), (2.18), and (3.12), this should translate directly into a constraint on the memory effect (see [74] for a statement along these lines), and it would be interesting to further explore the implications of this relation.Moreover, the uncertainty in the location of the shockwave is tied to the Goldstone mode canonically conjugate to the soft graviton.Goldstone modes have previously been shown to appear in gravitational dressings [50,[52][53][54][55], suggesting a relation between shocks and coherent clouds of gravitons.
However, this relation remains quite mysterious to us, and we hope to return to it in future work.
Throughout this paper we restricted our analysis to the semiclassical regime: gravitational memory is a classical effect and we are only implementing canonical quantization.
However, one can promote the stress tensor to a quantum operator.This implies that the matter stress tensor can have non-zero fluctuations even in the vacuum where the expectation value vanishes.That is, the two-point function of T M µν is nonzero while the one-point vanishes, and this may correspond to nontrivial two-point functions involving the soft graviton and Goldstone modes N and C, respectively [49,52,75,76].
Our results suggest that we may be able to probe quantum gravity by measuring the memory effect sourced by quantum fluctuations of spacetime rather than classical gravitational waves.This is the essence of the proposal to detect observable effects in quantum gravity in [25-27, 36, 37, 44, 77, 78].Of course, to have any observable effects from the quantum metric fluctuations, it is not enough for a single shockwave to be sourced by quantum fluctuations.Rather, the memory effect must accumulate sufficiently, over multiple quantum shocks, in order to be observable.Furthermore, to truly tie our results to an observational signature, we need to connect the time delay in the shockwave geometry and transverse displacement associated with the gravitational memory effect with a gauge invariant observable as measured by an interferometer.An idea for how this could happen was at the heart of the proposals in [25,27].Nevertheless, by treating shockwaves as a gravitational memory effect, we open up additional directions to understand the infrared behavior of shockwave geometries, and we leave this study for future work.which corresponds precisely to (2.20).This is the time delay acquired by the source when propagating in the background of the probe.As such the commutation relations (2.21) capture the quantum uncertainty in the positions of the source and probe introduced by their scattering [10].

B Diffeomorphism Transforming the Shockwave Metric
In this appendix, we apply the diffeomorphism generated by ξ from (2.49) with F given in (3.3) to the shockwave metric (2.1), and show that we recover the metric (3.5).For convenience, we record here that the vector generating the diffeomorphism is Using (3.2), we get Recall we chose F to obey (B.1) such that we eliminate the h uu component in the shockwave metric (2.1).We now work out how the rest of the components transform.Recalling and that the nonzero Christofel symbols of the unperturbed metric are In particular, after performing the diffeomorphism, we have Comparing with (2.39) , which is we can identify Acting on both sides of this expression with + 2, we can derive from (B.1) that Finally, we observe that we can rewrite (B.5) in terms of C and N using (B.8), such that In terms of these fields, the resulting metric is In the limit where r → ∞, this reduces precisely to the metric studied in [26].

C Symplectic Form of Transformed Metric
In this appendix, we compute the symplectic form of the transformed shockwave metric Following [79][80][81][82], the symplectic form for a metric in linearized gravity is given by (we where ω µ is the symplectic current, Σ a Cauchy slice, dΣ µ the volume 3-form labeled by its orthogonal vector x µ , and δh µν ≡ g µσ g νρ δh σρ .For the case we are interested in, we want to take the Cauchy slice Σ to be I + , which is defined as the u + 2r = v 0 hyperplane with v 0 → ∞.This means that n µ = ( 1 2 , 1, 0, 0), or equivalently that the normal vector is n µ = (−1, 1  2 , 0, 0).The relevant components of the symplectic current are therefore ω u and ω r .It is straightforward to show that ω u = O(r −3 ) and hence its contribution to the symplectic form vanishes in the limit r → ∞.We are left with evaluating ω r .
To this end, we start by computing the nonzero components of the linearized metric Since √ −g = r 2 γ z z, we will be interested in the terms of ω r that fall off slower than r −3 .
Keeping only terms quadratic in h, we first focus on the terms in ω r that could lead to divergences when multiplying √ −g.We find Upon integration by parts over the z, z coordinates and antisymmetrization, this term drops out.
We conclude that the divergent terms cancel out and the symplectic form in the large r limit reduces to Note that after integration by parts and dropping total derivatives on the sphere, the first four terms simplify and we are left with For σ = 0, this agrees with (2.41) as expected.For σ given in (C.2), which was obtained from the shockwave metric (2.1) via the diffeomorphism constructed in Section 3.1 and N AB entirely "soft", i.e. associated with the shear in (2.39), the symplectic form (C.9) vanishes.This is consistent with the vanishing of the symplectic form for shockwave metrics of the form (2.1), which can be straightforwardly verified.

Figure 1 :
Figure 1: We depict the spherical shockwave spacetime in lightcone coordinates x − = t − r = u ) implies the 't Hooft commutation relation (1.1) between P − and X − .Had we considered an ingoing shockwave instead of an outgoing one, we would have established the 't Hooft commutation relation between P + and X + .

. 42 )
As a result (quantum commutators [•, •] are identified with the Dirac brackets i{•, •}), [N(w, w), C(z, z)] = 8iG N S log |z − w| 2 , (2.43) nents of the vector fields preserve the relaxed Bondi boundary conditions where h uu = O(1).The time-independent component preserves the asymptotically flat large-r fall-offs in (2.26) while modifying the boundary value of the vacuum shear (2.28), and is considered a large gauge transformation.As such, for g = 0, (2.1) and (3.5) describe different theories.One could arrive at (3.5) without performing a large gauge transformation (3.7) by considering a shockwave metric in an arbitrary Minkowski vacuum.This can be achieved by turning on C vac zz in (2.1), as we shall see in Section 3.2.On the other hand, the time-dependent component violates the asymptotically flat boundary conditions at large r and is therefore not typically included in the asymptotic symmetry group of asymptotically flat spacetimes.Analogous gauge transformations are known to trade Coulombic for soft/memory degrees of freedom in gauge theories [68-70], where they were classified as residual gauge transformations.Indeed, we show in Appendix C that the symplectic form of (3.5) vanishes subject to (3.6), and hence remains invariant under such time-dependent diffeomorphisms.We leave a complete understanding of time-dependent diffeomorphisms and their associated charges (if any) to future work.

(C. 7 )
Defining NAB ≡ N AB − D A D B σ ,(C.8)the symplectic form becomesΩ I + [δ 1 h, δ 2 h] = 1 16πG N I + √ −gω r n r du dz dz = − 1 32πG N I + γ z z δ 1 C zz δ 2 Nzz + δ 1 C z z δ 2 Nzz − δ 1 Nzz δ 2 C zz − δ 1 Nzz δ 2 C zz du dz dz = − 1 16πG N I + γ z z δ 1 C zz δ 2 Nzz − δ 2 C zz δ 1 Nzz du dz dz + • • • , (C.9)where we integrated by parts to obtain the last equality, and • • • denote additional boundary terms.From this we can immediately read off the Poisson bracket involving the bulk degrees of freedom to be C zz (u, z, z), Nww (u ′ , w, w) = 16πG N γ z zδ(u − u ′ )δ (2) (z − w) .(C.10) by(2.49) that relates the shockwave metric (2.1) to a metric that resembles (2.26) up to contact terms.This diffeomorphism provides a relation between N(z, z) in (2.39) and the shock profile α(z, z).Both metrics are shown in Appendix C to have vanishing symplectic form, suggesting the diffemorphism is "small."InSection3.2, the same diffeomorphism is shown to relate another shockwave metric to the Bondi metric(2.26).Both metrics now have non-vanishing symplectic forms that agree, subject to the identification between the shockwave and memory variables α and N.In Section 3.3 we show that the canonical commutation relation between the Goldstone and soft graviton modes C 1) and Bondi metrics of the form(2.26)withCzzgiven by(2.39).In particular, we first demonstrate in Section 3.1 that there exists a diffeomorphism parametrized