Quantum Entanglement of Locally Excited States in Maxwell Theory

In 4 dimensional Maxwell gauge theory, we study the changes of (Renyi) entangle-ment entropy which are defined by subtracting the entropy for the ground state from the one for the locally excited states generated by acting with the gauge invariant local operators on the state. The changes for the operators which we consider in this paper reflect the electric-magnetic duality. The late-time value of changes can be interpreted in terms of electromagnetic quasi-particles. When the operator constructed of both electric and magnetic fields acts on the ground state, it shows that the operator acts on the late-time structure of quantum entanglement differently from free scalar fields.


Introduction and Summary
Quantum entanglement significantly distinguishes quantum states from classical states. It can characterize conformal field theories [1,2,3] and topological phases [4,5,6]. In Gauge/Gravity correspondence [7,8,9], the structure of quantum entanglement in quantum field theories (QFTs) living on the boundary is expected to be related to the gravity in the bulk [10,11]. There are a lot of works done to reveal how the structure of quantum entanglement on the boundary corresponds to the geometry in the bulk [12,13,14,15,16,17,18]. Therefore it is important to uncover the fundamental features which quantum entanglement possesses. (Rényi) entanglement entropy is one of the useful quantities to investigate them.
However the definition of (Rényi) entanglement entropy in gauge theories has subtleties [19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]. In gauge theories, physical states have to be gauge invariant. It obeys constraints which guarantee its gauge invariance. They make it difficult to divide the Hilbert space into subsystems A and B because the physical degrees of freedom in A depends on the freedom in B due to the constraints. Their boundary is ∂A. Then the definition of (Rényi) entanglement entropy needs the precise method of dividing Hilbert space and defining the reduced density matrix ρ A which is given by tracing out the degrees of freedom in B, ρ A = tr B ρ. (1.1) On the other hand, the entropy in QFTs depends on a UV cutoff (ultraviolet cutoff) δ because by definition it has the UV divergence. It is given by a series expansion in conformal field theories. The physical degrees of freedom around ∂A have the significant effect on the terms which depend on δ. The method of dividing the Hilbert space is expect to affect the degrees of freedom around ∂A in the direct fashion. In the present paper, we study the changes of (Rényi) entanglement entropy ∆S

(n)
A which is defined by subtracting the entropy for the ground state from the one for the locally excited state, which is defined by acting with a local operator on the ground state. Here we assume that the operator is located far from ∂A. We will explain it more in the next section. As in [35,36,37,38,39,40,41], their changes do not possess the UV divergence. More precisely, they measure how the local operator changes the structure of quantum entanglement. Therefore they are expected to avoid the subtleties which (Rényi) entanglement entropy has.
In this paper we study ∆S

(n)
A in 4d Maxwell gauge theory, which is a free CFT [42]. The previous works [37,38,39] show the time evolution of ∆S (n) A can be interpreted in terms of relativistic propagation of entangled quasi-particles which are created by local operators. In the free theories, the late-time value of ∆S (n) A is given by the constant, which depends on the operators. It comes from the quantum entanglement between quasi-particles. As in [39], the late-time entanglement structure depends on the kind of quasi-particles. The authors in [40] show that in the specific 2d CFTs, it is related to the quantum dimension of the operator which acts on the ground state. The Authors in [41,43] have shown that in holographic theories the late time value of ∆S (n) A logarithmically increases similarly to the behavior of entanglement entropy for the local quenches [44,45]. ∆S

(n)
A in the finite temperature system was investigated by the authors in [46]. There are many works done to study the fundamental properties of ∆S A depend on theories and the quasi-particles created by the local operator. Then we study how the structure of quantum entanglement is changed by gauge invariant local operators such as electric and magnetic fields. In particular, we study how the late-time structure of quantum entanglement depends on them. More precisely, we study how quasiparticles have the effect on the structure. We also study whether ∆S

(n)
A for gauge invariant locally excited state reflects electric-magnetic duality.

Summary
Here we briefly summarize our results in this work. We study how they change the structure of quantum entanglement by measuring the time evolution of ∆S (n) A for various gauge invariant local operators. We also study whether ∆S (n) A is invariant the electric-magnetic duality transformation.

