A Cardy formula for off-diagonal three-point coefficients; or, how the geometry behind the horizon gets disentangled

In the AdS/CFT correspondence eternal black holes can be viewed as a specific entanglement between two copies of the CFT: the thermofield double. The statistical CFT Wightman function can be computed from a geodesic between the two boundaries of the Kruskal extended black hole and therefore probes the geometry behind the horizon. We construct a kernel for the AdS3/CFT2 Wightman function that is independent of the entanglement. This kernel equals the average off-diagonal matrix element squared of a primary operator. This allows us to compute the Wightman function for an arbitrary entanglement between the double copies and probe the emergent geometry between a leftand right-CFT that are not thermally entangled.


Introduction
The holographic mapping between the two-sided black hole in Anti de Sitter space (AdS) and the thermofield double (TFD) state of a strongly-coupled conformal field theory [1] is one of the most useful settings for the understanding of quantum gravity. It allows for the reformulation of questions concerning the black hole interior in terms of unitary quantum mechanics involving both copies of the boundary field theory. In this context, the black hole entropy is accounted for by the entropy associated to the entanglement between the two boundaries. This can be understood in the bulk as a consequence of the Ryu-Takayanagi formula [2]: the minimal surface homologous to one of the boundaries is located at the horizon.
In [3], the idea of using the Ryu-Takayanagi formula to understand how entanglement "stitches up" the bulk geometry was introduced. More precisely, a thought experiment was proposed: suppose the two boundaries are progressively disentangled. What would happen to the bulk geometry? The Ryu-Takayanagi formula tells us that the area of the minimal surface separating the two boundaries (i.e. the black hole horizon) would decrease. The geodesic approximation, e −∆L(x 1 ,x 2 ) ∼ O(x 1 )O(x 2 ) (1.1)

JHEP09(2018)005
and the relationship between mutual information and correlation functions tells us that the proper distance between the two boundaries would increase. These two observations lead to the conclusion that the bulk would pinch off. This general picture can be checked by examining the geometry of black holes of different masses, since decreasing the mass leads to a decrease of entanglement. A bulk calculation was also presented in [4]. The idea of probing the geometry behind the horizon using two-point functions was also extensively studied in [5][6][7].
In this paper, we make the answer to the question of what happens to the geometry when entanglement is decreased more precise. It builds on the observation that the twosided black hole is more appropriately viewed as dual to thermal Schwinger-Keldysh pathintegral [8] and the correlation function between operators inserted on opposite boundaries dual to geodesics for large conformal dimension computes the Wightman function [7]. The Wightman function vanishes for T = 0, where there is no entanglement in the TFD description, and one directly infers that at T = 0 the two Schwinger-Keldysh copies of the CFT geometrically disconnect. Clearly Feynman Green's functions dual to operators inserted on one side remain well-defined, showing that it is only the geometry inside the horizon that changes drastically while the geometry outside the horizons stays smooth as one changes T . Computing the Wightman function in the TFD for an arbitrary entanglement between the two copies that interpolates between no entanglement and thermal thus allows one to probe in some detail the emergence of the geometry behind the horizon. This is the aim of this paper.
Focusing on AdS 3 /CFT 2 , we do so and explicitly calculate boundary two-point functions between scalar primary operators inserted on the two boundaries for states which have less entanglement than the TFD. A crucial step in this computation is our construction of an expression for the average matrix element squared of a scalar primary operator O of dimension ∆ at given energies: where e S(E) is the number of states at energy E. 1 This can be seen as a generalization of the expression found for average diagonal elements found in [9]. We obtain this formula by inverting the relationship between F (E a , E b , ∆) and the two-point function of primary operators inserted on both boundaries in the thermofield double (1.3) using the fact that the holographic expression for the left-hand side is known. As we were completing this article, we were alerted that Brehm, Das, Datta also computed F (E a , E b , ∆) by different means [10]. Similar work was also subsequently done by Hikida, Kusuki, Takayanagi [11]. Moreover, an average over heavy primaries, based on large-c JHEP09(2018)005 conformal blocks, has been found in [12]. The resulting expression for F (E a , E b , ∆) reflects properties of eigenstate thermalization. For high enough ∆ in non-integrable CFTs we expect the eigenstate thermalization hypothesis (ETH) to hold [13][14][15][16]: where O(E) is the average of the diagonal matrix elements at energy E, f (E a , E b ) is a smooth function and R ab is a random number with vanishing average and mean-squared equal to unity. In particular F (E a , E b , ∆) is closely related to the variance function f (E a , E b ), as we shall show. With our construction of an expression for F (E a , E b , ∆) and the density of states, we can compute Wightman functions for states with entanglement patterns different than that of the thermofield double. The simplest such pattern is one which brutally interpolates between effective thermal entanglement below an arbitrary hard cut-off Λ and zero entanglement above this cut-off.
For this set of non-thermally, lesser entangled "doubled" states we examine the dependence of the Wightman function on time separation and cutoff. They clearly illustrate that erasing entanglement destroys a conventional notion of a smooth connected geometric spacetime behind the horizon.
2 Constructing F (E a , E b , ∆): the average of the square of matrix elements

