## Entanglement Entropy of Disjoint Regions in Excited States : An Operator Method

## Noburo Shiba

Yukawa Institute for Theoretical Physics (YITP), Kyoto University, Kyoto 606-8502, Japan

Abstract: We develop the computational method of entanglement entropy based on the idea that Tr ρ n Ω is written as the expectation value of the local operator, where ρ Ω is a density matrix of the subsystem Ω . We apply it to consider the mutual Rényi information I ( n ) ( A,B ) = S ( n ) A + S ( n ) B -S ( n ) A B ∪ of disjoint compact spatial regions A and B in the locally excited states defined by acting the local operators at A and B on the vacuum of a ( d +1) -dimensional field theory, in the limit when the separation r between A and B is much greater than their sizes R A,B . For the general QFT which has a mass gap, we compute I ( n ) ( A,B ) explicitly and find that this result is interpreted in terms of an entangled state in quantum mechanics. For a free massless scalar field, we show that for some classes of excited states, I ( n ) ( A,B ) -I ( n ) ( A,B ) | r →∞ = C ( n ) AB /r α d ( -1) where α = 1 or 2 which is determined by the property of the local operators under the transformation φ → -φ and α = 2 for the vacuum state. We give a method to compute C (2) AB systematically.

## Contents

| 1 Introduction   | 1 Introduction                                                 | 1 Introduction                                                 |   1 |
|------------------|----------------------------------------------------------------|----------------------------------------------------------------|-----|
| 2                | Operator formalism                                             | Operator formalism                                             |   4 |
|                  | 2.1                                                            | Operator representation of Tr ρ n Ω                            |   4 |
|                  | 2.2                                                            | General properties of the glueing operator E Ω                 |   6 |
| 3                | Locally excited states in the general QFT which has a mass gap | Locally excited states in the general QFT which has a mass gap |   7 |
| 4                | Free scalar fields                                             | Free scalar fields                                             |  10 |
|                  | 4.1                                                            | Explicit calculation of the glueing operator E Ω               |  10 |
|                  |                                                                | 4.1.1 The case n = 2                                           |  12 |
|                  | 4.2                                                            | The mutual information of locally excited states               |  14 |
|                  |                                                                | 4.2.1 The case O i ′ A = O j ′ B = 0                           |  15 |
|                  |                                                                | 4.2.2 The case O iA , O jB , O i ′ A and O j ′ B are nonzero   |  16 |
| 5                | Conclusion and discussions                                     | Conclusion and discussions                                     |  17 |
| A                | Rényi entropy for the vacuum state in free scalar field theory | Rényi entropy for the vacuum state in free scalar field theory |  18 |

## 1 Introduction

The entanglement entropy in the quantum field theory (QFT) plays important roles in many fields of physics such as the string theory, condensed matter physics, and the physics of the black hole. The entanglement entropy is a useful quantity which characterize quantum properties of given states. For example, the entanglement entropy of ground states follows the area law [1-4] if we consider a local quantum field theory with a UV fixed point, while non-local field theories [5, 6] or QFTs with fermi surfaces [7] at UV cut off scale can violate the area law.

For a given density matrix ρ of the total system, the entanglement entropy of the subsystem Ω is defined as

$$S _ { \Omega } = - \text{Tr} \rho _ { \Omega } \ln \rho _ { \Omega },$$

where ρ Ω = Tr Ω c ρ is the reduced density matrix of the subsystem Ω and Ω c is the complement of Ω . The Rényi entropy S ( n ) Ω is defined as

$$S _ { \Omega } ^ { ( n ) } = \frac { 1 } { 1 - n } \ln \text{Tr} \rho _ { \Omega } ^ { n }.$$

The limit n → 1 coincides with the entanglement entropy lim n =1 S ( n ) Ω = S Ω .

In this paper we develop the computational method of Rényi entanglement entropy based on the idea that Tr ρ n Ω is written as the expectation value of the local operator at Ω . We apply this method to consider the mutual Rényi information I ( n ) ( A,B ) = S ( n ) A + S ( n ) B -S ( n ) A B ∪ of disjoint compact spatial regions A and B in the locally excited states defined by acting the local operators at A and B on the vacuum of a ( d +1) -dimensional field theory, in the limit when the separation r between A and B is much greater than their sizes R A,B .

Our method is based on the idea that Tr ρ n Ω is written as the expectation value of the local operator at Ω . This idea was originally used to compute I ( n ) ( A,B ) in the vacuum state by Cardy [8], Calabrese et al. [9] and Headrick [10]. We generalize this idea to an arbitrary density matrix ρ and construct explicitly the local operator. The density matrix of the total system ρ can be a mixed state and an excited state. We consider the general scalar field and do not specify its Hamiltonian. (Our method is applicable to QFT with interaction. ) We summarize our method. We consider n copies of the scalar fields and the j -th copy of the scalar field is denoted by { φ ( j ) } . Thus the total Hilbert space, H ( n ) , is the tensor product of the n copies of the Hilbert space, H ( n ) = H ⊗ H · · · ⊗ H where H is the Hilbert space of one scalar field. We define the density matrix ρ ( n ) in H ( n ) as

$$\rho ^ { ( n ) } \equiv \rho \otimes \rho \otimes \cdots \otimes \rho$$

where ρ is an arbitrary density matrix in H . We can express Tr ρ n Ω as

$$\text{Tr} \rho _ { \Omega } ^ { n } = \text{Tr} ( \rho ^ { ( n ) } E _ { \Omega } ),$$

where

$$\text{where} \\ E _ { \Omega } = \int \prod _ { j = 1 } ^ { n } \prod _ { x \in \Omega } D J ^ { ( j ) } ( x ) D K ^ { ( j ) } ( x ) \exp [ i \int d ^ { d } x \sum _ { l = 1 } ^ { n } J ^ { ( l + 1 ) } ( x ) \phi ^ { ( l ) } ( x ) ] \\ \times & \exp [ i \int d ^ { d } x \sum _ { l = 1 } ^ { n } K ^ { ( l ) } ( x ) \pi ^ { ( l ) } ( x ) ] \times \exp [ - i \int d ^ { d } x \sum _ { l = 1 } ^ { n } J ^ { ( l ) } \phi ^ { ( l ) } ], \\. & \quad. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot. & \cdot.$$

where π x ( ) is a conjugate momenta of φ x ( ) , [ φ x , π ( ) ( y )] = iδ d ( x -y ) , and J ( j ) ( x ) and K ( j ) ( x ) exist only in Ω and J ( n +1) = J (1) and we normalize the measure of the functional integral as ∫ DJ ( j ) exp[ i ∫ d xJ d ( j ) ( x f ) ( x )] = ∏ x ∈ Ω δ f ( ( x )) where f ( x ) is an arbitrary function. Notice that φ and π in (1.5) are operators and the ordering is important. We call this operator E Ω as the glueing operator . When ρ is a pure state, ρ = Ψ | 〉 〈 Ψ | , the equation (1.4) becomes

$$\text{Tr} \rho _ { \Omega } ^ { n } = \langle \Psi ^ { ( n ) } | \, E _ { \Omega } \, | \Psi ^ { ( n ) } \rangle$$

where

$$| \Psi ^ { ( n ) } \rangle = | \Psi \rangle \, | \Psi \rangle \dots | \Psi \rangle \,.$$

For a free scalar field, we can rewrite E Ω in (1.5) using the normal ordering. In the case n = 2 , we obtain a simple expression of E Ω and reproduce the result that I (2) ( A,B ) in the vacuum state is proportional to the product of the electrostatic capacitance of each regions obtained by Cardy [8]. Furthermore the simple expression of E Ω is useful for numerical calculation.

The advantages of this method are that we can use ordinary technique in QFT such as OPE and the cluster decomposition property and that we can use the general properties and the explicit expression of the glueing operator to compute systematically the Rényi entropy for an arbitrary state.

We apply this method to the mutual Rényi information I ( n ) ( A,B ) in the locally excited states. We consider the following locally excited state,

$$\left | \Psi \right \rangle = N ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) \left | 0 \right \rangle,$$

where N is a real normalization constant and O iA and O i A ′ ( O jB and O j B ′ ) are operators on A ( B ) and i and i ′ ( j and j ′ ) label a kind of operators. For the general QFT which has a mass gap, we compute I ( n ) ( A,B ) explicitly and find that this result is interpreted in terms of an entangled state in a quantum mechanical system which has finite degrees of freedom. For a free massless scalar field, we show that for some classes of excited states, I ( n ) ( A,B ) -I ( n ) ( A,B ) | r →∞ = C ( n ) AB /r α d ( -1) where α = 1 or 2 which is determined by the property of the local operators under the transformation φ →-φ and α = 2 for the vacuum state.

The mutual information for a free scalar field in higher dimensions has been studied in only a few papers [8, 11-15, 24]. In these papers, the authors considered only the mutual information for the vacuum state. The relation between the mutual information and the physics of the black hole was considered in [13]. Recently, it is proposed [16] that the mutual information is obtained by the quantum correction to the holographic entanglement entropy formula [17]. It would be interesting to use our results to check this proposition.

The entanglement entropy for an excited state defined by acting the local operator on the vacuum was considered in [18-21]. In [18], the subsystem is a half of the total space and the local operator exists at the complement of the subsystem and the time evolution of the entanglement entropy was considered. It was found that the entanglement entropy at late time is interpreted in terms of an entangled state in quantum mechanics [18]. This result is analogous to our result in the general QFT which has a mass gap.

The present paper is organized as follows. In section 2.1 we develop the computational method of the entanglement entropy. We derive the basic formula (1.4) and construct explicitly the glueing operator E Ω . In section 2.2 we investigate the general properties of E Ω . In section 3 we consider the mutual Rényi information I ( n ) ( A,B ) in the locally excited states in the general QFT which has a mass gap. We compute I ( n ) ( A,B ) explicitly and find that this result is interpreted in terms of an entangled state in quantum mechanics. In section 4 we consider free scalar fields. In section 4.1 we rewrite E Ω in (1.5) using the normal ordering. In the case n = 2 , we obtain a simple expression of E Ω and reproduce the result that I (2) ( A,B ) in the vacuum state is proportional to the product of the electrostatic capacitance of each regions obtained by Cardy. In section 4.2 we consider a free massless scalar field. We show that for some classes of excited states I ( n ) ( A,B ) -I ( n ) ( A,B ) | r →∞ = C ( n ) AB /r α d ( -1) where α = 1 or 2 which is determined by the property of the local operators under the transformation φ → -φ and α = 2 for the vacuum state. In section 5 we summarize our conclusion.

## 2 Operator formalism

## 2.1 Operator representation of Tr ρ n Ω

We represent the trace of the n th power of the reduced density matrix as the expectation value of the local operator. As a model amenable to unambiguous calculation we deal with the scalar field as a collection of coupled oscillators on a lattice of space points, labeled by capital Latin indices, the displacement at each point giving the value of the scalar field there. The local Hermitian variables ˆ q A and ˆ p B (coordinates and the conjugate momentum) obey the canonical commutation relations

$$[ \hat { q } _ { A }, \hat { p } _ { B } ] = i \delta _ { A B }, \quad [ \hat { q } _ { A }, \hat { q } _ { B } ] = [ \hat { p } _ { A }, \hat { p } _ { B } ] = 0.$$

The density matrix ρ of the total system in coordinate representation is

$$\rho ( q _ { A } ; q _ { B } ^ { \prime } ) = \langle \{ q _ { A } \} | \, \rho \, | \{ q _ { B } ^ { \prime } \} \rangle$$

