Correlation function with the insertion of zero modes of modular Hamiltonians

Zero modes of modular Hamiltonian of one interval are found in momentum space for two dimensional massless free scalar theory. Finite correlators are extracted from separate region connected correlation functions with the insertion of zero modes. Correlators of (n, 1)-type are claimed to be conformal block up to a set of theory dependent constants. We fix the correlators of (2, 1)-type with the knowledge of three point function coefficients in 2d CFTs. e-mail: longjiang@hust.edu.cn ar X iv :1 91 1. 11 48 7v 2 [ he pth ] 5 D ec 2 01 9


Introduction
The interest on a subregion is motivated by quantum information theory [1][2][3][4]. However, the domain of dependence D(A) of a spacelike subregion A is a well defined spacetime and one can study it as an independent system. Given the evolution equations, fields in D(A) are completely determined by the initial conditions on Cauchy surface A. Operators inside D(A) form a close algebra. The most important operator of subregion D(A) is modular HamiltonianĤ A which generates modular flow [5]. Modular Hamiltonian is the logarithmic of reduced density matrix ρ A ,Ĥ A = − logρ A . Usually, it is highly non-local, only in special cases it has analytic form [6][7][8][9]. For a conformal field theory in a state which has a gravitational dual, modular Hamiltonian is claimed [10,11] to be an area operator of Ryu-Takayanagi (RT ) surface [12] with quantum corrections.
Operators in D(A) could be classified according to their eigenvalues under modular transformation, From the point of view of quantum mechanics, operators whose commutator with modular Hamiltonian vanish (ω = 0) are rather intersting since they may provide new quantum numbers to classify operators. These operators are called zero modes of modular Hamiltonian [13,14]. In holography [15], bulk fields close to AdS boundary are identified with local CFT operators [16]. Correlation functions of local CFT operators are mapped to correlators of dual bulk fields according standard dictionary of AdS/CFT [17,18]. As a generalization, zero mode is claimed to have a gravitational dual [14] which is an integral of dual bulk field over RT surface. A nature question would be to understand correlators of zero modes from both CFT and gravity side.
The first step in this direction has been taken by [19] recently. In that paper, connected correlators of modular Hamiltonians are shown to be finite and universal for 2d CFTs. The universal property is a direct consequence of Ward identity of stress tensor [20]. As a result, correlators of modular Hamiltonians are free from quantum corrections in AdS 3 . In the present work, zero modes are inserted into correlation functions. Connected correlation functions with m zero modes inserted into region A and n zero modes in region B will be called (m, n)-type. We will mainly focus on (m, 1)-type correlators 2 . Any (m, 1)-type correlator may be a conformal block up to a theory dependent coefficient.
An outline of this paper is as follows. In section 2 we will review necessary background to study connected correlation functions in 2d massless free scalar theory. In section 3, zero modes of 2d massless free scalar theory are constructed in momentum space. Some of them are identified to OPE blocks. Then connected correlators with zero modes insertion are studied extensively in 2d massless free scalar theory in section 4. Correlators of (n, 1)-type have been discussed. Universal correlators of (2, 1)-type are found for general 2d CFTs in the same section. In section 5, we will briefly discuss operator product expansion(OPE) of deformed reduced density matrix.
We are able to interpret several universal correlators using OPE approach. We will discuss the results and point out some future directions in the last section.

Review
In this section we review basic results on modular Hamiltonian in momentum space. Conventions used in this paper are also introduced at the same time. We refer the reader to [19] for a detailed derivation. The system is in vacuum with Lagrangian in 2d Minkowski spacetime. The signature of spacetime is (−, +), spacetime coordinates are x µ = (t, z). The system can be factorized to left moving and right moving part, we will only focus on right moving modes, fields are functions of t − z only. The two point function is fixed where x, y are right moving spacetime coordinates. The spacelike region A (B) we will study is an interval with radius R A (R B ) whose center is located at z A (z B ), where the end points of the intervals are The theory can be quantized in D(A), the domain of dependence of A. More explicitly, field φ in D(A) can be expanded as where {g v } is a complete set of positive frequency modes. Annihilation and creation operators Vacuum in D(A) may be defined by