Setup
Our starting point is the global AdS 3 black hole. This is dual to a two-dimensional CFT defined on a spatial circle of size L at finite inverse temperature β. Thermal correlation functions in this CFT can be computed from the thermofield double state, where the latter equals by definition Here the sum is to be taken over all energy eigenstates |a of a single copy of the CFT. The one-sided Hamiltonian H is such that H|a = E a |a . Z(β) is the usual thermal partition function, which ensures that the state is properly normalized. We consider time-evolution in the thermofield double with the Hamiltonian 3)

JHEP09(2018)005
so that the whole state is stationary. Therefore, for an operator inserted on the left boundary we have matching the notion that "time runs backwards" on one of the boundaries. 2 The expectation value of the TFD two-point function -dual to operators inserted on opposite boundaries -is related to the thermal Wightman function by analytic continuation. Assuming that O is Hermitian where ρ β is the usual thermal density matrix. In what follows, we will simplify the calculation by taking both operators to be inserted at the same coordinate φ. 3 Since they are on opposite boundaries (or, if analytically continued to a one-sided calculation, separated in time by at least iβ/2), we do not have to worry about issues of ordering or contact divergences.
Substituting the definition of the TFD one has We can split the sum over states into a sum over energies and a sum of states at that energy: 2 Note that we will stick with this convention for the direction of time. Operators on the left boundary should always be thought of as being time-evolved with the two-sided Hamiltonian, so that if they are at a positive value of t they are inserted in the bottom half of the Penrose diagram. This is because, with our choice of Hamiltonian, real time on either boundary matches the real part of the Schwarzschild time coordinate in the bulk in the appropriate patch. In other words, the symmetry associated with the bulk Killing vector consists in shifting the time coordinate on both boundaries by the same constant. 3 In appendix A, we generalise the results of this section to the case of operators inserted at different spatial coordinates φ1 = φ2.

JHEP09(2018)005
Replacing the weighted sum over energy levels with an integral over the density of states ρ(E a,b ), 4 and making the following change of variables: with E L = − πc 6L being the Casimir energy, we obtain an expression for the TFD two-point function in terms of the average of the square of matrix elements F (E a , E b , ∆), as defined in (1.2).
where Θ is a step function enforcing the obvious bound χ ≤ 2E. We now use the fact that the l.h.s. is known from holography in the large c limit [1]. The thermal Wightman function in the deconfined state dual to a black hole equals 5 where C O is a normalization constant and the sum over images ensures spatial periodicity φ φ + 2πn -recall that both operators are inserted at the same spatial coordinate φ = 0. Note that as L → ∞, only the n = 0 term survives. This is also the case if we take ∆ → ∞. In holography, the n = 0 images correspond to winding geodesics, which should be suppressed.
Matching Equations (2.10) and (2.11) gives us: Notice that the left-hand side is nothing more than a Fourier transform in χ followed by a Laplace transform in E. Inverting these two transforms will give us an expression for F (E, χ, ∆). We now proceed to do so.

