Arxiv:1905.00652V2 [Quant-Ph] 16 Jul 2019 Nteoehn Eitoueauiaytransformation Unitary a Introduce We Hand Speciﬁcally, One Locality

Home , Lattice gauge theory

arXiv:1905.00652v2 [quant-ph] 16 Jul 2019 nteoehn eitoueauiaytransformation unitary a introduce we hand Specifically, one locality. the keeping eliminate on while completely matter to fermionic us the allows bosons, hard-core into ihorrcn ok[] hr esoe htfermionic that showed we in where matter [7], work that recent Combining our hard- configuration. with by staggered represented a is in matter bosons the core where LGTs to [3] Ref. [6]. fermions representation staggered Hamiltonian with the 5] in [4, (LGTs) theories gauge tice iia rcdr a edvsdas oeiiaethe eliminate a to that of show also context the will devised in we matter be work fermionic can well-known this the procedure In simply similar is fixing. it trans- gauge of unitary unitary care and The takes constraints, fields. where that local matter field, [3], formation the the gauge version of all the part removing lattice by eliminating mass, its a absorbed in it are as giving mech- bosons Higgs well Goldstone the the- as of the the part 2] of important [1, an locality is anism the this keeping fact, still In terms in ory. while fields others, the the of lo- some of the expressing solving and simply by constraints partially, cal least at lifted, redundancy be this could whether question This under the invariant physical. raises considered are immediately are that - fields transformations gauge the the of quanti- combinations invariant - gauge ties Only description. mathematical and matter con- the local freedom. both to of by rise degrees giving satisfied gauge space, laws) motion of (Gauss of point straints lo- constants each a many at is lo- involves It defined mediate that particles. - symmetry, matter The bosons cal the gauge between physics. - interactions particle fields cal of gauge model of standard inter- excitations and the forces in fundamental the actions describes it as physics, nodrt banti eut eetn h ehdof method the extend we result, this obtain to order In ag ymtyas mle eudnyi the in redundancy a implies also symmetry Gauge in important very is invariance gauge of principle The eoigSagrdFrincMte in Matter Fermionic Staggered Removing uihCne o unu cec n ehooy(CS) S (MCQST), Technology and Science Quantum for Center Munich aa nttt fPyis h erwUiest fJerusa of University Hebrew The Physics, of Institute Racah SU a-lnkIsiu ¨rQatnpi,Hans-Kopfermann für Quantenoptik, Max-Planck-Institut ere ffedm ttecs flclyetnigteinte finiteness the the to extending due locally constraints, of local fermion cost new work, staggered the introducing this with at I and In theories freedom, gauge freedom. of re calculation. lattice of degrees that any degrees in constraints in laws gauge local and Gauss addressed laws, matter and Gauss the the of relate through care phrased taken is be must that a-lnkIsiu ¨rQatnpi,Hans-Kopfermann für Quantenoptik, Max-Planck-Institut ( ag hois hog h oa ymtywihi ntheir in is which symmetry local the through theories, Gauge N .INTRODUCTION I. and ) U ( N uihCne o unu cec n ehooy(MCQST), Technology and Science Quantum for Center Munich Gscnb transformed be can LGTs ) celnsr ,D879Muce,Gray and Germany, München, D-80799 4, Schellingstr. U ( N and ) SU Dtd uy1,2019) 18, July (Dated: ( N .IncoCirac Ignacio J. lat- ) rzZohar Erez U ag edi o laspsil;te epoedt a to proceed we then the possible; the for arbi- always of it demonstration in not of fields is solution complete field matter a gauge for while solved - be dimensions trary could law Gauss the but invariance, gauge local constrains. local breaks other local, elim- introduces degrees is is gauge to contains freedom, latter only of converted the theory resulting secondly, is The and matter inated. one transfor- fermionic bosonic two the hard-core of firstly, consists ones. de- procedure mations: bosonic fermionic complete hard-core the by the broken replaced which Thus, not are in is freedom subgroup laws of symmetry Gauss finite grees gauge obtains a the one only group, and ex- uses gauge locally, procedure the but this a of slightly, As As is field. range tended. gauge interaction the the to result, statistics the ones, hand, transferring bosonic (hard-core) other while the by the freedom replaces of On [7] degrees satis- in fermionic arise. introduced constraints alone transformation local unitary fields other as the gauge [3], laws the case Gauss by field the fied scalar solving the explicitly sym- by in gauge original broken the is though space fields, even metry Hilbert Thus, gauge local bosons. the the the of for of the finiteness constraint the for from local state stems theory. new which product a the trivial the and of a in matter, locality enforces fields the it gauge Furthermore, keeping the from while matter Hamiltonian, the decouples that the egnrlz h rcdr to procedure demonstrations, these the After com- and generalize [7]. by we statistics) of fields results fermionic the matter with no bining fermionic but staggered (that to spaces matter eventually bosonic local hard-core finite to it has extend [3], of lattice eetdmti n etridcs nessae differ- stated unless indices, vector and matrix repeated fermions. staggered ag hois ihany with theories, gauge ( h oki raie sflos is,w ru that argue we First, follows. as organized is work The hogotti okw suetesmainof summation the assume we work this Throughout N U ) Srß ,878Grhn,Gray and Germany, Garching, 85748 1, -Straße 1 ntr ag o h ope clrfil nthe on field scalar complex the for gauge unitary (1) e 10,GvtRm euae,Israel. Jerusalem, Ram, Givat 91904, lem and Srß ,878Grhn,Germany, Garching, 85748 1, -Straße fteoiia oa atrspaces. matter local original the of ato ag,bekn h symmetry the breaking range, raction hligt.4 -09 München, Germany, D-80799 4, chellingstr. o opeeyrmvn h matter the removing completely for s tittepyia ibr pc and space Hilbert physical the strict oe xii aylclconstraints, local many exhibit core, epeetawyta ssalthe all uses that way a present we SU h aitna itr this picture Hamiltonian the n ( N ) U atc ag Theories Gauge Lattice 1 ae trigwt eiwof review a with starting case, (1) > N ,culdt fundamental to coupled 1, U ( N and ) SU ( N lattice ) 2 ently. show how to use that for completely eliminating matter fields in particular lattice gauge theories.

II. THE GAUSS LAW B. Lattice Gauge Theories A. The Classical Case The physical Hilbert space of a quantum gauge theory, Already in the context of classical electrodynamics, the phys, is contained in the product of the Hilbert spaces ofH the gauge field, , and the matter, : gauge fields and matter are related through the Gauss Hgauge Hmatter law, that has nothing to do with quantization. It is given (3) by the equation Hphys ⊂ Hgauge ⊗ Hmatter It is not equal to the product, as the gauge field and mat- E (x)= ρ (x) (1) ∇ · ter degrees of freedom are connected through the Gauss law, that becomes (for electrodynamics) the eigenvalue E x stating that the divergence of the electric field ( ) at equation any space point x is equal to the local charge density ρ (x) associated with the matter. It is one of Maxwell’s E (x) ψ = ρ (x) ψ (4) equations, that is obtained as an equation of motion in ∇ · | i | i the Euler-Lagrange formalism; nevertheless, it is a static From now on, we will focus on lattice gauge theories equation that includes no time derivatives - a set of local [4, 5]. In the Hamiltonian picture [5] time is continuous, constraints. Indeed, when one uses instead the Hamil- and the degrees of freedom reside, in d + 1 dimensions, ton formalism that does not treat time and space on an on a Zd lattice. We will begin our demonstration with equal footing, it does not appear as an equation of mo- U(1) gauge fields. tion anymore (but rather as a constraint that is added The matter, either bosonic or fermionic, resides on the to the Hamiltonian with a Lagrange multiplier after the vertices. Later on we will focus on particular types of Legendre transformation). matter, but now we it will be enough for us to define The question is, then, whether we could solve the con- charge operators Q (x) at each vertex x Zd. In the U(1) straints and use them for reducing the number of degrees case it will take integer eigenvalues, either∈ bounded or of freedom in our system by completely eliminating ei- not, depending on the type of matter used. Gauge fields, ther the gauge field or the matter. If it is the gauge field on the other hand, reside on the lattice’s links (x,i), ema- that we wish to eliminate, we should solve (1) for E (x). nating from the vertex x in the direction i =1, ..., d. The In this case it is a differential equation, and it can only local Hilbert space of a U(1) gauge field on a link is this be integrated in one space dimension where it becomes a of a particle on a ring: the role of the vector potential simple first order equation, ∂E/∂x = ρ (x), giving rise to will be played by a compact variable φ (x), canonically a nonlocal solution, conjugate to the electric field operator - the ”angular mo- mentum” E (x,i), with an unbounded integer spectrum. x Thanks to the canonical relation E (x)= dx′ρ (x′) (2) Z [φ (x,i) , E (y, j)] = iδij δ (x, y) (5) This is a well known solution that may be generalized (where both deltas are Kronecker’s) the group element to non-Abelian groups and quantum cases as well, but operator defined by not beyond 1+1d: in more dimensions the electric fields x has more components and there are just not enough con- U (x)= eiφ( ,i) (6) straints (constants of motion) that can be used for integration (the electric field is a rotation-vector, and the is simply an electric field raising operator, equations are rotation-scalars).

The solution for the matter, on the other hand, is [E (x,i) ,U (y, j)] = δij δ (x, y) U (y, j) (7) completely diﬀerent, as the charge density is a rotation- scalar, and (1) is just a simple algebraic equation for it, The pure-gauge part of the Hamiltonian of such theo- already explicitly and locally solved. In this work we will ries is usually the Kogut-Susskind Hamiltonian,

g2 1 H = E2 (x,i) cos(φ (x,i)+ φ (x + e , j) φ (x + e ,i) φ (x, j)) (8) KS 2 − g2 i − j − x,i x,i

where ei is a unit vector in the i direction. The interaction of the matter with the gauge field takes 3 the form physical state may be expanded in a more restrictive way. For that, we define the charge states of the matter - eigen- Hint = ǫ (x,i) a† (x) U (x,i) a (x + ei)+ h.c (9) states (not necessarily unique) of the charge operators x,i Q (x) , X { } where a† (x) is bosonic or fermionic, such that Q (x) q = q (x) q (17) |{ }i |{ }i Q (x) ,a† (y) = δ (x, y) a† (x) (10) Similarly, we can define (non-unique) eigenstates of the electric field divergence operators D (x) (not unique) as This can be added to some free matter part, HM , and altogether H = HKS + HM + Hint. D (x) d = d (x) d (18) Gauge invariance is the invariance under transforma- |{ }i |{ }i tions generated by the local generators The gauge invariant states may be expanded as

(x)= (E (x,i) E (x e ,i)) Q (x) (11) G − − i − i ψ = f ( q, d ) δ x x d q X | i { } d( ),q( ) |{ }igauge⊗|{ }imatter q,d " x # {X} Y These operators commute with the Hamiltonian, (19) The notations q and d are, as mentioned above, [H, (x)]=0 x Zd (12) |{ }i |{ }i G ∀ ∈ generally non-unique: different local matter configura- - a local symmetry, or a set of local constraints. A physi- tions, corresponding to different quantum states, may give rise to similar local q eigenvalues of the charge cal state ψ is gauge invariant, that is, it is an eigenstate { } of all the| locali generators (x), operators Q , and similarly with the divergences of the G electric field{ }D and their eigenvalues d (which is the { } { } (x) ψ = λ (x) ψ (13) usual case, in which more indices can be added); the no- G | i | i tations above are merely illustrative and their accurate and the commutation relations (12) imply that states details will be later discussed, here the important thing with different λ (x) are not connected by the dynamics, we wanted to emphasize is that the physical Hilbert space and give rise to{ disconnected} sectors, is not spanned by all the product states).

phys = phys ( λ (x) ) (14) H H { } D. Solving the Gauss Laws in a Lattice Gauge [ The eigenvalue equation (13) could be rewritten as Theory

(E (x,i) E (x e ,i)) ψ = (Q (x)+ λ (x)) ψ After having described the Hilbert space of a lattice − − i | i | i i gauge theory, and the implication of the local symme- X (15) tries on its structure, we are ready to see how to use - which we recognize as the Gauss law: the divergence of the Gauss laws, or their explicit solutions, for simplify- electric fields at a vertex equals the charge there, which is ing the description and reduction of redundant degrees composed of the dynamical charge - the operator Q (x), of freedom. and the eigenvalues λ (x) which we can now recognize as We wish to discuss particular types of unitary trans- static charges. From now on we shall choose λ (x)=0 formations, that take a gauge invariant state ψ as in everywhere - that is, restrict our physical Hilbert space (19), satisfying (16), and completely eliminate either| i the to the sector with no static charges, but the results may gauge field or matter degrees of freedom, while conserv- be generalized in a straightforward way also to any other ing the physical information - the amplitudes of elements charge sector. Defining the local electric field divergence in superposition. This will be done by using the gauge operator as D (x) = (E (x,i) E (x ei,i)), we can symmetry and will, eventually, break it. One can con- − − i sider, obviously, other types of transformations that leave now rewrite the GaussP law we will use from now on as some ingredients of the field and do not remove them D (x) ψ = Q (x) ψ (16) completely, but we will not discuss such transformations | i | i here. A complete removal of the gauge field would be done C. Gauge Invariant States by a transformation , as follows: W ψ = 0 f ( q ) q (20) A general state in gauge matter could be expanded W | i | igauge ⊗ { } |{ }imatter H ⊗ H q as A (g,m) g gauge m matter, using some arbitrary X{ } g,m | i ⊗ | i basesP g , m of the gauge field and matter Hilbert spaces where 0 gauge is some ”empty” gauge field state that is respectively.| i | i However, thanks to the Gauss law (16), a factored| i out. 4

This transformation, in the cases it can be defined, is and change the local charges as follows: a controlled operation, making use of the symmetry and uq q = 0 (25) the local set of constraints. One locally transforms the | i | i degrees of freedom that are to be decoupled, controlled for each integer q = 0. Since u is unitary, negative values by those that remain. This requires the solution of the involve hermitian6 conjugation. Gauss law that was previously discussed; for that rea- We still have to check in which cases such transforma- son, the transformation is only possible, in general, W tions exist; but when they do, in order to eliminate the for one spatial dimension, where the gauge field can be charge, we will have to act on each vertex with a power integrated, using the quantum, lattice analogue of (2). of this transformation, that is exactly equal to the initial Just like there the divergence of a vector quantity was amount of charge that was there. We can do it thanks replaced by a simple differential equation, here we ob- to the Gauss law; let us define a controlled, local unitary tain a simple difference equation, that can be inverted transformation (x) that lowers the charge at x by an non-locally, since the number of links (electric fields) is amount given byU D (x), the divergence of electric field equal to this of vertices (charges and constraints), while there. Thanks to the Gauss law (16), this is equal to the in more dimension there are not enough equations: in charge and hence, acting on gauge invariant states it will d + 1 dimensions, the links scale as d times the vertices. reduce Q (x) to zero. Let us see that explicitly; consider Therefore, Only in this case the gauge field is completely an expansion (19) of a gauge invariant state ψ . Then, redundant - in more dimensions it cannot be completely | i x removed. In 1 + 1d we define (for open boundary con- Q (x) u (x)D( ) ψ = ditions - periodic ones do not allow to remove the gauge | i D(x) field completely) f ( q ) δ x x u (x) d q = { } d( ),q( ) |{ }igauge ⊗ |{ }imatter q,d " x # {X} Y = exp i φ (x) Q (y) (21) W d(x) x y

A. Complex Scalar Matter Before we move on to other types of matter, one could ask what happens if we couple the same gauge field to First, we consider the case in which each vertex hosts a more matter fields, residing at the vertices and adding up complex scalar field, Φ (x), for which the relevant trans- to the local charges, Q (x) = Qi (x). In this case, in formation - the lattice version of the unitary gauge of the i general, a transformation of theP form will not be pos- Brout-Englert-Higgs mechanism [1, 2] was discussed by sible, since the spectrum of the localU charge operators Fradkin and Shenker in [3]. The field may be expanded Q (x) will be degenerate. Microscopical matter config- in a polar form, with a radial and angular part (phase); urations must be completely distinguishable in terms of we assume, for simplicity, that the radial part is frozen, their charges, if one wishes to decouple the matter in for example due to the Higgs mechanism. Following [3], the manner described above. The controlled operation x is based on the divergence of the electric fields which is x iθ( ) Φ ( )= R0e− (31) equal to the total fermionic charge at the vertex, with no At each vertex we have a non-bounded charge operator way to distinguish different charge contributions. Q (x) with an infinite spectrum of integers. It is raised by the operators eiθ(x), which mutually commute, and B. Staggered Hard-Core Bosonic matter therefore

iθ(x) u (x)= e− (32) Our next stop en route to fermions will be hard-core bosonic matter. In this case, the matter Hilbert space satisfies both (29) and (24), allowing us to obtain a well on each vertex is that of a spin-half particle. To be able defined using (27) and (28). Since u here is a uni- U to eventually generalize to staggered fermions, we will tary lowering operator, the Gauss laws will vanish for stagger the hard-core bosons. On each vertex we define the transformed state using (30). the ”number operator”, The complete transformation then takes the form

x x e x 1 i θ(x)D(x) i E( ,i)(θ( + i) θ( )) n (x)= (σz (x)+1) (37) − x x,i − = e P = e P (33) 2 U which, seen now as a transformation of the gauge ﬁelds and the staggered charge operators controlled by the matter, is recognized as the unitary Q (x)= n (x) s (x) (38) gauge [3]. In such theories, the interaction terms take − the form Φ (x) U (x,i) Φ (x + e ), which, after freezing † i where s (x) = 0 (1) on the even (odd) sublattice rep- the radial ﬁeld, become (if ǫ is real) resenting particles (anti-particles). The Gauss law (16) may then be rewritten as H =2R2 ǫ (x,i)cos(φ (x,i)+ θ (x) θ (x + e )) int 0 − i x,i X G (x) ψ = n (x) ψ (39) (34) | i | i Note that introducing

φ (x,i) † = φ (x,i)+ θ (x + e ) θ (x + e ) (35) U U i − i G (x)= D (x)+ s (x) (40) - therefore, after the transformation the matter field de- The parts of the Hamiltonian that depend on the mat- couples from the interaction terms, and we obtain mas- ter take the form sive gauge fields: s(x) 2 HM = M ( 1) n (x) (41) H˜ = H † =2R ǫ (x,i)cos(φ (x,i)) (36) − int U intU 0 x x,i X X and 2 thanks to the terms proportional to R0 cos(φ (x,i)). These terms explicitly break gauge invariance, that does Hint = (ǫ (x,i) σ+ (x) U (x,i) σ (x + ei)+ h.c) − not exist anymore as anticipated. The other parts of the x,i X Hamiltonian commute with the transformation and do (42) not transform (the rest are pure gauge terms, and the Unlike in the scalar case, now the charges are bounded, G (x) are generators of pure gauge transformations and because the operators n (x) are bounded. The Gauss law therefore commute with them). in the form (39) implies that the spectrum of the opera- Similar transformations are possible for other Higgs tors G (x) in the physical Hilbert space consists only of scenarios, with different groups. For example, see [14] 0, 1. Therefore we can write down extra local constraints for Z2. If the radial component of the field is not frozen, that, at this point, are completely redundant: it will not be eliminated, since it is not coupled to the gauge field and thus not subject to any local constraints. G (x) (G (x) 1) ψ = 0 (43) − | i 6

We define, on each vertex x, the operators Pg (x), which Hamiltonian terms commute with the operators G (x), it project to the subspaces where G (x) = g. The physical is sufficient to do it for the interaction part - the only part Hilbert space is contained within the subspace of g = that does not commute - and even there it is sufficient to 0, 1 everywhere, and therefore we can multiply the entire include only the most relevant local projectors. Finally Hamiltonian by the projectors P (x) = P1 (x)+ P0 (x) we obtain, using the Gauss law, that within the physical and have exactly the same spectrum and dynamics within Hilbert space, the physical Hilbert space. However, since most of the

phys Hint = (ǫ (x,i) P1 (x) σ+ (x) U (x,i) σ (x + ei) P0 (x + ei)+ h.c) (44) − x,i X The Gauss law also helps us to rewrite, within the physical Hilbert space, the mass part of the Hamiltonian as

x x Hphys = M ( 1)s( ) G (x)=2M ( 1)s( ) E (x,i) + const. (45) M − − x x,i X X

Another implication is that now, that we do not have an This implies that infinite ladder of charges, we cannot define a unitary u that lowers Q. Instead, we can define x (1 2G (x)) ψ˜ = ( 1)G( ) ψ˜ (50) − − u (x)= σx (x) (46) E E

which will satisfy (25) but not (29). and (48) simplifies to We would like to construct a local transformation (x) that eliminates the charges. Using the modified,U stag- σz (x) ψ˜ = ψ˜ (51) gered Gauss law (39), we construct the local transforma- − as expected, implying that allE the matterE degrees of free- tion dom are decoupled and in a spin-down state, or that x x (x)= u (x)G( ) = σ (x)G( ) (47) U x n (x) ψ˜ = 0 (52) That takes a spin up state (corresponding thanks to the E Gauss law to G = 1) to spin down, and leaves spin down invariant. The operators G (x) are left invariant, and us- How does the Hamiltonian transform? As before, the G(x) ing (x) σz (x) † (x) = ( 1) σz (x) we obtain that pure gauge part H does not transform at all, as it U U − KS the transformed state, ψ˜ = ψ satisfies commutes with the transformation. The mass part (45) U | i commutes as well and does not transform, and we are left E x with the transformation of the interaction part. However, σ (x) ψ˜ = ( 1)G( ) (1 2G (x)) ψ˜ (48) z − − − in this case we have to be more careful with the transfor- E E mation of the interaction part of the Hamiltonian that

From this equation, it looks as if there is still some local in the complex scalar case. The interaction term does coupling between the gauge ﬁeld and the matter. How- not commute with G (x) and thus changes their eigen- ever, this is not the case; recall the conditions (43) on values. Unlike in the scalar case, here we started with the spectrum of the G (x) operators that were redundant matter that resides in ﬁnite local Hilbert spaces, giving before the transformation; now they are not redundant rise to the local constraints (43,49), which have to be anymore, and in fact, they are invariant under the trans- incorporated into the interaction explicitly before transformation, that is forming it, otherwise we will get terms that breaks them. Thus, we transform the physical interaction Hamiltonian, G (x) (G (x) 1) ψ˜ = 0 (49) within the physical subspace (44), and obtain − E

phys phys H˜ = H † = ǫ (x,i) P1 (x) U (x,i) P0 (x + ei) x x e + h.c (53) int U int U ⊗ |↓i h↓| ⊗ |↓i h↓| + i x,i X We can now completely forget about the matter degrees of freedom as they are decoupled. The ﬁnal Hamiltonian 7 to be used in the transformed physical Hilbert space is, therefore,

phys s(x) x H † x = HKS +2M ( 1) E (x,i)+ (ǫ (x,i) P1 (x) U (x,i) P0 (x + ei)+ h.c) (54) ⊗x h↓| U U ⊗x |↓i − x,i x,i X X without any local Gauss laws but with the local con- One could be tempted to use a Majorana mode as the straints (49) that are taken care of by the projectors. local unitary, u (x)= ψ (x)+ψ† (x), and indeed it is uni- These projectors extend the range of interaction, and this tary and satisfies (25), but as it has an odd fermionic par- is the price we have to pay for having, originally, local ity (it changes the parity of the states it acts upon), the matter spaces that are finite. commutation property (24) does not hold and one can- Here, once again, we could not add multiple matter not define (28) with it, since the local terms will not have species coupled to the same gauge field, because this a well defined fermionic parity and hence their product would destroy the unique mapping between a charge con- will have to be in some fixed order, giving rise to nonlocal figuration and a matter state. This also gives a good strings. motivation for staggering: if, instead, we had two hard- In a recent work [7] we have shown that lattice gauge core bosonic species per site, coupled to the same gauge theories with staggered fermionic matter whose gauge field, we would not be able to perform such a decoupling group contains Z2 as a normal subgroup may be mapped transformation. to lattice gauge theories with hard-core bosonic matter (spins), using a local and unitary transformation that does not involve nonlocal strings. This was done as well C. Staggered Fermionic Matter thanks to the fact that in a gauge theory one has local constraints, that allow to transfer the statistics informa- Finally, we are ready to deal with staggered fermionic tion to the gauge fields. However, since in that case only matter [6]. In this case, at each vertex there is one the parity is discussed, and it has to do with the finite Z2 fermionic species, created by the operator ψ† (x). The group that is not continuous, the procedure carried out fermionic number operators, as usual, are n (x) = there did not break the symmetry and only allowed to ψ† (x) ψ (x), and with them one may define the charges replace the fermionic matter by hard-core bosons, with just like in the hard-core bosonic case (38). the only ”payment” of extra, but local, appearances of The parts of the Hamiltonian that involve the matter signs of electric fields in the Hamiltonian, that account take the form for the fermionic statistics.

s(x) Since U(1) contains Z2 as a normal subgroup, one can HM = M ( 1) n (x) (55) x − replace the fermions in a U(1) gauge theory by hard- X core bosons. Following [7], we obtain that our model is and equivalently described by the Hamiltonian H′ = HKS′ + Hint′ + HM′ , in which HM′ = HM (but with the spin n Hint = ǫ (x,i) ψ† (x) U (x,i) ψ (x + ei)+ h.c (56) operators), x,i X

2 g 2 1 H′ = E (x,i) ξp (x,i,j)cos(φ (x,i)+ φ (x + e , j) φ (x + e ,i) φ (x, j)) (57) KS 2 − g2 i − j − x,i x,i

Hint = (ξ (x,i) ǫ (x,i) σ+ (x) U (x,i) σ (x + ei)+ h.c) − x,i X (58) where ξp (x,i,j) , ξ (x,i) are local functions of the electric ﬁelds on links that belong to the plaquette/link they are associated with and its neighboring links [7]. Since these After having stated the procedure that allows one to functions commute with G (x), one may use the proce- eliminate the matter of U(1) lattice gauge theories, with dure introduced for staggered hard-core bosonic matter either hard-core bosonic or fermionic matter, we will now to eliminate the matter in this fermionic case as well, generalize it to non-Abelian groups. For that, let us once the fermions are converted to hard-core bosons. brieﬂy review the structure of those theories. 8

A. SU(N) and U(N) Lattice Gauge Theories there are N 1 such generators - the maximal number of N dimensional− diagonal traceless matrices. We shall N 1 1. The groups U(N) and SU(N) and their Cartan denote them by T − . Since the identity matrix triv- { µ}µ=1 subalgebra ially commutes with any other matrix, the Cartan subalgebra of U(N) will include the N 1 elements of that − We consider here lattice gauge theories whose gauge of SU(N), as well as T0 = τ0 - altogether N mutually group G is either U(N) or SU(N). We denote group el- commuting generators. ements by g G, and irreducible representations by j; unitary matrix∈ representations are given by the Wigner j matrices, Dmn (g). As these are Lie groups, the repre- 2. The gauge field sentation j is generated by a set of matrix generators, j τa , whose dimension is referred to as the representation’s Let us begin with the description of the gauge degrees dimension dim (j). For SU(N), there are N 2 1 such − of freedom. As in the Abelian case, they reside on the generators, satisfying the group’s Lie algebra links of the lattice. The gauge field on each link is described by a set of operators [5, 15]: the group element j j j j τa , τb = ifabcτc (59) operators Umn, matrices of gauge field operators that h i transform as group elements in the j representation (in all of which Hermitian and traceless matrices. We shall the fundamental representation, where we omit j, these discuss the fundamental representation, whose dimension are N N matrices), and the left and right transforma- × for SU(N) and U(N) is N, and in this case the represen- tion generators, La and Ra respectively. These are two tation index will be omitted. In this representation, one independent sets of operators that fulfill the group’s Lie usually chooses the Cartan-Weyl basis for the generators, algebra, and the normalization condition [La,Lb]= ifabcLc 1 − Tr (τ τ )= δ (60) [Ra, Rb]= ifabcRc (62) a b 2 ab [La, Rb]=0 follows. j In order to obtain the algebra of U(N) from that of that generate transformations of the Umn: SU(N), one only has to introduce one extra generator, j j j La,Umn = τa ′ U ′ 1 mm m n (63) τ0 = 1 (61) j j j R ,U = U ′ τ ′ √2N a mn mn a n n

2 The normalization is chosen in accordance with (60). These operators satisfy RaRa = LaLa J . We will con- SU(N) and U(N) (N > 1) are non-Abelian Lie groups, sider, from now on, only the fundamental≡ representation whose generators, in general, do not commute, as in (59). for the gauge ﬁeld. However, there exists a maximal subset of generators that The pure gauge (Kogut-Susskind) part of the Hamil- commute - forming the Cartan subalgebra. For SU(N) tonian takes the form [5]

2 g 2 1 H = J (x,i) Tr U (x,i) U (x + e , j) U † (x + e ,i) U † (x, j) + h.c. (64) KS 2 − 2g2 i j x,i x,i

3. The matter Out of the N fermionic species, one constructs the local charges Q (x). In the SU(N) case, there are N 2 1 such a − As usual, the matter resides on the vertices. We will charges, consider once again staggered matter, and will restrict our discussion to the fundamental representation, where on each vertex there is a spinor with N components. In normal cases it is fermionic, with creation operators ψm† (x) satisfying the fermionic algebra 2 Qa (x)= ψm† (x) (τa)mn ψn (x) a =1, ..., N 1 ψm (x) , ψ† (x) = δmnδ (x, y) − n (65) (66) ψ (x) , ψ (x) =0 In the U(N) case, one adds the extra U(1) charge which { m n } 9 is explicit, law is

D (x) ψ = Q (x) ψ (71) N a | i a | i Q0 (x)= ψm† (x) (τ0)mn ψn (x) s (x) − 2 - R ,L play the role of non-Abelian, right and left elec- r (67) a a N tric fields. = n (x) /N s (x) 2 m − r m ! X B. Eliminating the matter in the U(N) case The charges Qa satisfy, on each vertex, the group algebra. The parts of the Hamiltonian that involve the matter As in the U(1) case, for the elimination of matter we take the form will first convert it to a hard-core bosonic form, following [7]. The mapping to bosons is possible in the U(N) case, s(x) H = M ( 1) n (x) (68) as it includes Z2 as a normal subgroup. The fermions M − m x,m are mapped to spinors ηm† (x), whose components anti- X commute on-site, and ηm (x) , η† (x) = δmn n (72) ηm (x) , ηn (x) =0 Hint = ǫ (x,i) ψm† (x) Umn (x,i) ψn (x + ei)+ h.c { } x,i X (69) but commute between different sites,

ηm (x) , ηn† (y) = [ηm (x) , ηn (y)]=0 (73) for x = y. 4. Gauge invariance Number6 operators become

2 Gauge invariance now is non-Abelian: there are N 1 nm (x)= ηm† (x) ηm (x) (74) local generators of gauge symmetry at each vertex− for SU(N), and one more for U(N). As the groups are non- and the charges are

Abelian, these do not commute. One defines the non- 2 Q (x)= η† (x) (T ) η (x) a =1, ..., N 1 Abelian divergence of electric fields by a m a mn n − N Q (x)= n (x) /N s (x) D (x)= (L (x,i) R (x + e ,i)) (70) 0 2 m − a a − a i r m ! i X X (75) The Hamiltonian transforms to H′ = HKS′ + Hint′ + The Hamiltonian is invariant under transformations HM′ , where HM′ = HM (but with the new definition of generated by D (x) Q (x); the non-Abelian Gauss number operators with the hard-core bosons), a − a

2 g 2 1 H′ = J (x,i) ξp (x,i,j) Tr U (x,i) U (x + e , j) U † (x + e ,i) U † (x, j) + h.c. (76) KS 2 − 2g2 i j x,i x,i

N Gm (x)=2ΛµmDµ (x)+ s (x) (83) Qµ (x)=Λµmnm (x) δµ0s (x) (80) − r 2 analogously to (40). where Combining (83) with (63), we obtain the very simple Abelian transformation rules Λµm = (Tµ)mm (no summation.) (81)

(where µ =0, ..., N 1,m =1, ..., N). The normalization [Gk (x) ,Umn (y,i)] = δ (x, y) δkmUmn (y,i) − T (84) of the generators (60) implies that ΛΛ µν = δµν /2 and δ (x, y + ei) δknUmn (y,i) (no summation). 1 T − therefore Λ− = 2Λ . We now consider the N commuting Gauss laws of the Cartan subalgebra. Using (71) and This, along with (80), they can be rewritten as

Gm (x) (Gm (x) 1) ψ = 1 (85) N − | i Dµ (x) ψ = Λµmnm (x) δµ0s (x) ψ (82) | i − r 2 ! | i that follows directly from (79) as in the Abelian case, allows us to write Hint′ projected to the physical Hilbert 1 T Using Λ− = 2Λ with the above equation, we obtain space,

phys m n Hint′ = ξ (x,i) ǫ (x,i) P1 (x) ηm† (x) Umn (x,i) ηn (x + ei) P0 (x + ei)+ h.c (86) x,i,m,n X m where Pg (x) projects Gm (x) to g. The Gauss laws also enable us to rewrite, within the physical Hilbert space, the mass part of the Hamiltonian as x x Hphys = M ( 1)s( ) G (x)=2M ( 1)s( ) E (x,i) + const. (87) M − − x x,i X X where E is the Abelian, U(1) * SU(N) electric field (78). local, commuting transformations With all that in hand, we can finally construct the local building blocks of the transformation that removes x x the matter. On each vertex, we define the unitaries (x)= u (x)GN ( ) u (x)G1( ) (89) U N ··· 1

um (x)= ηm (x)+ ηm† (x) (88) Gm The order of the um in the product matters, since they They do not commute with each other on-site, but do do not commute, but different orders give rise, finally, to commute on different sites, which allows us to construct physically equivalent results, and as the local products the transformation from local, commuting pieces. We commute with one another it does not matter. Eventu- will have then a transformation = (x) with the ally, one obtains that U x U Q m n P (x) η† (x) U (x,i) η (x + e ) P (x + e ) † = U 1 m mn n i 0 i U m−1 n−1 (90) Gi(x) Gi(x+ei) =1 m n =1 ( 1) iP P (x) η (x) η† (x) U (x,i) η (x + e ) η† (x + e ) P (x + e ) ( 1) iP − 1 m m mn n i n i 0 i − One can see that the matter is completely decoupled here: the only instances of which is through projectors to nm = 0 everywhere. Therefore, we conclude that the original Hamiltonian is equivalent to this of a matter-less theory without local symmetries (all the Gauss laws are transformed to trivial 0 = 0 equations as in the Abelian case), whose Hamiltonian is phys s(x) nm =0 x H′ † nm =0 x = H′ +2M ( 1) E (x,i) ⊗x h | U U ⊗x | i KS − x,i X m−1 n−1 (91) Gi(x) Gi(x+ei) =1 m n =1 + ǫ (x,i) ξ (x,i) ( 1) iP P (x) U (x,i) P (x + e ) ( 1) iP + h.c  − 1 mn 0 i −  x,i,m,n X   11

The price we pay for having originally both cases. Then, on each link we introduce an extra Z2 fermionic matter is double: the factors Hilbert space - a two level space of a single spin - and m−1 n−1 Gi(x) Gi(x+ei) wherever Umn (x,i) appears in the Hamiltonian, we re- x i=1 i=1 ξ ( ,i) ( 1) P ( 1) P for the statis- place it by U ′ (x,i)= Umn (x,i) Z (x,i), where Z (x,i) − − m n mn tics, and the projectors P1 (x) and P0 (x + ei) for the is the Pauli z operator acting on the auxiliary field on finiteness of the local matter spaces. Both extend the that link. We do not include dynamics for this field, as interaction range to nearest neighbor links, but not explained in [7]. beyond. The auxiliary field, when coupled to the matter in the What happens if we wish to include a larger represen- above way, introduces, in the extended Hilbert space, an tation of the fermions? In this case, the method will not extra local Z2, given by work, for the simple reason that we have exactly N com- x muting Gauss laws in the Cartan subalgebra, that can nm( ) X (x) ψ [X (x,i) X (x ei,i)] ψ = ( 1)Pm ψ be inverted to define N different G operators. These | i≡ − | i − | i m i can only correspond to N fermionic number operators Y (92) per vertex - no more. Once again, we will also not be both before and after transforming to hard-core bosons. able to extend the method for non-staggered, or flavored Instead of the Cartan-Weyl basis we used in the U(N) fermions. case, we will now use another form for the traceless generators, replacing Λ by

C. The SU(N) case Λµm′ = δµm δNm; µ =1, ..., N 1,m =1, ..., m − − (93) This brings the N 1 Cartan Gauss laws (82) to the form Finally, we wish to discuss the case of another very rel- − evant gauge group - SU(N). Can we repeat the same pro- D (x) ψ = (n (x) n (x)) ψ (94) cedure there? First, of course, we need to map fermionic µ | i µ − N | i matter to hard-core bosonic ones. According to [7], this Note that is possible only for SU(2N) without extra ingredients. If one wishes to do it for SU(2N +1), an auxiliary Z2 gauge D (x) ψ Dµ (x) ψ = nm (x) NnN (x) ψ field has to be introduced, without dynamics, to absorb | i≡ | i − ! | i the parity of the fermions and enable the transformation. µ m X X (95) However, even if we do that, it will not allow us to and in particular repeat the procedure used for U(N), since in SU(N) we nm(x) only have N 1 generators in the Cartan subalgebra - D(x) NnN (x) − ( 1) ψ = ( 1) ( 1)Pm ψ (96) and the fundamental representation is N dimensional. − | i − − | i We only have N 1 Gauss laws of the form (82) per Using the auxiliary Z2 local symmetry we obtain vertex, while we need− N equations of the form (83). The x x inversion discussed in the U(N) case will not be possible ( 1)D( ) ψ = ( 1)NnN ( ) X (x) ψ (97) now (the matrix Λ is no longer square). As we shall show, − | i − | i for SU(2N) there is only one way to proceed, through and now the roads for even and odd values of N split. x the introduction of an auxiliary U(1) gauge field, but In the even case, ( 1)NnN ( ) = 1, and the above equa- for SU(2N + 1) one can also use another method, the tion shows us that− we gain nothing from introducing the auxiliary Z2 gauge field introduced in [7] is enough. auxiliary field, since a local Z2 exists in the SU(2N) case anyway, as it is a subgroup of the group’s center, Z2N . Therefore, the current discussion can only be valid for x x 1. The SU(2N + 1) case SU(2N + 1), where, since ( 1)NnN ( ) = ( 1)nN ( ) = − − 1 2nN (x), we obtain the desired equation We seek for another, independent equation, that will − 1 D(x) add up to the N 1 Gauss laws (82), allowing us to ex- nN (x) ψ = 1 X (x) ( 1) ψ (98) plicitly solve for− each n separately. One could try to | i 2 − − | i m think of using the center ZN symmetry of SU(N) - an and immediately define Abelian subgroup that is a subgroup of the missing U(1) 1 x component we had in U(N). However, the elements of G (x)= 1 X (x) ( 1)D( ) (99) the center are diagonal SU(N) elements that are gen- N 2 − − erated with the N 1 Cartan generators, so this will Combining it with (94), we can finally define introduce no further,− independent equation.

Suppose that we add an auxiliary Z2 ﬁeld, following Gµ (x)= Dµ (x)+ GN (x) (100) the procedure of [7]. While for SU(2N + 1) it is required for the conversion to hard-core bosons, it is not required as well, which completes a linearly independent set for SU(2N) - but we can try to add it nevertheless in of N mutually commuting operators Gm (x) satisfying 12

G (x) ψ = n (x) ψ as in the U(N) case - but con- that add the new U(1) component to the former U op- m | i m | i structed differently. One can then construct the desired erators. We modify the Hamiltonian HSU(N) to another transformation using these G (x) operators. one, H , by replacing any SU(N) operator U by the U m U(N) As before, one needs to constrain the operators Gµ extended U(N) operator U ′, without adding any dy- 2 to have only 0, 1 eigenvalues, which gives rise to local namics: the J terms in HKS are left only with the constraints in the final, transformed model. However, SU(N) generators. Then, the original SU(N) Hamilto- note that such a constraint is not required for GN now, nian may be obtained by projecting the new U(N) Hamil- since this operator only has in its spectrum 0, 1 anyway. tonian to a configuration with φ = 0 everywhere: since φ =0 U ′ φ =0 = U, and H completely commutes For the SU(2N) case, however, we will have to use h | | i U(N) another method, that will be discussed next. with all the φ operators, we get that φ =0 H φ =0 = H (102) h{ }| U(N) |{ }i SU(N)

Therefore, if ψ is an eigenstate of HSU(N), ψ 2. The SU(2N) case | i | i ⊗ φ =0 is an eigenstate of the extended HU(N) with |{the same}i energy. However, this state is not invariant In order to solve the SU(2N) case - where we simply do under the complete set of U(N) gauge transformations, not have enough commuting, linearly independent equa- since the Abelian phase is fixed on all the links. tions to invert - we will need to introduce an auxiliary We therefore define, for each state ψ in the SU(N) U(1) gauge field, and embed SU(N) in U(N). This ap- physical Hilbert space, the state Ψ as| follows:i plies to any N, and thus we will not restrict ourselves | i i√2N (D0(x) Q0(x))α(x) only to SU(2N) in the discussion. As U(N) satisfies 1/2 − x − Ψ = − αe ψ φ =0 the requirements for transforming fermions to hard-core | i N D P | i⊗|{ }i Z bosons, no Z2 field has to be introduced, and the U(1) (103) auxiliary field should be introduced before converting the where α dα (x), α (x) is a set of local phases, and D ≡ x { } fermions to bosons. is a normalizationQ constant. The state ψ satisfies the On each link of our SU(N) system we introduce an ad- N 2 1 Gauss laws (71) corresponding to| SUi (N). The ditional U(1) gauge field, with electric field E and phase state− Ψ constructed from it using (103) preserves this operator φ. We define the extended U(N) group element symmetry,| i and adds up the missing U(N) Gauss law as operators, well - it has the complete U(N) gauge invariance. Let us show that the mapping ψ Ψ is an isomor- iφ(x,i) | i → | i U ′ (x,i)= U (x,i) e (101) phism. We begin with the norm of Ψ . | i

i√2N Q(x)(α(x) β(x)) i√2N D0(x)(α(x) β(x)) 1 x − − x − Ψ Ψ = − α β ψ e P ψ φ =0 e P φ =0 (104) h | i N D D h | | i h{ }| |{ }i Z Changing variables of integration, one obtains

i√2N Q(x)α(x) i√2N D0(x)α(x) 1 x − x Ψ Ψ = − β α ψ e P ψ φ =0 e P φ =0 (105) h | i N D D h | | i h{ }| |{ }i Z Z

i√2N D0(x)α(x) − x Since e P generates a pure-gauge transformation on the U(1) part, we obtain that

i√2N D0(x)α(x) − x φ =0 e P φ =0 = φ =0 φ = α (x) α (x + e ) = δ (α (x) α (x + e )) (106) h{ }| |{ }i h | − i i − i x,i x,i Y Y

We have more delta functions than integrations: the number of delta functions is the number of links, L, and the number of integrations is the number of vertices, (both depend on the system size and topology). WeN choose a NV path of V links that goes through all the vertices and use it to perform all the integrations, and then we are left N L V with the same phase everywhere - α0, and an inﬁnite constant, δ (0)N −N . Finally we obtain for the norm

i√2N Q(x) α0 1 L V V x Ψ Ψ = − δ (0)N −N (2π)N ψ e P ψ (107) h | i N h | | i The state ψ does not have a local U(1) symmetry, but has a global one (conservation of total number of particles). The | i i√2N Q(x) α0 x transformation we are left with, for the computation is the norm, is global, and hence ψ e P ψ = ψ ψ . h | | i h | i 13

So if we set

L V V = δ (0)N −N (2π)N (108) N we get that the mapping ψ Ψ preserves the norm. What about inner products?| i →Using | i exactly the same arguments as above, we obtain that, in general,

Ψ Ψ = ψ ψ (109) h 1| 2i h 1| 2i Next, we show that it preserves the Hamiltonian matrix elements.

i√2N (D0(x) Q0(x))β(x) i√2N (D0(x) Q0(x))α(x) 1 x − − x − Ψ H Ψ = α β − ψ φ =0 e P H e P ψ φ =0 h 1| U(N) | 2i D D N h 1|⊗h{ }| U(N) | 2i⊗|{ }i Z (110) i√2N (D0(x) Q0(x))α(x) But since e− x − is a gauge transformation it commutes with H , and H ψ φ =0 = P U(N) U(N) | 2i⊗|{ }i HU(N), and HU(N) ψ2 φ =0 = HSU(N) ψ2 φ =0 . Therefore, and using the same methods for the integration etc., | i ⊗ |{ }i | i ⊗ |{ }i Ψ H Ψ = ψ H ψ (111) h 1| U(N) | 2i h 1| SU(N) | 2i

Indeed, the physical Hilbert spaces of HSU(N) and troduce other local constraints, accounting for the finite- HU(N) are isomorphic, and one may use HU(N) to study ness of the local Hilbert spaces of the original matter, HSU(N). This is true for open boundary conditions, and slightly extend the range of interactions due to the where there are no topological sectors. In the case of same fact, and also accounting for the fermionic statis- periodic boundaries, only states in the same topological tics, building up on the result of [7]. sector are connected unitarily to φ =0 , and hence the This opens the way to possibly easier variational and right sector must be sought. But|{ in general}i we see that numerical computations for Hamiltonian lattice gauge studying the extended theory where SU(N) is embedded theories, without having to deal with extra degrees of into U(N) is physically equivalent. Therefore, one can freedom for the matter - in particular with fermionic Fock use the above prescription to embed an SU(N) lattice spaces, as well as for quantum simulation of lattice gauge gauge theory into a U(N) one, and then eliminate the theories [16–20] with fermionic matter, without the ac- fermions as explained in the previous section. tual use of fermionic degrees of freedom in the simulator. After the completion of this work, we became aware of another work [21] discussing the elimination of the matter V. SUMMARY AND CONCLUSIONS degrees of freedom, in 1+1 compact QED (the Schwinger model). In this paper we have shown how to completely remove staggered fundamental fermionic matter from U(N) and SU(N) lattice gauge theories, by solving the Gauss law ACKNOWLEDGMENTS for the matter, which, unlike the solution for the electric field, is independent of the dimension. We extended EZ would like to thank David B. Kaplan, Mar- the well known procedure for complex scalar fields, aris- tin J. Savage and John Preskill for insightful discus- ing from the Higgs mechanism, to fermionic matter and sions. JIC is partially supported by the EU, ERC grant showed that it can be eliminated as well, by making use of QUENOCOBA 742102. This work was supported by the the gauge symmetry and breaking it. However, unlike in EU-QUANTERA project QTFLAG (BMBF grant No. the scalar case, when the matter is fermionic one must in- 13N14780).

[1] F. Englert and R. Brout, [6] L. Susskind, Phys. Rev. D 16, 3031 (1977). Phys. Rev. Lett. 13, 321 (1964). [7] E. Zohar and J. I. Cirac, [2] P. W. Higgs, Phys. Rev. Lett. 13, 508 (1964). Phys. Rev. B 98, 075119 (2018). [3] E. Fradkin and S. H. Shenker, [8] C. J. Hamer, Z. Weihong, and J. Oitmaa, Phys. Rev. D 19, 3682 (1979). Phys. Rev. D 56, 55 (1997). [4] K. G. Wilson, Phys. Rev. D 10, 2445 (1974). [9] B. Bringoltz, Phys. Rev. D 79, 105021 (2009). [5] J. Kogut and L. Susskind, Phys. Rev. D 11, 395 (1975). [10] M. C. Ba˜nuls, K. Cichy, J. I. Cirac, and K. Jansen, 14

Journal of High Energy Physics 2013, 158 (2013), 10.1007/JHEP11(2013)158.[16] L. Tagliacozzo, A. Celi, A. Zamora, and M. Lewenstein, [11] E. A. Martinez, C. A. Muschik, P. Schindler, D. Nigg, Annals of Physics 330, 160 (2013). A. Erhard, M. Heyl, P. Hauke, M. Dalmonte, T. Monz, [17] U.-J. Wiese, Annalen der Physik 525, 777 (2013). P. Zoller, and R. Blatt, Nature 534, 516 (2016). [18] E. Zohar, J. I. Cirac, and B. Reznik, Rep. Prog. Phys. [12] P. Sala, T. Shi, S. K¨uhn, M. C. Ba˜nuls, E. Demler, and 79, 1 (2015). J. I. Cirac, Phys. Rev. D 98, 034505 (2018). [19] M. Dalmonte and S. Montangero, Cont. Phys. , 1 (2016). [13] E. Zohar, arXiv:1807.01294 [quant-ph] (2018). [20] J. Preskill, arXiv:1811.10085 [hep-lat] (2018). [14] J. Haegeman, K. Van Acoleyen, N. Schuch, J. I. Cirac, [21] F. Surace, P. P. Mazza, G. Giudici, A. Lerose, A. Gam- and F. Verstraete, Phys. Rev. X 5, 011024 (2015). bassi, and M. Dalmonte, arXiv:1902.09551 [cond- [15] E. Zohar and M. Burrello, mat.quant-gas] (2019). Phys. Rev. D 91, 054506 (2015).