Gauge choices and Entanglement Entropy of two dimensional lattice gauge fields

In this paper, we explore the question of how different gauge choices in a gauge theory affect the tensor product structure of the Hilbert space in configuration space. In particular, we study the Coulomb gauge and observe that the naive gauge potential degrees of freedom cease to be local operators as soon as we impose the Dirac brackets. We construct new local set of operators and compute the entanglement entropy according to this algebra in $2+1$ dimensions. We find that our proposal would lead to an entanglement entropy that behave very similar to a single scalar degree of freedom if we do not include further centers, but approaches that of a gauge field if we include non-trivial centers. We explore also the situation where the gauge field is Higgsed, and construct a local operator algebra that again requires some deformation. This should give us some insight into interpreting the entanglement entropy in generic gauge theories and perhaps also in gravitational theories.


Introduction
The seminal paper [1] initiated the important question of defining the notion of entanglement entropy of a gauge theory. The definition of the entanglement entropy has been heavily based on a tensor product structure in the Hilbert space. To study entanglement in configuration space, it requires a tensor product structure in configuration space. Naively, this is a natural feature in the Hilbert space, since the world we experienced around us is local, and that it gives the impression that the Hilbert space naturally factorizes as a tensor product of spaces defined locally at each point in space. There are clear subtlties when we work with gauge theories, where it is well known that gauge theories are by construction made local by including gauge degrees of freedom. The gauge constraints such as the Gauss law implies that degrees of freedom at different locations are not entirely independent, and thus a naive factorization of the Hilbert space is not possible. This is clearly a significant issue both for gauge theories and gravitational theories [2,3]. A clear understanding is thus crucial, also in ultimately formulating a theory of quantum gravity.
Until very recently, the replica trick has been the main tool employed to computing the entanglement enropy or the Renyi entropy in field theories, which in turn can be formulated as a path-integral in a conical space. This allows one to momentarily brush off issues of the Hilbert space and obtain some results -until it is realized that the issue in fact re-emerge as some ambiguities with edge modes that are localized at the entangling surface that is only recently understood [4,5]. See also [6,7].
Since the introduction of the notion of center in [1] into the discussion of entanglement entropy, it has made sense of entanglement even as the Hilbert space does not admit a factorization. Instead of directly considering the Hilbert space, one formulates the question in terms of the choice of an algebra attached to some region. There is some ambiguity in the selection of the algebra, and they can characterized by different centers. A large amount of work is inspired to understand the physical significance of these difference choices [8,9,10,11,12,13,14,15,16,17,18], and to demonstrate that these arbitrary choices can approach the same value as the continuous limit is taken in the mutual information for example [1].
The discussion in [19] is mainly phrased directly in terms of gauge invariant degrees of freedom. However, most lattice gauge theories or the field theoretic studies of gauge theories are formulated in terms of the gauge potential, and a Fock space is constructed for the gauge potential. Most of the discussions of entanglement in lattice gauge theories proceed by picking the temporal gauge A 0 = 0, and that the gauge potential lives on the links of the lattice, with gauge invariance imposed at the vertex. For example in the original paper [1] and most other references, this is essentially the choice. This suggests a natural question: have we exhausted all the subtleties posed by the non-locality of gauge theories? Do we understand the operator algebra and how they are attached to local regions in an arbitrary gauge? Notwithstanding the introduction of centers, it is still necessary that there is an approximate choice of local operators that can be associated to some region for entanglement in configuration space to be meaningful.
We therefore take the first step in this direction, and explore the quantization of a U(1) gauge theory in Coulomb gauge. It is well known that by imposing the Coulomb gauge, standard Poisson brackets have to be modified by Dirac brackets to take into account the gauge condition which in this case are second class constraints. Perhaps unsurprisingly, the gauge potential and its conjugate momenta cease to be local operators. Their Dirac commutators are non-vanishing even as the operators are separated by very large distances.
To make sense of the entanglement, it becomes necessary to construct local operators. We found a suitable construction making use of the duality relation between a vector and a scalar in 2d. We proposed a construction of local operators in terms of the now nonlocal gauge potentials, and obtained the entanglement entropy accordingly. It appears that the proposal is robust -it is insensitive to various local prescriptions at the boundary. However the log term associated to corners seem to approach that of a local scalar field in the case of a trivial center, but approaches that of a gauge field when we pick a center mimicking the electric/magnetic center.
Then we consider a U(1) theory in the Higgsed phase, and apply our prescription still picking the Coulomb gauge. We recover a mutual information that decays exponentially according to the ratio of the mass and the lattice scale. A gauge theory coupled to matter in 1+1 d has been considered in [20], although the gauge theory has no dynamical degrees of freedom and the only remnant is Gauss's law. Here we have an example in which matter interacts with dynamical gauge fields.
Our paper is organized as follows. In section 2 we review [19] briefly the computation of entanglement entropy of a quadratic theory making use of the values of the correlators and commutators. In section 3, we discuss the quantization of the U(1) gauge field in Coulomb gauge and the associated Dirac brackets. We discretize the theory in section 4, and introduce dual scalar variables, which allowed to construct truely local operators. In setion 5, we computed the entanglement entropy based on our prescription. In section 6-8, we generalize our method to the Higgsed gauge theory and computed the entanglement entropy accordingly. We look into mutual information in section 9 and finally conclude in section 10.