Inverting the Fourier transform
Inverting the Fourier transform on the r.h.s. of eq. (2.12) is straightforward. For the l.h.s. we use the following integral identity This gives us Note that F n is a smooth function of β and χ and is even in χ. Furthermore, F n (β = 0, χ, ∆) = 0 and F n (β, χ = 0, ∆) ∝ β. For the special case of n = 0, we have (2.17) The identity eq. (2.12) after inverting the Fourier transform becomes:

Inverting the Laplace transform
We now use the Mellin inversion theorem to invert the Laplace transform. This inversion theorem states that, given a function H(z) analytic in C\{z 1 , . . . , z n } (the z i being isolated singularities), its inverse Laplace transform is Res(e zw H(z); z j ) , γ 0 > Re(z i ) . (2.19)

JHEP09(2018)005
The integral is taken over a vertical line in the complex z plane which is set to the right of all the singularities of H(z). We then use the residue theorem by closing the contour to the left (the so-called " Bromwich contour"). It turns out to be simplest to split the inversion by finding the inverse Laplace transform of e βE L Z(β) and of 2π β 2∆ F n (β, χ, ∆) separately and convolving them. Going from the β domain to the E domain, we have: The identity determining F (E, χ, ∆) in terms of these functions is (2.22) Restricting ourselves to integer ∆ > 1, we computeF n (E, χ, ∆) as defined in eq. (2.21) with the inverse Mellin transform. We observe that the combination 2π β 2∆ F n (β, χ, ∆) has simple poles at: and a pole of order 2∆ − 1 at β = 0. Summing over the residues at these poles, we find that Next we compute ρ(E +E L ) defined in eq. (2.20). This function is the density of states. Recall that E L = − πc 6L is the Casimir energy. By definition, Shifting the range of integration one has

JHEP09(2018)005
Thus ρ(E + E L ) indeed equals the density of states. To compute it, we need to inverse Laplace transform the partition function. In [17], it was shown that, to leading order in c, the partition function in the deconfined phase is given by Thus, using the Mellin inversion formula (for E + E L > 0), we have since the only singularity is the essential singularity at β = 0. We can write the Laurent series to evaluate the residue: The residue at β = 0 corresponds to the k = m − 1 term, hence where I 1 is a modified Bessel function. In terms of the usual energy coordinateẼ = E +E L ∈ − πc 6L , ∞ and substituting the explicit expression for the Casimir energy E L < 0, this equals This inherently c 1 expression is also only valid forẼ > 0, i.e. E > −E L = |E L | = πc 6L , since we can only close the Bromwich contour to the left when e β(E+E L ) vanishes for Re(β) = −∞. For values −E L > E > E L the Mellin transform as it stands appears illdefined, but this must be cured by all lower-order in c corrections to the partition function, which are not universal. 6 We shall only consider E |E L |. Then, keeping c fixed and takingẼL → ∞, we recover the Cardy formula: The limit L → ∞: decompactifying the CFT from the circle to the line Equations (2.25) and (2.31) give a formal solution to F (E, χ, ∆). However, the fact that this solution is expressed in terms of two infinite sums -one over images from the periodic boundary condition on the interval L, see eq. (2.22), and one over infinite poles, eq. (2.25) -is an obstacle to explicitly evaluating correlation functions in our modified state. We will therefore study the system with the size L of the system to be large. This way we JHEP09(2018)005 discard one of the sums and restrict ourselves to only the n = 0 image in eq. (2.22). In that regime E L → 0, so we are justified in ignoring the non-universal part of the spectrum mentioned above. Restricting to only the n = 0 term, it turns out the remaining sum over poles can be formally done. One can show that the identity for F (E, χ, ∆) reduces to and ρ(E) as in eq. (2.32). Here, 2∆ F 2∆−1 and (a) n ≡ a(a + 1) . . . (a + n − 1) = Γ(a+n) Γ(a) being the Pochhammer symbol. We should comment that strictly speaking neither side has a well-defined decompactification limit L → ∞, but diverges. However, it is easy to see that under a scale transfor- Thus, as long as we keep the dimensionless variables EL and χL fixed and divide by the appropriate power of L, we can think of our expressions as being valid in a decompactifying limit L → ∞. In that limit, the n = 0 winding terms in the Wightman function are suppressed. Since the transforms we are doing are linear, we can neglect the images and focus on the n = 0 term as L → ∞. 7 3 Analysis of matrix elements and implications for ETH