Electric-Magnetic Duality
As it will be explained later, ∆S A reflects the electric-magnetic duality, the entropy for B i is equal to that for E i , which is the electric field along the same direction as that of magnetic field. ∆S

(n)
A for the electric and magnetic fields along the direction vertical to the entangling surface increases slower than those for fields along the directions parallel to the surface. However there are no difference between the effects of electromagnetic field and that of scalar one on the entanglement structure at the late time 1

Composite Operators
If the composite operator such as B 2 acts on the ground state, they lead to the late-time structure of quantum entanglement in the same manner as a specific scalar operator. Then the late-time value of ∆S

(n)
A for that can be interpreted in terms of quasi-particles created by the scalar operator. However ∆S

(n)
A for some specific operators (e.g, E 2 B 3 ) constructed of both electric and magnetic fields can be interpreted in terms of not the scalar quasi-particles but electromagnetic one, which is explained in section 4. Here B 3 (E 3 ) and B 2 (E 2 ) are the magnetic (electric) fields along the direction perpendicular to ∂A as we will explain it later.

Late-time Algebra
We interpret the late-time values of ∆S (n) A in terms of electromagnetic quasi-particles created by an electromagnetic field, and derive a late-time algebra which they obey. There are commutation relations between the particles of the same kinds of fields . As we will mention later, there are also additional relations between E 2 (E 3 ) and B 3 (B 2 ), which are parallel to the entangling surface. They make the effect of electromagnetic fields different from that of scalar fields on the late-time structure of quantum entanglement.

Organization
This paper is organized as follows. In section 2, we will explain locally excited states and how to compute ∆S A in terms of entangled quasi-particles in section 4. We study how they have the effect on the late-time structure of quantum entanglement. We finish with the conclusion, future problems and the detail of propagators is included in appendices.

How to compute Excesses of (Rényi) Entanglement Entropy
By measuring the excess of (Rényi) entanglement entropies ∆S

(n)
A , we study how local gaugeinvariant operators changes the structure of quantum entanglement in the 4d Maxwell gauge theory: In this section, we explain the definition of locally excited state and how to compute ∆S (n) A in the replica method.

The Definition of Locally Excited States
The locally excited state is defined by acting with a gauge invariant local operator O such as F µν on the ground state: where N is a normalization constant and |0 is a gauge invariant state. As in Figure.1, O is located at t = −t, x 1 = −l and x = (x 2 , x 3 ).

(n)
A approaches to a constant, which comes form quantum entanglement between entangled quasi-particles. In this paper we would like to study how the constant depends on gauge invariant operators. Therefore the region in Figure.1 is chosen as A.

Excesses of (Rényi) Entanglement Entropy
Here we explain more about the definition of ∆S (n) A . (Rényi) entanglement entropy for the ground state is a static quantity, which does not depend on time. Then we define the excesses of (Rényi) entanglement entropy ∆S A for the ground state from those for locally excited states, ∆S , are (Rényi) entanglement entropies for the excited states in (2.2) and the ground state |0 , respectively. In the sense that ∆S

(n)
A does not depend on δ, it is a "renormalized" (Rényi) entanglement entropy.

The Replica Trick
We would like to study the time evolution of ∆S (n) A in 4d Minkowski spacetime. However in this paper we do not directly study the changes of entanglement structure in the spacetime. Without doing so, we compute ∆S (n) A in Euclidean space by the replica trick. After that we perform the analytic continuation, which we will explain later. Then we compute the real time evolution of ∆S (n) A . As in [37,38,39,40], a reduced density matrix in Euclidean space is given by where τ is Euclidean time. By introducing a polar coordinate, (τ l,e , −l) is mapped to (r 1,2 , θ 1,2 ) as in Figure.2. In the replica trick, (Rényi) entanglement entropies for (2.4) and the ground state are respectively given by 2 where θ k 1,2 = θ 1,2 + 2(k − 1)π. The actions S n and S 1 are defined on n-sheeted geometry Σ n (see Figure.3) and the flat space Σ 1 , respectively. By substituting (Rényi) entanglement We only need to compute propagators on Σ n in order to compute ∆S If we choose a specific gauge 3 , their green functions obey the same equation of motion as that for 4d free massless scalar field theory, The solution of the equation is given by where t 0 is defined by (2.8) has been obtained by the authors in [37,38,39,41,48,49].