Entropies of Gaussian States in terms of correlation functions
In this section, we briefly review the methods described in [19]. It is demonstrated that expressions for entanglement entropy can be readily expressed in terms of commutators and correlations of some (canonical) variables in a quadratic theory. We consider the general commutation relations and correlation functions In the trivial center case, we have the entanglement entropy of region V where Θ = C −1 X(C −1 ) T P . X and P are matrices of correlation functions and C are matrix of commutators. The method can be generalized to the case with center. We consider the algebra generated by q i , p j , with i ∈ V = 1, . . . n and j ∈ B = k + 1, . . . n with B ⊂ V . We assume [q i , p j ] = 0 for i ∈ A = 1, . . . k,j ∈ B, such that q i , i ∈ A span the center of the algebra. In the case with center, the entanglement entropy is defined as [21] where S Q (B) is an average of quantum contributions and H(A) is the classical Shannon entropy. With the expressions of correlation functions and commutators, we have [19] S Q (B) = tr((Θ + 1/2) log(Θ + 1/2) − (Θ − 1/2) log(Θ − 1/2)), (2.7) Here C V B is the commutation matrix (2.1) between q i with i ∈ V and p j with j ∈ B. The classical part has the form The case of center formed by p i with i ∈ A = 1, . . . k can be analyzed in the same way, interchanging X ↔ P .

Local Operators of U (1) Gauge Fields with Coulomb Gauge
We consider the U (1) gauge fields with Coulomb gauge in 2 + 1 dimensions. The Lagrangian is and the Coulomb gauge fixing is The temporal component is non dynamical and it needs only to satisfy a constraint following from the Gauss law and the gauge constraint, relating it to the total charge. Because we consider free gauge fields with no charged matter, we can set the temporal component A 0 to be 0. We have the canonical momentum The Gauss law and the gauge constraint also imply To impose the two constraints (3.2) and (3.4), we have to consider the Dirac bracket. A detailed discussion can be found in [22]. The commutators and Since we consider a 2+1 dimensional theory, there is only one degree of physical freedom in Maxwell fields, also only one polarization. The mode expansions for gauge fields and their canonical momenta are given by and π j ( x, t) = 1 2π To satisfy the commutator (3.5), the polarization e i ( k, σ) have to satisfy Interestingly, by modifying the brackets by the Dirac method, we have made some rather drastic change to the tensor product structure of the Hilbert space. The gauge potentials and their conjugate momenta, even before applying the Gauss's constraint, can no longer be considered as a local degree of freedom. i.e. The operators A i and π j are not local i.e. From (3.5), one can see that their commutators are not local. They remain non-vanishing even though the fields are separated by large distances. To discuss entanglement entropy, we need to recover a basis of local operators. To our knowledge, we are not aware of a standard method of defining a suitable set of basis in such a situation. We therefore propose the following. Consider the operatorŝ The two constraints of the new variablesÂ i and π j are and We find that the commutators of new operatorsÂ i and π j are and We can see that the operatorsÂ i and π j are local. We will consider the duality of them in the lattice and calculate the entanglement entropy.