Asymptotic formulae
In the previous section, we have given an exact analytical expression for F (E, χ, ∆) for any integer ∆, an arbitrary value of χL at large EL and large c. However, there are a few ap- 7 Additionally, if we take ∆ to be large then the n = 0 term dominates even at finite L.

JHEP09(2018)005
proximations we can do to get asymptotic expressions that are more transparent. To do so it is convenient to take one step back before we perform the Mellin transform. Formally inverting eq. (2.18) in the decompactification limit where we only keep the term n = 0, we have where the integral is taken over the Bromwich contour. The most obvious thing to do is a saddle-point approximation of the integral, in the regime where ELc 1. In that regime, e βE Z(β) is strongly peaked around β 0 = πcL 6E . Including the Gaussian corrections, this leads to This approximation converges to the correct value in the large E limit for any χ/E, but does so at different rates since the width of peak in the integrand depends on χ. It can therefore be used to get a good approximation of F (E, χ = aE, ∆) as E → ∞ for fixed χ/E = a, but should not be used to examine the χ/E = a dependence of that expression. Note that doing this same saddle-point approximation for dβe βE Z(β) recovers the Cardy formula with the correct prefactor, given in equation (2.33), as seen in [18]. This approach is similar to the one used in [9] to obtain the average of the diagonal matrix elements of O. Alternatively, we can expand F 0 (β, χ, ∆) in χ first, i.e. βχ, χL, χ E 1 and then do the inverse Laplace transform exactly, term by term. For χ where ψ (1) (∆) is the Polygamma function of order one. Performing then the contour integral using residue calculus in the identity eq. (3.1) gives: where I n is the modified Bessel function of order n. Note that all the subleading terms

JHEP09(2018)005
This is a good approximation to examine the χ dependence for a fixed value of χ/E 1. Furthermore, as EL → ∞, this approximation overlaps with the saddlepoint approximation, provided χ is held fixed and does not scale with E. Indeed, when χ = 0 = E a − E b , i.e. for the average of the squared matrix elements of the same energy, this approximation becomes exact, and we have 3.2 Does F (E, χ, ∆) reflect eigenstate thermalization?
It is instructive to compare the asymptotic results for F (E, χ, ∆) obtained in the previous subsection to the results predicted by eigenstate thermalization: the near thermal form of probe correlation functions in a closed many-body quantum system. There are more and more hints that the actual physics underlying black hole thermodynamics is this eigenstate thermalization, rather than conventional dynamical classical thermalization (see e.g. [19,20]). Suppose that ETH is indeed at work. Then the mean-square of an operator expectation value is well approximated by where the sums over a and b are restricted to states with energy E a and E b respectively. The term linear in R ab does not appear since it averages to zero. At χ = 0, ETH predicts that (3.8) A priori, this includes equal contributions from diagonal matrix elements (the O(E) term) and off-diagonal matrix elements of the same energy (the f (E, χ = 0) term). In particular, note that while the off-diagonal terms were suppressed by a factor of the density of states in eq. (3.6), they overcame that suppression in the average by their sheer number. Furthermore, we expect that O(E) is in fact small compared to f (E, χ = 0) since it is dual to the thermal one-point function of a free field (which vanishes at tree level). In [24], the one-loop JHEP09(2018)005 correction was obtained. Using that result and the exact expression for F (E, χ = 0, ∆), eq. (3.5), it is easy to see that, so long as E/c is finite as c → ∞, where |µ is the lightest state making the right-hand side non-zero and A depends on c only through E/c. µ|O|µ decays as c −1 or faster, as we would expect when comparing a oneloop effect to a tree-level result [24]. The suppression of any discontinuities at χ = 0 (the result at χ = 0 does not include a contribution from O(E)) implies that the natural way to view F (E, χ = 0, ∆) in holographic theories is as the variance of O(E). Using equation (3.5) it is easy to see that F (E, χ = 0, ∆) is peaked around Ec ∼ ∆ 2 . This is shown in figure 1.