Analytic Continuation
After computing green functions on Σ n in Euclidean space, we perform the following analytic continuation, where ǫ is a smearing parameter which is introduced to keep the norm of the excited state finite. Analytic-continued green functions depend on ǫ. We are interested in the behavior of ∆S

Excesses of (Rényi) Entanglement Entropy
In this section, we study the time evolution and the late time value of ∆S A does not depend on the operator. It can be interpreted in terms of the quasi-particle created by a scalar operator φ.
(ii) Composite operators which act on the ground state are constructed of only electric or magnetic fields such as E 2 and B 2 . ∆S

(n)
A for E 2 is equivalent to the entropy for B 2 . Then ∆S

(n)
A for them invariant under the electric-magnetic duality transformation. There are no differences between their effect on the (Rényi) entanglement entropy. Its late-time values can be interpreted in terms of quasi-particles, which are created by the operator constructed of massless free scalar fields 3 a=1 (φ a ) 2 . Here a denotes the kinds of fields. (iii) Local operators are constructed of both electric and magnetic fields such as E 2 A shows that there is a significant difference between the effect of E 1 (B 1 ) and E 2,3 (B 2,3 ) on the late-time entanglement structure. Here E 1 (B 1 ) is the electric (magnetic) field along the direction vertical to the entangling surface. On the other hand, E 2,3 (B 2,3 ) is the electric (magnetic) field along the direction parallel to the entangling surface. As it will be explained in the next section, the difference comes from the commutation relation between electromagnetic quasi-particles created by E 2 (E 3 ) and the particles created by B 3 (B 2 ).
Here locally excited states are defined by acting with only E i or B i on the ground state. ∆S

(n)
A is given by (2.6) in the replica method with Euclidean signature. After performing the analytic continuation in (2.10) and taking the limit ǫ → 0, their time evolution is given as follows. ∆S

(n)
A vanishes before t = l(> 0), but after t = l, they increase. The detail of their time evolution is summarized in Table.1. After taking the late time limit (0 < l ≪ t), they are given by ∆S Their late time value is the same as that for φ in free massless scalar field theories with any spacetime dimensions. It can be interpreted as (Rényi) entanglement entropy for maximally entangled state in 2 qubit system. Therefore they do not show the difference between the effect of electromagnetic fields and that of free scalar one on the late-time structure of quantum entanglement. However time evolution of ∆S where · · · are contributions from the higher · · · are contributions from the higher order O t−l l 3 . Their time evolution shows that quasi-particles created by E 2,3 (B 2,3 ) enter the region A faster than those generated by E 1 (B 1 ). These behaviors seem to be natural since particles created by E 1 (B 1 ) do not propagate along the direction parallel to x 1 . ∆S

(n)
A in Table.1 shows that they are invariant under the transformation, where ǫ µνρσ is an anti-symmetric tensor. Under the transformation in (3.4), the local operator E i (B i ) changes to −B i (E i ). Therefore this duality changes a locally excited state to a different one.

Composite Operators Constructed of Only Electric or Magnetic Fields
The excited states which we consider here are generated by acting with the following operators: (a) E i E j or B i B j , (b) E 2 or B 2 . We study the time evolution and the late-time value of ∆S (n) A for them.
A for E 1 (B 1 ) and E 2,3 (B 2,3 ). The horizontal and vertical axes correspond to time t and ∆S (2) A , respectively. The red and blue lines correspond to ∆S (2) A for E 1 (B 1 ) and E 2,3 (B 2,3 ), respectively.  It is the same as that of ∆S

(n)
A for φ 2 in the massless free scalar field theories as in [37,38]. Therefore the late-time value of ∆S which can be interpreted as maximum (Rényi) entanglement entropy for ρ A = 1 4 diag (1, 1, 1, 1). It is the same as ∆S