Modular Hamiltonian of region A in real space iŝ
where T tt is the stress tensor of massless free scalar and is evaluated at t = 0 timeslice. By transforming to momentum space, modular Hamiltonian becomes rather simplê The constant term may be fixed by requiring the normalization of reduced density matrix The field can also be quantized in Minkowski spacetime, where {f ω } forms a complete set of positive frequency modes in Minkowski spacetime. Annihilation and creation operators a ω , a † ω obey standard commutation relations. The Minkowski vacuum is defined by We will denote it as | to simplify notation. Annihilation and creation operators in region D(A) are related to those in Minkowski spacetime by Bogoliubov transformation where the Bogoliubov coefficients are They could be regarded as matrix elements of Bogoliubov matrices α, β. Any bounded operator in D(A) may always be rewritten as a function of annihilation and creation operators b v , b † v . Field in D(B) can be studied in a similar way. The corresponding annihilation and creation operators are denoted as bṽ, b † v . A quadratic operator has been studied extensively in [19]. The subscript I denotes possible quantum number. It is v in region D(A) andṽ in region D(B). ThereforeĤ could be a multi-region operator. Vacuum expectation value of the exponential ofĤ is shown to be where the matrix T is . We extract finite correlators of modular Hamiltonians from above results. However, the results actually don't requireĤ being linear superposition of modular Hamiltonians. More explicitly, x I could be any reasonable function of quantum number I. This sheds light on the study of zero modes in this paper. In modular fourier space,

Zero modes of modular Hamiltonian
one finds the commutation relation In the limit ω → 0, operator O 0 belongs to zero modes of modular Hamiltonian since it commutes with modular Hamiltonian. Zero mode is claimed to have a gravitional dual [14] O 0 (x) = 4π in the large N limit. Roughly speaking, zero mode can be expressed as an integral over RT surface [12] of the bulk dual field where the integral kernel is the corresponding bulk-to-boundary propagator. In the following subsections, we will construct zero modes of modular Hamiltonian of interval A in massless free scalar theory in momentum space. They have rather simple expressions and form an infinite dimensional algebra. Some zero modes are actually OPE blocks. We establish the relationship between zero modes in momentum space and OPE blocks at the end of this section.

Two dimensional massless free scalar theory
A basis of operators in D(A) could be chosen to be The commutator between modular Hamiltonian and b Therefore the commutator vanishes if and only if Hence a complete set of zero modes is found for massless free scalar. We list first few zero modes according to the degree of b and b † .
3. m + n = 2, zero modes are any linear superposition of The subscript "zm" denotes zero modes of modular Hamiltonian. An interesting fact is that (3.9) is exactly the form of (2.16), therefore all the conclusions in previous section can be applied to zero modes.
Zero modes with higher degree (≥ 3) may be useful in other context, however we will not study them in this paper. We also mention that all zero modes form an infinite dimensional algebra. Actually, from Jacobi identity Therefore, ifQ 1 andQ 2 are two zero modes and their commutator is non-zero, then [Q 1 ,Q 2 ] is also a zero mode.

Two dimensional conformal field theory
In CFT 2 , one example of zero mode is the following OPE block, The correlation function ofQ where G h (η) is two dimensional conformal block [24,25]

Zero modes in momentum space
There are infinite many chiral primary operators [26] with spin s and conformal weight h = s for two dimensional massless free boson. The coefficients A s k are fixed to be up to a k independent non-zero constant. The partial derivative operator is Since we are considering real boson, φ * = φ, odd spin currents J s vanish. The first few currents are is proportional to stress tensor.

s = 4,
For each even spin current, the corresponding zero mode iŝ All zero modes (3.26) are quadratic in terms of b, b † using the quantization (2.6). Therefore, they should be a linear superposition of quadratic zero modeŝ up to a constant. The constant term is proportional to identity which is also commute with modular Hamiltonian, it may appear due to commutator To determine the function Using the identity in Appendix A and summation in Appendix B, we find The integral in (3.26) leads to Dirac delta function, then C s (v, v ) has no contribution and one can read for even spin. We will choose the normalization coefficients c A [J s ] such that The first few functions f s (v) and normalization coefficients c A [J s ] are (3.37)

Correlation functions with the insertion of zero modes
In [19], connected correlation functions of modular Hamiltonians are finite. This is still true by inserting zero modes into connected correlation functions. In this section, we will establish the general framework to extract these finite correlators. Then we will compute several examples in massless free scalar theory. (n, 1)-type correlators are always conformal blocks in these examples.

Generator of connected correlation functions
The first step is to generalize the reduced density matrixρ We will call it deformed reduced density matrix since it is deformed by zero modeQ A . We still useρ A though it is already not the usual reduced density matrix. The subscript denotes the region. For simplicity, we just include one zero mode, the generalization to multiple zero modes is straightforward 4 . The generator of connected correlator is Connected correlators are defined as This is the connected correlators of (m+n, k + )-type with the insertion of zero modes. Since A and B are disjoint, all the operators are commute with each other, the order is not important. The generator is