U (1) Gauge Fields with Coulomb Gauge duality in the lattice
In (2+1) dimensional U (1) gauge fields, the polarization constraint (3.9) implies a solution where ω = | k| in this section. From this solution, we can see that the Maxwell fields are dual to a scalar field χ in a fixed time slice in the two dimensional spatial slice. The duality is written as giving the following identifications 3) From (3.7) and (3.9), we have the mode expansion of χ χ( x, t) = 1 2π Therefore, it follows that Let us definê For the local operatorsÂ i and π j , we have similar relationŝ and π j = jk ∂ kχ . (4.12) Now we discretize the model in a square lattice. We define the operatorsÂ 1 and π 1 associated to horizontal links,Â 2 and π 2 to vertical links, as shown in figure (1). For example, we haveÂ 1(ij,i+1j) associated to a horizontal link andÂ 2(ij,ij+1) to vertical link, where (ij, i j ) are coordinates of the initial and final points of the links. For simplicity, we label them with the initial vertex of the vector,

13)
A 2ij ≡Â 2(ij,ij+1) . (4.14) The discrete version of (4.11) and (4.12) is also shown in figure (1). The operatorsÂ and π are related to the differences of the scalar field operatorsχ andχ in the orthogonal direction in the dual lattice respectively, such aŝ and Because there are redundant degrees of freedoms in gauge fields, we have two constraints (3.11) and (3.12) with operatorsÂ i and π j . In the discrete lattice, the two constraints become where the sum is over all the links (ab) with the common vertex a. In the above equations, it is assumed that the field component is the corresponding one to the link direction. The links have orientations, which changes the field attached to it when changing the orientation, such asÂ ab = −Â ba . With the above dualities, the non-zero commutators of the discrete version of operatorŝ A and π are We can see that the discrete version of operatorsÂ and π are almost local. We use these operators to calculate the entanglement entropy in section 5. From (4.8) and (4.7), the vacuum correlation functions of operatorsχ andχ are found to be The vacuum correlation functions of the discrete version are The scalar field operatorχ is in the center of the plaquette, and the gauge field operatorÂ in some link is equal to a difference of scalar fields across the link in the dual lattice which is perpendicular to the one corresponding tô A. The duality of gauge field momentum operator π is in the same way. (4.28) The vacuum correlation functions of discrete variablesÂ ij and π ij can be expressed with the above correlation functions, such as

Entanglement Entropy of two dimensional lattice gauge fields with Coulomb gauge
To calculate the entanglement entropy of some "region" in the lattice, we have to choose an algebra of local operators to define the "region". In the case of gauge fields, the gauge fields operators are associated to the links. We study four possible choices of algebras, which are shown in figure (2), figure (3) and figure (4). The four choices of algebras are full trivial center A, trivial center with one physical degree of freedom removedÂ,Â center AÂ and π center A π . In figure (2), we illustrate graphically the different algebras also in terms of the gauge potential and also their dual scalar variables.
Due to the redundancy of degrees of freedom in the gauge potentials, we have to remove any remaining unphysical degrees of freedom. From the two constraints (4.19) and (4.20), we have to remove one degree of freedom with a vertex. That is, for every vertex, we have to remove one link connected to it. Note that it is not possible to keep all the external links of any regions, since there is an overall constraint. They are not independent.
In the trivial center choices, as shown in figure (3), both operatorsÂ and π are associated to every link in the figure. We keep the same number ofÂ and π. In the left figure, we keep all the physical degrees of freedom, which is obtained by removing the unphysical degrees of freedom of gauge potential in the left panel of figure (2), while in the right figure, we remove a physical degree of freedom.
To get theÂ center choice, we remove the unphysical degrees of freedom of the gauge field in the right panel of figure (2), in which we remove all the operators π associated to the boundary links, and then remove a degree of freedom, as shown in figure (4). In such a choice, all the operators associated to the boundary links commute with the rest of the operators on the algebra. Hence, they form a center. In the π center choice, we do it in the same way, interchangingÂ ↔ π. Now let us see the results of four algebra choices. We expect the entropy to take the following form as a function of the square region size n, The entanglement entropy with algebra A is shown in figure (5). We have the coefficients The entanglement entropy with algebraÂ is shown in figure (6). We have the coefficients We can see that the entanglement entropy of the two different trivial centers are quite close. The effect of one degree of freedom is very small and can be neglected when the region becomes large. To calculate the entanglement entropy of gauge fields with nontrivial centers, we have to remove one physical degree of freedom after gauge fixing. The results are as follows. c)Â center − algebra AÂ Links with solid lines mean both the corresponding operatorsÂ and π belong to the algebra. Links with dashed lines mean the corresponding operatorÂ or π does not belong to the algebra. Marked dots correspond to both the scalar field operatorsχ andχ. Circle dots mean the scalar field operatorχ orχ does not belong to the algebra. The left panel shows the trivial center choice, while the right panel is related to theÂ center or π center choice (but not exactly), according to the meaning of dashed lines.
The entanglement entropy with algebra AÂ is shown in figure (7). We have the coefficients The entanglement entropy with algebra A π is shown in figure (8). We have the coefficients c 0 = 2.59095, c 1 = 3.2852, c log = 0.460412, c 2 = 0.503424, c 3 = −0.0732796.
(5.5) We see that, in the two trivial center choices, the logarithmic coefficients are very close, and they are also close to the logarithmic coefficient of entanglement entropy of massless scalar field in [19]. While, in the two non-trivial centers, the logarithmic coefficients of their entanglement entropy are also very close. They are close to the logarithmic coefficient of entanglement entropy of gauge field in [19]. In [19], they calculate the entropy from the gauge invariant electric and magnetic fields, while we from the perspective ofÂ and π. The final results are very close.

