Bulk locality from modular flow

We study the reconstruction of bulk operators in the entanglement wedge in terms of low energy operators localized in the respective boundary region. To leading order in $N$, the dual boundary operators are constructed from the modular flow of single trace operators in the boundary subregion. The appearance of modular evolved boundary operators can be understood due to the equality between bulk and boundary modular flows and explicit formulas for bulk operators can be found with a complete understanding of the action of bulk modular flow, a difficult but in principle solvable task. We also obtain an expression when the bulk operator is located on the Ryu-Takayanagi surface which only depends on the bulk to boundary correlator and does not require the explicit use of bulk modular flow. This expression generalizes the geodesic operator/OPE block dictionary to general states and boundary regions.


Introduction
Bulk locality, operators commuting at spacelike separation, is a striking feature of holographic CFTs. In the N → ∞ limit, low energy fields satisfy the free wave equation in a fixed gravitational background. The extrapolate dictionary [1] identifies bulk fields close to the boundary with local CFT operators: scalar fields Φ with a given mass are dual to single traces O with a fixed conformal weight, ∆. In [1,2,3], it was observed that the extrapolate dictionary together with the free wave equation determines the local bulk fields in terms of boundary operators: where the boundary integration region consists of all points that are space like separated from Φ(X). At leading order in 1/N bulk locality is then simply a consequence of large-N factorization and the specific two point function of O determined in this background. 1/N corrections can be systematically included by solving the wave equation including interactions. True holography enters when we realize that O(x) can be written in terms of Heisenberg operators on a constant boundary time slice Σ by explicitly evolving the operators with the boundary Hamiltonian. Then Φ(X) encodes bulk locality in a highly non-trivial, yet necessarily approximate, way.
On the other hand the emergence of bulk locality is somewhat obscured by the fact (1.1) is a completely non-local mapping: a bulk field at some given point depends on all boundary operators which are spacelike separated from X. This is particularly disturbing when trying to relate the notions of bulk and boundary locality, where for example (1.1) does not smoothly reduce to the extrapolate dictionary as X limits close to the boundary. Similarly, subregion subregion duality [4,5], an attempt to relate the information contained in boundary sub-regions to that of a dual bulk sub-region, is not usefully constrained by (1.1).
In [6], a generalization of (1.1) has been postulated which is better suited for the purpose of having a more local mapping between the bulk and the boundary. They proposed that bulk fields in a certain bulk subregion r can be written purely in terms of (Heisenberg) CFT operators in a given boundary subregion R where we will take r, R to be spacelike regions of a bulk Cauchy slice Σ such that R is the intersection of r with the boundary Cauchy slice, ∂Σ. A natural mapping of subregions r(R) is suggested by the holographic entanglement entropy formula [7,8,9]: A ext (∂r) 4G N + S bulk (r) (1.2) which tells us that the entanglement entropy of the boundary subregion R is given by the area of the extremal Ryu-Takayanagi (RT) surface ∂r anchored to ∂R at the boundary 1 , plus the entanglement entropy of bulk QFT reduced to r. More generally it is expected that bulk fields in the entanglement wedge D(r), the domain of dependence of r, can be reconstructed in this way. A piece of evidence in favor of the entanglement wedge reconstruction proposal is based on the following special case. For the vacuum state in a CFT and R a ball shaped region, one can go to AdS-Rindler coordinates which cover D(r) and use Rindler mode functions to write the bulk operator in r in terms of the R operators Φ(X ∈ D(r)) = where D(R) is the boundary causal domain of R, which in this case is simply a double light cone. Again, with Hamiltonian evolution, we can write this as a non-local operator acting at R. This is the simplest expression that one could have hoped for, however it is easy to see that such a simple expression can't be correct if one considers more general regions and states. This can be traced to the fact that more generally the entanglement wedge contains a spacetime subregion which is entirely space-like separated from D(R) (the so called causal shadow region of [5] ). This means that bulk operators in that region would commute with all the local operators in D(R). Then a reconstruction formula analogous to (1.3) would imply they trivially commute with one another which is inconsistent.
In [10], it was proposed that the simplest expression that would take the previous complication into account should read: As we will explain later, this conjugation by the density matrix is a natural operation in the field theory called modular flow. The operators O s (x R ) are non-local and can't all commute with Φ(X r ). If R is a sphere in the vacuum, this expression reduces to the Rindler expression (1.3).
The modular hamiltonian is defined as the logarithm of the density matrix. In theories with a holographic dual, one can think of the modular hamiltonian as an operator in bulk perturbation theory and it is given by [11,10]: This expression implies that the commutator of bulk operators (which are spacelike separated to ∂r) with the boundary modular hamiltonian is equal to the commutator with the bulk modular hamiltonian. So, an operator in r commutes with the modular hamiltonian ofR. This property lead to the conjecture (1.4) and was argued in [12] to be equivalent to quantum correctability. In this paper, we will derive the expression (1.4) and write a formula for the smearing function which only depends on boundary information and is in principle computable. In the limit where the bulk operator lives in ∂r, this formula simplifies significantly, allowing us to compute complicated boundary quantities (which depend on the modular hamiltonian) in terms of simple bulk calculations. The formula we would like to advertise relates the operators on the RT surface ∂r to the modular average of the boundary operator: The outline of the paper is as follows. We start in section 2 by introducing various properties of modular hamiltonians and modular flows. In section 3, we derive and explore the expression (1.4). Section 4 deals with entanglement wedge reconstruction when the bulk operator sits at the RT surface, which is much simpler. In section 5, we comment on state dependence and the inclusion of interactions . We conclude with some closing thoughts in section 6 .
2 Modular evolution, fourier space and free fields Given our state of interest, ρ, and its associated algebra of operator A, the modular hamiltonian 2πK = − log ρ generates an automorphism: it sends operators in the algebra to operators in the algebra. This is called modular flow [13] 2 : This modular flow is interesting for various reasons. Formally, understanding its properties was important to develop Tomita-Takesaki theory, which made type III algebras tractable. In some very symmetric cases, such as half-space in the vacuum (or a sphere in the vacuum of a CFT), modular flow is generated by (local) Hamiltonian evolution [14,15,16] : Given that in general this flow is non-local it will be useful to illustrate some of its properties by going to modular fourier space, a basis where the action of the modular hamiltonian is simple: From now on, we are going to focus on a pure state and we will consider the density matrices of subregions in that state. If we consider a subregion R, purity implies that (K R − KR)|Ψ = 0. From this, it follows that where A, B are two operators in the algebra of bounded operators in R.
Furthermore, these operators satisfy the KMS conditions (β = 2π) with respect to the modular flow since correlators are analytic in the strip 0 < Im(s) < 2π: The main observation is that the KMS condition allows us to express the commutator 2 Note that our 2πs = −s Haag in order to identify s with Rindler time in the local case.
at finite frequency in terms of the correlator. This is an important property: the commutator of any operator and any other operator at finite modular frequency will be zero iff they are not correlated. In the following, we are going to consider local operators, O(x) in region R. We will denote their modular frequency two point function G w : By definition, correlators in modular time are time translation invariant. If we had a translation invariant state, we could go to momentum space and canonically normalize the modes. In general, G w (x, y) is not translation invariant and one can't simplify the modes further. As a side note, the generator of the full modular flow U (t) = e −i(K R −KR)t is well defined in the continuum limit: it acts like the usual modular flow when acting in an operator in R orR but it is also an good operator, which doesn't depend on details near the entangling surface ∂R. In this way, sometimes is useful to think of the previous correlators in terms of: which is equivalent to the previous discussion when A, B are in the same subregion, but it is well defined more generally.