Correlators in two dimensional massless free scalar theory
We will setQ A toQ A [J s ] for massless free scalar theory. Therefore similar computation as [19] leads to The explicit form of A, B, C, D are The function P is We also replace v and v to x and x , respectively. Expanding T A∪B (a, µ; b, ν) as (4.13) The first few orders of T n are They are where P (x, y; a, µ; b, ν) = xy sinh πf 1 sinh πf 2 sinh πx sinh πy sinh π(f 1 + x) sinh π(f 2 + y) , (4.18) The function F(x, x , y) is We define Several examples are shown in the following. We will only include spin 4 operator and restrict to m+n+k + ≤ 4 in this work. Therefore all zero modesQ A [J 4 ] is rewritten asQ A to simplify notation.

Two zero modes
We evaluate three correlators.
The correlator of two zero modes is proportional to conformal block which has been mentioned in (3.14). The first correlator has been studied in [19]. The second correlator is zero since the correlator of stress tensor and spin 4 primary operator is zero.

Three zero modes
Without losing generality, we insert only one zero mode into region B while other two zero modes are inserted into region A. Therefore we could evaluate the following six (2, 1)-type correlators.
The correlators are proportional to conformal block. We will reproduce this property later.

Four zero modes
There are two cases according to the number of zero modes in region B.
1. (3, 1)-type correlators. We could find eight different correlators. As correlators of (2, 1)-type, four zero modes correlators of (3, 1)-type are also proportional conformal block. This may indicate general properties of OPE of zero modes. We will discuss this point later.
2. (2, 2)-type correlators. The correlators are not proportional to a single conformal block, we don't find exact result. However, one can evaluate them in the small cross ratio limit order by order. We will present the leading order and leave the discussion on higher orders to another seperate paper.  wherec O is a function of a, µ. It should be the generator of coefficients c O . For massless free scalar theory, We check numerically that (4.51) and (4.52) indeed factorize for general positive a, µ, η. Therefore we match the lowest order of η in (4.51) and (4.52), , (4.53) .(4.54) They can be simplified tõ , (4.55) . (4.56) The integrals are rather hard to evaluate for general a, µ, however, we can expand them for small a and µ,c where Riemann Zeta function is The details of the integral can be found in Appendix C. The deformed reduced density matrix is still non-vanishing even when a → ∞. Usually, an exponential function e −aI(x) vanishes when I(x) > 0 in the limit a → ∞. In free boson case, modular Hamiltonian operatorĤ (4.72) All the modes with v > 0 should not contribute to deformed reduced density matrix in the limit a → ∞. The exception is the soft mode v → 0. We claim that soft mode contributes this non-zero effect. We expect to return to this problem in the near future.

Universal correlators
In previous subsection, three zero modes correlators are always proportional to conformal block. We will reproduce these results for general CFT 2 and find universal correlators with three zero modes in this subsection. (2, 1)-type correlators are where C 123 is the three point function coefficient (4.74) (4.75) where ": · · · :" means one should regularize the integral since it has poles near z 1 = z 2 . The parameter z is related to cross ratio by (4.76) In small cross ratio limit, We refer reader to reference [19] for more details on the parameter z and how to regularize the integral. We check that the regularized integral is indeed proportional to conformal block G h 3 (η) for various integer h 1 , h 2 , h 3 > 0. This is also supported by the explicit examples in previous subsection. Therefore one may assume the regularized integral is always To fix the coefficient λ(h 1 , h 2 , h 3 ) we expand both sides in the small cross ratio limit and match the coefficient before η h 3 , It matches the regularized integration for general positive integer h 1 , h 2 , h 3 . We discuss the regularization of this integral in Appendix D. Therefore,

OPE of deformed reduced density matrix
One of the intriguing result in previous sections is the connected correlation function between deformed reduced density matrix and another OPE block is proportional to conformal block This property may be true for general 2d CFTs. In this section, we will try to understand this result from the point of view of OPE. Similar discussion can be found in [19]. The normalized deformed reduced density matrixρ A (a, µ) is an exponential operator, it is a summation of different operators of corresponding CFT. They should be possible to reorganized as a summation of complete orthogonal operators defined in region D(A) with proper coefficients, we single out the operator O since this is the unique primary operator which has nonvanishing two point function withQ B [O]. The terms in "decendants" are fixed by conformal symmetry. We also note (3.14) is proportional to conformal block exactly, therefore operator product expansion of deformed reduced density matrix may bê One immediate consequence is that the correlator We present the contribution of conformal block G h O (η). The formula (5.5) is conformal block expansion of correlator ρ AρB . It has similar structure to conformal block approach to mutual information [27][28][29]. This is nature since conformal block is the object associated with four points and accommodate with conformal invariance. Several properties are collected blow: 1. In the limit a → 0, µ → 0,ρ A (a, µ) should be identity operator, therefore d O can be expanded around a = µ = 0 as Then we find As a polynomial of four variables a, µ, b, ν, it is at least degree of 2. For 2d massless free scalar, the first few s O (m, n, k, ) are listed in the Table 1 and 2.