Local Operators of U (1) Gauge Fields coupling with matter
Now we consider U (1) gauge fields coupling with matter. We consider the situation of a Higgsed U(1) theory, where the gauge fields gain a mass. The Lagrangian considered is We take the gauge fixing ∂ i A i = mσ. We denote such an algebra as AÂ. When the dashed links correspond to operators π, notÂ, we get the algebra of π center. we denote it as A π . From the gauge fixing, Gauss's law reduces to  With the gauge fixing (6.2), we have the equations of motion  and π σ = ∂L ∂σ =σ − mA 0 =σ. (6.8)

4)
With A 0 = 0, we get two second class constraints and χ 2 = ∂ i π i − mσ = 0. (6.10) Because they are the second class constraints, we have to consider the Dirac bracket. The Poisson bracket is defined as .
For the two constraints, we have The Dirac bracket is defined as where the matric C is defined as For other Poisson brackets, we have {σ(x), χ 1 (y)} P = 0, (6.15) {σ(x), χ 2 (y)} P = −m{σ(x),σ(y)} P = −mδ( x − y), (6.16) {σ(x), χ 2 (y)} P = 0, (6.18) From the above equations and eq(6.13), we have and {σ( x, t), σ( y, t)} D = {π σ ( x, t), π σ ( y, t)} D = {A i ( x, t), A j ( y, t)} D = {π i ( x, t), π j ( y, t)} D = 0. (6.27) To quantize the fields, we canonically quantize, and impose the Dirac brackets to obtain the following commutators and [σ( x, t), σ( y, t)] = [π σ ( x, t), π σ ( y, t)] = [A i ( x, t), A j ( y, t)] = [π i ( x, t), π j ( y, t)] = 0. (6.32) Here we know, there are 2 degrees of physical freedom in the total fields in 2+1 dimension. For simplicity, we treat σ as A 0 . We do the mode expansion and we get (6.34) Here ω 2 = | k| 2 + m 2 . From the above Dirac brackets, we have the constraints for the polarizations However, from the commutators above, we find that the operators σ, A i and their canonical momentums are again not local. As before, we need to construct a set of local operators in this model. We defineÃ and consider the operatorsσ,Ã i and π σ ,π j . The two constraints of the new variables are The commutators of the new operators are and We can see that the new operators are local. We will discretize them in the lattice and calculate the entanglement entropy. In the next section, we will find that the new operators are very useful when we consider the lattice.