Gaussian states in free field theory
Gaussian states in free theories have been explored extensively in the literature (see for example [17] and references therein). Since all correlators are fixed by the two point function, the density matrix will be gaussian and the modular hamiltonian bilinear. We will be focusing on scalar fields, the operators in the algebra of the region R 3 are {Φ(X), Π(Y ); ∀x ∈ R} and the respective modular hamiltonian will be: where n µ is the unit term to the Cauchy slice containing R.
In general, K is unknown unless the system is specially symmetric: one can write it formally in terms of matrix inversions and log's of correlators, but these are hard to obtain in the continuum limit (and only understandable in a distributional sense) although one can obtain it numerically using a lattice. Free fields are interesting in the holographic setting because to leading order in 1/N , the bulk quantum theory consists of free fields. Subleading orders introduce weak interactions that can be included perturbatively in the previous discussion: for example if we have a Φ 3 interaction, there will be a contribution to the modular hamiltonian which will be trilinear on the fields.
The modular fourier modes will then be a linear combination of Φ, Π: where the coefficients are determined by the canonical commutation relations: These modular frequency modes can also be used as a basis of operators. This is how the entanglement entropy is computed in [17]. The idea is that instead of labeling the modes by (d − 1) coordinates x, we can label them by w and (d − 2) coordinates X S , which correspond to a codimension 2 surface S ⊂ R. If we think of the c w (X S , Y ) as a matrix C Y ;Y ∝ G w (X S , Y ) with index Y = (w, X S ), then different Y, Y have different values and C Y ;Y is always different from zero (we don't expect the correlators to vanish). Furthermore, no linear combination of Φ,Π can give zero, so this matrix is invertible. The invertibility of this matrix guarantees that Φ w (X S ) is a good basis, ie there is a Bogoliubov transformation that sends This basis is useful because the modular hamiltonian is simpler and the entropy can be computed easily.
In terms of the modes Φ w (X S ), the modular hamiltonian will read: This hamiltonian is slightly simpler compared with (2.13): it has one less integral. The commutation relations determine K in terms of G w : where we should think of K, the modular hamiltonian kernel as a distribution, which can be thought as the matrix inverse of G w . Of course, this discussion seems very formal because in order to define the modular frequency modes, one has to know the modular hamiltonian in the first place. However, as a pratical tool, it seems like it might be simpler to compute G w (X, Y ) than evaluating the modular hamiltonian (2.13) in terms of the Φ, Π directly. The reason is that this two point function has a very particular analytic structure: it picks up a phase when going around the entangling surface ∂R . Solving the wave equation with these branch cut boundary condition seems like an easier way to handle the problem and this is how [18] dealt with the entanglement entropy massive scalar field and a spherical region. This approach to the entanglement entropy of free fields seems to be practically useful also for less generic regions [19].

A smearing function in the entanglement wedge
The extrapolate dictionary relates boundary and bulk operators in the boundary of AdS: In the future, we will denote this as Φ(x, z = 0) = O(x). Given that the bulk field is a free field, one can use the equations of motion to move this field to the bulk, at the expense of including more boundary operators: One normally presents f (X|y) as the Green's function for the scalar field with support at spacelike separation. However, f (X|x) can be understood by formally inverting the bulk-to-boundary correlator 4 : which should be understood as a distribution. For example, if the system under consideration has spatial translation symmetry and time translation symmetry (as in the case of the vacuum, finite temperature or Rindler), these kernels are quite simple in fourier space [20]. We can fourier transform all bulk directions but one, which we will denote generally as z. Given that the two point function is diagonal in momentum space,K w,k = G −1 w,k . Φ(z, x) is normally expanded in terms of canonically normalized annihilation and creation modes Φ w,k (z) = ψ w,k (z)a w,k , the extrapolate dictionary relates a w,k = O w,k G −1/2 w,k , so f w,k (z) = ψ w,k (z)G −1/2 w,k . This summarizes the explicitly known cases of HKLL and, as emphasized in [20,22], makes it clear that the kernel f (X|x) is defined only as a distribution.

Entanglement wedge reconstruction
In [10], the modular hamiltonian of holographic theories was studied to leading order in 1/N . They concluded that the commutator of a bulk field in the entanglement wedge with the modular hamiltonian was given by the commutator with the bulk modular hamiltonian: That is, in holographic theories, modular flow is bulk modular flow 5 : This equation is suggestive for subregion reconstruction and [10] conjectured that there is a smearing function in terms of these modular evolved modes: In this section we are going to derive and explore this formula. To do this, we will modular evolve the extrapolate dictionary Φ(x, z = 0) = O(x). If we apply (3.22) to the extrapolate dictionary, we get: This new version of the extrapolate dictionary seems quite shocking at first: it maps a non-local boundary operator to a smeared bulk operator with support in all the entanglement wedge. When the modular hamiltonian is local 6 , this seems standard since we can go to the Schrodinger picture, but, in the Heisenberg picture, this should seem surprising at first. Only if the modular hamiltonian is local can O s (x) be understood as a local operator at some other point in the causal domain of R. This will be the key property behind entanglement wedge reconstruction. It is important to note that, while the modular flow of a single trace will be in general a very complicated operator in the region, (3.24) implies that they have a simple expression in terms of single traces at all boundary times, because one can just do global HKLL on the modular flow of the bulk operator. 5 At this point one could complain about the fact that [10] only proved this relation for [K, Φ], but not the exponentiation. However, it is easy to argue that [K, Φ] is a low energy operator, and it keeps us within the "code subspace" of low energy bulk perturbation theory. To make this argument one should think of computing these commutators by analytically continuing the (derivatives of the) Renyi entropies: it is clear that a low energy operator in the original geometry will be low energy in the replicated geometry. We thank Don Marolf for asking this question. 6 Because of the extrapolate dictionary, the bulk modular hamiltonian is local if and only if the boundary modular hamiltonian is local.

A change of basis for free fields
We can now apply the free field theory discussion of the previous section when the subregion is R = r, a space-like section of the bulk entanglement wedge. As we discussed, the expression In order to use the extrapolate dictionary, we should consider S to be the intersection of r with the boundary that is we want to set S to be the boundary region R. The extrapolate dictionary then implies that: As we argued before, we can invert this matrix and its inverse gives the smearing function for entanglement wedge reconstruction. The key point is that, in the bulk, we have been able to write any bulk operator in the entanglement wedge in a basis where the modular extrapolate dictionary acts naturally. The fact that entanglement wedge reconstruction could be implemented by a Bogoliubov transformation was mentioned in [23], however [23] didn't identify what set of modes in the entanglement wedge has a localized boundary representation.

Symmetry based argument
Alternatively, given that the bulk and boundary modular flow are the same, one could look directly for a boundary representation of the bulk operator that satisfies the following requirements: The latter requirement implies that this operator is in R and the other two conditions force this operator to have fixed modular frequency and have a fixed conformal dimension ∆. The only single trace operators that satisfy these requirements are the operators of region R with fixed modular frequency: O w (x R ) 7 . This way of finding a boundary representation of the bulk operator is reminiscent of [24].
To leading order, the bulk field should be expandable in single trace with fixed modular frequency. This symmetry based approach seems more useful to justify the existence of a smearing function to higher order in matter interactions (or changes in the background metric), since this would require higher traces coming from bulk interactions to appear with fixed modular frequency.

The smearing function
The smearing function is determined by requiring it reproduces the correct bulk-toboundary correlator: Note that f w (X|x) satisfies the free wave equation in the bulk variable and K w (x, y) can be thought of as being defined implicitly in terms of the boundary correlators (3.28) or alternatively as the bulk modular hamiltonian of (2.15) 8 . All equalities are to be understood as expansions around a fixed state in the G N → 0 limit. We expect that one can introduce corrections due to interacting matter by just demanding this expression to give the right bulk-boundary-boundary three point function with f satisfying the corrected wave equation in the presence of interactions. The symmetry based approach of the previous section might be more convenient for this.
One often thinks of the HKLL kernel f w (X|x) as the solution to the wave equation with certain boundary conditions in a fixed background. Our approach is a combination of a change of basis in the bulk and the modular extrapolate dictionary. However, note that we could have also obtained f w (X|x) following the usual approach: requiring it that it has δ function support when X → x, that it satisfies the wave equation and that it adquires a e −2πw monodromy when going around RT surface (this condition appears naturally when considering the wave equation that Φ w (X)Φ(Y ) satisfies).
Note that f w (X|x) is an expression purely in the boundary: can be computed using global HKLL and the knowledge of the boundary modular flow. Alternatively, one can compute this bulkto-boundary modular correlator in terms of the bulk modular flow which might be easier to implement explicitly.
To make contact with (3.23), we can also write the explicit f in modular time: with K s the fourier transform of K w . As was explained in [22,20], in the local case we don't expect f to be a proper function, but a distribution, which when evaluated within correlators gives a definite answer, this is clear from our approach since the bulk modular hamiltonian, K w , is a distribution.

Different boundary representations for the same bulk operator
As we have seen, it seems like there are different equivalent representations for the bulk fields: they can be written in terms of boundary fields smeared over different subregions. In [6], it was suggested that this is realization of quantum error correction. In more simple terms, these two boundary representations for the bulk operator have to differ by an operator that is zero within low energy correlators. As was explained in [25,26], in the vacuum, these operators correspond to single trace operators with space-like momenta.
More generally, we expect that the reason why there are different representations for bulk operators has to do with the fact that [K R − K r,bulk , Φ(X r )] = 0 when inserted in low energy correlators. It is this approximate equality between the action of bulk and boundary modular hamiltonians that allows us to localize the field: from the global HKLL expression it is not clear whether [KR, Φ(X r )] = 0, but since this is equal to the commutator with the bulk modular hamiltonian this is trivially zero.

The effective generalized free fields modular hamiltonian
From the boundary point of view, the theory satisfies large N factorization. To leading order in 1/N we can think of the single trace operators as generalized free fields (GFF). Generalized free fields are defined as fields whose commutator is a c-number. By themselves they are not very nice quantum field theories, see for example [27], but they often arise in an approximate manner. GFF look like free fields but they don't satisfy any equation of motion, which means that operators at different times are independent. In holographic theories, O s (x) is a single trace operator (because of the extrapolate dictionary) and independent of O(x), so it seems that if one wanted to define the algebra of generalized free fields constrained in a subregion 9 , one should include all operators of the form O s (x).
One could construct the effective modular hamiltonian of the subregion, and because to leading order in the large N , we just have a gaussian theory, K GF F,R will be bilinear. That is, if O(x) is a generalized free field, K GF F is the bilinear operator that effectively (within two point functions) generates the same modular evolution as the modular hamiltonian: . This condition is enough to get a explicit expression for K GF F : Given the previous section, we don't expect this to be an independent object and, of course, K w (x, y) is just the boundary limit of the bulk modular hamiltonian (2.15). One can understand the 1/N perturbation theory by adding higher order terms, as one does in the bulk. From this expression, one could equivalently derive the smearing function by demanding the right commutation relations with the modular hamiltonian , which expresses Φ w (X) in terms of O w because the commutator of two single trace operators is just a c-number to leading order.

The zero mode
In the previous section, we demonstrated how one should think about entanglement wedge reconstruction if one has access to the boundary (or bulk) modular hamiltonian. The resulting formula were however not very explicit and involved formal inversions of integral kernels. It would be nice if one can get some more explicit expressions to work with. In this section we show how to do this in the case of the modular zero modes.
Consider the bulk operators Φ w (X), for small modular frequency w ∼ 0. These should approximately commute with the modular hamiltonian. From the bulk point of view, this operator has to be roughly localized near the RT surface. In the next subsection, we will show that the previous expression (3.26) simplifies dramatically in the w → 0 limit, to give: When the modular hamiltonian is local, such as for a spherical region in a CFT, this expression generalizes the OPE block story of [28,29,30]. In those papers, instead of having the zero mode of a local operator on the left hand side, one smears the zero mode over the boundary region R. Because a sphere in the vacuum is conformally equivalent to a the hyperboloid at temperature 1/2π, they consider uniformly integrating the field over the hyperboloid, which is a particularly symmetric operator. This allows one to use symmetry arguments to show that the result is the bulk operator integrated over the RT surface. Our expression is consistent with this -one can explicitly check, in the hyperboloid conformal frame, that integrating (4.32) over the boundary region R gives exactly the geodesic (minimal surface) bulk operator on the right hand side. When the modular Hamiltonian is non-local (4.32) is a highly non-trivial and surprising generalization of this OPE block story. Note that due to the reduced symmetry in this case there is no natural way to integrate (4.32) over the boundary region R. At this point we will address a subtlety that might have bothered some readers. One might worry that we don't really want the operator to sit exactly on the RT surface -the operators in the algebra have to be properly smeared into the region D(r). In particular we don't expect the bulk scalar field theory to have a center -so the notion of a well defined operator in the algebra commuting with the modular Hamiltonian K R and KR will necessarily be approximate. Or perhaps a more constructive way to think about this is that both sides of (4.32) should really be thought of as operator valued distributions, just as is the case for local quantum fields. So the zero mode may diverge when inserted into correlation functions. This divergence could arise either because of the ∞ −∞ ds integral defining the zero mode, or equivalently because the right hand side of (4.32) is not property smeared over a co-dimension 0 region in the bulk spacetime. Of course (4.32) will also often be finite when inserted in correlation functions, so it is perfectly reasonable to talk about this operator. We will see some of these delicate issues in the derivation below.
If the reader is not interested in the derivation, she can jump ahead to the next subsection, where we discuss the implications of (4.32).

Derivation of the zero mode formula
The derivation we present here works from a slightly different starting point compared to our discussion in Section 3. In particular we have found several ways to derive this zero mode formula, we present the method here that we find most instructive. For a derivation of (4.32) which starts with (3.28) see Appendix B.2.
It is convenient to consider a more general formula to (4.32) where we take any local operator in the entanglement wedge and compute it's zero mode using the bulk modular Hamiltonian Φ 0 (X). This operator also commutes with the modular Hamiltonian so one might expect it can be reconstructed on the RT surface also.
Firstly note that since the bulk is a free theory, for Gaussian states, the resulting modular zero mode will necessarily be an integral of Φ and Π on a fixed Cauchy slice: where Π = n µ ∂ µ Φ, n µ is the unit normal to Σ and h is the induced metric on Σ. Note that we can compute the unknown functions f Π,Φ by commuting both sides with the field operators and using the canonical commutation relations. We find: where we have used the fact that the bulk theory is quadratic to replace the c-number commutator with the expectation value. Note that this expression for f Π is independent of the Cauchy slice Σ on which it is evaluated. At this point we simply have to constrain f Π . Since the zero mode that appears in these commutators has very special properties we will easily be able to fix f Π . For example we can apply the KMS condition: where the later contour C lies along the lines Ims = 0, 2π at the boundaries of the complex s strip. Since via very general considerations the modular flow correlator in (4.36) should be analytic in the strip we can complete the contour, assuming the large |s| behavior is sufficiently decaying, and (4.36) will yield f Π = 0. This will be correct for Y any bulk point contained in the entanglement wedge D(r). For points in D(r), the original commutator vanishes (4.35) because the two operators are contained in commuting algebras. However there is a question of exaclty what happens at the RT surface ∂r. In particular our Cauchy slice will be chosen to pass through the RT surface. Since this is where the zero mode is expected to live we must analyze f Π carefully here.
In order to smooth out the behaviour of f Π (Y ) near the RT surface, let us approach ∂r as follows. We will pick the Cauchy slice Σ = Σ H to lie along the future part of the "horizon" ∂D(r) close to the RT surface where it will be a null light sheet. We also pick Σ H to continue smoothly onto a segment of the light sheet contained in ∂D(r) as pictured in Figure 1. Further away from the RT surface we then take Σ H to smoothly match onto a spacelike surface. Figure 1: (left) We reconstruct the zero mode using the "Cauchy slice" Σ H with a small null segment close to the RT surface. (right) We should really take a limit of space-like slices that approaches Σ H .
Since Σ H has a null segment we imagine approaching it as a limit of a sequence of spacelike Cauchy slices Σ → Σ H . This however brings about some additional complications. We now have to worry about an order of limits issue -consider a regulated version of the zero mode: where we will eventually send T → ∞. The associated function f Π can be computed for this regulated zero mode as above. We want to first push Σ to Σ H and then take the zero modes, so the appropriate order of limits we should take is: and equivalently for n µ ∂ µ f Π (Y H ∈ Σ H ). If we had taken the other order of limits we would have again found f Π = 0 away from the RT surface, and although this would still likely give rise to the correct answer for Φ 0 , no value is gained compared to simply taking Σ space-like.
Keeping the order of limits (4.38) in mind we now argue for a different behavior for f Π for points along the null segments of Σ H . Firstly note that f Π still vanishes on Σ H ∩∂D(r) since the operators within D(r) always commute with Φ s (X) which means the first limit in (4.38) will yield zero even for finite T . Now let us address the points on Σ H ∩∂D(r). Since we are free to act with modular flow in (4.38) to the left onto the operator Φ(Y ) one might guess that sending s → T → −∞ will now yield an operator that gets "stuck" on the RT surface and so the correlator will not decay for large −s. Thus we are no longer able to complete the s contour as we did when Y was a bulk point in D(r). Let us assume this is true, and we will confirm this a posteriori. Now the contour argument we used in (4.36) will give a non-zero contribution from the vertical piece of the contour: Note that the large T → +∞ limit will still give a correlator that vanishes sufficiently fast so we can complete the contour there. For example as we take Y → Y RT the modular flow of this operator will be independent of T and θ so we will have: We will now show that f Π (Y ) is constant along null generators of the future lightsheet in ∂D(r) and so it will take the value given by the limit (4.40) everywhere on the future part of the null generator 10 . We will only need to do this analysis in the neighborhood of the RT surface where we pick local Rindler coordinates: Consider the surface Σ defined as u = 0 for vr < v < v 1 where vr < 0 lies on the appropriate part of ∂D(r) and u = u (v) < 0 for v > v 1 . That is Σ H leaves the null sheet beyond v > v 1 . We evaluate Φ 0 (X) here with what we so far know about f Π : We have split the contribution into two parts -the first from the null segment and the second from close to where Σ H leaves the null segment. For the later we have dropped f Π n µ ∂ µ Φ since this will be smooth and thus vanish as we send → 0. All fields are evaluated at (u = u (v), v) where we take u (v) = 0 for v < v 1 . Note that the appropriate (rescaled) normal vector is: We can evaluate the final term in (4.42) by noting that: where the first term will be smooth as a function of v so we can drop it and the later term will have a delta function since as we leave the null sheet f Π goes from non-zero to zero 11 . Putting this together we have: where we have integrated by parts. The boundary term vanishes at v = vr since f Π vanishes there and at v = v 1 − it cancels the delta function term in (4.44). Let us now imagine picking a slightly different surface Σ H,2 which leaves the null segment at a slightly later null coordinate v = v 2 . The answer should be independent of 10 This suggests that f Π is a highly singular function especially when any of these null geodesics leaves D(r) along seams that emanate from a caustic (see [5] for some nice pictures of how this happens in a relevant context). Since in this case two different null generators, with different values of f Π intersect. This is to be expected since modular flow will have some singular behavior for operators close to these seams, however we can avoid these by picking Σ to leave the lightsheet before we ever encounter any caustics. Note that f Π is already a distribution since it vanishes everywhere inside D(r) and can be anyway regulated by picking T large but finite. 11 This delta function will be naturally resolved using the T regulator with a width = O(e −T ).
this choice, so we find that: and since this is an operator statement the only way we can satisfy this is if ∂ v f Π = 0 for v > 0. Note that we can pick the null coordinates differently on each generator v 1,2 (y) .
We have thus shown that the equations of motion necessitate that f Π is constant along null generators of ∂D(r). We made the argument near the RT surface using local Rindler coordinates, but that was just for simplicity and this argument should be applicable away from this coordinate patch (although our zero mode derivation does not rely on this fact.) Across the horizon now f Π jumps to zero so (4.45) gives: (4.47) where we have used (4.40). This was the claimed result. Since this derivation might seem delicate we give another derivation, making use of a different limit towards the horizon, in Appendix B.2 and find the same answer. This new derivation relies on some reasonable assumptions about modular flow for free QFTs close to the entanglement cut/RT surface.

Bulk reconstruction from zero modes
In the context of bulk reconstruction, we can invert the previous formula to write This formula gives a simple formula for reconstructing bulk fields in the entanglement wedge if we can foliate r by RT surfaces of smaller regions. This requires inverting the bulk-to-boundary correlator over the RT surface, which is easier than having to solve the full modular flow.

More on smearing and conserved charges
In quantum field theories, one normally works with smeared operators: Φ f = d d xf (x)Φ(x), which have finite two point functions. These operators are normally smeared in codimension 0 surfaces, but one could wonder if they could be defined in higher codimensional surfaces. Let's consider the case of a free scalar field, smeared around a codimension n surface (see for example [31]): So one can't smear a scalar field in a codimension 2 surface. This means that the operator localized in the RT is not a well defined operator and it is not contained in the algebra of operators of the entanglement wedge.
For this reason, whenever we discuss these zero modes, we should think of smearing it inside the entanglement wedge. In this way, this becomes an operator of the respective region R (orR) and not an operator which is seemingly in the center. Equivalently, since we are working to leading order in 1/N , we can't resolve frequencies which are O(1/N ) which means that we should put some O(N ) cutoff in the boundary integral over modular flow.
To some extent, this is a subtlety, but it is important because we don't expect to have a center for this simple example with bulk scalar fields.
An exception to this analysis would be that of conserved charges. Conserved charges, such as the electric field in the RT surface or the area operator can be localized in a codimension two surface, because Gauss' law allows it. From the algebraic point of view, this is no mystery, since we expect these algebras of operators to have a center, see [32].
Understanding conserved charges is really interesting but we are not going to explore these in detail. We expect that their boundary smearing function should simplify and it is likely that the vector, ξ, with which they are contracted is also simple. For the case of gravity, we expect that graviton contribution to the linearized area operator should correspond to the zero mode of the stress tensor: In the electromagnetic case, we expect that: The electric flux along the RT surface was studied before in [33]. For Rindler, these conserved charges were studied in [29] and since the modular hamiltonian is local there, the s integral drops out. Note that we only expect these integrals of conserved charges to be well defined without further smearing when Gauss' law applies. We don't expect this to be true for general bulk smearings over the RT surface.

Dressing
In all our discussion, we have ignored the contribution from the area to K bdy ∼ K bulk . As explained in the previous section, the zero modes shouldn't be exactly localized in the RT surface, but even if they were, they would commute with the area operator, since the area operator just generates a time translation. The only operators with which the area operator does not commute are operators in its past/future 12 . In gravitational theories, operators can't be exactly local, they have to be defined in a gauge invariant way. This means that they can't commute with the stress tensor everywhere in the boundary because Gauss's law relates the energy in the bulk with the energy flux at infinity. One normally thinks of this in terms of gravitational dressing: the bulk operator has its own gravitational field which gives a non-trivial commutator in the boundary, see [34] for more details.
From the bulk point of view, there are many ways that a bulk field could be gravitationally dressed. In analogy with the wilson line in electromagnetism, a geodesically dressed case is often considered, where the gravitational flux is localized at a point in the boundary [35,36,37,38], but more non-local dressings seem also natural from the boundary perspective [39] .
Our zero modes (slightly smeared into R) commute with all the operators inR , (approximately) with the modular hamiltonian K R and within correlators they commute with all operators in the region: [A 0 , B] ≈ 0, for an arbitrary operator B ∈ R. If the modular hamiltonian is local, [H R , A 0 ] = [K R , A 0 ] ∼ 0 would mean that Gauss' law is properly satisfied and that the gravitational field is localized around ∂R. More generally, we expect that, given that the operator approximately commutes with the modular hamiltonian, this is enough to guarantee that this operator should be localized in ∂R. If we had a nice quantum field theory, we certainly expect that these zero modes are basically operators localized in ∂R, but because of large N , these could be non-trivial operators.
While we leave a more detailed analysis for the future, we expect that these zero modes are dressed to the boundary codimension 2 region, ∂R. It is unclear to what extent a zero commutator with the modular hamiltonian is enough to conclude that Gauss' law is satisfied but our operators satisfy [Φ 0 , T µν ] = 0. We expect that, after taking care of the gravitational interactions, [Φ 0 , T µν (x R )] = 0 inside correlators.

Zero modes and G 0 (x, y)
Equation (4.32) also gives us a simple bulk method to compute G 0 (x, y): This quantity is very non-trivial to compute in the field theory, but in gravity it can be computed explicitly given the RT surface and the bulk-to-boundary correlator in a particular state. For example, if the operators have a large external dimension where L(X, y) is the geodesic distance between a bulk and a boundary point and Y RT,min is the point in the RT surface which minimizes the previous expression.
In appendix A, we prove this formula explicitly in the local case, using the methods of [40]. In the next subsection, we compute G 0 (x, y) for some simple holographic examples and compare it with the field theory expectations.

Two interval example
A simple, illustrating example is that of two intervals in two dimensions in the vacuum. Consider the t = 0 spacelike slice and put the two intervals at (−1, −y) ∪ (y, 1). G 2−int 0 (x 1 , x 2 ) will just be the sum over two one interval modular correlators G In this way, the two intervals zero mode correlators are given by: , if y > y 0 (4.55) In Figure 2, we plot G 2−int (x 1 , x 2 ). We can compare it with the modular flow of the large interval (−y, y) and the modular flow of the single interval (y, 1).

Comparison with 2d free fermions
The only analytic example of non-local modular flow was studied in [41,42] for 2d free fermions and n intervals. Given that this is a usual quantum field theory, we expect these zero modes to be roughly localized in ∂R.
In this case, one can compute the modular flow for two intervals, at (−1, −y), (y, 1) and the modular flow will be non-local but it will only mix two trajectories, one in each interval. Furthermore, it doesn't mix chiralities. The modular evolution of the chiral fermionic field operators undergoes mixing: Ψ s (x 1 ) = cos θ(x 1 , s)Ψ(x 1 (s)) + sin θ(x 1 , s)Ψ(x 2 (s)) (4.56) with x 1 being the initial location of the operator, inside one of the intervals (of course, Ψ(x 1 , s = 0) = Ψ(x 1 )). The explicit expression for θ(x 1 , s) can be found in [41,42] and x 1 (s), x 2 (s) are the corresponding local modular trajectories , which are given explicitly ( 1−y 2 , y + z(1 − y)) as a function of z ∈ (0, 1) for ∆ = 1 and y = 0.1 (which corresponds to a cross ratio of x = 0.65). The large interval flow corresponds G (−y,y) 0 (x 1 , x 2 ) and the one interval flow to that of G (y,1) 0 (x 1 , x 2 ). We see that when x 2 is close to the endpoints of the interval, the flow is roughly that of the smaller interval, as we expect from locality. For z > 1/2, the modular flow is well approximated by that of the large interval.
by the two solutions to: In [42] they compute the modular evolved correlator explicitly to make sure that it satifies the KMS condition The zero mode correlator will then be After integrating over modular flow the pole at x = z disappears and one instead gets two branch cuts when X(x)Z(z) = 0, that is, when x or z are light-like separated from the left boundaries of the intervals. This is what one might have expected for chiral fermions.

Interactions, background dependence and backreaction
It is interesting to understand how these formulas for entanglement wedge reconstruction, which depend on the modular flow of the state under consideration, are compatible with the usual background dependent description of HKLL. Similarly it is interesting to understand how to include interactions. Of course, our expressions are valid for any state as long as we don't need to shift the entanglement wedge (and thus the algebra of operators stays the same). However if the background fields change or one considers higher orders in bulk perturbation theory, this will be no longer be true.

Interactions
In the presence of 1/N corrections, the HKLL dictionary between bulk and boundary fields is not diagonal anymore, let us illustrate this with a λΦ 3 bulk interaction. In the vacuum, these corrections to HKLL read [3]: where the new term is a sum over double trace operators. Of course, the 1/N expansion encodes the expression of the bulk fields in other states which are related by a shift in the one point functions. It contains more information since it reproduces higher order correlators, which can't be understood as two point functions in a different state. From the argument of the previous sections, we expect that our modular fourier mode expression has the same expansion: That is, bulk operators of fixed modular frequency can be understood in terms of a single trace operator plus double trace operators of the same frequency. Gravitational or gauge interactions can be taken into account by solving the wave equation perturbatively and will add corrections for the boundary field which include the boundary stress tensor and currents.
In the context of zero modes, one can imagine accounting for the interactions in two different ways. One could stick to the zero mode of a single trace operator O 0 and thus the bulk field integrated in the RT surface would get 1/N corrections, or keep fixed the bulk field Φ integrated in the RT surface and study what boundary zero mode it corresponds to. While the first approach was explored in the Rindler case in [29], we expect the second approach to be simpler, since it is aligned with the usual ideas of HKLL, but we reserve a more careful analysis to the future.

Backreaction
The discussion of the previous sections referred to matter interactions or changes in the background fields, but didn't include backreaction. In this work, we have been using heavily the fact that the bulk and boundary modular flow are equivalent. However, if we want to account for higher G N corrections, the results of [10] will receive order G N corrections and the two modular flows won't be exactly the same [43]. Despite this, [44] have shown that when difference between bulk and boundary relative entropies is small (as in the case of G N corrections), one can still reconstruct the bulk operator. Even if the position of the surface changes by a small amount, the modular hamiltonian will change everywhere [45], but the previous arguments seems to suggest that one can still work in perturbation theory. In this way, we expect then that [44] together with [43] allow our story to be generalizable in the presence of backreaction.

Background dependence
As discussed before, in order to define bulk operators in a gauge invariant way, one can do it by throwing geodesics from the boundary. In this way, one can define the operator Φ(x, z) by throwing a spacelike geodesic from x in the boundary up to some renormalized proper distance z. As discussed in [46,47,6], one can define the same operator Φ(x, z) for the family of states which contain the same geodesic. However, since Φ(x, z) will satisfy different wave equations in different backgrounds, when writing a boundary expression for this operator, it would seem to depend on the state. This is often called "background dependence".
Of course, whether two of these state independent (within the family of states with the same geodesic) operators are spacelike separated is going to depend on the state where the commutator is evaluated. Similarly, whether an operator Φ(x R , z) is in the entanglement wedge of R is state dependent. In this way, a gauge invariant description of the entanglement wedge in terms of boundary geodesics (x R , z) will be state dependent, so we don't have any reason to expect a state independent boundary representation for a bulk field (geodesically dressed to the boundary) in the entanglement wedge.
However, it could be that there is a more state independent definition of bulk operators in the entanglement wedge which doesn't rely on throwing geodesics from the boundary. For example, one can define the state independent bulk operator integrated over the extremal surface anchored in ∂R. In any given state, this corresponds to a lin-ear combination of our bulk zero modes 13 and, while in the boundary the zero modes depend on the modular hamiltonian explicitly, it is possible that this is just "background dependence". This could be understood from the fact that we can think of background changing operator as as a bulk coherent state, so |ψ = U |ψ . Since a bulk coherent state is linear in the fields, the unitary factorizes U = U r Ur and the bulk modular hamiltonian changes by K bulk,r = U r K bulk,r U † r . However, at this moment it is not clear to us if this argument is enough to show that these operators are state independent.

Discussion
In this paper, we have studied how to reconstruct bulk operators in the entanglement wedge in terms of boundary operators. Because of the equivalence between bulk and boundary modular flow, the natural boundary operators dual to bulk operators in the entanglement wedge have fixed modular frequency. The smearing function f w (X|x) is simply a convolution of the modular hamiltonian with the correlator between the bulk operator and a boundary operator with frequency w. This is well defined, although it might be hard to compute explicitly. In the limit where the bulk operator goes to the RT surface, the expression simplifies and the smearing function is just the inverse of a bulk-to-boundary correlator.
[ADDED] As we have explained before, in [44] a different connection between modular flow and entanglement wedge reconstruction was discussed. We would like to emphasize that while eq. (7) of [44] and our expression (3.23) both have modular flow, these formulas are not related at all. The reason why there is modular flow in (3.23) in the semiclassical limit where we are working (ie G N = 0), the bulk and boundary modular flows are equivalent and this provides a convenient way to parametrize operators. In contrast, [44] use the modular hamiltonian to make the discussion of [12] stable for small but finite G N , by smearing it with a fixed kernel in modular time, if one sets G N = 0 in their formula, one morally gets O R = trRΦ(X r ).
We would like to conclude with some further observations.

Bulk emergence
It would be nice to understand to what extent one can use the expression for the zero modes (4.32) to probe the bulk locality directly from the boundary, by considering the zero modes of different regions. A key ingredient for bulk locality is the non-commutativity of [Π(X), Φ(X)], but the zero modes seem to be naturally probing Φ rather than Π. In other words, the equal time commutator of zero modes of different spatial regions (lying on the same spacelike slice in the bulk) is zero. Of course, ∂ t O 0 would naturally give us the Π operator, but this doesn't seem a natural operator in R. One could think that in principle, bulk locality can be studied by understanding the singularities of the correlators of zero modes (4.52) in the spirit of [48,49]. Similar ideas have been pursued recently in [50,51].
Note also that, even if modular flow is very complicated in general, (3.24) implies that, in holographic theories, it is simpler than we would have expected, since it can be effectively described in terms of a single trace operator smeared over the whole boundary.
Locality and the "algebra" of generalized free fields If we consider a given bulk t = 0 slice, the usual HKLL discussion suggests that the algebra of bulk operators in that Cauchy slice is given in terms of 1/N corrected generalized free fields at the boundary t = 0 slice combined with hamiltonian evolution 14  Even if it is clear that there should be such factorization at finite N , it doesn't seem like this factorization would have a smooth large N limit: trying to take trace overR of a single trace operator (GFF) which has support across ∂R seems rather ill defined. If one could do this, then by taking the trace over R of the HKLL global representation of an operator in the entanglement wedge of R, one should get the identity inR (plus operators which annihilate any low energy state).
What we have observed is that instead, we can think of the bulk operators at t = 0 as being generated by boundary operators at t = 0 and U K = e i(K R −KR)t : {Φ(X, t = 0)} = {O(x, t = 0), e i(K R −KR)t }, for an arbitrary choice of R. For any given R, this representation of the bulk algebra explicitly preserves subregion locality.
It is unclear to us how the algebra generated by the 1/N corrected generalized free fields and U K is (in the large N sense) isomorphic to the algebra generated by the GFF's and U H (isomophic in the sense that the respectives GNS Hilbert spaces are the same). Generalized free fields at a given time can't be expanded in terms of generalized free fields at earlier times because they don't satisfy any equation of motion (this is the reason why they violate the time slice axiom), but we don't need the whole U H for all times to get all linearly independent operators, somehow U K is enough. This is reminiscent of the interesting idea [53] that modular time might play and interesting role in quantum gravity.

Quantum error correction
In [6,12], an expression for bulk reconstruction was given which depends on the Schmidt decomposition into R,R of states created by bulk excitations (see [44] for a formula which generalizes the previous when the relative entropies are not exactly the same). It is not clear to what extent these formulas are different from writing the global HKLL representation of the bulk operator and tracing overR , to get a localized operator in R plus something that annihilates all simple states. And, as discussed in the previous few paragraphs, we don't expect this to have a smooth large N description.
Our discussion didn't use any of the insights from quantum error correction, we have just followed the usual HKLL approach with the knowledge that the bulk and boundary modular flows are equivalent. However, if we wanted to use the formulas of [6,12] in this bulk perturbative setting, their interpretation would be that they make the extrapolate dictionary more precise. That is, we would first write the bulk fields in a basis labeled by a position in the boundary AdS and then map them to boundary operators where it is now clear how to take the trace. Given that the bulk and boundary operators related by the extrapolate dictionary live in different Hilbert spaces, one has to be careful when taking the trace over the complement of the region and this is what the formulas in the appendix of [6] achieve in this setting: they map the bulk operator Φ r (x, z = 0) to the boundary operator O R (x).
Finally, note that in order to write the bulk operators in a basis where it is easy to take the boundary trace, we had to use an insight of bulk low energy subspaces: the modular flow of a simple operator is a low energy operator and thus acts within the code subspace (see footnote 5). We don't expect this to be a generic property of quantum error correcting codes. of local modular flow to fields near the entanglement cut where locally the cut looks like a Rindler cut.

B.1 Local modular hamiltonians
For certain symmetric cases, the modular hamiltonian is local: the flow is geometric with respect to a killing vector ξ(x) : That is, O(x, s) = O(x ξ (s)) and the modular hamiltonian in position space will simply be the integral of the stress tensor. This means that we can write K = R dΣ µ ξ ν T µν .

Ultralocality
Consider the half space/Rindler cut in vacuum. The entangling surface in lightcone coordinates dudv + dy 2 corresponds to u = v = 0. We can consider the horizon u = 0 as the Cauchy slice where we consider our free fields. The Rindler cut corresponds to the future horizon v > 0. Fields in the future horizon are ultralocal: the algebra of operators is a direct product of the two dimensional algebra of operators at every (y, v), that is correlators between different lightrays vanish: In this case, modular evolution just amounts to rescaling v → ve s . Given that correlators in the vacuum state are ultralocal, the modular hamiltonian acts independently on each light ray. Locality of the modular hamiltonian and ultralocality are equivalent statements: if a state has ultralocal correlators, then the modular hamiltonian necessarily acts independently in each lightray. If the modular hamiltonian is local, given that the "bilocal" kernel in modular Fourier space defined in (2.15) is the inverse of the commutator G w , the commutator has to be ultralocal which then implies the correlator is ultralocal.
The algebra of operators in each lightray (y, u) is basically that of the (chiral half of a) 2d massless scalar CFT. So for example the correlator Φ(v, y)Φ(0, y) ∼ (4π) −1 log(v)δ(y− y ) has a strong IR divergence which we would usually try to avoid. In the next section these ln divergences naturally arise for more general entangling cuts due to the approximate locality of modular flow close to the cut.

Locality of the modular flow close to the entangling surface
One would like to understand to what extent the modular flow is local close to more general entangling surfaces. For free fields, one can explore this by studying the extent to which correlators are ultralocality in this more general case. This is because for free fields (and a gaussian state satisfying Wick's theorem) knowledge of the two point correllator of a state is enough to reproduce the modular Hamiltonian.
For general surfaces , it will be convenient to write the free fields in terms of adapted lightfront coordinates. For a given codimension two surface we will denote coordinates along this surface as y. By shooting lightrays from this surface, we can define the future/past horizons H ± . We can introduce coordinates u, v which denote some fixed affine distance along these lightrays. In general, the future/past horizons will have caustics, so these coordinates will breakdown at some finite affine parameter. We would always like to consider a family of space-like Cauchy slices which limits to the horizon H + -at least for some portion of H + avoiding the caustics. So in the limit we will choose to parameterize our Cauchy slice with the null coordinate v along H + for 0 < v < v after which v becomes a space-like coordinate and the Cauchy slice moves into the interior of the domain of dependence of the entanglement region. We will denote this smoothed out Cauchy-slice as H + smooth . For u ∼ 0, v 1, we can use adapted coordinates to write the metric of this surface: The presence of this extrinsic curvature implies that for neighbourging light rays (y ∼ y ), ultralocality will be preserved as long as vK + 1. This can be seen explicitly from the correlator: This quantifies what one means by the modular hamiltonian being approximately local -it will be approximately local if we act on operators with v 1/K + . Thus we can approximate the flow via v → ve s in this case.
Note that for y far from y ultralocality will certainly break down since and the nonlocality of the modular hamiltonian will depend on more global geometric properties. For example, if we have two disconnected surfaces with K + ij = 0, it is clear that the correlators between two lightrays emanating from the different surfaces at y, y won't be ultralocal, rather ∂ v Φ(v, y)Φ(v , y ) → C(y, y ) for small v, v . However this correlator will be suppresed in the distance between the lightrays and comparing to (B.72) we expect modular flow to produce non-local mixing with distant light rays, but this effect is supressed by the correlator C(y, y ) between these two distant points and positive powers of v, v .
To sum up, for general surfaces, H + has to be smoothed out because it contains caustics. Generally the presence of caustics spoils ultralocality and thus the locality of the modular hamiltonian. However, for fields close to the entangling surface, we expect that as long as v is smaller than any of the other length scales of the system, we still expect the modular hamiltonian to be local.
While this analysis seemed to rely explicitly on the ultralocality of free fields, the lesson is more general: we expect the modular flow to act locally close to the boundary of a region, as long as we don't probe any other scale of the system. One can argue for this by considering the expectation value of operators in the presence of conical singularity of strength 2π n (which would compute the Renyi entropies), inserted along the entangling surface. In conformal field theories, if we consider a field close to the conical singularity (denote ρ the polar distance to the singularity), we expect: That is, we expect these one point functions to diverge with their anomalous dimension as they get close to the singularity, and for a general state and region, this UV conformal behaviour will be leading as long as we consider ρ L, where L is the shortest length in the system. Upon analytic continuation in n, we expect that this notion of approximate locality also applies to modular evolution. in the limit where the bulk operator goes to the RT surface, X → X RT . In this limit the commutator [Φ s (X), y] becomes a constant for a long period of modular time − log ρ < s < log ρ, where ρ is the approximate Rindler radial coordinate, the constant is 0 since the RT surface and the boundary point are space-like separated (given than the modular flow is local as long as ρ is smaller than any other scale). For most of this discussion we assume the fourier transform is well defined, so the commutator vanishes sufficiently fast for large s. The correlator similarly is constant for a long period of modular time − log ρ < s < log ρ and the fourier transform evaluates to a resolved delta function:

B.2 Zero mode from local modular flow
where the resolution scale is related to log(ρ) and in particular δ log ρ −1 (0) ≈ 2 log ρ.
In this way, in the limit where the bulk operator goes to the RT surface, (3.28) reduces to Φ(X RT ) = R dxF (X RT |x)O 0 (x) , F (X RT |x) = dyK 0 (x, y) Φ(X RT )O(y) (B.76) Now the idea is to take correlators of (B.76) with respect to Φ(Y ) for Y close to the entangling surface. Doing this carefully and using the same as above to take Y → Y RT we have: The log ρ divergence in the right hand side above could seem surprising at first. Of course, as discussed in the previous section, we expect this logarithmically divergent term from the short distance behaviour of the correlator: where ρ is a regulator which separates the two operators away from the RT surface and x RT and y RT are d − 2 flat coordinates along the RT surface in a patch around X RT and Y RT . One can work out that N −1 d = 2πS d−3 where S d−3 is the area of a d − 3 sphere. As a distribution this behaves as a delta function: log( ρ L )δ (d−2) (X RT , Y RT ) + (finite as ρ → 0) (B.79) where L is some IR scale associated to the curvature scale or the mass of the scalar and these formulas are consistent with our previous discussion around (B.72) if we set ρ = √ uv. Comparing the log ρ divergent terms we see that: Let's write this equation as 1 = 4πF · Φ where integration has become a matrix product and we assume all these matrices are invertible (both left and right inverse exist), which means that F −1 = 4πΦ. Likely these are anyway only formal manipulations since we already know that F is only defined as a distribution acting on CFT correlation functions. The reconstruction of operators on the RT surface is defined as Φ = F · O 0 so that we can write 4πΦ · Φ = O 0 , which when written out looks like: which is the same formula as we have derived in the main text.