where { q A } denotes the collection of all q A 's. We consider the arbitrary density matrix ρ . Now consider a subsystem (or subregion) Ω in the space. The oscillators in this region will be specified by lowercase Latin letters, and those in its complement Ω c will be specified by Greek letters. We can obtain a reduced density matrix ρ Ω for Ω by integrating out over q α ∈ R for each of the oscillators in Ω c , and then we have

$$\rho _ { \Omega } ( q _ { a } ; q _ { b } ^ { \prime } ) = \int \prod _ { \alpha } d q _ { \alpha } \rho ( q _ { a }, q _ { \alpha } ; q _ { b } ^ { \prime }, q _ { \alpha } )$$

We obtain the trace of the n th power of the reduced density matrix ρ Ω as

$$\text{Tr} \rho _ { \Omega } ^ { n } = \int \prod _ { j = 1 } ^ { n } \prod _ { a } d q _ { a } ^ { ( j ) } \rho _ { \Omega } ( q _ { a } ^ { ( 1 ) } ; q _ { a } ^ { ( 2 ) } ) \rho _ { \Omega } ( q _ { a } ^ { ( 2 ) } ; q _ { a } ^ { ( 3 ) } ) \dots \rho _ { \Omega } ( q _ { a } ^ { ( n ) } ; q _ { a } ^ { ( 1 ) } )$$

$$= \int \prod _ { j = 1 } ^ { n } \prod _ { a } ^ { k } d q _ { a } ^ { ( j ) } \prod _ { \alpha } d q _ { \alpha } ^ { ( j ) } \rho ( q _ { a } ^ { ( 1 ) }, q _ { \alpha } ^ { ( 1 ) } ; q _ { a } ^ { ( 2 ) }, q _ { \alpha } ^ { ( 1 ) } ) \rho ( q _ { a } ^ { ( 2 ) }, q _ { \alpha } ^ { ( 2 ) } ; q _ { a } ^ { ( 3 ) }, q _ { \alpha } ^ { ( 2 ) } ) \dots \rho ( q _ { a } ^ { ( n ) }, q _ { \alpha } ^ { ( n ) } ; q _ { a } ^ { ( 1 ) }, q _ { \alpha } ^ { ( 1 ) } ).$$

(2.4)

We consider n copies of the oscillators and the j -th copy of the oscillators is denoted by { q ( j ) A } . Thus the total Hilbert space, H ( n ) , is the tensor product of the n copies of the Hilbert space, H ( n ) = H ⊗ H · · · ⊗ H . We define the following density matrix,

$$\rho ^ { ( n ) } ( q _ { A } ^ { ( 1 ) }, \dots, q _ { A } ^ { ( n ) } ; q _ { B } ^ { ( 1 ) ^ { \prime } }, \dots, q _ { B } ^ { ( n ) ^ { \prime } } ) \equiv \rho ( q _ { A } ^ { ( 1 ) } ; q _ { B } ^ { ( 1 ) ^ { \prime } } ) \rho ( q _ { A } ^ { ( 2 ) } ; q _ { B } ^ { ( 2 ) ^ { \prime } } ) \dots \rho ( q _ { A } ^ { ( n ) } ; q _ { B } ^ { ( n ) ^ { \prime } } ), \quad ( 2. 5 )$$

i.e.

$$\rho ^ { ( n ) } \equiv \rho \otimes \rho \otimes \cdots \otimes \rho.$$

Then we can rewrite Tr ρ n Ω as

$$\mathbb { T } _ { \prime }$$

$$\text{Then we can rewrite } & \text{Tr} \rho _ { \Omega } ^ { n } \text{ as } \\ & \text{Tr} \rho _ { \Omega } ^ { n } = \int \prod _ { j = 1 } ^ { n } \prod _ { a } d q _ { a } ^ { ( j ) } \prod _ { \alpha } d q _ { \alpha } ^ { ( j ) } \prod _ { b } d q _ { b } ^ { ( j ) ^ { \prime } } \prod _ { \beta } d q _ { \beta } ^ { ( j ) ^ { \prime } } \\ & \times \rho ( q _ { a } ^ { ( 1 ) }, q _ { \alpha } ^ { ( 1 ) } ; q _ { b } ^ { ( 1 ) ^ { \prime } }, q _ { b } ^ { ( 1 ) ^ { \prime } } ) \rho ( q _ { a } ^ { ( 2 ) }, q _ { a } ^ { ( 2 ) } ; q _ { b } ^ { ( 2 ) ^ { \prime } }, q _ { \beta } ^ { ( 2 ) ^ { \prime } } ) \dots \rho ( q _ { a } ^ { ( n ) }, q _ { a } ^ { ( n ) } ; q _ { b } ^ { ( n ) ^ { \prime } } ) \\ & \times \delta ( q _ { a } ^ { ( 1 ) } - q _ { \beta } ^ { ( 1 ) ^ { \prime } } ) \delta ( q _ { a } ^ { ( 2 ) } - q _ { \beta } ^ { ( 2 ) ^ { \prime } } ) \dots \delta ( q _ { b } ^ { ( n ) ^ { \prime } } - q _ { a } ^ { ( 1 ) } ) \\ & \times \delta ( q _ { b } ^ { ( 1 ) ^ { \prime } } - q _ { a } ^ { ( 2 ) } ) \delta ( q _ { b } ^ { ( 2 ) ^ { \prime } } - q _ { a } ^ { ( 3 ) } ) \dots \delta ( q _ { b } ^ { ( n ) ^ { \prime } } - q _ { a } ^ { ( 1 ) } ) \\ & = \int \prod _ { j = 1 } ^ { n } \prod _ { a } d q _ { A } ^ { ( j ) } \prod _ { B } d q _ { B } ^ { ( j ) ^ { \prime } } \rho ^ { ( n ) } ( q _ { A } ^ { ( 1 ) }, \dots, q _ { A } ^ { ( n ) } ; q _ { B } ^ { \prime }, \dots, q _ { B } ^ { \prime } ) \\ & \times E _ { \Omega } ( q _ { B } ^ { ( 1 ) ^ { \prime } }, \dots, q _ { B } ^ { ( n ) ^ { \prime } } ; q _ { A } ^ { ( 1 ) }, \dots, q _ { A } ^ { ( n ) } ) \\ & = \text{Tr} ( \rho ^ { ( n ) } E _ { \Omega } ), \\ & \text{where}$$

where

$$& E _ { \Omega } ( q _ { B } ^ { ( 1 ) ^ { \prime } }, \dots, q _ { B } ^ { ( n ) ^ { \prime } } ; q _ { A } ^ { ( 1 ) }, \dots, q _ { A } ^ { ( n ) } ) \equiv \langle \{ q _ { B } ^ { ( 1 ) ^ { \prime } } \}, \dots, \{ q _ { B } ^ { ( n ) ^ { \prime } } \} | E _ { \Omega } \, | \{ q _ { A } ^ { ( 1 ) } \}, \dots, \{ q _ { A } ^ { ( n ) } \} \rangle \\ & \equiv \prod \prod _ { a } \delta ( q _ { \alpha } ^ { ( 1 ) } - q _ { \beta } ^ { ( 1 ) ^ { \prime } } ) \delta ( q _ { \alpha } ^ { ( 2 ) } - q _ { \beta } ^ { ( 2 ) ^ { \prime } } ) \dots \delta ( q _ { \alpha } ^ { ( n ) } - q _ { \beta } ^ { ( n ) ^ { \prime } } ) \times \delta ( q _ { b } ^ { ( 1 ) ^ { \prime } } - q _ { a } ^ { ( 2 ) } ) \delta ( q _ { b } ^ { ( 2 ) ^ { \prime } } - q _ { a } ^ { ( 3 ) } ) \dots \delta ( q _ { b } ^ { ( n ) ^ { \prime } } - q _ { a } ^ { ( 1 ) } ).$$

We call E Ω as the glueing operator . When ρ is a pure state, ρ = Ψ | 〉 〈 Ψ | , the equation (2.7) becomes

$$\text{Tr} \rho _ { \Omega } ^ { n } = \langle \Psi ^ { ( n ) } | \, E _ { \Omega } \, | \Psi ^ { ( n ) } \rangle$$

where

$$| \Psi ^ { ( n ) } \rangle = | \Psi \rangle \, | \Psi \rangle \dots | \Psi \rangle \,.$$

We represent E Ω as a function of ˆ q and ˆ p . We can rewrite E Ω as

$$& \text{We represent } E _ { \Omega } \text{ as a function of }$$

where we have written the delta functions as the Fourier integrals. Note that ˆ q ( ) l a is a operator and not a integral variable. The middle term in (2.11) is the tensor product of the following operator,

$$\Phi _ { \Phi }, \quad \int d q ^ { \prime } \$$

where we have omitted the subscripts for simplicity. We can check easily that (2.12) is correct just by taking the matrix elements of its both sides. Thus we can rewrite the middle term in (2.11) as