U (1) Gauge Fields coupling with matter on the lattice
For U (1) gauge fields coupling with matter σ, we have to modify the duality relation (4.2). and which satisfies the constraints (6.9) and (6.10) automatically. We have the mode expansions of the dual scalar fields and the time derivative of them (7.8) Similarly, for the new operatorsσ,Ã i and π σ ,π j , we have the duality relatioñ where we definẽ (7.14) The non-vanishing commutators ofχ,χ,φ andφ are The non-vanishing vacuum correlation functions are 17) Figure 9: Dual lattice: The scalar field operatorχ is in the center of the plaquette and operatorφ is in the vertex. The gauge field operatorÃ in some link is equal to a difference of scalar field operatorsχ across the link in the dual lattice which is perpendicular to the one corresponding toÃ plus a difference of scalar field operatorsφ along the same link in the dual lattice. The duality of gauge field momentum operator π is in the same way.
Now we discretize the model in a square lattice. We define the operatorsÃ 1 and π 1 associated to horizontal links,Ã 2 and π 2 to vertical links,σ and π σ to the vertices. Because of the redundant degrees of freedom in this model, we have two constraints (6.41) and (6.42). In the discrete lattice, the two constraints become bÃ ab = mσ a (7.21) and bπ ab = mπ σa , (7.22) where the sum is over all the links (ab) with the common vertex a. In the above equations, it is assumed that the gauge field component is the corresponding one to the link direction and matter field component is associated to the vertex. The links are oriented. The field attached to changes sign if the orientation is flipped i.e.Ã ab = −Ã ba .
Because of the above constraints, we fix the scalar field variablesσ and π σ to make gauge fields the physical degrees of freedom. As shown in figure (9), we don't show the scalar field variableσ there because we have fixed it. The way to label the gauge field operators on the lattice is the same as the pure gauge theory in section 4. We label them with the coordinates of the initial vertices of the vectors, such as The discrete version of (7.9) and (7.11) is also shown in figure (9). The operatorÃ is related to the difference of the scalar field operatorsχ in the orthogonal direction in the dual lattice and the difference of the scalar field operatorsφ at the two end vertices of the same link asÃ. The dual of operator π is defined in the same way, interchangingÃ ↔ π, χ ↔ċhi andφ ↔φ. For example, With the above dualities, the non-vanishing commutators of the discrete version of operatorsχ,χ andφ,φ are dk y (4 sin 2 k x 2 +4 sin 2 k y 2 ) cos(k x (i−k)) cos(k y (j−l)), (7.29) and [φ ij ,φ kl ] = im 2 δ ik δ jl . (7.30) Note that in the second term on the rhs of (7.29), we have left the result inside an integral. It might at first sight lead one to question the locality ofχ. We note however that these terms are precisely corresponding to the derivatives of the delta function. They are thus local on the lattice. We leave them in the current form for simpler manipulation in the numerics. The commutators of the discrete version of operatorsÃ, π can be expressed by the above commutators. The non-vanishing ones are and We can see that the discrete version of operatorsÃ and π are local. We use them to calculate the entanglement entropy in section 8. The non-vanishing vacuum correlation functions of operatorsχ,χ andφ,φ are (4 sin 2 k x 2 + 4 sin 2 k y 2 ) cos(k x (i − k)) cos(k y (j − l)), (4 sin 2 k x 2 + 4 sin 2 k y 2 ) cos(k x (i − k)) cos(k y (j − l)), (7.38) The vacuum correlation functions of discrete variablesÃ ij and π ij can be expressed with the above correlation functions, such as