Matrix elements for different energies
We now turn our attention to the matrix elements for states of different energies. Note that all our expressions are even in χ, as everything is symmetric under interchange of a and b.
is what characteristically distinguishes ETH from random matrix theory (see e.g. [25]), and such exponential behavior in particular has been seen in numerical simulations of ETH systems [26,27].

JHEP09(2018)005
where Ψ(E a ) is an arbitrary function of E a , subject only to the constraint that the state have finite norm: The Wightman function is then simply The function Ψ(E) gives us control over the entanglement between the two copies of the CFT. In particular Ψ(E) = e −βE/2 corresponds to the usual Boltzman distribution that leads to the TFD state, and Ψ(E) = δ E,0 corresponds to two decoupled CFTs in the ground state.
Our ability to compute Wightman functions in arbitrary entangled doubled-CFTs gives us a unique probe to study the spacetime behind the horizon, and its emergence/dissolution due to (dis)entanglement. In the thermally entangled TFD state the (logarithm of the) thermo-field double Wightman function for large ∆ computes the length of geodesics between the two sides of the Kruskal extended BTZ black hole. For small deviations from the thermal state, a geodesic interpretation should still hold, but clearly for large deviations any geometrical notion should break down. This is clear both from the gedankenexperiment by Van Raamsdonk [3] quoted in the introduction and the simple observation that at T = 0 the two copies of the CFT should be disentangled and hence disconnectedno geodesic exists. According to the AdS/CFT correspondence, the quantum extension of this family of geodesics is precisely given by the Wightman function above for arbitrary entanglement, i.e. for arbitrary double-sided quantum geometries. Choosing different entanglement functions Ψ(E a ) thus allows us to study the disentanglement of the spacetime behind the horizon, see figure 3. Note that the Feynman two-point function for large ∆ is dual to geodesics ending on the same boundary. These should be smooth, no matter the function Ψ(E a ), as is clear from taking Ψ(E a ) = δ Ea,0 .
We will choose here a rather drastic function Ψ(E) = e −βE/2 Θ(Λ−E) that interpolates discontinuously between the entangled thermofield double state and un-entangled vacuum beyond the cutoff Λ. Thus we consider the family of doubled states By construction these states have less entanglement between the two sides than the thermofield double at inverse temperature β. To make this precise, we can compute the entanglement entropy between the two sides. It is given by  Comparing to the entropy of the thermal state (the BTZ black hole) this gives This ratio S/S BTZ starts at zero and grows linearly as a function of the cutoff Λ until it saturates to 1, as seen in figure 4. Thus, varying Λ lets us directly dial the entanglement between the two copies of the CFT. The linearity and the saturation point can easily be understood using a saddle-point approximation. In the Cardy regime, ρ(E)e −βE ∼ e −βE+2 √ πc 6 EL is strongly peaked around E = πc 6 L β 2 . The width of this peak is 1 β πcL 3 . Therefore, we expect the average energy to match closely the thermal answer when the range of integration includes most of the peak, i.e when Λ ≥ πc 6 L β 2 + 1 β πcL 3 . This is indeed what we observe. The linearity can be explained by the fact that for E πcL 6β 2 the integral is dominated by the value of the integrand at E = Λ, and so E Λ ∼ Λ. Note that in the large c regime the relative width of the peak goes to zero since it scales as c −1/2 .