Discussion and conclusion
The correlator with the insertion of zero modes of modular Hamiltonians in 2d massless free scalar theory has been studied. Connected correlators are still finite in the presence of zero modes, this is a direct generalization of the observation in [19]. There are several universal structures of these correlators.
1. Any (n, 1)-type correlator may be conformal block (4.49) up to a theory dependent coefficient. We check this point numerically for 2d massless free scalar theory. Analytical results are obtained up to n ≤ 3. We find the same conclusion using an OPE argument for general CFT 2 . Then we derive the leading order result for general correlators, see (5.12). Interestingly, leading order result is completely fixed using the knowledge of (n, 1)-type zero modes correlators. We could check (5.12) for 2d massless free scalar theory up to m + n + k + = 4.
2. (2, 1)-type correlator is universal, at least for positive integer conformal weight h 1 , h 2 , h 3 , see (4.80). We could fix the coefficient using the knowledge of three point function coefficients C ijk . The function λ(h 1 , h 2 , h 3 ) are obtained by regularization of a divergent integral. We could check (4.80) for massless free scalar theory. A rigorous proof of (4.80) is still lacking. It would be better to study (2, 1)-type correlator with non-integer conformal weight in the future.
To derive connected correlators, we used an exponential operator which is a deformation of reduced density matrix. It may be meaningful by itself. We briefly discussed its OPE. It would be better understand it to higher order. The correlator of zero modes is a rather interesting topic, there are many problems to be solved in this direction: 1. Gravitational dual. Zero mode of modular Hamiltonian in conformal field theory are claimed to be dual to integral over RT surface in large N limit in the gravitatational side, see (3.4). Correlator of zero modes may provide rich details on this conjecture. In AdS/CF T , the gravitational dual of correlator of primary operators has been well established in [17,18]. It would be interesting to understand how to generalize the story to zero modes correlators. Technically, the integral in the bulk becomes messy and suffers divergent problem when there are more than two zero modes in the correlator. It is great to see how to fix these problems.
2. Quantum corrections. Connected correlation functions of modular Hamiltonians are classical in the gravitional side [19]. This conclusion is not true when there are zero modes inserted. In the theory with W symmetry [30], the four point correlator of spin 3 current has O(1/c) correction [31]. Then it may provide an explicit example of O(1/N ) correction in the conjecture [14].
3. Zero modes which are not OPE blocks. In this case, the property (4.49) is not likely to be true. From 2d massless free scalar theory, their connected correlation functions are still finite. It would be better to generalize this point to other 2d CFTs.
4. It would be nice to generalize the discussion to higher dimensions.

B Summation
Assume k is a positive integer, then we will show The summation is independent of x. To prove this point, we will assume Therefore, At the first line, we used the identity of differentiation of Hypergeometric function At the second line, we relabel k to k + 1 in the summation. At the last step, we used the identity Therefore S(s; x) is independent of x, we can evaluate it at x = 0, This is actually a polynomial function of a, b, we redefine S(s; a, b) = −a(s − 1) The first few functions are It is even under the exchange of a, b for even spin and odd for odd spin S(s; a, b) = (−1) s S(s; b, a). (B.13)

C Integrals
In this work, we will use the following integrals, We use two identities The integral can be evaluated iteratively.

D Regularization
In this section, we will discuss on how to regularize the integral We will tackle this problem by trial and error using examples. The method has been used in [19]. We just use one example, h 1 = 2, h 2 = 3, h 3 = 3, to illustrate the regularization. Then we will list other results below.
J[2, 3, 3] = : At the first step, we integrate z 1 as if it is an indefinite integral. Since −1 ≤ z 2 ≤ 1, the integral is actually divergent around z 1 = z 2 . However, at the second step, we just ignore the divergence and taking the integration bound directly. This is the regularization used in this work. At the third step, the integral is well defined as usual. where κ = (−1) h + 1 which is zero when h is odd. These integrals are divergent superficially when h being some odd integers. However, they are actually zero since κ is zero in these cases. One could guess the general result to be . (D.11) Note to compute the integration of z 2 , we assume z 1 > 1 and then try to continue it to −1 ≤ z 1 ≤ 1. To calculate the integration of z 1 , we assume h 1 + h 2 − h 3 − 1 < 0 and then continue it to other region of h 3 . At the third line, one should keep the integral to be real since our original integration is real. To get the final result, we used the value of hypergeometric 2) is obtained using the regularization method of [19], this indicates that the result may be indepedent of regularization method. In the text, we also test this general result by matching it with free scalar theory .

Tables
This is a list of coefficients s O (m, n, k, ) for 2d massless free scalar. In the table, i is an integer.