(n)
A for the excited state given by acting with the operator φ a φ b on the ground state. Here a, b denote the kind of scalar fields, and a = b. They are two kinds of massless free scalar fields. The time evolution of ∆S

(n)
A for them is summarized in Table.2. Table.2 shows ∆S In order to study whether E 1 acts on the late-time structure of quantum entanglement differently from E 2,3 5 , we study the late-time value of ∆S

(n)
A for the given locally excited state: Before studying its late-time value, we comment on its time evolution. Before t = l, ∆S (n) A for the state in (3.7) vanishes and after t = l, it increases. Its time evolution is summarized in Table.3. After t = l, as in Table.3, ∆S

(n)
A is given by where D, N i and P i are defined in Table.3. If we take the late time limit (0 < l ≪ t), the ratios of P i and N i to D reduce to constant numbers [50], where we ignore the higher order contribution O l t . Amazingly, The sum of them is 1. Therefore if the effective reduced density matrix is defined by which can be interpreted in terms of quasi-particles created by (φ 1 ) is constructed of three kinds of free scalar fields. Therefore, there are no differences between the effect of E 1 and that of E 2,3 on the late-time structure of quantum entanglement. As in the Table.2, ∆S (n) A for E 2 is equivalent to that for B 2 . Therefore they is the electric-magnetic duality invariant. .

Composite Operators Constructed of Both Electric and Magnetic Fields
In the previous two subsection we study how the entanglement structure changes at the late time if either electric or magnetic fields act on the ground state. However we do not uncover how it changes at the late time when both of them act on the ground state. Here we study ∆S

(n)
A for (a) E 2 1 + B 2 1 , (b) E i B j and (c) F µν F µν and B · E, which can show that E i and B i act on the late-time structure of quantum entanglement differently from scalar fields such as φ a .

E
Here in order to study whether there are differences between the effects of electric and magnetic fields on the late-time structure of quantum entanglement, we study ∆S

(n)
A for the following excited state: Before investing the late time value of ∆S

(n)
A , let's study the time evolution of ∆S (n)

A . ∆S (n)
A vanishes before t = l. After t = l, its time evolution is summarized in Table.3. If you take the late time limit t → ∞, the late time value of ∆S (n) A reduces to the (Rényi) entanglement entropy whose effective reduced density matrix is given by diag (1, 1, 1, 1) . (3.14) Its entropies are given by ∆S (n) It shows there are no differences between the effects of electric and magnetic fields on the late-time structure.

E i B j
Here let's find out how the operators constructed of both electric and magnetic fields, E i B j , affect the late-time structure of quantum entanglement. The late-time values of ∆S (n) A for E i B j except for E 2 B 3 and E 3 B 2 are the same as (3.6). Their time evolution is summarized in Table.2.
On the other hand, after t = l the time evolution of ∆S (n) A for E 2 B 3 or E 3 B 2 is summarized in Table.3. We can see that it has the electric-magnetic duality from the Table.3. The latetime value of ∆S (n) A is given by (Rényi) entanglement entropy whose reduced density matrix is given by diag (25,7,7,25).  It shows how different the effect of E 1 (B 1 ) is from that of E 2,3 (B 2,3 ) on the structure. The value can not be interpreted in terms of quasi-particles created by scalar fields such as φ a φ b . As we will explain later, in the entangled quasi-particle interpretation, there is a commutation relation between the quasi-particle created by E 2 (B 2 ) and that by B 3 (E 3 ).

B · E and F µν
We finally study ∆S

(n)
A for more complicated operators, B · E, F µν F µν and B 2 E 3 − B 3 E 2 . Before t = l, ∆S (n) A for them vanish, but after t = l, they increases. The detail of them is summarized in Table.3 6 . It shows that ∆S