JHEP09(2018)005
We now use the Wightman function in the modified state |Λ , based on the entanglement function Ψ(E) = e −βE/2 Θ(Λ − E), to study what happens to the geometry as we vary the entanglement between the two boundaries. The leading behavior for large ∆ is related to the length of geodesics connecting the points (φ, t/2) on both boundaries 9 in the bulk dual of our modified state. The distance between the two boundaries is of course formally infinite. To regulate we compare to the geodesics computed in the thermal state dual to the BTZ black hole, i.e. we compute (4.8) We start our analysis at t = 0. Figure 5 shows that as Λ → ∞, we smoothly recover the black hole. The t = 0 geodesics are special in that they do not cross the horizon: they connect the two boundaries through the bifurcation surface. This indicates that as we decrease the entanglement between the two boundaries, the two patches move apart (as expected) but are not otherwise affected. In particular, we can use the Killing symmetry to move the endpoints of the geodesic up on one side and down on the other without affecting the answer. This corresponds to sweeping the entire geometry outside of the horizon. To see this, consider a point outside the horizon. This point always lies on a spacelike geodesic connecting a boundary point at some time (say t 0 ) and the bifurcation surface. By symmetry, that geodesic will always intersect the other boundary at Schwarzschild time −t 0 −i β 2 , i.e boundary time t 0 . As also remarked before, one sees that the bulk geometry dual to our modified state is not dramatically affected outside the horizon by the decrease in entanglement. In particular, we know that as Λ → 0, we simply recover two disconnected patches of AdS.
As we look at finite t = 0, we start probing the geometry behind the horizon. Now we experience the drastic effects of our severe change in entanglement. Figure 6 shows what happens for geodesics corresponding to early times. Instead of the ratio of the two-point functions increasing monotonically as a function of cutoff, it increases past one and then decreases again, eventually saturating to one. The intuition that the geometry behind the horizon simply becomes more connected in a smooth fashion as Λ → ∞ -which appears to be supported by the t = 0 geodesics through the bifurcation surface -is clearly wrong! This becomes even clearer as we increase the time, which corresponds to geodesics in the black hole background having a greater fraction of their length behind the horizon. Figure 7 shows that the changes in the ratio of the Wightman functions as a function of ΛL increase for larger t. The Wightman function G Λ (t) actually becomes negative for some values of 9 There is a symmetry that lets us shift the time of the insertions on both boundaries by the same amount in opposite direction. This corresponds to the action of the usual timelike Killing vector in the black hole geometry. Note that this symmetry is preserved in our modified state. Λ. 10 This signals a complete breakdown of the geodesic approximation. At best, one could argue that there are multiple geodesics that interfere destructively, with the path integral whose saddlepoints give the geodesics very sensitive to small variations in Λ. This cannot be expected to describe a classical space-time. Note that since the time shift symmetry is manifestly preserved in the expression for the Wightman function for any value of Λ, we cannot think of the geometric breakdown behind the horizon as being localized in one particular region.
Rather than the cut-off dependence for fixed t, one can also look at the time dependence for fixed cutoff. 11 Doing so dispels the notion that the non-geometric behaviour described previously only occurs for some values of ∆. For example, note that figure 7 shows that for t/L = 1/2π, ΛL ∼ 76 and ∆ = 4, the Wightman function is positive. However, we see from figure 8 that for the same cutoff and 0.20 t/L 0.35 the Wightman function for ∆ = 4 does become negative. 10 There is nothing from the field theory perspective that would prevent this: at finite t, the only constraint is that the Wightman function be real. 11 Note that there is of course no expectation that the difference in lengths saturate as t → ∞. Clearly, when the two copies of the CFT are too disentangled, no conventional notion of spacetime survives. Vice versa, by dialing the entanglement we can now literally see classical spacetime emerge/dissolve in a rather precise way. A geometric notion may survive for very large but finite Λ. Indeed, for all values of t (including t = 0), there is a value of Λ past which the Wightman function essentially matches the thermal answer. This value of the cutoff increases with c. This was to be expected, since the total energy of the black hole also increases in c and we certainly do not expect a classical bulk for Λ E BTZ . Given the strongly peaked nature of ρ(E)e −βE around E = πcL 6β 2 and of F (E, χ, ∆) around χ = 0, we expect large fluctuations around an energy of order E ∼ cL β 2 , within a window of width scaling as √ cL β . Our numerical results are consistent with this.

