-
The quantum field theory entanglement is widely studied in the literature [1-3]. Recently, the source of entanglement entropy was likewise considered in [4]. It expresses the vacuum state in terms of the superposition of particle pair states in the configuration space by the Bogoliubov transformation, which transforms the original degrees of freedom to actual local degrees of freedom. Interestingly, it tells us that the entanglement entropy originates from the particle pairs across the boundary. The area law [1, 5-7] emerges naturally by this method. In [4], only the free scalar field is discussed. In this study, we generalize it to 1+1 and 2+1 dimensional Majorana (spin-
12 ) field, 2+1 dimensional U(1) gauge (free spin-1) field, and 3+1 dimensional weak gravitational (free massless spin-2) field. Because there is a Hilbert space with a tensor product structure for the spin-12 field, the generalization is straight forward. As for higher spin fields, more treatments are required because of the gauge redundancy.The definition of entanglement entropy in gauge theory was discussed in [8, 9]. The definition of entanglement entropy is heavily dependent on the tensor product structure of the Hilbert space. However, because of the gauge redundancy, there is no natural structure as such in the Hilbert space of the gauge field. Some studies addressed the entanglement entropy of the gauge field with various prescriptions [10, 11], which give up the direct tensor product structure and introduced a "center". It is clear that this issue is also present in the gravitational theories [12, 13], as gravitational theories likewise are gauge theories. The tensor product structure problem of gravitational theories is still under discussion [14]. In this study, we deal with the gauge redundancy by the light-cone gauge for both the free U(1) gauge field and the gravitational field (or strictly speaking, a free spin-2 field). We demonstrate that this gauge-fixing provides a natural tensor product structure in the Hilbert spaces belonging to them, while surrendering the explicit Lorentz invariance.
Imposing the light-cone gauge, we find that the Lagrangian of the U(1) gauge field and gravitational field simplifies significantly and behaves like a massless scalar field. Subsequently, we quantize them and obtain their Hilbert spaces with the tensor product structure. We perform the Bogoliubov transformation of the theories and express the ground state in terms of a local basis of degrees of freedom. Finally, we calculate the second Rényi entropy, which is a good approximation of the entanglement entropy. We find that in such a prescription, the area law is also naturally preserved.
Our paper is organized as follows. In Section 2, we make a Bogoliubov transformation of the 1+1 dimensional Majorana field and consider its second Rényi entropy. In Section 3, we perform this procedure for the 2+1 dimensional Majorana field. In Section 4, we derive the Lagrangian of the 2+1 dimensional free U(1) gauge field in the light-cone gauge and obtain its Hilbert space with the natural tensor product structure. Then we perform a Bogoliubov transformation and calculate the second Rényi entropy. In Section 5, we generalize our method to the 3+1 dimensional gravitational (free spin-2) field in the light-cone gauge. In Section 6, we provide conclusions of our work.
-
The Majorana field
ψ is a field sharing the same Lagrangian with the Dirac field, but satisfying the conditionψ∗=ψ , which means that the anti-particle and particle of the Majorana field are the same [15]. In 1+1 dimensions at t = 0, a Majorana field withγμ=(σy,iσx) has the form [16]ψ(x)=∫dk√2π1√2ω[1√ω+k(ω+kiμ)bk+1√ω−k(ω−k−iμ)b†−k]eikx,
(1) where the annihilation and creation operators bk and
b†p satisfy the anti-commutation relation{bk,b†p}=δ(k−p) . We define the vacuum|0⟩ as the vacuum of this set of operators, i.e.bk|0⟩=0 .Let us consider the following Bogoliubov transformation [17]
ψ(x)=∫dk√2π(a1ka2k)eikx,
(2) which transforms the operators
(bk,b†−k) into the new set of operators(a1k,a2k) . For the Majorana field, we have the relation{ψα(x),ψβ(y)}=δαβδ(x−y).
(3) where
α,β=0,1 represent the two components of the Majorana field. To satisfy (3), the operatorsa1k anda2p should have the relation{aik,ajp}=δijδ(k+p),
(4) where
i,j=1,2 . We choosea1k=1√2(ak+a†−k),
(5) a2k=−i1√2(ak−a†−k),
(6) where
(ak,a†−k) are new annihilation and creation operators that satisfy the anti-commutation relation{ak,a†p}= δ(k−p) . We can verify that (5) and (6) satisfy the constraint (4). We define the vacuum|Ω⟩ , which is annihilated by ak, i.e.ak|Ω⟩=0 . Hence, we obtain the Bogoliubov transformation between(ak,a†−k) and(bk,b†−k) (bkb†−k)=12√ω(√ω+k−√ω−k√ω+k+√ω−k√ω+k+√ω−k√ω−k−√ω+k)×(aka†−k).
(7) Certainly, we have
bk=12√ω((√ω+k−√ω−k)ak+(√ω+k+√ω−k)a†−k).
(8) Because we have
bk|0⟩=0 , the new set of operators satisfy((√ω+k−√ω−k)ak+(√ω+k+√ω−k)a†−k)|0⟩=0.
(9) The vacuum
|0⟩ can be expressed as a state constructed from the new set of operators|0⟩=1γe−∑kCka†ka†−k|Ω⟩
(10) =1γ∏k(1−Cka†ka†−k)|Ω⟩,
(11) where
γ is the normalization factor,|Ω⟩ is the vacuum defined byak|Ω⟩=0 and the coefficient Ck is fixed byCk=√ω+k+√ω−k√ω+k−√ω−k=ω+μk.
(12) Following [4], let us consider a system with a finite number of sites. The inverse Fourier transform of the operator ak is
ak=∑NaNe−ikN,
(13) where N is the site label. The vacuum
|Ψ⟩≡|0⟩ of the operator bk can thus be written as|Ψ⟩=1γ∏k(1−∑N,LCkeik(N−L)a†Na†L)|Ω⟩≃1γ(1−∑k∑N,LCkfkNLa†Na†L)|Ω⟩,
(14) where
fkNL=eik(N−L) . To compute the entanglement entropy, the configuration space is divided into two regions A andˉA , where the sites in each region are labelled by small letters (n) and small letters with bar (ˉn ), respectively. The region A we choose can be any subregion of the total system. The state|Ψ⟩ is|Ψ⟩≃1γ(1−∑k∑nlCkfknla†na†l−∑k∑ˉnˉlCkfkˉnˉla†ˉna†ˉl−12∑k∑ˉnlCk˜fkˉnla†ˉna†l)|Ω⟩,
(15) where
˜fkˉnl=fkˉnl+fklˉn .A reduced density matrix of
ρA is constructed by tracing out the degrees of freedom in regionˉA .ρA=TrˉA(|Ψ⟩⟨Ψ|)≃1γ2[(1−∑k∑nlCkfknla†na†l)|Ω⟩AA⟨Ω|(1−∑k′∑n′l′C∗k′f∗k′n′l′an′al′)]ˉA⟨Ω|Ω⟩ˉA+|Ω⟩AA⟨Ω|[ˉA⟨Ω|(∑k′∑ˉn′ˉl′C∗k′f∗k′ˉn′ˉl′aˉl′aˉn′)(∑k∑ˉnˉlCkfkˉnˉla†ˉna†ˉl)|Ω⟩ˉA]+a†l|Ω⟩AA⟨Ω|al′[ˉA⟨Ω|(∑k′∑ˉn′l′C∗k′˜f∗k′ˉn′l′aˉn′)(∑k∑ˉnlCk˜fkˉnla†ˉn)|Ω⟩ˉA].
(16) This expression is already quite revealing. We can see that the first and second line actually lead to no entanglement. They contain operators acting only on A or
ˉA . The third line explains the large amount of entanglement between regions A andˉA - entanglement is created when pairs of particles are created, one inside A and the other inˉA . Moreover, as we shall see below, these particles that are created cannot be separated by too large of a distance, because the amplitude decreases rapidly with separation.Tracing out the degrees of freedom in region
ˉA gives the following coefficients inρA κ0=ˉA⟨Ω|(∑k′∑ˉn′ˉl′C∗k′f∗k′ˉn′ˉl′aˉl′aˉn′)(∑k∑ˉnˉlCkfkˉnˉla†ˉna†ˉl)|Ω⟩ˉA=∑kk′C∗k′Ck(∑ˉnˉl(f∗k′ˉnˉlfkˉnˉl+f∗k′ˉlˉnfkˉnˉl)+∑ˉn2),ˉA⟨Ω|(∑k′∑ˉn′l′C∗k′˜f∗k′ˉn′l′aˉn′)(∑k∑ˉnlCk˜fkˉnla†ˉn)|Ω⟩ˉA=∑kk′∑ˉnll′C∗k′Ck˜f∗k′ˉnl′˜fkˉnl.
(17) The reduced density matrix
ρA is thenρA≃1γ2[(1+κ0)|Ω⟩⟨Ω|−|2⟩⟨Ω|−|Ω⟩⟨2|+|1⟩⟨1|+|2⟩⟨2|],
(18) where the subscript label A for the vacuum in region A is dropped,
|2⟩=∑k∑nlCkfknla†na†l|Ω⟩,
(19) |1⟩⟨1|=∑kk′∑ˉnll′C∗k′Ck˜f∗k′ˉnl′˜fkˉnla†l|Ω⟩⟨Ω|al′.
(20) Let us consider the second Rényi entropy S2. The square of the reduced density matrix is
ρ2A=1γ4[((1+κ0)2+κ2)|Ω⟩⟨Ω|−(1+κ0+κ2)|Ω⟩⟨2|−(1+κ0+κ2)|2⟩⟨Ω|+κ1|1⟩⟨1|+(1+κ2)|2⟩⟨2|].
(21) The coefficients
κ1 andκ2 are given byκ1|1⟩⟨1|=|1⟩⟨1|1⟩⟨1|=(∑kk′∑ˉnll′C∗k′Ck˜f∗k′ˉnl′˜fkˉnla†l|Ω⟩⟨Ω|al′)×(∑k″k‴∑ˉpqq′C∗k‴Ck″˜f∗k‴ˉpq′˜fk″ˉpqa†q|Ω⟩⟨Ω|aq′)=∑kk′k″k‴∑ˉnll′ˉpq′C∗k′CkC∗k‴Ck″˜f∗k′ˉnl′˜fkˉnl˜f∗k‴ˉpq′˜fk″ˉpl′a†l|Ω⟩⟨Ω|aq′,
κ2=⟨2|2⟩=⟨Ω|(∑k′∑n′l′C∗k′f∗k′n′l′al′an′)(∑k∑nlCkfknla†na†l)|Ω⟩=∑kk′C∗k′Ck(∑nl(f∗k′nlfknl+f∗k′lnfknl)+∑n2).
S2 is computed by taking the trace of (21)
S2=−lnTrρ2A=−ln((1+κ0)2+κ2+|κ1|+(1+κ2)κ2),
(22) where
|κ1| is given by|κ1|=∑kk′k″k‴∑ˉnˉpll′C∗k′CkC∗k‴Ck″˜f∗k′ˉnl′˜fkˉnl˜f∗k‴ˉpl˜fk″ˉpl′.
(23) In the 1+1 dimensional Majorana fermion case, the coefficients Ck control the range of the interaction. The factor
∑kCkfknl can be converted to an integral by taking a continuum limit with an IR cut-offϵ and UV cut-off KFnl=∑kCkfknl≃∫Kϵdkω+μkeik(n−l)a
(24) ≈F(x)=∫Kϵdkω+μkeikx.
(25) This is evaluated numerically and shown in Fig. 1.
Figure 1. (color online) Red and blue lines represent the real and imaginary part of the integral (25), respectively. Both the real and imaginary parts of the integral decrease quickly with distance, and the long range contribution is very small.
ϵ=0.1 and K = 10 are used in this plot.The figure shows that the interaction between sites falls over large distances. The sites near the boundary make the main contribution to the entanglement entropy.
-
In 2+1 dimensions at t = 0, the Majorana field with
γμ=(σy,−iσz,iσx) has the form [16]ψ(→x)=∫d2→k2π1√2ω[1√ω+ky(ω+kykx+iμ)bk+1√ω−ky(ω−ky−kx−iμ)b†−k]ei→k⋅→x,
(26) where the annihilation and creation operators bk and
b†p satisfy the anti-commutation relation{bk,b†p}=δ(→k−→p) . We define the vacuum|0⟩ as the vacuum of this set of operators, i.e.bk|0⟩=0 .Let us consider the following Bogoliubov transformation
ψ(→x)=∫d2→k2π(a1ka2k)ei→k⋅→x,
(27) which transforms the operators
(bk,b†−k) to the new set of operators(a1k,a2k) . For the 2+1 dimensional Majorana field, we have the relation{ψα(→x),ψβ(→y)}=δαβδ(→x−→y),
(28) where
α,β=0,1 represent two components of the Majorana field. To satisfy (28), the operatora1k anda2p should have the relation{aik,ajp}=δijδ(→k+→p),
(29) where
i,j=1,2 . We choosea1k=1√2(ak+a†−k),
(30) a2k=−i1√2(ak−a†−k),
(31) where
(ak,a†−k) are new annihilation and creation operators that satisfy the anti-commutation relation{ak,a†p}=δ(→k−→p) . We can verify that (30) and (31) satisfy the constraint (29). We define the vacuum|Ω⟩ , which is annihilated by ak, i.e.ak|Ω⟩=0 . Hence, we obtain the Bogoliubov transformation between(ak,a†−k) and(bk,b†−k) (bkb†−k)=−12√ω√kx−iμkx+iμ(−kx+iμ√ω−ky+i√ω−ky−kx+iμ√ω−ky−i√ω−ky−kx+iμ√ω−ky−i√ω+ky−kx+iμ√ω+ky+i√ω+ky)(aka†−k).
(32) Certainly, we have
bk=−12√ω√kx−iμkx+iμ((−kx+iμ√ω−ky+i√ω−ky)ak−(kx+iμ√ω−ky+i√ω−ky)a†−k).
(33) Because we have
bk|0⟩=0 , the new set of operators satisfy((−kx+iμ√ω−ky+i√ω−ky)ak−(kx+iμ√ω−ky+i√ω−ky)a†−k)|0⟩=0.
(34) The vacuum
|0⟩ can be expressed as a state constructed from the new set of operators|0⟩=1γe−∑kCka†ka†−k|Ω⟩
(35) =1γ∏k(1−Cka†ka†−k)|Ω⟩,
(36) where
γ is the normalization factor,|Ω⟩ is the vacuum defined byak|Ω⟩=0 , and the coefficient Ck is fixed byCk=kx+iμ+iω−ipykx+iμ−iω+ipy.
(37) Following [4], let us consider a system with a finite number of sites. The inverse Fourier transform of the operator ak is
ak=∑→Na→Ne−i→k⋅→N,
(38) where
→N is the site label. The vacuum|Ψ⟩≡|0⟩ of the operator bk can now be written as|Ψ⟩≃1γ∏→k(1−∑→N,→LCkei→k⋅(→N−→L)a†→Na†→L)|Ω⟩≃1γ(1−∑→k∑→N,→LCkfkNLa†→Na†→L)|Ω⟩,
(39) where
fkNL=ei→k⋅(→N−→L) . For a simplification of notations, we neglect the vector nation and denote→n as n. To compute the entanglement entropy, the configuration space is divided into two regions A andˉA , where the sites in each region are labelled by small letters (n) and small letters with bars (ˉn ), respectively. The region A we choose can be any subregion of the total system. The state|Ψ⟩ is|Ψ⟩≃1γ(1−∑k∑nlCkfknla†na†l−∑k∑ˉnˉlCkfkˉnˉla†ˉna†ˉl−12∑k∑ˉnlCk˜fkˉnla†ˉna†l)|Ω⟩,
(40) where
˜fkˉnl=fkˉnl+fklˉn .Following the same procedure of the case in 1+1 dimensional Majorana field, we can obtain the reduced density matrix and the second Rényi entropy S2, which have the same form as (18) and (22), respectively. The only difference is the expression of Ck.
In the 2+1 dimensional Majorana fermion case, the coefficients Ck also control the range of the interaction. The factor
∑kCkfknl can be converted to an integral by taking a continuum limit with an IR cut-offϵ and UV cut-off KFnl=∑kCkfknl≃∫Kϵ∫Kϵdkxdkykx+iμ+iω−ikykx+iμ−iω+ikyei→k⋅(→n−→l)a
(41) ≈F(x,y)=∫Kϵ∫Kϵdkxdkykx+iμ+iω−ikykx+iμ−iω+ikyei(kxx+kyy).
(42) This is evaluated numerically and shown in Fig. 2.
Figure 2. (color online) Plots of integral (42) in 2+1 dimensional Majorana fermion model. The left and right figures represent the real and imaginary parts of the integral, respectively. The amplitude of the surface decreases quickly with distance, and the long range contribution can be neglected.
ϵ=0.1 and K = 10 are used in these plots.The figures show that the interaction between sites diminishes over large distances. The sites near the boundary contribute most significantly to the entanglement entropy. Hence, an area law is expected in this model.
-
In this section, we consider the 2+1 dimensional free U(1) gauge field in the light-cone gauge. To this end, we derive the Lagrangian, and verify that it is the same as a massless scalar field. Subsequently, we perform the Bogoliubov transformation to calculate the entropy.
-
We consider the entanglement entropy of 2+1 dimensional Maxwell fields with light-cone gauge-fixing. In this study, we use the Minkowski metric
ημν=diag(−1,1,1) with the Minkowski coordinate(x0,x1,x2) . We also introduce the light-cone coordinatex+,x−,x2 withx±≡1√2(x0±x1).
(43) With regard to the light-cone coordinate, interested readers can read the book [18]. In the light-cone gauge, the vector is defined as
a±≡1√2(a0±a1),
(44) and the metric becomes
ˆημν=(0−10−100001).
(45) The vector in the light-cone gauge has the properties
a+=−a− anda−=−a+ .Following is the calculation of the 2+1 dimensional Maxwell fields. We start with the Lagrangian
L=−14FμνFμν . BecauseL is a scalar,μ ;ν can be 0, 1, 2 or +, −, 2. We start with the light-cone coordinate. Because of the asymmetry ofFμν ,F++=F−−=F22=0 , the Lagrangian becomesL=−14FμνFμν=−14(2F+−F+−+2F+2F+2+2F−2F−2).
(46) For the light-cone gauge,
A+=0 andA−=−A+=0 . By imposing the gauge-fixing, we find that the Lagrangian becomesL=12(∂−A−)2+∂−A2(∂+A2+∂2A−).
(47) With the light-cone gauge fixing
A+(p)=0 , we haveA−=1p+(p2A2) [18]. We perform the Fourier transformationAμ(x)=∫d3p(2π)3eipxAμ(p).
(48) The Fourier transformation of
L is˜L=−12(p−A−)2−p−A2(p+A2+p2A−)=−12(p−1p+p2A2)2−p−A2(p+A2+p21p+p2A2)=−12(p2A2)2+p−A2p−A2+p2A2p2A2=12(p2A2)2+p−A2p−A2.
(49) When written in the momentum space of Minkowski spacetime, it becomes
˜L=12(p2A2)2−12((p0)2−(p1)2)(A2)2=−12(p0A2)2+12((p1A2)2+(p2A2)2).
(50) When we perform the inverse Fourier transformation and return to Minkowski coordinates, the expression becomes
L=12(∂0A2)2−12[(∂1A2)2+(∂2A2)2].
(51) This is exactly the same Lagrangian as the massless scalar field, and the corresponding Hamiltonian is
H=12(∂0A2)2+12[(∂1A2)2+(∂2A2)2].
(52) To simplify the notation in the next subsection, we depict
A2 asAy . The Lagrangian and Hamiltonian can be written asL=12(∂tAy)2−12[(∂xAy)2+(∂yAy)2],
(53) H=12(∂tAy)2+12[(∂xAy)2+(∂yA2)2].
(54) Here, there is no gauge redundancy in the Lagrangian (53). The Hilbert space of
Ay should have the tensor product structure, so we can consider its entanglement entropy. It behaves like a free scalar field, which coincides with the trivial center case of [11]. -
In the 2+1 dimensional U(1) gauge field with the light-cone gauge, from (53), we have the equation of motion
◻Ay=0,
(55) where
◻=−∂2∂t2+∂2∂x2+∂2∂y2 . We can find that there is only one physical degree of freedom in the 2+1 dimensional U(1) gauge field, and the component with physical freedomAy satisfies the equation for the massless scalar field. For the non-zero components of the gauge field, we have the solution(Ay˙Ay)=∫d2→k2π1√2ω(1−iω)bkei→k⋅→x+(1iω)b†−kei→k⋅→x.
(56) Let us consider a Bogoliubov transformation, which transforms the set of operators
(bk,b†−k) to the following set of operators(Ay˙Ay)=∫d2→k2π(a1ka2k)ei→k⋅→x.
(57) With operator
a1k anda2k , the commutation relations ofAy and˙Ay are expressed as below[Ay(→x,t),˙Ay(→y,t)]=∫d2→k2π∫d2→p2π[a1k,a2p]ei(→k⋅→x+→p⋅→y)=iδ(→x−→y),
(58) [Ay(→x,t),Ay(→y,t)]=[˙Ay(→x,t),˙Ay(→y,t)]=0.
(59) We need to have
[a1k,a2p]=iδ(→k+→p),
(60) [a1k,a1p]=[a2k,a2p]=0.
(61) We choose
a1k=1√2α(ak+a†−k),
(62) a2k=−i√α2(ak−a†−k),
(63) where
α is a real parameter. The operator ak anda†p have the commutation relations[ak,a†p]=δ(→k−→p),
(64) [ak,ap]=[a†k,a†p]=0.
(65) We define
ak|Ω⟩=0 , where|Ω⟩ is the vacuum of the new set of operators. We find that (62) and (63) satisfy (60) and (61). Thus, from (56), (57), (62), and (63), we have the Bogoliubov transformation(bkb†−k)=1i√2ω(iω√2α+i√α2iω√2α−i√α2iω√2α−i√α2iω√2α+i√α2)(aka†−k),
(66) with
ω=√k2x+k2y . From the above Bogoliubov transformation, we havebk=1i√2ω((iω√2α+i√α2)ak+(iω√2α−i√α2)a†−k).
(67) Here
(ak,a†−k) are the annihilation and creation operators of the new modes. We also have the vacuum|0⟩ , which is annihilated by bkbk|0⟩=0.
(68) From (67), we find that the new set of operators satisfy
((iω√2α+i√α2)ak+(iω√2α−i√α2)a†−k)|0⟩=0.
(69) The vacuum
|0⟩ can be expressed in terms of a state constructed from the new set of operators|0⟩=1γe−∑kCka†ka†−k|Ω⟩
(70) ≃1γ∏k(1−Cka†ka†−k)|Ω⟩,
(71) where
γ is the normalization factor,|Ω⟩ is the vacuum defined byak|Ω⟩=0 , and the coefficient Ck is fixed byCk=ω−αω+α.
(72) We find that the Bogoliubov transformation of the light-cone gauge field is very similar to that of the free scalar field [4].
Let us consider a system with a finite number of sites. The inverse Fourier transform of the operator ak is
ak=∑→Na→Ne−i→k⋅→N,
(73) where
→N is the site label. The vacuum|Ψ⟩≡|0⟩ of the operator bk can now be written as|Ψ⟩≃1γ∏→k(1−∑→N,→LCkei→k⋅(→N−→L)a†→Na†→L)|Ω⟩≃1γ(1−∑→k∑→N,→LCkfkNLa†→Na†→L)|Ω⟩,
(74) where
fkNL=ei→k⋅(→N−→L) . For a simplification of notations, we neglect the vector nation and denote→n as n. To compute the entanglement entropy, the configuration space is divided into two regions A andˉA , where the sites in each region are labelled by small letters (n) and small letters with bars (ˉn ), respectively. The region A can be chosen as any subregion of the total system. The state|Ψ⟩ is|Ψ⟩≃1γ(1−∑k∑nlCkfknla†na†l−∑k∑ˉnˉlCkfkˉnˉla†ˉna†ˉl−12∑k∑ˉnlCk˜fkˉnla†ˉna†l)|Ω⟩,
(75) where
˜fkˉnl=fkˉnl+fklˉn .Following the same procedure of previous cases, we can obtain the reduced density matrix and second Rényi entropy S2, which have the same form as (18) and (22), respectively. The only difference is the expression of Ck.
In 2+1 dimensional U(1) gauge field theory in the light-cone gauge, the coefficients Ck also control the range of the interaction. The factor
∑kCkfknl can be converted to an integral by taking a continuum limit with an IR cut-offϵ and UV cut-off KFnl=∑kCkfknl≃∫Kϵ∫Kϵdkxdkyω−αω+αei→k⋅(→n−→l)a
(76) ≈F(x,y)=∫Kϵ∫Kϵdkxdkyω−αω+αei(kxx+kyy),
(77) with
ω=√k2x+k2y . This is evaluated numerically and shown in Fig. 3.Figure 3. (color online) Plots of integral (77) in 2+1 dimensional U(1) gauge field with light-cone gauge. The left and right figures represent the real and imaginary parts of the integral, respectively. The amplitude of the surface decreases quickly with distance, and the long range contribution can be neglected.
ϵ=0.1 , K = 10 andα=1 are used in these plots.The figures show that the interaction between sites diminishes over large distances. The sites near the boundary contribute most significantly to the entanglement entropy. Hence, an area law is expected in this model.
-
In this section, we consider the 3+1 dimensional free spin-2 field theory in the light-cone gauge. To this end, we derive the Lagrangian and find that it is also the same as the massless scalar field. Subsequently, we perform the Bogoliubov transformation to calculate the entropy.
-
We consider the entanglement entropy of the 3+1 dimensional weak gravitational field
gμν=ημν+hμν with light-cone gauge-fixing. Here,hμν is a spin-2 field. In this study, we use the Minkowski metricημν=diag(−1,1,1,1) with the Minkowski coordinate(x0,x1,x2,x3) . We also introduce the light-cone coordinatex+,x−,x2,x3 withx±≡1√2(x0±x1).
(78) In the light-cone gauge, the vector is defined as
a±≡1√2(a0±a1),
(79) and the metric becomes
ˆημν=ˆημν=(0−100−100000100001).
(80) The vector in the light-cone gauge has the properties
a+=−a− anda−=−a+ .Now we come to the calculation of weak gravitational (spin-2) field. We start with the Ricci scalar R. Because R is a scalar,
μ ;ν can be 0, 1, 2, 3 or +, −, 2, 3. We begin with the light-cone coordinate. We impose the gauge-fixingh++=h+−=h+I=0 , whereI=2,3 . For the components with subscript indexes, the gauge fixing ish−−=h−+=h−I=0 . For calculating the Ricci scalar R, we useΓμαβ=gμν2[∂gαν∂xβ+∂gβν∂xα−∂gαβ∂xν]
(81) and
Rμν=∂αΓαμν−∂νΓαμα+ΓαβαΓβμν−ΓαβνΓβμα.
(82) The Christoffel symbols are
Γ+++=12∂h++∂x−,
(83) Γ++I=12∂h+I∂x−,
(84) Γ+IJ=12∂hIJ∂x−,
(85) Γ++−=Γ+−−=Γ+−I=0,
(86) Γ−++=−12∂h++∂x+,
(87) Γ−+−=−12∂h++∂x−,
(88) Γ−−−=0,
(89) Γ−+I=−12∂h++∂xI,
(90) Γ−−I=−12∂hI+∂x−,
(91) Γ−IJ=−12(∂hI+∂xJ+∂hJ+∂xI−∂hIJ∂x+),
(92) ΓI++=−12(2∂h+I∂x+−∂h++∂xI),
(93) ΓI+−=12∂h+I∂x−,
(94) ΓI−−=0,
(95) ΓI+J=12(∂h+I∂xJ+∂hJI∂x+−∂h+J∂xI),
(96) ΓI−J=12∂hJI∂x−,
(97) ΓIJK=12(∂hJI∂xK+∂hKI∂xJ−∂hJK∂xI).
(98) The Ricci scalar is given by
R=ˆημνRμν=ˆη+−R+−+ˆη−+R−++ˆηIJRIJ=−2R+−+RII,
(99) where the repeated I, J, and K are summed, and we take this convention below. We have
R+−=−12∂2h++∂x−2+12∂2h+I∂xI∂x−−14∂hIJ∂x−(∂h+I∂xJ+∂hJI∂x+−∂h+J∂xI)
(100) and
\setcounterequation101RII=∂2hI+∂x−∂xI+∂2hIJ∂xJ∂xI−12(∂h+I∂x−)2−∂hIJ∂x−(∂hIJ∂x+−∂hI+∂xJ)−14[(∂hJK∂xI)2−(∂hIK∂xJ−∂hJI∂xK)2].(101)
The Ricci scalar can be written as
\setcounterequation102R=∂2h−−∂x−2+∂hIJ∂xI∂xJ+12∂hIJ∂x−(−∂h−I∂xJ+∂hJI∂x++∂h−J∂xI)−12(∂h−I∂x−)2−∂hIJ∂x−(∂hIJ∂x++∂hI−∂xJ)−14[(∂hJK∂xI)2−(∂hIK∂xJ−∂hJI∂xK)2].(102)
When expressed in Fourier space and considering the light-cone gauge-fixing, we can obtain [18]
hI−=1p+pJhIJ
(103) and
h−−=1p+pIh−I=pIpJ(p+)2hIJ.
(104) With the above relations, we have
R=−2pIpJhIJ+12p−p+hIJhIJ−12pJpKhIJhIK+14[(pIhJK)2−(pJhIK−pKhJI)2].
(105) Because
hIJ is symmetric and traceless, there are only two degrees of freedomh22 andh23 . We expand the above expression withI,J,K=2,3 , and obtainR=−2pIpJhIJ+(p+p−−12p22−12p23)[(h22)2+(h23)2]=−2pIpJhIJ+12((p0)2−(p1)2−(p2)2−(p3)2)×[(h22)2+(h23)2].
(106) Apart from the total derivative term
∂I∂JhIJ , the Lagrangian of the gravitational field can be written asL=12((∂thyy)2−(∇hyy)2+(∂thyz)2−(∇hyz)2),
(107) where
∇=∂x→i+∂y→j+∂z→k . We have replacedh22 andh23 withhyy andhyz respectively.There is no gauge redundancy in the Lagrangian (107). We can expect that the Hilbert space ofhyy andhyz should have the tensor product structure, so we can consider the entanglement entropy of this model. Hence, we provide a prescription of the entanglement entropy of a free spin-2 field. -
In the 3+1 dimensional gravitational (spin-2) field with the light-cone gauge, we have two independent degrees of freedom. From (107), we have the equations of motion
◻hyy=0,◻hyz=0.
(108) The equations of motion are the same as the free massless scalar field. Their solutions are
hyy(x)=∫d3→k√2ω(bke−i→k⋅→x+b†kei→k⋅→x)
(109) and
hyz(x)=∫d3→k√2ω(b′ke−i→k⋅→x+b′†kei→k⋅→x)
(110) respectively, where
(bk,b†k) and(b′k,b′†k) are the creation and annihilation operators of two physical degrees of freedom.Let us consider the mode
(bk,b†k) . For the operatorhyy , its canonical momentum is˙hyy . We obtain the solution(hyy˙hyy)=∫d3→k(2π)321√2ω(1−iω)bkei→k⋅→x+(1iω)b†−kei→k⋅→x.
(111) Let us consider a Bogoliubov transformation, which transforms the set of operators
(bk,b†−k) to the following set of operators(hyy˙hyy)=∫d3→k(2π)32(a1ka2k)ei→k⋅→x.
(112) With operator
a1k anda2k , the commutation relations ofhyy and˙hyy are expressed as below[hyy(→x,t),˙hyy(→y,t)]=∫d3→k(2π)32∫d3→p(2π)32[a1k,a2p]ei(→k⋅→x+→p⋅→y)=iδ(→x−→y),
(113) [hyy(→x,t),hyy(→y,t)]=[˙hyy(→x,t),˙hyy(→y,t)]=0.
(114) We find that apart from the dimension, the Bogoliubov 3+1 dimensional spin-2 field with the light-cone gauge is the same as that of the 2+1 dimensional U(1) gauge field. Moreover, they are both the same with the free scalar field. The form of the Bogoliubov transformation is the same as (66), with
ω=√k2x+k2y+k2z .We define the vacuum
|0⟩ bybk|0⟩=0 , and the vacuum|Ω⟩ byak|Ω⟩=0 in the 3+1 dimensional spin-2 field with the light-cone gauge. Because of the Bogoliubov transformation, the vacuum|0⟩ can be expressed by the vacuum|Ω⟩ and operator(ak,a†p) as|0⟩=1γe−∑kCka†ka†−k|Ω⟩
(115) ≃1γ∏k(1−Cka†ka†−k)|Ω⟩,
(116) where
γ is the normalization factor, and the coefficient Ck is fixed byCk=ω−αω+α,
(117) with
ω=√k2x+k2y+k2z .We can consider a system with a finite number of sites. We can divide the system into two parts, A and
ˉA . The region A we choose can be any subregion of the total system. Considering the state|Ψ⟩=|0⟩ , we can obtain the reduced density matrixρA and the second Rényi entropy S2, which have the same form of the cases in previous sections. The only difference is the expression of Ck.In the 3+1 dimensional spin-2 field with the light-cone gauge, the coefficients Ck also control the range of the interaction. The factor
∑kCkfknl can be converted to an integral by taking a continuum limit with an IR cut-offϵ and UV cut-off KFnl=∑kCkfknl≃∫Kϵ∫Kϵ∫Kϵdkxdkydkzω−αω+αei→k⋅(→n−→l)a
(118) ≈F(x,y,z)=∫Kϵ∫Kϵ∫Kϵdkxdkydkzω−αω+αei(kxx+kyy+kzz),
(119) with
ω=√k2x+k2y+k2z . Numerical evaluation in the special case of x = y = z gives the result shown in Fig. 4.Figure 4. (color online) Plots of integral (119) in 3+1 dimensional gravity with the light-cone gauge in the special case of x = y = z. Both the real and imaginary parts of the integral decrease quickly with distance, and the long range contribution can be neglected.
ϵ=0.1 , K = 10 andα=1 are used in this plot.The figure shows that the interaction between sites once again diminishes over large distances. The sites near the boundary contribute most significantly to the entanglement entropy. Hence, an area law is expected in this model as well.
-
In this study, we explore the entanglement of free spin-
12 , spin-1 and spin-2 fields. First, we consider the 1+1 dimensional Majorana field, which is just a pair of left and right moving fermions, and the 2+1 dimensional Majorana field. We perform the Bogoliubov transformation of their modes and express the vacuum with a particle pair state in the configuration space. Subsequently, we calculate the second Rényi entropies in the finite systems. Let us emphasize that while a Majorana Weyl fermion is well known to be non-local, a local Hilbert space can be defined when both chiralities are present. This is demonstrated explicitly in the current note. After that, we generalize the method to the 2+1 dimensional free U(1) spin-1 gauge field and the 3+1 dimensional gravitational (free spin-2) field. Because of the gauge redundancy of the higher spin field, there is no Hilbert space with a natural tensor product structure. We take the light-cone gauge for both fields and find that their Lagrangians behave like a free massless scalar field. The light-cone gauge allows simple quantization, while surrendering explicit Lorentz invariance. Nonetheless, it provides a candidate tensor product structure. The definition of entanglement entropy is dependent on both the state and the operator algebra. If the operator algebra is gauge-invariant [10], the corresponding entanglement entropy is likewise gauge-invariant. In this work and in [11], the operator algebras implicitly chosen are not gauge-invariant, such that the corresponding results follow trend. In our past study [11], we explored several different algebras and demonstrated which of those would reproduce the universal log terms found in Casini [10]. It is, however, expected that generic non-gauge invariant algebra choices, such as those considered in the current note, lead to a result that is gauge-dependent. For the U(1) gauge field, the Lagrangian behaves like a scalar field, which coincides with the trivial center case of [11]. As for the gravitational (free spin-2) field, we provide a prescription to observe the tensor product structure of the Hilbert space. Before doing so, there is no prescription of the gravitational (free spin-2) field. After we obtain their Hilbert spaces with a tensor product structure, we calculate the second Rényi entropies. This method can be helpful in dealing with the Hilbert space and entanglement of the perturbative gravitational field, i.e. weak gravitational field. In the non-perturbative regime, the structure of the Hilbert space is still not clear. In all the cases studied, we find that the entropy originates from the particle pairs across the boundary, and the area law emerges naturally.We would like to thank Prof. Yong-Shi Wu for critical and meticulous reading of our manuscript. LYH acknowledges the Thousands Young Talents Program.
![](/fileZGWLC/journal/article/zgwlc/2019/5/PIC/CPC-2019-0031-1.jpg)