(n)
A for B · E is the same as that for F µν F µν . The effective reduced density matrices for B · E (F µν F µν ) and B 2 E 3 − B 3 E 2 are given by diag (30,30,16,16,49,49,1,1) , diag (50,50,7,7,7,7) ,  which are for O = B 2 E 3 − B 3 E 2 . As we will explain in the next section, they can be reproduced by using a late-time algebra which electromagnetic quasi-particles obey.

A Late-time Algebra
We interpret the late-time value of ∆S (n) A in terms of quasi-particles. More precisely, let's interpret the effective reduced density matrix in (3.10) in terms of quasi-particles. The effective reduced density matrix for the excited state generated by a composite operator O(−t, −l, x) is defined by whereN is a normalization constant. The operator O is assumed to be constructed of electric and magnetic fields 7 . As in [37,38,39,41], these fields can be decomposed into left moving and right moving electromagnetic quasi-particles as follows, where since we take x 1 ≥ 0 as A in this paper, left-moving and right-moving quasi-particles correspond to particles included in B and A at late time, respectively. The ground state for them is defined by The late-time algebra which quasi-particles obey is given by which is obtained so that the results by the replica trick are reproduced. Here C is a real number 8 . In the gauge theory in addition to (4.2), we need the following commutation relation for different particles: where X R,L and Y R,L are given by Here X R,L and Y R,L are real numbers 9 . The commutation relation between electric (magnetic) quasi-particles is determined so that the effective density matrices computed by (4.4) are consistent with (3.11) respectively. The relation for the quasi-particles by E 1 should be the same as that for B 1 so that the effective density matrix in (4.1) reproduces ∆S

(n)
A for the matrix in (3.14). That between quasi-particles generated by E 2 (E 3 ) and those by B 3 (B 2 ) reproduces the matrix in (3.16). We also check that ∆S are reproduced by using the commutation relation in (4.4) and (4.5). 7 Here ρ e A is not the same as the reduced density matrix for the locally excited state. It is for a "effective" stateN O |0 . It is different form the "original" locally excited state. 8 The redefinition of quasi-particles can absorb the constant C. 9 We find the correspondence between propagators and commutation relations. The commutations can be defined by the late time limit of propagators. We will discuss the detail of the correspondence in [50]. When we use this correspondence, X L = − 3 4 C.
The relation in (4.5) shows that the effect of fields along the direction vertical to ∂A is significantly different from that along the direction parallel to ∂A on the late-time structure. It makes the effects of electromagnetic fields different from that of free scalar fields on the late-time structure of quantum entanglement.

Conclusion and Future Problems
We also studied how gauge invariant operators such as E i , B i and the composite operators constructed of them changes the structure of quantum entanglement by studying ∆S A , which we studied in this paper, is invariant under the duality transformation. If only E i or B i acts on the ground state, without taking the late time limit, the time evolution of ∆S (n) A depends on them. Due to the duality, ∆S A for E 2,3 (B 2,3 ) increases slower than that for E 1 (B 1 ). However they can not show the difference between the effects of electromagnetic fields and that of scalar fields on the late-time structure because the late-time values of ∆S (n) A for them can be interpreted in terms of quasi-particle created by scalar fields.
On the other hand, the late-time values of ∆S (n) A for the specific operators constructed of both electric and magnetic fields can not be interpreted in terms of quasi-particles by scalar fields. They show that there are differences between the effects of electromagnetic and that of scalar fields on the late-time structure of quantum entanglement. If their late-time values are interpreted in terms of electromagnetic quasi-particles in (4.2), there are commutation relations between E 2 (E 3 ) and B 3 (B 2 ), which make the effect of electromagnetic field significantly different from that of scalar fields on the late-time structure. The effect of E 1 and B 1 on the late-time structure is different from that of E 2,3 and B 2,3 .

Future Problem
We finish with comments on a few of future problems: • In this paper we only consider 4d Maxwell gauge theory which has conformal symmetry. D( = 4) dimensional Maxwell gauge theory is not a CFT. Therefore it is interesting to generalize the analysis in 4d Maxwell theory to that in the theories with general dimensions.
• We expect that the structure of the late-time algebra depends on the spacetime dimension D. Then it is also interesting to study it in general dimensions.
comments on this work.