Summary and discussion
In this article, we have achieved two things. First, we have derived an expression for the average of the square of matrix elements of a scalar primary of integer dimension taken be-

JHEP09(2018)005
tween states of different energies. This is the function we call F (E, χ, ∆). We have obtained an exact, formal expression for this function valid for any holographic CFT at finite temperature and radius (in the deconfined phase). A more useful, Cardy-like formula, was also found. If interpreted in the context of the eigenstate thermalization hypothesis, F (E, χ, ∆) computes the square of the universal part of the variance, f (E, χ) for holographic CFTs. Secondly, we have demonstrated that the knowledge of this function can be used to make strong statements about the dual bulk geometry. More precisely, it allows us to compute the Wightman function for states with an arbitrary energy-dependent entanglement pattern between the two copies of the theory. We have focused on the most brutal possible modification to the thermofield double state, chopping off the contribution of all energy levels above some cutoff. This allowed us to dial the entanglement from a thermal state to two disentangled copies of the vacuum. We have seen explicitly that below some amount of entanglement (which scales with c), no clear notion of a semi-classical geometry behind the horizon can exist.
There are several possible extensions to this work. First, let us note that the space of entanglement patterns we can explore is now essentially the space of square-normalizable functions. An obvious generalization of our state |Λ is to remove states within an energy window. Another interesting avenue to explore may be how this function can help understand traversable wormholes [28] or dynamical space-times.

A Spatial dependence
Here we generalize our formal result from eq. (2.14) for operators inserted at different spatial coordinate φ. In this case, eq. (2.6) is where φ 12 ≡ φ 1 − φ 2 and t 12 ≡ t 1 − t 2 . As in the case of operators inserted at the same φ, we have split the sum into a sum over energies and defined F (E a , E b , ∆) in the same way, except that it contains an additional spatial depencende in φ 12 . The calculation follows in JHEP09(2018)005 exactly the same way as explained in section 2. Namely, we compute F (E a , E b , φ 12 , ∆) by inverting [cosh (2π/β(φ 12 + nL)) + cosh (2π(t 12 )/β)] ∆ . (A.4) The inversion of the Fourier transform has already been given in section 2.2, eqs. (2.14)-(2.16), where now b n (φ 12 ) = cosh(2π/β(φ 12 + nL)) instead. Similarly, the inversion of the Laplace transform carries through as explained in section 2.3. The final result is therefore, whereF n (E − E , χ, ∆), given in eq. (2.25), now depends on φ 12 through b n (φ 12 ) = cosh(2π/β(φ 12 + nL)). Note that F (E, χ, φ 12 , ∆) contains information about more than the spectrum of the CFT, as is necessary to reconstruct Wightman functions for operators inserted at different points in space.

B Comments on corrections to the density of sates
In section 2.3 we have obtained the density of states using the leading-in c partition function. We restricted to the region E − πc 6L > 0 so that we could ignore the subleading corrections in c to the partition function. Here, we explain how these non-universal corrections determine the density for E − πc 6L < 0. The partition function to leading order in c is [17] log Z = , where f contains all the subleading corrections. This parametrisation is chosen so that modular invariance is evident. We are now in a position to see that the density of states for E < πc where, as before, γ 0 is fixed so that the integration line is set to the right of all singularities of the integrand. For E > πc 6L the first integral is obtained by closing the contour to the left (Bromwich contour), while the second integral is obtained by closing the contour to the right giving zero contribution to the density of states. On the other hand for E < πc 6L , the situation is reversed and the density of states is obtained from the second integral with the contour closed to the left. Notice that, due to modular invariance, the corrections f (β) are needed non-perturbatively in β.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.