$$\begin{array} {$$

We substitute (2.13) into (2.11) and obtain

$$& E _ { \Omega } = \int \prod _ { j = 1 } ^ {$$

From (2.14), for the ( d +1) dimensional scalar field theory, we obtain

$$\dots \dots \dots \dots \dots \dots \dots \$$

where π x ( ) is a conjugate momenta of φ x ( ) , [ φ x , π ( ) ( y )] = iδ d ( x -y ) , and J ( j ) ( x ) and K ( j ) ( x ) exist only in Ω and J ( n +1) = J (1) and we normalize the measure of the functional integral as ∫ DJ ( j ) exp[ i ∫ d xJ d ( j ) ( x f ) ( x )] = ∏ x ∈ Ω δ f ( ( x )) where f ( x ) is an arbitrary function.

## 2.2 General properties of the glueing operator E Ω

We investigate some general properties of E Ω .

(1) Symmetry: From (2.15), E Ω is invariant under the sign changing transformation φ, π →-φ, -π , i.e.

$$E _ { \Omega } ( \phi ^ { ( 1 ) }, \dots, \phi ^ { ( n ) }, \pi ^ { ( 1 ) }, \dots, \pi ^ { ( n ) } ) = E _ { \Omega } ( - \phi ^ { ( 1 ) }, \dots, - \phi ^ { ( n ) }, - \pi ^ { ( 1 ) }, \dots, - \pi ^ { ( n ) } ). \quad \ ( 2. 1 6 )$$

- (2) Locality: When Ω = A ∪ B and A ∩ B = ∅ ,

$$E _ { A \cup B } = E _ { A } E _ { B }.$$

- (3) For n arbitrary operators F j ( j = 1 2 , , . . . , n ) on H ,

$$\text{Tr} ( F _ { 1 } \otimes F _ { 2 } \otimes \dots \otimes F _ { n } \cdot E _ { \Omega } ) = \text{Tr} ( F _ { 1 \Omega } F _ { 2 \Omega } \dots F _ { n \Omega } ),$$

where F j Ω ≡ Tr Ω c F j . This property is the generalization of (2.7). We can prove easily (2.18) by using the matrix elements of E Ω in (2.8).

(4) The cyclic property: From the cyclic property of the trace in the right hand side in (2.18), we obtain

$$\text{Tr} ( F _ { 1 } \otimes F _ { 2 } \otimes \dots \otimes F _ { n } \cdot E _ { \Omega } ) = \text{Tr} ( F _ { 2 } \otimes F _ { 3 } \otimes \dots \otimes F _ { n } \otimes F _ { 1 } \cdot E _ { \Omega } ). \quad \ ( 2. 1 9 )$$

(5) The relation between E Ω and E Ω c for pure states: We consider two pure states | φ 1 〉 | φ 2 〉 . . . | φ n 〉 and | ψ 1 〉 | ψ 2 〉 . . . | ψ n 〉 in H ( n ) where | φ j 〉 and | ψ j 〉 ( j = 1 2 , , . . . , n ) are arbitrary pure states. We can prove the following equation:

$$\langle \psi _ { 1 } | \, \langle \psi _ { 2 } | \dots \langle \psi _ { n } | \, E _ { \Omega } \, | \phi _ { 1 } \rangle \, | \phi _ { 2 } \rangle \dots | \phi _ { n } \rangle = [ \langle \phi _ { 2 } | \, \langle \phi _ { 3 } | \dots \langle \phi _ { n } | \, \langle \phi _ { 1 } | \, E _ { \Omega ^ { c } } \, | \psi _ { 1 } \rangle \, | \psi _ { 2 } \rangle \dots | \psi _ { n } \rangle ] ^ { * }.$$

We prove (2.20) as follows: From (2.8) we obtain

$$\langle \psi _ { 1 } | \, \langle \psi _ { 2 } | \dots \langle \psi _ { n } | \, E _ { \Omega } \, | \phi _ { 1 } \rangle \, | \phi _ { 2 } \rangle \dots | \phi _ { n } \rangle = \int \prod _ { \substack { j = 1 \\ j = 1 \\ k } } ^ { n } \prod _ { a } d q _ { a } ^ { ( j ) } \prod _ { \alpha } d q _ { \alpha } ^ { ( j ) }$$

× ψ 1 ( q (1) a , q (1) α ) ∗ ψ 2 ( q (2) a , q (2) α ) ∗ . . . ψ n ( q ( n ) a , q ( n ) α ) ∗ φ 1 ( q ( n ) a , q (1) α ) φ 2 ( q (1) a , q (2) α ) . . . φ n ( q ( n -1) a , q ( n ) α )

$$= [ \int \prod _ { j = 1 } ^ { n } \prod _ { a } d q _ { a } ^ { ( j ) } \prod _ { \alpha } d q _ { \alpha } ^ { ( j ) } \phi _ { 2 } ( q _ { a } ^ { ( 1 ) }, q _ { \alpha } ^ { ( 2 ) } ) ^ { * } \phi _ { 3 } ( q _ { a } ^ { ( 2 ) }, q _ { \alpha } ^ { ( 3 ) } ) ^ { * } \dots \phi _ { n } ( q _ { a } ^ { ( n - 1 ) }, q _ { \alpha } ^ { ( n ) } ) ^ { * } \phi _ { 1 } ( q _ { a } ^ { ( n ) }, q _ { \alpha } ^ { ( 1 ) } ) ^ { * }$$

$$& \stackrel { \nu } { \psi _ { 1 } } ( q _ { a } ^ { ( 1 ) }, q _ { \alpha } ^ { ( 1 ) } ) \psi _ { 2 } ( q _ { a } ^ { ( 2 ) }, q _ { \alpha } ^ { ( 2 ) } ) \dots \psi _ { n } ( q _ { a } ^ { ( n ) }, q _ { \alpha } ^ { ( n ) } ) ] ^ { * } \\ & = [ \langle \phi _ { 2 } | \, \langle \phi _ { 3 } | \dots \langle \phi _ { n } | \, \langle \phi _ { 1 } | \, E _ { \Omega ^ { e } } \, | \psi _ { 1 } \rangle \, | \psi _ { 2 } \rangle \dots | \psi _ { n } \rangle ] ^ { * }. \quad \square$$

By using (2.20), we can obtain the following basic property (see e.g. [22]) for a pure state ρ = Ψ | 〉 〈 Ψ | ;

$$\text{Tr} \rho _ { \Omega } ^ { n } = \langle \Psi ^ { ( n ) } | \, E _ { \Omega } \, | \Psi ^ { ( n ) } \rangle = [ \langle \Psi ^ { ( n ) } | \, E _ { \Omega ^ { c } } \, | \Psi ^ { ( n ) } \rangle ] ^ { * } = \text{Tr} \rho _ { \Omega ^ { c } } ^ { n }.$$

where we have used the fact that Tr ρ n Ω is real. So (2.20) is the generalization of Tr ρ n Ω = Tr ρ n Ω c for a pure state ρ = Ψ | 〉 〈 Ψ | .

## 3 Locally excited states in the general QFT which has a mass gap

We consider the mutual Rényi information I ( n ) ( A,B ) ≡ S ( n ) A + S ( n ) B -S ( n ) A B ∪ of disjoint compact spatial regions A and B in the locally excited states of the general QFT which has the mass gap m in the limit when the separation r between A and B is much greater than their sizes R A,B and 1 /m r ( glyph[greatermuch] R A,B , 1 /m ). We consider the following locally excited state,

$$\left | \Psi \right \rangle = N ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) \left | 0 \right \rangle,$$

where N is a real normalization constant and O iA and O i A ′ ( O jB and O j B ′ ) are operators on A B ( ) and i and i ′ ( j and j ′ ) label a kind of operators. We assume the distance between the positions of O A ( O B ) and the boundary of A B ( ) is much greater than 1 /m R ( A,B glyph[greatermuch] 1 /m )

(we have omitted the labels i s for simplicity.) We impose following orthogonal conditions for simplicity,

$$\left \langle 0 \right | O _ { i A } ^ { \dagger } O _ { i ^ { \prime } A } \left | 0 \right \rangle = \left \langle 0 \right | O _ { j B } ^ { \dagger } O _ { j ^ { \prime } B } \left | 0 \right \rangle = 0.$$

This state is similar to the EPR state. We compute the mutual Rényi information of this state. In this calculation, the general properties (2), (5) and the cluster decomposition property in the QFT play important roles.

From the normalization condition 〈 Ψ Ψ | 〉 = 1 we obtain

$$N ^ { - 2 } \simeq \left \langle 0 \right | O _ { i A } ^ { \dagger } O _ { i A } \left | 0 \right \rangle \left \langle 0 \right | O _ { j B } ^ { \dagger } O _ { j B } \left | 0 \right \rangle + \left \langle 0 \right | O _ { i ^ { \prime } A } ^ { \dagger } O _ { i ^ { \prime } A } \left | 0 \right \rangle \left \langle 0 \right | O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } \left | 0 \right \rangle, \quad \left ( 3. 3 \right )$$

where we have used the cluster decomposition property 〈 0 | O O O O † A A † B B | 0 〉 glyph[similarequal] 〈 0 | O O † A A | 0 〉 〈 0 | O O † B B | 0 〉 and the orthogonal conditions (3.2). By using (2.9), we obtain Tr ρ n Ω for ρ = Ψ | 〉 〈 Ψ | as

$$& \text{Tr} \rho _ { \Omega } ^ { n } = \langle \Psi ^ { ( n ) } | \, E _ { \Omega } \, | \Psi ^ { ( n ) } \rangle \\ & = N ^ { 2 n } \, \langle 0 ^ { ( n ) } | \, ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { j ^ { \prime } B } ^ { \dagger } ) ^ { ( 1 ) } \dots ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { j ^ { \prime } B } ^ { \dagger } ) ^ { ( n ) } \\ & \times E _ { \Omega } ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) ^ { ( 1 ) } \dots ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) ^ { ( n ) } \, | 0 ^ { ( n ) } \rangle \,,$$

where O ( ) l are operators in the Hilbert space of the l -th copy and the Ω is A,B or A ∪ B . First, we consider Tr ρ n A B ∪ . From (2.20) and (3.4) we obtain

$$& \text{Tr} \rho _ { A \cup B } ^ { n } = N ^ { 2 n } [ \langle 0 ^ { ( n ) } | \, ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { j ^ { \prime } B } ^ { \dagger } ) ^ { ( 1 ) } \dots ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { j ^ { \prime } B } ^ { \dagger } ) ^ { ( n ) } \\ & \times E _ { ( A \cup B ) ^ { * } } ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) ^ { ( 1 ) } \dots ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) ^ { ( n ) } \, | 0 ^ { ( n ) } \rangle ] ^ { * } \\ & \simeq \langle \Psi | \Psi \rangle \, [ \langle 0 ^ { ( n ) } | \, E _ { ( A \cup B ) ^ { c } } \, | 0 ^ { ( n ) } \rangle ] ^ { * } = \langle 0 ^ { ( n ) } | \, E _ { A \cup B } \, | 0 ^ { ( n ) } \rangle = \langle 0 ^ { ( n ) } | \, E _ { A } E _ { B } \, | 0 ^ { ( n ) } \rangle \\ & \simeq \langle 0 ^ { ( n ) } | \, E _ { A } \, | 0 ^ { ( n ) } \rangle \, \langle 0 ^ { ( n ) } | \, E _ { B } \, | 0 ^ { ( n ) } \rangle = \text{Tr} \rho _ { 0 A } ^ { n } \text{Tr} \rho _ { 0 B } ^ { n },$$

where we have used the cluster decomposition property and the conditions r, R A , R B glyph[greatermuch] 1 /m and ρ 0 A B ( ) is the reduced density matrix of the vacuum state.

Next we consider Tr ρ n A . We expand the product in (3.4). The terms in the expansion in (3.4) have the following form,