Entanglement Entropy of two dimensional lattice gauge fields in the Higgs phase
We calculate the entanglement entropy of gauge fields coupling with matter in a square region now. We consider the square region with four different algebra choices, which are shown in figure (11). The four choices of algebras are full trivial center B, trivial center with some physical degrees of freedom removedB,Ã center BÃ and π center B π .
In figure (10), we illustrate the different algebras in the different duality frames. In the top figures, we show the gauge field operators and matter field operators on the lattice, while in the bottom figures, we show the corresponding dual scalar field operators. The figure (10) is without gauge fixing. We remove the redundant degrees of freedom by fixing the variablesσ and π σ . By gauge fixing, we get the full trivial center algebra from the left panel of figure (10). After fixing the matter degrees of freedom in the right panel of figure (10) and remove some operators, we can get theÃ center algebra or theπ center algebra, as shown in figure (11).
In the trivial center choices, as shown in the first and second figure of figure (11), both operatorsÃ and π are associated to every link. We have the same number of operatorsÃ and π. In the first figure, we keep all the physical degrees of freedom, while in the second figure, we remove some links to get another algebra.
To get theÃ center choice, we remove operators π on the boundary links of the second figure in figure (11) which leads to the third figure. We find that all the operatorsÃ associated to the boundary links commute with the rest of the operators in the algebra. Hence, they form a center. In the π center choice, we do it in the same way, interchangingÃ ↔ π. Now let us see the results of fours algebra choices. We also expect the entropy has the following form as a function of the square region size n, S n = c 0 + c 1 n + c log log n + c 2 n + c 3 n 2 .  (14). We can see that the leading term of entanglement entropy is different for two trivial centers, but from the fitting coefficients, we find that the logarithmic term is very close. In the middle figure, we remove some links. In the remaining links, both operatorsÃ and π are there. It also forms an algebra of trivial center. We denote it asB. In the right figure, links with solid lines mean both the corresponding operatorsÃ and π belong to the algebra and links with dashed lines mean the corresponding operatorÃ orπ does not belong to the algebra. Without loss of generality, let us assume operatorsÃ to be on the dashed links, not π. The operatorsÃ on dashed links commute with all the operators π on the solid links and they form a center. We denote such an algebra as BÃ. When the dashed links correspond to operators π, notÃ, we get the algebra of π center. we denote it as B π .  We can see that the logarithmic term is independent of the algebra choices.  The mutual information is monotonously increasing with the algebra. Hence, we expect to have where A and B are the two square regions.
In figure (17), we show the numerical calculation of mutual information of gauge fields with Coulomb gauge between two square regions of the same area n 2 and separated by a distance n for different algebra choices. The figure shows that the relation (9.3) holds, as we expect. The mutual information for algebra A andÂ is difficult to distinguish when n ≥ 4. When the regions become large, the effect of one pair of operators on the entropy and mutual information is negligible. This is because the algebraÂ gets closer and closer to algebra A, as the regions become larger and larger. The mutual information of gauge fields with Coulomb gauge between two square regions of size n and separation distance n for four different algebra choices, from top to bottom A andÂ (purple curve), AÂ (red curve) and A π center (green curve). There are two curves in the purple curve. One is mutual information of algebra A and the other is of algebraÂ. At n = 3, the mutual information of algebra A is a little larger than algebraÂ. For other values of n, they are hard to distinguish.
10 Mutual Information of 2D lattice gauge field coupling with matter Now we consider the mutual information of gauge fields coupling with matter between two squares as in the previous section, again for the four different algebra choices in figure (11). In figure (11), we have the four algebras The mutual information is monotonously increasing with the algebra. Hence, we expect to have where V and W are the two square regions.
In figure (18), we take m = 1 in the calculation. We show the numerical calculation of mutual information of gauge fields coupling with matter between two square regions of the same size n and separated by a distance n for different algebra choices. The figure shows that the relation (10.2) holds, as we expect. The mutual information of all four different algebra choices decays to 0 very quickly. This is one of the properties of gauge fields in the Higgs phase, the decay reflecting the mass gained. We take m = 1 in this figure. In the limit n → ∞, the mutual information for all four algebra choices goes to 0 very fast.
Now we compare the mutual information of gauge fields coupling with matter of different masses. We consider the same algebra choices and regions as the above case of m = 1. As the mass of matter increases, the vacuum correlation functions ofÃ and π decays faster with distance. We expect that the larger the mass, the smaller the mutual information. We compare the results of m = 1 2 , 1, 2 for four different algebra choices in figure (19), figure (20), figure (21) and figure (22) respectively.      (22) show the mutual information with different mass for four algebra choices B,B, BÃ and B π respectively. We can see that, for the four different algebra choices, the mutual information all decreases as the mass of matter increases.
In this paper, we have explored the effect of a different gauge choice on the tensor product structure of the Hilbert space constructed for the gauge potential in a U(1) gauge theory in 2+1 dimensions. In particular, departing from the usual choice of temporal gauge, we studied the Coulomb gauge and explicitly demonstrates that by imposing the constraints as second class constraints and obtaining Dirac brackets, the gauge potentials attain commutators that correspond to highly non-local operators-commutators are non-vanishing all the way to infinite separation. To make sense of entanglement entropy associated to a local region, it is therefore necessary to construct a local operator algebra. We proposed a way making use of the duality with a scalar degree of freedom in the two dimensional spatial slice, in the cases of the massless theory and also in the Higgs phase. In the massless case, we recover the entanglement entropy of a scalar field (the log term in particular) and that of a U(1) theory when we pick a center. We also studied the case of the Higgs phase in the Coulomb gauge and construct a modified set of local operators accordingly. In this case, the choice of algebra has little detectable effect on the log term of the entanglement entropy.
We have shown, at least in 2+ 1 dimensions how to recover known results in terms of the gauge potential when the gauge potentials are non-local. This gives us some insights into the structure of the Hilbert space in a gauge theory. In terms of the gauge potentials, the Hilbert space is in some sense fluid -local operators could turn into non-local ones based on different gauge choices. However, suitable construction of local degrees of freedom appear to recover the expected result. This should have implications in gravitational theories as well. We also note that our construction of a set of local operators is based on a duality relation with scalars in 2 spatial dimensions. It would therefore be interesting to generalize our construction to higher dimensions.