$$& \ln \{ 0, 4 \} \text{ have the following form}, \\ & \langle 0 ^ { ( n ) } | O ^ { ( 1 ) } _ { i _ { 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } E _ { A } O ^ { ( 1 ) } _ { i _ { n + 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } \cdot O ^ { ( 1 ) } _ { j _ { 1 } B } \dots O ^ { ( n ) } _ { j _ { n } B } O ^ { ( 1 ) } _ { j _ { n + 1 } B } \dots O ^ { ( n ) } _ { j _ { n } B } | 0 ^ { ( n ) } \rangle \\ & \simeq ( 0 ^ { ( n ) } | O ^ { ( 1 ) } _ { i _ { 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } E _ { A } O ^ { ( 1 ) } _ { i _ { n + 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } | 0 ^ { ( n ) } \rangle \, ( 0 ^ { ( n ) } | O ^ { ( 1 ) } _ { j _ { 1 } B } \dots O ^ { ( n ) } _ { j _ { n } B } O ^ { ( 1 ) } _ { j _ { n + 1 } B } \dots O ^ { ( n ) } _ { j _ { n } B } | 0 ^ { ( n ) } \rangle \\ & \quad = \langle 0 ^ { ( n ) } | O ^ { ( 1 ) } _ { i _ { 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } E _ { A } O ^ { ( 1 ) } _ { i _ { n + 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } | 0 ^ { ( n ) } \rangle \prod _ { l = 1 } ^ { n } \langle 0 | O ^ { \dagger } _ { j _ { l } B } O _ { j _ { l + n } B } | 0 \rangle \\ & \quad = ( 0 ^ { ( n ) } | O ^ { ( 1 ) } _ { i _ { 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } E _ { A } O ^ { ( 1 ) } _ { i _ { n + 1 } A } \dots O ^ { ( n ) } _ { i _ { n } A } | 0 ^ { ( n ) } \rangle \prod _ { l = 1 } ^ { n } \delta _ { j _ { l } j _ { l + n } } \langle 0 | O ^ { \dagger } _ { j _ { l } B } O _ { j _ { l } B } | 0 \rangle \\ & \quad \text{where} \, & \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \, \iota \end{cases}$$

where ( i , j l l ) is ( i, j ) or ( i , j ′ ′ ) ( l = 1 2 , , . . . , 2 n ) and we have used the cluster decomposition property and the condition r glyph[greatermuch] 1 /m in the second line and used the orthogonal conditions

in (3.2) in the last line. From (3.6), the number of the nonzero terms in the the expansion in (3.4) is 2 n . From (3.4) and (3.6), we obtain

$$& \cdots ^ { \prime } \cdots ^ { \prime } \cdots ^ { \prime } \\ & T r p _ { A } ^ { n } = N ^ { 2 n } [ ( 0 | O _ { j B } ^ { \dagger } O _ { j B } | 0 ) ^ { n } \, ( 0 ^ { ( n ) } | O _ { i A } ^ { ( 1 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i A } ^ { ( 1 ) } \dots O _ { i A } ^ { ( n ) } | 0 ^ { ( n ) } \rangle \\ & + ( 0 | O _ { j B } ^ { \dagger } O _ { j B } | 0 ) ^ { n - 1 } \, ( 0 | O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } | 0 ) \, ( ( 0 ^ { ( n ) } | O _ { i ^ { \prime } A } ^ { ( 1 ) \dagger } O _ { i A } ^ { ( 2 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i ^ { \prime } A } ^ { ( 1 ) } O _ { i A } ^ { ( 2 ) } \dots O _ { i A } ^ { ( n ) } | 0 ^ { ( n ) } \rangle + \dots ) \\ & + ( 0 | O _ { j B } ^ { \dagger } O _ { j B } | 0 ) ^ { n - 2 } \, ( 0 | O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } | 0 ) ^ { 2 } \, ( ( 0 ^ { ( n ) } | O _ { i ^ { \prime } A } ^ { ( 1 ) \dagger } O _ { i ^ { \prime } A } ^ { ( 2 ) \dagger } O _ { i A } ^ { ( 3 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i ^ { \prime } A } ^ { ( 1 ) } O _ { i ^ { \prime } A } ^ { ( 2 ) } O _ { i A } ^ { ( 3 ) } \dots O _ { i A } ^ { ( n ) } | 0 ^ { ( n ) } \rangle + \dots ) \\ & + \dots + ( 0 | O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } | 0 ) ^ { n } \, ( 0 ^ { ( n ) } | O _ { i ^ { \prime } A } ^ { ( 1 ) \dagger } \dots O _ { i ^ { \prime } A } ^ { ( n ) \dagger } E _ { A } O _ { i ^ { \prime } A } ^ { ( 1 ) } \dots O _ { i ^ { \prime } A } ^ { ( n ) } | 0 ^ { ( n ) } \rangle ]. \\ & \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \ \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \.$$

The number of the terms which are proportional to 〈 0 | O † jB O jB | 0 〉 n -l 〈 0 | O † j B ′ O j B ′ | 0 〉 l in (3.7) is n C l . By using (2.20), we obtain

$$( 3. 7 ) & \text{ is } n C _ { l }. \text{ By using } ( 2. 2 0 ), \, \text{ we obtain} \\ & \langle 0 ^ { ( n ) } | O _ { i _ { 1 } A } ^ { ( 1 ) \dagger } \dots O _ { i _ { n } A } ^ { ( n ) \dagger } E _ { A } O _ { i _ { 1 } A } ^ { ( 1 ) } \dots O _ { i _ { n } A } ^ { ( n ) } | 0 ^ { ( n ) } \rangle \\ & = [ ( 0 ^ { ( n ) } | O _ { i _ { 2 } A } ^ { ( 1 ) \dagger } \dots O _ { i _ { n } A } ^ { ( n - 1 ) \dagger } O _ { i _ { 1 } A } ^ { ( n ) \dagger } E _ { A ^ { c } } O _ { i _ { 1 } A } ^ { ( 1 ) } \dots O _ { i _ { n } A } ^ { ( n ) } | 0 ^ { ( n ) } \rangle ] ^ { * } \\ & \simeq [ ( 0 ^ { ( n ) } | E _ { A ^ { c } } \left | 0 ^ { ( n ) } \right \rangle \left \langle 0 ^ { ( n ) } | O _ { i _ { 2 } A } ^ { ( 1 ) \dagger } \dots O _ { i _ { n } A } ^ { ( n - 1 ) \dagger } O _ { i _ { 1 } A } ^ { ( n ) \dagger } O _ { i _ { 1 } A } ^ { ( 1 ) } \dots O _ { i _ { n } A } ^ { ( n ) } | 0 ^ { ( n ) } \right \rangle ] ^ { * } \\ & = \langle 0 ^ { ( n ) } | E _ { A } \left | 0 ^ { ( n ) } \right \rangle \prod _ { l = 1 } ^ { n } \langle 0 | O _ { i _ { l } A } ^ { \dagger } O _ { i _ { l + 1 } A } \left | 0 \right \rangle \\ & = \begin{cases} \text{Tr} \rho _ { 0 A } ^ { n } \left ( 0 | O _ { i _ { l } A } ^ { \dagger } O _ { i _ { l } A } \left | 0 \right \rangle ^ { n } & \text{ for } i = i _ { 1 } = i _ { 2 } = \dots = i _ { n } \\ \text{Tr} \rho _ { 0 A } ^ { n } \left \langle 0 | O _ { i _ { l } ^ { \prime } A } ^ { \dagger } O _ { i ^ { \prime } A } \left | 0 \right \rangle ^ { n } & \text{ for } i ^ { \prime } = i _ { 1 } = i _ { 2 } = \dots = i _ { n } \\ 0 & \text{ otherwise} \end{cases}$$

where i n +1 = i 1 and we have used the cluster decomposition property and the condition R A glyph[greatermuch] 1 /m in the third line and used the orthogonal conditions in (3.2) in the last line. From (3.8) the nonzero terms in (3.7) are only the first term and the last term and we obtain

$$\text{Tr} \rho _ { A } ^ { n } = \text{Tr} \rho _ { 0 A } ^ { n } \cdot N ^ { 2 n } [ \langle 0 | O _ { i A } ^ { \dagger } O _ { i A } \, | 0 \rangle ^ { n } \, \langle 0 | O _ { j B } ^ { \dagger } O _ { j B } \, | 0 \rangle ^ { n } + \langle 0 | O _ { i ^ { \prime } A } ^ { \dagger } O _ { i ^ { \prime } A } \, | 0 \rangle ^ { n } \, \langle 0 | O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } \, | 0 \rangle ^ { n } ].$$

By the same way, we obtain Tr ρ n B as

$$\text{Tr} \rho _ { B } ^ { n } = \text{Tr} \rho _ { 0 B } ^ { n } \cdot N ^ { 2 n } [ \langle 0 | O _ { i A } ^ { \dagger } O _ { i A } \, | 0 \rangle ^ { n } \, \langle 0 | O _ { j B } ^ { \dagger } O _ { j B } \, | 0 \rangle ^ { n } + \langle 0 | O _ { i ^ { \prime } A } ^ { \dagger } O _ { i ^ { \prime } A } \, | 0 \rangle ^ { n } \, \langle 0 | O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } \, | 0 \rangle ^ { n } ].$$

From (3.3), (3.5), (3.9) and (3.10), we obtain the mutual Rényi information I ( n ) ( A,B ) = ( n -1) -1 ln Tr ρ n A B n ∪ Tr ρ n A Tr ρ B as

$$I ^ { ( n ) } ( A, B ) = \frac { 2 } { n - 1 } \ln \frac { ( x + y ) ^ { n } } { x ^ { n } + y ^ { n } }$$

where

$$x \equiv \langle 0 | \, O _ { i A } ^ { \dagger } O _ { i A } \, | 0 \rangle \, \langle 0 | \, O _ { j B } ^ { \dagger } O _ { j B } \, | 0 \rangle \,, \ \ y \equiv \langle 0 | \, O _ { i ^ { \prime } A } ^ { \dagger } O _ { i ^ { \prime } A } \, | 0 \rangle \, \langle 0 | \, O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } \, | 0 \rangle \,.$$

Taking the limit n → 1 leads to the the mutual information

$$I ( A, B ) = 2 [ \ln ( x + y ) - \frac { 1 } { x + y } ( x \ln x + y \ln y ) ].$$

We can reproduce these results from the quantum mechanics. Let us consider the following state,

$$\left | \Psi \right \rangle _ { q m } = N ( \left | i \right \rangle _ { A } \left | j \right \rangle _ { B } + \left | i ^ { \prime } \right \rangle _ { A } \left | j ^ { \prime } \right \rangle _ { B } ),$$

where N is the real normalization constant, and | i ( i ′ ) 〉 A and | j ( j ′ ) 〉 B are the pure state of the subsystem A and B and

$$\langle i | i ^ { \prime } \rangle _ { A } = \langle j | j ^ { \prime } \rangle _ { B } = 0.$$

From (3.14), (3.15) and the normalization condition 〈 Ψ Ψ | 〉 qm = 1 , we obtain

$$\text{Tr} \rho _ { q m A } ^ { n } = \text{Tr} \rho _ { q m B } ^ { n } = N ^ { 2 n } [ \langle i | i \rangle _ { A } ^ { n } \, \langle j | j \rangle _ { B } ^ { n } + \langle i ^ { \prime } | i ^ { \prime } \rangle _ { A } ^ { n } \, \langle j ^ { \prime } | j ^ { \prime } \rangle _ { B } ^ { n } ]$$

where

$$N ^ { - 2 } = \langle i | i \rangle _ { A } \, \langle j | j \rangle _ { B } + \langle i ^ { \prime } | i ^ { \prime } \rangle _ { A } \, \langle j ^ { \prime } | j ^ { \prime } \rangle _ { B } \,.$$

Because the total system is a pure state, Tr ρ n qmA B ∪ = 1 . From (3.16) and (3.17) we obtain the mutual Rényi information as

$$I _ { q m } ^ { ( n ) } ( A, B ) = \frac { 2 } { n - 1 } \ln \frac { ( x _ { q m } + y _ { q m } ) ^ { n } } { x _ { q m } ^ { n } + y _ { q m } ^ { n } },$$

where

$$x _ { q m } \equiv \langle i | i \rangle _ { A } \, \langle j | j \rangle _ { B } \,, \ y _ { q m } \equiv \langle i ^ { \prime } | i ^ { \prime } \rangle _ { A } \, \langle j ^ { \prime } | j ^ { \prime } \rangle _ { B } \,.$$

(3.11) is the same as (3.18) when we replace the states as follows:

$$O _ { i ( i ^ { \prime } ) A } \left | 0 \right \rangle \rightarrow \left | i ( i ^ { \prime } ) \right \rangle _ { A }, \ \ O _ { j ( j ^ { \prime } ) B } \left | 0 \right \rangle \rightarrow \left | j ( j ^ { \prime } ) \right \rangle _ { B }.$$

Interestingly, the mutual information in the QFT measures only the quantum entanglement in the limit r → ∞ although the mutual information measures generally the total of the quantum entanglement and the classical one [22]. So, in this limit, the mutual information is a good measure of quantum entanglement in this sense.

## 4 Free scalar fields

## 4.1 Explicit calculation of the glueing operator E Ω

We consider ( d +1) dimensional free scalar field theory. For free scalar fields, it is useful to represent the glueing operator E Ω in (2.15) as the normal ordered operator. We decompose φ and π into the creation and annihilation parts,

$$\phi ( x ) = \phi ^ { + } ( x ) + \phi ^ { - } ( x ), \quad \pi ( x ) = \pi ^ { + } ( x ) + \pi ^ { - } ( x ),$$

where

$$\lambda & \phi ^ { + } ( x ) = \int \frac { d ^ { d } p } { ( 2 \pi ) ^ { d } } \frac { 1 } { \sqrt { 2 E _ { p } } } a _ { p } e ^ { i p x }, \ \phi ^ { - } ( x ) = ( \phi ^ { + } ( x ) ) ^ { \dagger }, \\ & \pi ^ { + } ( x ) = \int \frac { d ^ { d } p } { ( 2 \pi ) ^ { d } } ( - i ) \sqrt { \frac { E _ { p } } { 2 } } a _ { p } e ^ { i p x }, \ \pi ^ { - } ( x ) = ( \pi ^ { + } ( x ) ) ^ { \dagger }, \\ &.$$

here E p is the energy and [ a , a p † p ′ ] = (2 π ) d δ d ( p -p ′ ) . The commutators of these operators are

$$& \text{are} \\ & \left [ \phi ^ { + } ( x ), \phi ^ { - } ( y ) \right ] = \langle 0 | \, \phi ( x ) \phi ( y ) \left | 0 \rangle = \int \frac { d ^ { d } p } { ( 2 \pi ) ^ { d } } \frac { 1 } { 2 E _ { p } } e ^ { i p ( x - y ) } \equiv \frac { 1 } { 2 } W ^ { - 1 } ( x - y ), \\ & \left [ \pi ^ { + } ( x ), \pi ^ { - } ( y ) \right ] = \langle 0 | \, \pi ( x ) \pi ( y ) \left | 0 \rangle = \int \frac { d ^ { d } p } { ( 2 \pi ) ^ { d } } \frac { E _ { p } } { 2 } e ^ { i p ( x - y ) } \equiv \frac { 1 } { 2 } W ( x - y ), \\ & \left [ \pi ^ { + } ( x ), \phi ^ { - } ( y ) \right ] = \left [ \pi ^ { - } ( x ), \phi ^ { + } ( y ) \right ] = - \frac { i } { 2 } \delta ^ { d } ( x - y ), \\ & \text{where we have defined the matrix $W$ which has continuous indices $x$} \, x \, u \, \text{in} \, ( 4. 3 ) \text{ and } W ^ { - 1 }$$

where we have defined the matrix W which has continuous indices x, y in (4.3) and W -1 is the inverse of W . By using (4.3) and the Baker-Campbell-Hausdorff (BCH) formula e X Y e = e [ X,Y ] e Y e X , e X + Y = e -1 2 [ X,Y ] e X Y e , for [[ X,Y ,X ] ] = [[ X,Y ,Y ] ] = 0 , we obtain

$$e ^ { \lambda } & e ^ { r } = e ^ { | \lambda, r | } e ^ { r } e ^ { \lambda }, \ e ^ { \lambda + r } = e ^ { - \bar { 2 } | \Lambda, r | } e ^ { \lambda } e ^ { r } \,, \text{ for } [ [ X, Y ], X ] = [ [ X, Y ], Y ] = 0, \, \text{ we obtain} \\ & \exp [ i \int d ^ { d } x J ^ { \prime } \phi ] \exp [ i \int d ^ { d } x K \pi ] \exp [ - i \int d ^ { d } x J \phi ] \\ & = \colon \exp [ i \int d ^ { d } x ( K \pi + ( J ^ { \prime } - J ) \phi ) ] \, \colon \\ & \times \exp [ \int d ^ { d } x d ^ { d } y ( - \frac { 1 } { 4 } K ( x ) W ( x - y ) K ( y ) - \frac { 1 } { 4 } ( J ^ { \prime } - J ) ( x ) W ^ { - 1 } ( x - y ) ( J ^ { \prime } - J ) ( y ) ) \\ & - \int d ^ { d } x \frac { i } { 2 } K ( x ) ( J ^ { \prime } + J ) ( x ) ], \\ & \text{where } \colon O \, \colon \text{ means the normal ordered operator of } O. \, \text{ From } ( 4. 4 ) \, \text{ we can rewrite } E _ { \Omega } \text{ in}$$

where : O : means the normal ordered operator of O . From (4.4) we can rewrite E Ω in (2.15) as the normal ordered operator,

$$E _ { \Omega } = \int \prod _ { j = 1 } ^ { n } \prod _ { x \in \Omega } D J ^ { ( j ) } ( x ) D K ^ { ( j ) } ( x ) \, \colon \exp [ i \sum _ { l = 1 } ^ { n } \int d ^ { d } x ( ( J ^ { ( l + 1 ) } - J ^ { ( l ) } ) \phi ^ { ( l ) } + K ^ { ( l ) } \pi ^ { ( l ) } ) ] \, \colon \exp [ - \tilde { S } ],$$

where J ( n +1) = J (1) and

$$& \text{where $\varrho \cdot \cdot \cdot -\varrho \cdot u$} \\ & \tilde { S } \equiv \sum _ { l = 1 } ^ { n } [ \int d ^ { d } x d ^ { d } y [ \frac { 1 } { 4 } K ^ { ( l ) } ( x ) W ( x - y ) K ^ { ( l ) } ( y ) + \frac { 1 } { 4 } ( J ^ { ( l + 1 ) } - J ^ { ( l ) } ) ( x ) W ^ { - 1 } ( x - y ) ( J ^ { ( l + 1 ) } - J ^ { ( l ) } ) ( y ) ] \\ & + \frac { i } { 2 } \int d ^ { d } x K ^ { ( l ) } ( x ) ( J ^ { ( l + 1 ) } + J ^ { ( l ) } ) ( x ) ].$$

(4.6)

For the vacuum state ρ 0 = 0 | 〉 〈 0 | , from (4.5) we obtain

$$\text{Tr} \rho _ { 0 \Omega } ^ { n } = \langle 0 ^ { ( n ) } | \, E _ { \Omega } \, | 0 ^ { ( n ) } \rangle = \int \prod _ { j = 1 } ^ { n } \prod _ { x \in \Omega } D J ^ { ( j ) } ( x ) D K ^ { ( j ) } ( x ) \exp [ - \tilde { S } ]. \text{ \quad \ \ } ( 4. 7 )$$

We can show that (4.7) reproduces the same result as that of the real time approach [1, 2, 23] which is based on the wave functional calculation (see Appendix A). By expanding the exponential in the normal ordered product in (4.5) and performing the Gauss integral of J and K , we can rewrite the E Ω as a series of operators. Note that the odd powers of the expansion of (4.5) vanish from the symmetric property (2.16).

## 4.1.1 The case n = 2

In the case n = 2 , it is useful to define the following linear combinations,

$$\phi _ { \pm } & = \frac { 1 } { \sqrt { 2 } } ( \phi ^ { ( 1 ) } \pm \phi ^ { ( 2 ) } ), \ \pi _ { \pm } = \frac { 1 } { \sqrt { 2 } } ( \pi ^ { ( 1 ) } \pm \pi ^ { ( 2 ) } ), \ J _ { \pm } = \frac { 1 } { \sqrt { 2 } } ( J ^ { ( 1 ) } \pm J ^ { ( 2 ) } ), \ K _ { \pm } = \frac { 1 } { \sqrt { 2 } } ( K ^ { ( 1 ) } \pm K ^ { ( 2 ) } ). \\ - & \quad \dots \quad \dots \quad \dots$$

From (4.8) we rewrite (4.5) as

$$E _ { \Omega } = \int D J _ { + } D J _ { - } D K _ { + } D K _ { - } \, \colon \exp [ i \int d ^ { d } x ( - 2 J _ { - } \phi _ { - } + K _ { + } \pi _ { + } + K _ { - } \pi _ { - } ) ] \, \colon$$

$$& \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ & E _ { \Omega } = \int D J _ { + } D J _ { - } D K _ { + } D K _ { - } \, \colon \exp [ i \int d ^ { d } x ( - 2 J _ { - } \phi _ { - } \, + \, K _ { + } \pi _ { + } \, + \, K _ { - } \pi _ { - } ) ] \, \colon \\ & \quad \times \exp [ \int d ^ { d } x d ^ { d } y ( - \frac { 1 } { 4 } K _ { + } ( x ) W ( x - y ) K _ { + } ( y ) \, - \, J _ { - } ( x ) W ^ { - 1 } ( x - y ) J _ { - } ( y ) ) ], \\ & = \int D J _ { - } D K _ { - } \, \colon \exp [ i \int d ^ { d } x ( - 2 J _ { - } \phi _ { - } \, + \, K _ { - } \pi _ { - } ) ] \, \colon \\ & \quad \times \exp [ \int d ^ { d } x d ^ { d } y ( - \frac { 1 } { 4 } K _ { - } ( x ) W ( x - y ) K _ { - } ( y ) \, - \, J _ { - } ( x ) W ^ { - 1 } ( x - y ) J _ { - } ( y ) ) ],$$

$$( 4. 9 )$$

where we have performed J + and K + integrals. In order to represent the Gauss integrals of K -and J -, we will use the following matrix notation,

$$W ( x - y ) = \begin{pmatrix} W ( x _ { \Omega } - y _ { \Omega } ) & W ( x _ { \Omega } - y _ { \Omega ^ { c } } ) \\ W ( x _ { \Omega ^ { c } } - y _ { \Omega } ) & W ( x _ { \Omega ^ { c } } - y _ { \Omega ^ { c } } ) \end{pmatrix} \equiv \begin{pmatrix} A & B \\ B ^ { T } & C \end{pmatrix}$$

$$W ^ { - 1 } ( x - y ) = \begin{pmatrix} \overset { \cdot } { W ^ { - 1 } ( x _ { \Omega } - y _ { \Omega } ) } & W ^ { - 1 } ( x _ { \Omega } - y _ { \Omega ^ { c } } ) \\ W ^ { - 1 } ( x _ { \Omega ^ { c } } - y _ { \Omega } ) & W ^ { - 1 } ( x _ { \Omega ^ { c } } - y _ { \Omega ^ { c } } ) \end{pmatrix} \equiv \begin{pmatrix} D & E \\ E ^ { T } & F \end{pmatrix}$$

where x Ω(Ω ) c and y Ω(Ω ) c are the coordinates in Ω(Ω ) c . Thus the propagators of J -and K -are

$$< J _ { - } ( x ) J _ { - } ( y ) > \equiv \frac { \int D J _ { - } J _ { - } ( x ) J _ { - } ( y ) e ^ { - \int d ^ { d } x d ^ { d } y J _ { - } ( x ) W ^ { - 1 } ( x - y ) J _ { - } ( y ) } } { \int D J _ { - } e ^ { - \int d ^ { d } x d ^ { d } y J _ { - } ( x ) W ^ { - 1 } ( x - y ) J _ { - } ( y ) } } = \frac { 1 } { 2 } D ^ { - 1 } ( x - y ) \ ( 4. 1 2 )$$

and

$$< K _ { - } ( x ) K _ { - } ( y ) > \equiv \frac { \int D K _ { - } K _ { - } ( x ) K _ { - } ( y ) e ^ { - \int d ^ { d } x d ^ { d } y \frac { 1 } { 4 } K _ { - } ( x ) W ( x - y ) K _ { - } ( y ) } } { \int D K _ { - } e ^ { - \int d ^ { d } x d ^ { d } y \frac { 1 } { 4 } K _ { - } ( x ) W ( x - y ) K _ { - } ( y ) } } = 2 A ^ { - 1 } ( x - y ).$$

Then we can expand the E Ω for n = 2 as

$$& E _ { \Omega } = \text{Tr} \rho _ { 0 \Omega } ^ { 2 } [ 1 - 2 \int d ^ { d } x d ^ { d } y < J _ { - } ( x ) J _ { - } ( y ) > \colon \phi _ { - } ( x ) \phi _ { - } ( y ) \, \colon \\ & - \frac { 1 } { 2 } \int d ^ { d } x d ^ { d } y < K _ { - } ( x ) K _ { - } ( y ) > \colon \pi _ { - } ( x ) \pi _ { - } ( y ) \, \colon + \dots ].$$

Next let us apply above results to the mutual Rényi information I (2) ( A,B ) of disjoint compact spatial regions A and B in the vacuum states of the massless free scalar field. We express E Ω as a sum of the local operators at some conventionally chosen points ( r A , r B ) inside A and B . From (4.14), we have

$$E _ { A ( B ) } = \text{Tr} \rho _ { 0 A ( B ) } ^ { 2 } [ 1 - 2 \, \colon \phi _ { - } ^ { 2 } ( r _ { A ( B ) } ) \, \colon C _ { A ( B ) } + \dots ]$$

where

$$C _ { A ( B ) } = \int d ^ { d } x d ^ { d } y < J _ { - } ( x ) J _ { - } ( y ) >.$$

From the local property (2.17), E A B ∪ = E E A B , and (4.15) we obtain

$$\frac { \text{Tr} \rho _ { 0 A \cup B } ^ { 2 } } { \text{Tr} \rho _ { 0 A } ^ { 2 } \text{Tr} \rho _ { 0 B } ^ { 2 } } = \frac { \langle 0 ^ { ( 2 ) } | \, E _ { A } E _ { B } \, | 0 ^ { ( 2 ) } \rangle } { \langle 0 ^ { ( 2 ) } | \, E _ { A } \, | 0 ^ { ( 2 ) } \rangle \, \langle 0 ^ { ( 2 ) } | \, E _ { B } \, | 0 ^ { ( 2 ) } \rangle } \simeq 1 + \frac { 1 } { 2 } C _ { A } C _ { B } ( W ^ { - 1 } ( r ) ) ^ { 2 }, \quad ( 4. 1 7 )$$

where we have used 〈 0 (2) | : φ 2 -( x ) :: φ 2 -( y ) : | 0 (2) 〉 = 2( 0 〈 | φ x φ y ( ) ( ) | 0 ) 〉 2 = 1 2 ( W -1 ( x -y )) 2 . We use the notation E p = | p | and

$$W ^ { - 1 } ( x - y ) = \int \frac { d ^ { d } p } { ( 2 \pi ) ^ { d } } \frac { 1 } { | p | } e ^ { i p ( x - y ) } = \frac { B _ { d } } { | x - y | ^ { d - 1 } },$$

where B d is a numerical constant. Thus the the mutual Rényi information I (2) ( A,B ) is

$$I ^ { ( 2 ) } ( A, B ) \simeq \frac { 1 } { 2 } C _ { A } C _ { B } \frac { B _ { d } ^ { 2 } } { r ^ { 2 d - 2 } }.$$

We can compute C A B ( ) numerically at least. The expression (4.16) of C A B ( ) is useful for numerical computation. Furthermore, we can obtain the alternative expression of C A B ( ) as follows. Let us consider the following generating function,

$$G ( \sigma ) = \, \int D J _ { - } \exp [ - \, \int d ^ { d } x d ^ { d } y J _ { - } ( x ) W ^ { - 1 } ( x - y ) J _ { - } ( y ) + i \sigma \, \int d ^ { d } x J _ { - } ( x ) ], \quad ( 4. 2 0 )$$

where σ is a c number. We can obtain the coefficients of : φ 2 m -( r A ) : in the expansion of E A from G σ ( ) . For example, we have

$$C _ { A } = \frac { - 1 } { G ( 0 ) } \ \frac { \partial ^ { 2 } G } { \partial \sigma ^ { 2 } } \Big | _ { \sigma = 0 } \,.$$

We can rewrite G σ ( ) by using the ( d +1) dim Euclidean path integral as

$$G ( \sigma ) = \frac { 1 } { Z ( 0 ) } \int D J _ { - } Z ( J _ { - } ) \exp [ i \sigma \int d ^ { d } x J _ { - } ( x ) ],$$

where

$$Z ( J _ { - } ) = \int D \varphi \exp [ - \frac { 1 } { 2 } \int d ^ { d + 1 } x ( \partial _ { \mu } \varphi ) ^ { 2 } - 2 i \int d ^ { d } x J _ { - } \varphi | _ { \tau = 0 } ], \quad \ ( 4. 2 3 )$$

where τ is the Euclidean time coordinate and ϕ τ, x ( ) is a scalar field in ( d +1) dimensional Euclidean space. First, we perform the J -integral and obtain

$$G ( \sigma ) = \frac { 1 } { Z ( 0 ) } \int D \varphi \exp [ - \frac { 1 } { 2 } \int d ^ { d + 1 } x ( \partial _ { \mu } \varphi ) ^ { 2 } ] \prod _ { x \in \Omega } \delta ( \sigma - 2 \varphi ( \tau = 0, x ) ).$$

We decompose ϕ into the quantum part ϕ q and the classical part ϕ cl . The ϕ cl is the solution of ∂ ϕ 2 = 0 and satisfy the boundary condition,

$$\varphi _ { c l } ( \tau = 0, x ) = \frac { \sigma } { 2 } \ \ ( x \in A ), \text{ and } \ \varphi _ { c l } ( \tau, x ) = 0 \ \ ( | \tau ^ { 2 } + x ^ { 2 } | \to \infty ). \quad ( 4. 2 5 )$$

Thus the region ( τ = 0 , x ∈ A ) acts like a conductor where electrostatic potential is σ/ 2 and ϕ cl ( τ, x ) is electrostatic potential at ( τ, x ) . We perform the ϕ q integral and obtain

$$G ( \sigma ) = G ( 0 ) \exp [ - \frac { 1 } { 2 } \int d ^ { d + 1 } x ( \partial _ { \mu } \varphi _ { d } ) ^ { 2 } ] = G ( 0 ) \exp [ - \frac { \sigma ^ { 2 } } { 8 } \mathbf C _ { A } ],$$

where 1 2 ∫ d d +1 x ∂ ( µ ϕ cl ) 2 is the electrostatic energy and we have rewritten it by using C A which is electrostatic capacitance of the conductor ( τ = 0 , x ∈ A ) . From (4.21) and (4.26) we have

$$C _ { A } = \frac { 1 } { 4 } \mathbf C _ { A }.$$

Thus, from (4.19) and (4.27), we reproduced the result that I (2) ( A,B ) in the vacuum state is proportional to the product of the electrostatic capacitance of each regions obtained by Cardy [8].

## 4.2 The mutual information of locally excited states

We consider the mutual Rényi information I ( n ) ( A,B ) = S ( n ) A + S ( n ) B -S ( n ) A B ∪ of disjoint compact spatial regions A and B in the locally excited states of the (d+1) dimensional free massless scalar field theory in the limit when the separation r between A and B is much greater than their sizes R A,B .

We consider the following locally excited state,

$$\left | \Psi \right \rangle = N ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) \left | 0 \right \rangle,$$

where N is a real normalization constant and O iA and O i A ′ ( O jB and O j B ′ ) are operators on A B ( ) and i and i ′ ( j and j ′ ) label a kind of operators. We impose following orthogonal conditions,

$$\left \langle 0 \right | O _ { i A } ^ { \dagger } O _ { i ^ { \prime } A } \left | 0 \right \rangle = \left \langle 0 \right | O _ { j B } ^ { \dagger } O _ { j ^ { \prime } B } \left | 0 \right \rangle = 0.$$

This state is similar to the EPR state. We compute the mutual Rényi information of this state. When r →∞ , the mutual information I ( n ) ( A,B ) is the same as that of the general

QFT which have a mass gap in (3.11). In the mass gap case there is a correction which is O e ( -mR A B ( ) ) for (3.11). In the free massless case there is a correction which is O (1 /R A B ( ) ) .

We consider the leading term of I ( n ) ( A,B ) which depends on r for r glyph[greatermuch] R ,R A B . Because the symmetric property (2.16) is important to determine the r dependence of I ( n ) ( A,B ) as we will show it later, from now on, we impose the condition that under the sign changing transformation ( φ, π ) → -( φ, -π ) the operators O in (4.28) is transformed as

$$O \rightarrow ( - 1 ) ^ { | O | } O,$$

where | O | = 0 or 1 . From the normalization condition 〈 Ψ Ψ | 〉 = 1 we have

$$N ^ { - 2 } = \langle 0 | \left [ O _ { i A } ^ { \dagger } O _ { i A } O _ { j B } ^ { \dagger } O _ { j B } + O _ { i A } ^ { \dagger } O _ { i ^ { \prime } A } O _ { j B } ^ { \dagger } O _ { j ^ { \prime } B } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { i A } O _ { j ^ { \prime } B } ^ { \dagger } O _ { j B } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { i ^ { \prime } A } O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ^ { \dagger } O _ { j ^ { \prime } B } \right ] \left | 0 \right ).$$

By using (2.9), we obtain Tr ρ n Ω for ρ = Ψ | 〉 〈 Ψ | as

$$& \text{Tr} \rho _ { \Omega } ^ { n } = \langle \Psi ^ { ( n ) } | \, E _ { \Omega } \, | \Psi ^ { ( n ) } \rangle \\ & = N ^ { 2 n } \, \langle 0 ^ { ( n ) } | \, ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { j ^ { \prime } B } ^ { \dagger } ) ^ { ( 1 ) } \dots ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } + O _ { i ^ { \prime } A } ^ { \dagger } O _ { j ^ { \prime } B } ^ { \dagger } ) ^ { ( n ) } \\ & \times E _ { \Omega } ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) ^ { ( 1 ) } \dots ( O _ { i A } O _ { j B } + O _ { i ^ { \prime } A } O _ { j ^ { \prime } B } ) ^ { ( n ) } \, | 0 ^ { ( n ) } \rangle \,,$$

where O ( ) l are operators in the Hilbert space of the l -th copy and the Ω is A,B or A ∪ B .

## 4.2.1 The case O i A ′ = O j B ′ = 0

First we consider the case O i A ′ = O j B ′ = 0 and O iA and O jB are nonzero

$$\left | \Psi _ { 1 } \right \rangle = N _ { 1 } O _ { i A } O _ { j B } \left | 0 \right \rangle.$$

From the normalization condition we have

$$N _ { 1 } ^ { - 2 } = \langle 0 | \, O _ { i A } ^ { \dagger } O _ { i A } O _ { j B } ^ { \dagger } O _ { j B } \, | 0 \rangle \,.$$

From the condition (4.30), | O † iA O iA | = 0 and the OPE of O † iA O iA has the form

$$O _ { i A } ^ { \dagger } O _ { i A } = \langle 0 | \, O _ { i A } ^ { \dagger } O _ { i A } \, | 0 \rangle \, [ 1 + C _ { i A } ^ { i \phi ^ { 2 } \colon } \colon \phi ^ { 2 } \colon ( r _ { A } ) + \dots ],$$

where ( r A , r B ) are some conventionally chosen points inside A and B . The key point is that there is not φ in the OPE (4.35) because of | O † iA O iA | = 0 . The OPE of O † jB O jB has the same form as that of O † iA O iA because of | O † jB O jB | = 0 . Thus we have

$$N _ { 1 } ^ { - 2 } = \langle 0 | \, O _ { i A } ^ { \dagger } O _ { i A } \, | 0 \rangle \, \langle 0 | \, O _ { j B } ^ { \dagger } O _ { j B } \, | 0 \rangle \, [ 1 + O ( 1 / r ^ { 2 d - 2 } ) ],$$

where we have used 〈 0 | : φ 2 ( r A ) :: φ 2 ( r B ) : | 0 〉 ∝ 1 /r 2 d -2 .

We consider Tr ρ n A . From (4.32) we have

$$& \text{Tr} \rho _ { A } ^ { n } = \langle \Psi _ { 1 } ^ { ( n ) } | \, E _ { \Omega } \, | \Psi _ { 1 } ^ { ( n ) } \rangle \\ & = N _ { 1 } ^ { 2 n } \, \langle 0 ^ { ( n ) } | \, ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } ) ^ { ( 1 ) } \dots ( O _ { i A } ^ { \dagger } O _ { j B } ^ { \dagger } ) ^ { ( n ) } E _ { A } ( O _ { i A } O _ { j B } ) ^ { ( 1 ) } \dots ( O _ { i A } O _ { j B } ) ^ { ( n ) } \, | 0 ^ { ( n ) } \rangle \quad ( 4. 3 7 ) \\ & = N _ { 1 } ^ { 2 n } \, \langle 0 ^ { ( n ) } | \, O _ { i A } ^ { ( 1 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i A } ^ { ( 1 ) } \dots O _ { i A } ^ { ( n ) } \cdot O _ { j B } ^ { ( 1 ) \dagger } \dots O _ { j B } ^ { ( n ) \dagger } O _ { j B } ^ { ( 1 ) } \dots O _ { j B } ^ { ( n ) } \, | 0 ^ { ( n ) } \rangle \,.$$

From (2.16) and (4.30), O (1) † iA . . . O ( n ) † iA E O A (1) iA . . . O ( n ) iA and O (1) † jB . . . O ( n ) † jB O (1) jB . . . O ( n ) jB are even under the sign changing transformation

$$( \phi ^ { ( 1 ) }, \dots, \phi ^ { ( n ) }, \pi ^ { ( 1 ) }, \dots, \pi ^ { ( n ) } ) \to ( - \phi ^ { ( 1 ) }, \dots, - \phi ^ { ( n ) }, - \pi ^ { ( 1 ) }, \dots, - \pi ^ { ( n ) } ).$$

So the OPEs of these operators are

$$O _ { i A } ^ { ( 1 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i A } ^ { ( 1 ) } \dots O _ { i A } ^ { ( n ) } = \langle 0 | O _ { i A } ^ { ( 1 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i A } ^ { ( 1 ) } \dots O _ { i A } ^ { ( n ) } \, | 0 \rangle \, [ 1 + \sum _ { l, m } \tilde { C } _ { i A } ^ { \phi ( l ) \phi ( m ) } \colon \phi ^ { ( l ) } \, \colon \langle r _ { A } \rangle + \dots ]$$

and

$$O _ { j B } ^ { ( 1 ) \dagger } \dots O _ { j B } ^ { ( n ) \dagger } O _ { j B } ^ { ( 1 ) } \dots O _ { j B } ^ { ( n ) } = \langle 0 | O _ { j B } ^ { ( 1 ) \dagger } \dots O _ { j B } ^ { ( n ) \dagger } O _ { j B } ^ { ( 1 ) } \dots O _ { j B } ^ { ( n ) } \, | 0 \rangle \, [ 1 + \sum _ { l } \tilde { D } _ { j B } ^ { \cdot \phi ( l ) \cdot \phi ( l ) } \colon \phi ^ { ( l ) } \phi ^ { ( l ) } \colon ( r _ { B } ) + \dots ].$$

Thus, we substitute (4.36), (4.39) and (4.40) into (4.37) and obtain the power of r of the subleading term of Tr ρ n A

$$\text{Tr} \rho _ { A } ^ { n } & = \text{Tr} \rho _ { A } ^ { n } | _ { r \rightarrow \infty } [ 1 + O ( 1 / r ^ { 2 d - 2 } ) ].$$

In the same way, we can obtain the power of r of the subleading terms of Tr ρ n B and Tr ρ n A B ∪

$$\text{Tr} \rho _ { B } ^ { n } = \text{Tr} \rho _ { B } ^ { n } | _ { r \rightarrow \infty } [ 1 + O ( 1 / r ^ { 2 d - 2 } ) ], \ \ \text{Tr} \rho _ { A \cup B } ^ { n } = \text{Tr} \rho _ { A \cup B } ^ { n } | _ { r \rightarrow \infty } [ 1 + O ( 1 / r ^ { 2 d - 2 } ) ]. \quad ( 4. 4 2 )$$

Thus the mutual Rényi information is

$$I ^ { ( n ) } ( A, B ) = O ( 1 / r ^ { 2 d - 2 } ),$$

where we have used I ( n ) ( A,B ) | r →∞ = 0 for O i A ′ = O j B ′ = 0 from (3.11).

## 4.2.2 The case O iA , O jB , O i A ′ and O j B ′ are nonzero

There are two cases which are different powers of r of the subleading term of I ( n ) ( A,B )

(i)The case | O | = | O ′

. iA i A |

In this case the operators at A in the expansions of N -2 , Tr ρ n Ω (Ω = A,B,A ∪ B ) in (4.31) and (4.32) are even under the sign changing transformation ( 4 38 . ). So the subleading terms of them come from the operators : φ ( ) l φ ( m ) : ( r A ) . Thus the mutual Rényi information is

$$I ^ { ( n ) } ( A, B ) = I ^ { ( n ) } ( A, B ) | _ { r \to \infty } + O ( 1 / r ^ { 2 d - 2 } ),$$

where I ( n ) ( A,B ) | r →∞ is the same form as (3.11).

glyph[negationslash]

(ii)The case | O iA | = | O i A ′ | and | O jB | = | O j B ′ |

glyph[negationslash]

There are some terms which are odd under the sign changing transformation ( 4 38 . ). As an example, in the expansion of Tr ρ n A B ∪ in (4.32), we consider the following term,

$$\langle 0 ^ { ( n ) } | \, O _ { i ^ { \prime } A } ^ { ( 1 ) \dagger } O _ { i A } ^ { ( 2 ) \dagger } \dots O _ { i A } ^ { ( n ) \dagger } E _ { A } O _ { i A } ^ { ( 1 ) } O _ { i A } ^ { ( 2 ) } \dots O _ { i A } ^ { ( n ) } \cdot O _ { j ^ { \prime } B } ^ { ( 1 ) \dagger } O _ { j B } ^ { ( 2 ) \dagger } \dots O _ { j B } ^ { ( n ) \dagger } E _ { B } O _ { j B } ^ { ( 1 ) } O _ { j B } ^ { ( 2 ) } \dots O _ { j B } ^ { ( n ) } \, | 0 ^ { ( n ) } \, \rangle.$$

In (4.45), the operators at A and at B are both odd under the sign changing transformation ( 4 38 . ). So the leading r dependent term of it come from the operators φ ( ) l . Thus the mutual Rényi information is

$$I ^ { ( n ) } ( A, B ) = I ^ { ( n ) } ( A, B ) | _ { r \to \infty } + O ( 1 / r ^ { d - 1 } ),$$

where I ( n ) ( A,B ) | r →∞ is the same form as (3.11).

## 5 Conclusion and discussions

We developed the computational method of Rényi entanglement entropy based on the idea that Tr ρ n Ω is written as the expectation value of the local operator at Ω . We expressed Tr ρ n Ω as the expectation value of the glueing operator E Ω , Tr ρ n Ω = Tr( ρ ( n ) E Ω ) . We constructed explicitly E Ω and investigated its general properties. For a free scalar field, we rewrote E Ω in (2.15) using the normal ordering. In the case n = 2 , we obtained a simple expression of E Ω and reproduced the result that I (2) ( A,B ) in the vacuum state is proportional to the product of the electrostatic capacitance of each regions obtained by Cardy [8]. The coefficients of the expansion of E Ω is obtained by the propagators of J -and K -in (4.12) and (4.13). We can compute these propagators numerically at least and the expression (4.9) is useful for numerical calculation.

The advantages of this methods are that we can use ordinary technique in QFT such as OPE and the cluster decomposition property and that we can use the general properties and the explicit expression of the glueing operator to compute systematically the Rényi entropy for an arbitrary state.

We applied this method to consider the mutual Rényi information I ( n ) ( A,B ) of disjoint compact spatial regions A and B in the locally excited states defined by acting the local operators at A and B on the vacuum of a ( d + 1) -dimensional field theory, in the limit when the separation r between A and B is much greater than their sizes R A,B . For the general QFT which has a mass gap, we computed I ( n ) ( A,B ) explicitly and find that this result is interpreted in terms of an entangled state in quantum mechanics. Interestingly, the mutual information in the QFT measures only the quantum entanglement in the limit r → ∞ although the mutual information measures generally the total of the quantum entanglement and the classical one [22]. So, in this limit, the mutual information is a good measure of quantum entanglement in this sense. For a free massless scalar field, we showed that for some classes of excited states I ( n ) ( A,B ) -I ( n ) ( A,B ) | r →∞ = C ( n ) AB /r α d ( -1) where α = 1 or 2 which is determined by the property of the local operators under the transformation φ →-φ and α = 2 for the vacuum state.

Finally we discuss the generalization of our method. Although we considered only the locally excited states, we can apply our method to more general excited states, for example, many particle states and thermal states. We might be able to generalize our method to fermionic fields. We could apply our method to perturbative calculation in an interacting field theory.

## Acknowledgments

I am grateful to Tadashi Takayanagi for a careful reading of this manuscript and useful comments and discussions. I also would like to thank Pawel Caputa, John Cardy, Song He, Masahiro Nozaki, Tokiro Numasawa, and Kento Watanabe for useful discussions. This work is supported by JSPS Grant-in-Aid for Scientific Research (B) No.25287058 and JSPS Grant-in-Aid for Challenging Exploratory Research No.24654057.

## A Rényi entropy for the vacuum state in free scalar field theory

In this appendix we show that (4.7) reproduces the same result as that of the real time approach [1, 2, 23] which is based on the wave functional calculation.

From (4.6) and (4.7), we perform the K integral in (4.7) and obtain

$$\text{Trp} _ { 0 \Omega } ^ { n } & = \left ( D e t \left ( \frac { A } { 4 \pi } \right ) \right ) ^ { - n / 2 } \int \prod _ { j = 1 } ^ { n } \prod _ { x \in \Omega } D J ^ { ( j ) } ( x ) \exp [ - \int d ^ { d } x d ^ { d } y ( J ^ { ( 1 ) } ( x ), \cdots, J ^ { ( n ) } ( x ) ) M _ { n } ( x, y ) \left ( \begin{matrix} J ^ { ( 1 ) } ( y ) \\ \vdots \\ J ^ { ( n ) } ( y ) \end{matrix} \right ), \\ \text{throw} \quad &$$

where

$$M _ { n } = \begin{pmatrix} X \ Y & 0 & \cdots & 0 & Y \\ Y \ X & Y & \cdots & 0 & 0 \\ 0 & Y \ X & \cdots & 0 & 0 \\ \vdots & \vdots & \vdots & & \vdots & \vdots \\ 0 & 0 & 0 & \cdots & X \ Y \\ Y & 0 & 0 & \cdots & Y \ X \end{pmatrix},$$

here

$$X = \frac { 1 } { 2 } ( A ^ { - 1 } + D ), \quad Y = \frac { 1 } { 4 } ( A ^ { - 1 } - D )$$

and we have used the matrix notation in (4.10) and (4.11).

The following calculation is analogous to that in [23]. From (A.1) and (A.2) we obtain

$$\text{Tr} \rho _ { 0 \Omega } ^ { n } = \left ( D e t \left ( \frac { A } { 4 \pi } \right ) \right ) ^ { - n / 2 } ( D e t ( 4 \pi M _ { n } ) ) ^ { - 1 / 2 },$$

where we have used the normalization condition of the J integral ∫ DJ ( j ) exp[ -∫ d xd yJ d d ( j ) ( x M x,y J ) ( ) ( j ) ( y )] = ( Det (4 πM )) -1 2 / .

We rewrite M n as

$$M _ { n } & = \frac { X } { 2 } \tilde { M } _ { n }, & ( \text{A.5} )$$

where

$$\tilde { M } _ { n } = \begin{pmatrix} 2 & - Z & 0 & \cdots & 0 & - Z \\ - Z & 2 & - Z & \cdots & 0 & 0 \\ 0 & - Z & 2 & \cdots & 0 & 0 \\ \vdots & \vdots & \vdots & & \vdots & \vdots \\ 0 & 0 & 0 & \cdots & 2 & - Z \\ - Z & 0 & 0 & \cdots & - Z & 2 \end{pmatrix},$$

here

$$Z = - 2 X ^ { - 1 } Y = ( 1 + A D ) ^ { - 1 } ( A D - 1 ).$$

We diagonalize Z and denote the eigenvalues of Z as z i . And we can diagonalize ˜ M n by Fourier transformation and obtain

$$D e t \tilde { M } _ { n } = \prod _ { i } \prod _ { r = 1 } ^ { n } [ 2 - 2 z _ { i } \cos ( \frac { 2 \pi r } { n } ) ] = \prod _ { i } 2 ^ { n } \frac { ( 1 - \xi _ { i } ^ { n } ) ^ { 2 } } { ( 1 + \xi _ { i } ^ { 2 } ) ^ { n } },$$

where ξ i is defined as

$$z _ { i } & = 2 \xi _ { i } / ( \xi _ { i } ^ { 2 } + 1 ). & ( A. 9 )$$

From (A.4) and (A.8) we obtain

$$\text{Tr} \rho _ { 0 \Omega } ^ { n } = \prod _ { i } \frac { ( 1 - \xi _ { i } ) ^ { n } } { ( 1 - \xi _ { i } ^ { n } ) }.$$

Thus we obtain the Rényi entropies S ( n ) 0Ω = (1 -n ) -1 ln Tr ρ n 0Ω and the entanglement entropy S 0Ω = -Tr ρ 0Ω ln ρ 0Ω = -∂ ∂n ln Tr ρ n 0Ω | n =1 as follows:

$$S _ { 0 \Omega } ^ { ( n ) } & = \sum _ { i } ( 1 - n ) ^ { - 1 } [ n \ln ( 1 - \xi _ { i } ) - \ln ( 1 - \xi _ { i } ^ { n } ) ], & ( A. 1 1 )$$

$$S _ { 0 \Omega } = \sum _ { i } [ - \ln ( 1 - \xi _ { i } ) - \frac { \xi _ { i } } { 1 - \xi _ { i } } \ln \xi _ { i } ].$$

In order to show that (A.11) and (A.12) are the same results as those of the real time approach, we rewrite (A.11) and(A.12). We define the matrix

$$\tilde { C } = \frac { 1 } { 2 } ( D A ) ^ { 1 / 2 },$$

and rewrite the matrix Z in (A.7) as

$$Z & = ( 1 + 4 ( \tilde { C } ^ { T } ) ^ { 2 } ) ^ { - 1 } ( 4 ( \tilde { C } ^ { T } ) ^ { 2 } - 1 ). & & ( \text{A.14} )$$

From (A.9) and (A.14) we obtain

$$\tilde { C } _ { i } ^ { T } & = \frac { 1 } { 2 } ( 1 + \xi _ { i } ) ( 1 - \xi _ { i } ) ^ { - 1 }.$$

From (A.15), we rewrite (A.11) and(A.12) as

$$S _ { 0 \Omega } ^ { ( n ) } & = \sum _ { i } ( n - 1 ) ^ { - 1 } \ln [ ( \tilde { C } _ { i } ^ { T } + 1 / 2 ) ^ { n } - ( \tilde { C } _ { i } ^ { T } - 1 / 2 ) ^ { n } ] \\ & = ( n - 1 ) ^ { - 1 } \text{tr} \ln [ ( \tilde { C } ^ { T } + 1 / 2 ) ^ { n } - ( \tilde { C } ^ { T } - 1 / 2 ) ^ { n } ] \\ & = ( n - 1 ) ^ { - 1 } \text{tr} \ln [ ( \tilde { C } + 1 / 2 ) ^ { n } - ( \tilde { C } - 1 / 2 ) ^ { n } ],$$

$$S _ { 0 \Omega } = \text{tr} [ ( \tilde { C } + 1 / 2 ) \ln ( \tilde { C } + 1 / 2 ) - ( \tilde { C } - 1 / 2 ) \ln ( \tilde { C } - 1 / 2 ) ].$$

(A.16) and (A.17) are the same results as those of the real time approach (see e.g. [12]).

## References

| [1]   | L. Bombelli, R. K. Koul, J. Lee and R. D. Sorkin, 'A Quantum Source of Entropy for Black Holes,' Phys. Rev. D 34 (1986) 373.                                                                                                                                                 |
|-------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [2]   | M. Srednicki, 'Entropy and area,' Phys. Rev. Lett. 71 (1993) 666 [hep-th/9303048].                                                                                                                                                                                           |
| [3]   | J. Eisert, M. Cramer and M. B. Plenio, 'Area laws for the entanglement entropy - a review,' Rev. Mod. Phys. 82 (2010) 277 [arXiv:0808.3773 [quant-ph]].                                                                                                                      |
| [4]   | J.I.Latorre and A.Riera, 'A short review on entanglement in quantum spin systems,' J. Phys. A 42 (2009) 4002 [arXiv:0906.1499 [hep-th]]                                                                                                                                      |
| [5]   | N. Shiba and T. Takayanagi, 'Volume Law for the Entanglement Entropy in Non-local QFTs,' arXiv:1311.1643 [hep-th].                                                                                                                                                           |
| [6]   | J. L. Karczmarek and P. Sabella-Garnier, 'Entanglement entropy on the fuzzy sphere,' arXiv:1310.8345 [hep-th].                                                                                                                                                               |
| [7]   | M. M. Wolf, 'Violation of the entropic area law for fermions', Phys. Rev. Lett. 96 (2006) 010404 [quant-ph/0503219]; D. Gioev, I. Klich, 'Entanglement entropy of fermions in any dimension and the Widom conjecture,' Phys. Rev. Lett. 96 (2006) 100503 [quant-ph/0504151]. |
| [8]   | J. Cardy, 'Some results on the mutual information of disjoint regions in higher dimensions,' J. Phys. A 46 (2013) 285402 [arXiv:1304.7985]                                                                                                                                   |
| [9]   | P. Calabrese, J. Cardy and E. Tonni, 'Entanglement entropy of two disjoint intervals in conformal field theory II,' J. Stat. Mech. 1101 (2011) P01021 [arXiv:1011.5482]                                                                                                      |
| [10]  | M. Headrick, 'Entanglement Renyi entropies in holographic theories,' Phys. Rev. D 82 (2010) 126010 [arXiv:1006.0047]                                                                                                                                                         |
| [11]  | H. Casini and M. Huerta, 'Remarks on the entanglement entropy for disconnected regions,' JHEP 03 (2009) 048 [arXiv:0812.1773]                                                                                                                                                |
| [12]  | H. Casini and M. Huerta, 'Entanglement entropy in free quantum field theory,' J.Phys. A42 (2009) 504007, arXiv:0905.2562 [hep-th].                                                                                                                                           |
| [13]  | N. Shiba, 'Entanglement Entropy of Two Black Holes and Entanglement Entropic Force,' Phys.Rev. D83 (2011) 065002, arXiv:1011.3760 [hep-th].                                                                                                                                  |
| [14]  | N. Shiba, 'Entanglement Entropy of Two Spheres,' JHEP 1207 (2012) 100, arXiv:1201.4865 [hep-th].                                                                                                                                                                             |
| [15]  | H.J.Schnitzer, "Mutual Renyi information for two disjoint compound systems", arXiv:1406.1161[hep-th].                                                                                                                                                                        |
| [16]  | T. Faulkner, A. Lewkowycz and J. Maldacena, 'Quantum corrections to holographic entanglement entropy," JHEP 1311, 074 (2013) [arXiv:1307.2892].                                                                                                                              |
| [17]  | S. Ryu and T. Takayanagi, 'Holographic derivation of entanglement entropy from AdS/CFT,' Phys. Rev. Lett. 96 (2006) 181602; 'Aspects of holographic entanglement entropy,' JHEP 0608 (2006) 045.                                                                             |
| [18]  | M. Nozaki, T. Numasawa and T. Takayanagi, "Quantum Entanglement of Local Operators in Conformal Field Theories," Phys. Rev. Lett. 112, 111602 (2014) [arXiv:1401.0539 [hep-th]].                                                                                             |
| [19]  | M. Nozaki, "Notes on Quantum Entanglement of Local Operators," arXiv:1405.5875 [hep-th].                                                                                                                                                                                     |

- [20] S. He, T. Numasawa, T. Takayanagi and K. Watanabe, "Quantum Dimension as Entanglement Entropy in 2D CFTs," arXiv:1403.0702 [hep-th].
- [21] P. Caputa, M. Nozaki and T. Takayanagi, "Entanglement of Local Operators in large N CFTs," arXiv:1405.5946 [hep-th]
- [22] M. Nielsen and I. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, Cambridge, England (2000), pg. 9.
- [23] C. G. Callan, Jr. and F. Wilczek, "On geometric entropy," Phys. Lett. B333,55 (1994), arXiv:hep-th/9401072 .
- [24] C.P.Herzog, "Universal Thermal Corrections to Entanglement Entropy for Conformal Field Theories on Spheres," arXiv:1407.1358 [hep-th] .
Noburo Shiba
Entanglement Entropy of Disjoint Regions in Excited States : An Operator Method
We develop the computational method of entanglement entropy based on the idea
that $Tr\rho_{\Omega}^n$ is written as the expectation value of the local
operator, where $\rho_{\Omega}$ is a density matrix of the subsystem $\Omega$.
We apply it to consider the mutual Renyi information
$I^{(n)}(A,B)=S^{(n)}_A+S^{(n)}_B-S^{(n)}_{A\cup B}$ of disjoint compact
spatial regions $A$ and $B$ in the locally excited states defined by acting the
local operators at $A$ and $B$ on the vacuum of a $(d+1)$-dimensional field
theory, in the limit when the separation $r$ between $A$ and $B$ is much
greater than their sizes $R_{A,B}$. For the general QFT which has a mass gap,
we compute $I^{(n)}(A,B)$ explicitly and find that this result is interpreted
in terms of an entangled state in quantum mechanics. For a free massless scalar
field, we show that for some classes of excited states,
$I^{(n)}(A,B)-I^{(n)}(A,B)|_{r \rightarrow \infty} =C^{(n)}_{AB}/r^{\alpha
(d-1)}$ where $\alpha=1$ or 2 which is determined by the property of the local
operators under the transformation $\phi \rightarrow -\phi$ and $\alpha=2$ for
the vacuum state. We give a method to compute $C^{(2)}_{AB}$ systematically.