Higgs Phases and Boundary Criticality
in memoriam Peter Higgs
Abstract
Motivated by recent work connecting Higgs phases to symmetry protected topological (SPT) phases, we investigate the interplay of gauge redundancy and global symmetry in lattice gauge theories with Higgs fields in the presence of a boundary. The core conceptual point is that a global symmetry associated to a Higgs field, which is pure-gauge in a closed system, acts physically at the boundary under boundary conditions which allow electric flux to escape the system. We demonstrate in both Abelian and non-Abelian models that this symmetry is spontaneously broken in the Higgs regime, implying the presence of gapless edge modes. Starting with the U(1) Abelian Higgs model in 4D, we demonstrate a boundary phase transition in the 3D XY universality class separating the bulk Higgs and confining regimes. Varying the boundary coupling while preserving the symmetries shifts the location of the boundary phase transition. We then consider non-Abelian gauge theories with fundamental and group-valued Higgs matter, and identify the analogous non-Abelian global symmetry acting on the boundary generated by the total color charge. For SU() gauge theory with fundamental Higgs matter we argue for a boundary phase transition in the O() universality class, verified numerically for . For group-valued Higgs matter, the boundary theory is a principal chiral model exhibiting chiral symmetry breaking. We further demonstrate this mechanism in theories with higher-form Higgs fields. We show how the higher-form matter symmetry acts at the boundary and can spontaneously break, exhibiting a boundary confinement-deconfinement transition. We also study the electric-magnetic dual theory, demonstrating a dual magnetic defect condensation transition at the boundary. We discuss some implications and extensions of these findings and what they may imply for the relation between Higgs and SPT phases.
Contents
- I Introduction
- II Boundary Symmetry Breaking in the Abelian Higgs Model
-
III Boundary Symmetry Breaking in Non-Abelian Higgs Models
- III.1 Preliminaries: Lattice Formulation
- III.2 Boundary Symmetry Breaking in the SU(2) Higgs Model
- III.3 Boundary Symmetry Breaking in Non-Abelian Group-Valued Higgs Models
- III.4 Boundary Symmetry Breaking in Non-Abelian Fundamental-Higgs Models
- III.5 Boundary Phase Transition Order and Universality Class
- III.6 Summary and Discussion: Non-Abelian Case
- IV Boundary Symmetry Breaking in Higher-Form Abelian Higgs Models
- V Discussion and Conclusion
- A Discrete Differential Calculus for Abelian Fields
- B U(1) Bulk Phase Diagram: Limits, Symmetries and Monte Carlo
I Introduction
Gauge fields and gauge invariance have a long and complex history in theoretical physics, deeply interwoven with the advent of quantum field theory, the formulation of the Standard Model of particle physics, and firmly embedded in the modern theory of quantum many body systems. In fundamental physics, gauge theories arise naturally in Lorentz covariant theories of massless particles where they resolve a mismatch between the physical degrees of freedom admitted by the Wigner little group—such as photon polarizations—and the vector potential used to describe the particle states. This is accomplished by rendering the surplus degrees of freedom redundant. As such they arise naturally in field theories of gravity, nuclear forces and electromagnetism. In condensed matter physics, gauge fields play a key role in describing a plethora of physical phenomena including highly-entangled emergent states of matter such as spin liquids and fractional quantum Hall fluids.
The redundancy inherent to gauge theories leads to important subtleties. In particular, while the theories are local, the physical gauge-invariant objects are non-local Wegner-Wilson string loops [1, 2]. This implies a tension in describing the spontaneous breaking of symmetries in the presence of dynamical gauge fields. Whereas many condensed matter systems demonstrate symmetry lowering phase transitions governed by spontaneous breaking of a global symmetry, gauge redundancy, being unphysical, cannot be broken. This fact, enshrined in Elitzur’s theorem, belies a rich landscape of different phases separated by phase transitions whose study began systematically in the 1970s. While Landau theory successfully accounts for a broad range of phase transitions in correlated many-body systems, the order parameters for theories with dynamical gauge fields are generically non-local, raising the question of how to understand the nature of the phase transitions in such theories. Recent advances generalizing notions of symmetries to higher-dimensional charged objects [3] have allowed for the extension of the Landau paradigm to describe such phase transitions [4, 5], for a review see [6].
Condensation of charged matter in gauge theories, governed by the Anderson-Higgs mechanism [7, 8, 9, 10], is one of cornerstones in physics. It plays a key role in our understanding of superconductivity phenomenon and clarifies the nature of electroweak interactions within the Standard Model of particle physics. One question which has stood the test of time is the relationship between the Higgs and confined phases of a gauge theory. Let us illustrate this with the 4D compact U(1) gauge theory, described in terms of a U(1) vector potential , i.e. Maxwell theory with magnetic monopoles. In the pure gauge theory at weak coupling there is a deconfined phase in which static electric charges interact via a Coulomb potential mediated by a gapless photon, while at strong coupling the theory enters a confined phase, via the proliferation of magnetic monopoles, where static electric charges interact via a linearly rising potential. If we couple the gauge field to a charged scalar Higgs field, then the Higgs field may condense, leading to a Higgs regime where the photon becomes massive. Can the Higgs regime be sharply distinguished from the confined regime? For a Higgs field in the fundamental representation of the gauge group, i.e. one carrying elementary charge, it is generally understood that the Higgs and confined regimes are actually the same phase, i.e. they are not separated by any thermodynamic bulk phase transition, as shown in [11, 12]. This Higgs-confinement continuity is believed to be true in generic models with a gauge group coupled to a scalar Higgs field in the fundamental representation, including both discrete and continuous Abelian and non-Abelian gauge groups [12].
Nonetheless, the question of whether there is a qualitative difference between these two regimes has been revisited time after time from many different perspectives [13]. To survey briefly the history of the endeavor to delineate these two regimes, one approach has been to perform a partial gauge fixing and observe symmetry breaking in an unfixed global subgroup, which shows a phase transition separating them, though the location of the transition line is gauge-dependent and thus lacks a clear physical meaning [14]. Other proposals seek to delineate them in the presence of global symmetries (whose realization is unaltered between the two regimes), see e.g. [15, 16, 17]. Yet another approach, advanced partially by one of the authors, emphasizes that (in a certain limit) Abelian Higgs phases with fundamental matter exhibit symmetry-protected topological (SPT) order [18, 19, 20]. This observation motivates the investigation of the Higgs mechanism in open geometries, and zooms in on low-energy excitations localized near boundaries. In contrast to the confined regime, where the ground state is unique, in the Higgs regime previous studies uncovered energy spectrum degeneracies, see for example [18, 19, 20]. The robustness of these degeneracies originating from boundary-localized modes originates from the interplay of the protecting (generalized) symmetries that depends on the gauge group and dimensionality of the problem. In summary, the presence of a boundary introduces a criterion by which one can delineate the Higgs and confined regimes of Abelian gauge theories with fundamental matter—they are separated by a boundary phase transition.
In this paper we explore in detail boundary symmetry breaking in Wilson lattice gauge theories. We find that the Higgs-confinement boundary criticality mechanism is in fact ubiquitous. We begin in Section II by showing the presence of a boundary phase transition in the Abelian Higgs model in four dimensions, where the magnetic one-form symmetry is broken explicitly. We discuss how, in the presence of boundaries which allow flux but not charge to exit the system, there is a bulk U(1) global matter symmetry which, by the Gauss law, is equivalent to an electric flux symmetry acting on the boundary. We show that in a particular limit of the theory, in which the action reduces to a 3D XY model on the boundary, this boundary U(1) symmetry can be broken spontaneously. We provide numerical evidence that there is a corresponding boundary phase transition in the presence of bulk fluctuations, and we trace the phase boundary in the bulk phase diagram. Next we turn to non-Abelian Higgs theories in Section III. We consider two types of Higgs models, those with group-valued Higgs fields and those with fundamental representation (vector-valued) Higgs fields, which coincide for gauge group but differ for other gauge groups. We demonstrate that these models also have a global charge symmetry which is realized at the boundary. Using large-scale lattice simulation, we first show that 4D SU(2) Higgs theory has a boundary phase transition in the 3D O(4) universality class, verifying our theoretical prediction. We show that this boundary symmetry breaking is expected to be generic in group-valued Higgs models, and provide a general argument that fundamental-Higgs models with gauge group SU() and SO() exhibit O() and O() boundary criticalities, respectively. We verify this prediction numerically for the case of the 4D SU(3) fundamental-Higgs. Lastly, in Section IV, we consider generalizing to higher-form Abelian-Higgs models, with a -form gauge field coupled to a -form Higgs field. We discuss how the higher-form matter symmetry is realized at the boundary through the Gauss law, and show that the action reduces in a limiting case to a boundary -form gauge theory which may exhibit a confinement-deconfinement phase transition in which the matter symmetry is spontaneously broken. In the same section, we perform a duality transformation and discuss how this symmetry breaking can be viewed from the perspective of magnetic defects which live at the boundary. Finally, we provide an overview of our findings and discuss how our work may connect to symmetry protected topology.
II Boundary Symmetry Breaking in the Abelian Higgs Model
The prototypical theory for the interaction of charges with a gauge field is the U(1) Abelian-Higgs model—i.e. scalar QED, electrodynamics in (3+1)D Lorentzian or 4D Euclidean dimensions coupled to scalar matter. The continuum action for this theory is
(1) |
where is the field strength tensor, is the vector potential, is a complex scalar Higgs field, is the covariant derivative, and is a potential for the Higgs field. This theory is known to exhibit two phases: a deconfined or Coulomb phase, where charged particles interact via a Coulomb potential mediated by massless photons and a gapped confinement-Higgs phase. The latter phase has two distinct regimes, a confined regime at strong coupling and a Higgs regime at weak coupling, which are continuously connected without any thermodynamic phase transition between them [12].
This theory is invariant under local gauge transformations of the form
(2) |
where is an arbitrary 0-form. By construction, gauge-non-invariant operators cannot exhibit a non-zero vacuum expectation which is formalized by Elitzur’s theorem. Thus, while it is commonly stated that this theory exhibits spontaneous breaking of the U(1) gauge symmetry leading to the Higgs phase, and thus cannot serve as a local order parameter for such a phase transition. Rather the gauge-invariant observables are non-local string operators such as Wilson lines, for example
(3) |
where is a curve from to . When properly normalized, such observables allow one to distinguish quantiatively the deconfined phase from the Higgs-confined phase as discussed in some detail by Fredenhagen and Marcu [21] and [22].
We demonstrate that, in the presence of open boundaries of certain type, there is a second-order boundary phase transition distinguishing the Higgs and confined phases. We derive an explicit boundary theory in a limit where the bulk is completely frozen, which we show is a 3D XY model, and demonstrate with Monte Carlo that the critical exponents of the boundary transition do not change when we restore bulk fluctuations. Surprisingly, the line of boundary transitions appears to merge with the bulk critical endpoint, see Fig. 1. We then discuss deformations of the model which tune the location of the boundary transition.
II.1 Preliminaries: Lattice Formulation
We begin our discussion with a quick summary of the formulation of the discretized lattice theory, the importance of the Gauss law, the role of magnetic monopoles, before introducing the open boundary problem and presenting our numerical results.
II.1.1 Action Formulation
To study this theory in more depth we consider regulating it by imposing a UV lattice cutoff. We undertake our exploration of Higgs phases in gauge theories within the Wilson-Fradkin-Shenker lattice formulation. We work on a 4D hypercubic lattice with linear dimension and periodic boundaries, a discretization of four-dimensional Euclidean spacetime. We consider a a complex 0-form Higgs field taking values at each vertex . Expanding the Higgs field at site as , we freeze the radial mode by fixing the radius , which does not affect the qualitative physics.111 This freezing corresponds to the limit of infinite bare Higgs self-coupling [12, 23], and the radial mode will be restored upon coarse-graining. Equivalently it may be regarded as a Stückelberg field, and the model can be viewed as a lattice discretization of a gauged nonlinear sigma model with target space U(1). We use these two perspectives to give two different generalizations to non-Abelian gauge groups in Section III. Thus we work with the compact -valued 0-form , i.e. the phase of the Higgs field, which is minimally coupled to the dynamical U(1) gauge field. We consider a compact 1-form gauge potential taking values on each oriented link , . Denoting the reversed orientation by , the gauge potential satisfies . We will often denote link variables by their endpoints, i.e. . See Appendix A for a more detailed explanation of lattice differential forms and the notation used here.
We demand the theory to be invariant under gauge transformations of the form of Eq. 2, which on the lattice become
(4) |
where is an arbitrary -valued 0-form, and is the discrete exterior derivative, . The minimal gauge-invariant building blocks are the Wilson links , defined on oriented links ,
(5) |
and the minimal Wilson loops , defined on oriented plaquettes ,
(6) |
where . The minimal gauge-invariant Euclidean lattice theory is the governed by what we will refer to as the Fradkin-Shenker action,
(7) |
which reduces upon substituting in Eqs. 5 and 6 to the Abelian-Higgs model, the lattice equivalent of Eq. 1 in the limit where the radial mode of the Higgs field is frozen,
(8) |
Note that may be interpreted as the squared length of the Higgs field. The generating function of the model is
(9) |
II.1.2 Hamiltonian Formulation and Gauss Law
It will also serve us to consider the Hamiltonian formulation of the model on a 3D cubic lattice with continuous time. This may be obtained from the action by fixing to temporal gauge ( on all timelike links) and taking the continuum limit in the time direction, expressing the partition sum in terms of transfer matrices [24]. The Higgs field phase and gauge connection become operators, denoted and respectively, acting on a local Hilbert space on each vertex or link. They each have a canonically conjugate operator, denoted and respectively, both with integer eigenvalues, satisfying and . Thus is the raising/lowering operator for , while is the raising/lowering operator for . The Hamiltonian may then be expressed as222Note that the and couplings in the Hamiltonian formulation cannot be quantitatively compared to their values in the Lagrangian formulation, as they are renormalized when taking the continuum limit in the timelike direction [25].
(10) |
The operator counts the amount of charge on site , while counts the number of electric field lines on oriented link .
Gauge transformations, Eq. 4, are implemented by the operators
(11) |
Here we have used the coexterior derivative (see Appendix A)
(12) |
where indicates the coboundary, the set of oriented links ending at site . Demanding that acts as the identity on physical states for arbitrary implies that physical states satisfy the Gauss law constraint
(13) |
at each site , where we used to rewrite the constraint in terms of the lattice divergence. Gauge invariant states satisfying the Gauss law are then created by Wilson line operators,
(14) |
where is a 1-dimensional contour in the lattice. Acting on the trivial vacuum state with everywhere, this operator creates a unit electric charge/anti-charge pair at the ends of the contour connected by a string of unit electric flux. If is a closed contour, this inserts a closed string of electric flux.
II.1.3 Magnetic Monopoles
In addition to the electric sector, there is also the magnetic sector, though it is not readily seen in this formulation, instead being exposed by duality transformations [26, 27, 28, 29, 30] (see Section IV.2. In the Hamiltonian formulation, the magnetic excitations are sources of divergence of the magnetic field, , i.e. magnetic monopoles. In the action formulation, they may be viewed as U(1) vortex defects of the gauge field, characterized by , which are allowed because the identity is only enforced modulo . In 4D, these homotopy defects form 1-dimensional closed strings in the dual lattice, which we refer to as ’t Hooft loops, the worldlines of magnetic monopoles. They are necessarily included in the Euclidean lattice gauge theory partition sum due to the compactness of the U(1) link variables. In the limit of Eq. 8, fluctuations of the gauge field are completely suppressed and no monopoles are present. Correspondingly, we may associate this absence of monopoles to a 1-form symmetry, called magnetic symmetry, indicated in Fig. 1.
II.2 Open Boundaries and Global Symmetry
We address now a well-known, but subtle point: by taking in (4) to be a constant function, it appears at first sight that this theory has a global 0-form U(1) symmetry, shifting , leaving unchanged since . In the Hamiltonian picture, this transformation is generated by the total charge,
(15) |
thus such a symmetry corresponds to global conservation of electric charge. Note, however, that in the absence of boundaries this “global symmetry” is pure gauge, because all physical quantum states belong to the zero-charge sector and thus carry the same quantum number. By the Gauss law, Eq. 13, is exactly zero for a system with periodic boundary conditions. In other words, since all electric flux lines must end somewhere inside the system, the system must be globally charge neutral. By construction, such a “symmetry” is therefore trivial and cannot be explicitly or spontaneously broken. However, in presence of specific boundary conditions, these global U(1) transformations actually generate a physical global symmetry that acts on the boundaries which can be spontaneously broken [19, 20].
We consider a lattice with open boundaries in the form illustrated in Fig. 2. The bulk of the system is a (hyper)cubic lattice with sites indicated by white spheres, links by black lines, and plaquettes by gray faces. At the boundary, we include a layer of cubic cells (green) which separate the bulk from the vacuum (blue sites, edges, and plaquettes). In particular, there is a set of links bridging between the bulk and the vacuum which carry dynamical gauge degrees of freedom and can therefore support electric flux lines which effectively exit the system. The vacuum side (blue) does not contain any dynamical degrees of freedom.
Key to our choice of boundary conditions is that we demand that the Gauss law, Eq. 13, is respected at every bulk site (white sphere), including those at the ends of the boundary links. No Gauss law constraints are imposed at vacuum sites. Let denote the gauge potential on the boundary link touching site , oriented “in” from the vacuum to the bulk. For the Gauss law to be respected at , we must have that under the gauge transformation Eq. 4, , i.e. it is “uncompensated” at the vacuum end of the link.
With this choice of boundary conditions, the global part of the gauge symmetry becomes physical—charge can pass in and out of the system, meaning there are different gauge-invariant charge sectors. More precisely, there are gauge-invariant half-open Wilson string operators with one end in the bulk and the other passing through the boundary,
(16) |
where is a contour starting in the vacuum and ending at bulk site . These operators create an isolated charge in the bulk, attached to an electric flux line that exits through the boundary, which is not possible in a closed system. These half-open string operators are charged under the global transformations generated by , Eq. 15. By the Gauss law, is equivalent to the net electric flux through the boundary,
(17) |
where the sum is over all boundary links (green in Fig. 2) oriented out. If the Hamiltonian contains no half-open string operators, then the bulk charge is conserved, and there is a global U(1) symmetry. Eq. 17 defines a “bulk-boundary correspondence” between charge and flux, and generates a global symmetry which may either be seen as acting on the bulk matter degrees of freedom or on the boundary gauge degrees of freedom.
II.3 Boundary Symmetry Breaking in the Abelian Higgs Model
The theory can now be chosen such that the Hamiltonian commutes with the charge . This now-physical global symmetry corresponds to global charge conservation or, equivalently, conservation of flux through the boundary. Since the open system has a global symmetry on the boundary, it may be spontaneously broken, a scenario we now study. With the boundary conditions shown in Fig. 2, the Euclidean action that we will study is given by
(18) |
with the boundary portion given by the Wilson plaquette loops on the boundary plaquettes (light green in Fig. 2),
(19) |
The absent links on the vacuum side are excluded, so that
(20) |
where indicates the value of on the boundary link touching site , oriented inwards from the vacuum to the bulk.333We may equivalently consider the exterior vacuum (blue in Fig. 2) to have trivial gauge field on all vacuum links and zero Higgs field on all vacuum sites.
II.3.1 Boundary XY Model at Infinite
To begin, we consider the behavior of this theory in the limit, i.e. deep in the Higgs regime. From the bulk action, Eq. 8, in this limit the bulk satisfies the constraint , i.e. the bulk gauge field is exact and thus pure gauge. Indeed there are no physical degrees of freedom left in the bulk—rotating to unitary gauge, , we end up with on all bulk links and a vanishing matter field. However, no such constraint is enforced on the boundary links bridging between the bulk and vacuum, and the gauge field on these links is free to fluctuate. Thus we obtain a dynamical 3D theory on the boundary of the system governed by the boundary action in Eq. 19.
In this limit, referring to Fig. 2, each boundary plaquette (green) has one edge in the bulk (black) with , and two edges straddling between the bulk and vacuum (green) which remain dynamical degrees of freedom. Substituting this into Eq. 20, we can recombine terms into the gauge-invariant variables
(21) |
corresponding to a half-open Wilson line coming from the vacuum and ending at site . The boundary action can then be written in the gauge-invariant form
(22) |
which is a 3D XY model at inverse temperature . This must exhibit a continuous phase transition from a paramagnet at small to a spontaneously broken phase at large . Thus we infer that along the line in the phase diagram there is a boundary phase transition in the 3D XY universality class indicated at the top of Fig. 1. We note that in the case of the Abelian-Higgs model, the same mechanism generates a boundary Ising model at [19].
II.3.2 Boundary Phase Transition at Finite
Next we consider to be large, , but finite. The constraint is no longer enforced exactly, so we expand in small fluctuations as . Assuming that the bulk action can be expanded in terms of (i.e. that topological defects are negligible), the bulk action becomes a Proca-type action,
(23) |
where , which describes a massive 1-form field. The boundary action is then
(24) |
While at infinite the theory reduces to an XY model on the boundary, at finite the XY model is minimally coupled to the weakly fluctuating massive bulk gauge field. While the bulk field lives in a higher dimension, the boundary remains quasi-3D, exponentially localized with a length scale determined by the mass of the bulk photon, . We therefore expect the symmetry breaking phase transition at the boundary to persist at large but finite .
We may ask where the boundary transition line may run from , finite . As it is a spontaneously broken symmetry it must end either on a boundary of the phase diagram or on a bulk transition line. The former case is ruled out as follows: It cannot end on the line because this is trivial from the point of view of both bulk and boundary variables. It also cannot end on the line because matter decouples on this line. The only remaining possibility is that the line ends at but there we understand the bulk theory as being pure gauge and an XY model so the physical degrees of freedom on the boundary drop out. We conclude that the boundary transition line must end on a bulk transition line.
These arguments are suggestive of the picture illustrated in Fig. 1, with a boundary phase transition between the Higgs and confinement regimes of the bulk phase diagram. To test this assertion, we have carried out Monte Carlo simulations of the full 4D lattice gauge theory with boundary. We compute the local XY order parameter on the boundary as well as gauge invariant bulk observables and . The results are summarized in Fig. 3. We find clear signs of a boundary phase transition, with Fig. 3(a) showing the boundary order parameter behavior as a function of while holding fixed showing behavior consistent with a continuous boundary phase transition. Fig. 3(b) shows the Binder cumulant for the order parameter taken along a cut at for different system sizes ranging from to , showing crossing behavior consistent with a second-order transition. Lastly, Fig. 3(c) shows the Binder parameter with scaled by using the 3D XY critical exponent [31], showing excellent scaling collapse, confirming a second-order phase transition on the boundary even for only moderately large .
Monte Carlo simulations of the 3D XY model in the literature put the critical point at [31]. We find that tends towards this value in the large limit. On general grounds we should expect that bulk fluctuations will serve to disorder the boundary. Therefore, by lowering we expect the transition shifts to larger (lower effective temperature in the statistical model). This is indeed what we observe numerically. For even smaller we find that the boundary transition line appears to intercept the bulk critical endpoint, as illustrated in Fig. 1. Below, by changing the boundary theory, we demonstrate that this result is not universal.
II.3.3 Tuning the Boundary Coupling
We now consider tuning the boundary coupling in Eq. 19 relative to the bulk, parameterized by the dimensionless ratio
(25) |
Such a change is allowed by gauge symmetry and does not affect either the electric or magnetic 1-form symmetries. The reason for this modification is that by tuning we can shift the location of the critical in the limit, as . This implies that the location of the boundary transition line must shift in the phase diagram in order to meet the location of the transition in the limit.
Indeed, this is precisely what we find numerically with the resulting transition lines shown in Fig. 4. We find, for all , that the transition is present and that it drifts to larger as is reduced, consistent with having to lower the temperature more to suppress the enhanced matter field fluctuations. For the line shifts to larger values of , and appears to separate from the tricritical point, instead terminating on the first order transition line separating the Higgs phase from the deconfined phase. For , the boundary transition line shifts to smaller values of , and appears to continue to terminate on the tricritical point.
II.4 Summary and Discussion: Abelian Case
In this section we have explored the case of the U(1) Abelian-Higgs model in spacetime dimensions with open boundary conditions of the form shown in Fig. 2, governed by the action defined by Eqs. 8, 19 and 18. In the case of open boundaries there is a global bulk charge symmetry for this model, because charge is not allowed to enter or leave the system. By the Gauss law, this is physically equivalent to a symmetry acting on the boundary of the system, yielding conservation of total electric flux through the boundary. We investigated the possibility of spontaneously breaking this symmetry.
In the limit, the bulk degrees of freedom are fully frozen while the boundary degrees of freedom remain fluctuating. The boundary theory in this limit can be written in the gauge-invariant form of a -dimensional 0-form U(1) model, i.e. an XY model, Eq. 22, which therefore exhibits a boundary phase transition at a critical value of in the appropriate XY universality class. By performing explicit numerical simulation in Euclidean spacetime dimensions and measuring the gauge-invariant order parameter, Fig. 3, we found that this boundary phase transition persists away from the limit, and appears to stay in the 3D XY universality class. The transition line appears to terminate at a critical point in the bulk phase diagram and can be tuned by adjusting the boundary coupling relative to the bulk coupling, as shown in Fig. 4.
This boundary transition line appears to delineate between the Higgs and confining regimes of the phase diagram. The idea that the Higgs phase can be characterized by boundary symmetry breaking was first raised in [19, 20]. This gives physical meaning to the non-gauge-invariant adage that the Higgs regime is a charge condensate by imposing open boundaries that make the symmetry physical. We note, however, that by tuning the boundary coupling relative to the bulk coupling, the location of the transition line in the phase diagram can be moved. This raises questions as to what the boundary phase transition implies as a probe of bulk physics, which we defer to the discussion at the end of the paper (Section V.2). First we broaden our perspective by studying more complex non-Abelian models exhibiting similar physics in Section III.
A few comments on extensions and exceptional cases are in order before we proceed. These results naturally carry over to the discrete Abelian gauge groups by restricting the Higgs field and gauge field to take discrete values in multiples of . The boundary theory Eq. 22 then becomes an -state clock model. The case was discussed in [19] where the boundary theory was shown to be an Ising model. They also can be generalized to higher-form extensions of the Abelian-Higgs model, which we discuss further in Section IV.
The U(1) Abelian-Higgs model is interesting because it reduces in the limit to the XY model on the boundary, which exhibits a BKT transition. It would be interesting to know if this BKT transition persists to finite . We complete our discussion of special cases by highlighting here gauge theory in four dimensions coupled to a charge Higgs field. For this case, a sharp distinction can be made in the bulk between confining and Higgs phases with a transition between them. The distinguishing feature of this theory is the partial Higgsing of the gauge group down to [32]. For sufficiently large the bulk transition is the confinement-deconfinement transition of the residual gauge theory. On the boundary, however, one still expects the emergence of the XY model we identified for the case. It would be quite interesting to investigate the interplay of boundary U(1) 0-form symmetry breaking with the bulk 1-form symmetry breaking and the resulting topological order present at large and .
III Boundary Symmetry Breaking in Non-Abelian Higgs Models
We now turn to extend the results of the previous section regarding Abelian Higgs models and boundary criticality to non-Abelian Higgs models. We will show that the general picture of boundary symmetry breaking persists, albeit with a richer structure owing to a set of non-commuting gauge transformations. We discuss two types of non-Abelian Higgs models: those with group-valued Higgs fields and those with fundamental representation vector-valued Higgs fields with fixed length. The two classes of models are equivalent for gauge group , and are distinct for other gauge groups. The group-valued case is a relatively straightforward extension of the Abelian case, because the fixed-length XY-rotor Higgs field considered previously is naturally a U(1) group element. The vector-valued case is more subtle because the limit does not trivialize all bulk degrees of freedom. We present numerical results for both and vector-valued cases and extract corresponding boundary criticalities.
III.1 Preliminaries: Lattice Formulation
Let be a compact connected Lie group, e.g. or . The lattice action is formulated in terms of group-valued link variables satisfying , which may be viewed as the exponentiated Lie-algebra-valued gauge field, , where indicates path-ordering. Under a gauge transformation, these transform as
(26) |
where are arbitrary group elements associated to each site . The minimal gauge-invariant quantity is the Wilson plaquette-loop (compare to Eq. 6),
(27) |
where the superscript on the product indicates path-ordering, and the trace may be taken in a representation . Normalizing by the dimension of the representation ensures that the trivial Wilson loop has unit magnitude. We focus primarily on the cases and , taking the trace in the fundamental representation as matrices.
III.1.1 Higgs Fields
For the Higgs field, many different models can be considered by putting the Higgs field in different representations of the gauge group. We consider two different types here for concreteness: vector-valued (fundamental representation) Higgs, and group-valued Higgs. These are both possible extensions of the Abelian U(1) rotor model considered in Section II, since a rotor may be viewed either as a fixed-length vector, or as a U(1) group element. In either case, the action is given by the Fradkin-Shenker form, Eq. 7, the only difference being the definition of the Wilson link . The generating function for the quantum theory is given by the Euclidean path integral, where the integration over the group-valued variables is performed with respect to the Haar measure.
The familiar model is the fundamental-Higgs, where the Higgs field is an -component vector, as in the Standard Model and analogous to Eq. 1. Denoting the Higgs vector at site by , we freeze the radial mode as in the Abelian Fradkin-Shenker model. Gauge transformations rotate the Higgs field as where is a group element in the fundamental matrix representation. The group-valued link variables define parallel-transport maps for the Higgs field, , i.e. they related the color frames at neighboring sites, and the generalization of the gauge-invariant Wilson link observable, Eq. 5, is
(28) |
where is the canonical inner product on , indicates complex conjugation, and we enforce the fixed-length constraint .
The second type of model we consider takes the Higgs field to be group-valued, like the link variables, denoted . In this case, the Higgs field transforms as under gauge transformations, and we can define a gauge-invariant Wilson link by
(29) |
where we take the trace in the fundamental representation. This is a lattice regularization of a gauged principal chiral model [33], a non-linear sigma model whose target space is the group manifold.
III.1.2 Hamiltonian Formulation and Gauss Law
The Hamiltonian formulation of the non-Abelian lattice gauge theory has a similar form to the Abelian case, but the electric field of the non-Abelian theory carries color indices and the different components do not commute. Fixing to temporal gauge and reformulating the partition function using transfer matrices, taking the continuum limit in the time direction one obtains a Hamiltonian for the time evolution [34, 25]. The basic ingredients are the group-valued link operators with eigenstates , such that and , along with a set of translation operators
(30) |
where left and right translations correspond to the two orientations of the link. Each group element may be expressed as , where are real numbers and a basis for the Lie algebra of gauge group . The serve as coordinates on the group manifold, and may be thought of as generalized Euler angles. The link operators can then be expressed as
(31) |
The operators are position operators on the group manifold, while are the color-electric fields, which serve as the conjugate momenta and can be expressed as derivatives with respect to the . The electric fields satisfy the same commutation relations as the group generators,
(32) |
where are the structure constants of .
For the Higgs field in the group-valued representation, we define the group-valued operators (analogous to ) and left- and right-translation generators and (analogous to ). For the Higgs field in the fundamental vector representation with frozen radial mode, the classical configuration space is that of a rigid rotor, and we define corresponding angular momentum operators . The Hamiltonian is then given by the Kogut-Susskind form [34, 35]
(33) |
where all sites, links, and plaquettes are purely spatial. Here, are the matter charge operators,
(34) |
and and are the corresponding quadratic Casimir operators, which do not depend on whether we use left- or right-generators. The operator is the operator analog of Eq. 27, the trace of the oriented product of on the links of spatial plaquette . Similarly, is the operator analog of Eq. 29.
The eigenstates of correspond to the irreducible representations of the gauge group, with the acting as raising and lowering operators [34, 35]. The same is true for the group-valued Higgs, with acting as the raising and lowering operators, while for the fundamental vector-valued Higgs, the charge eigenstates are angular momentum eigenstates of a rigid rotor.
Gauge transformations are performed by the operators
(35) |
where the lattice divergence is defined as
(36) |
with the sum taken over the links emanating from site oriented out. acts as the identity on physical, gauge-invariant states, which therefore satisfy the color-electric Gauss laws,
(37) |
one for each color index.
Note that this is similar but subtly distinct from the Abelian case, Eqs. 11, 12 and 13, where we used . This relationship is not true in non-Abelian gauge theory. Instead, left- and right-translations of the gauge field are related by
(38) |
which implies that the electric field in the two directions along a link are related by
(39) |
where is the adjoint representation of . As such the gauge field itself is charged in the adjoint representation with respect to color rotations, generated by the charge operators
(40) |
which are manifestly orientation-independent. The classic (though heuristic) way to think of this is that the gauge bosons (gluons) carry a distinct charge and anti-charge in the two directions along the link. In the Abelian case, in Eq. 39, and the link charge Eq. 40 is exactly zero.
III.1.3 Open Boundaries and Global Symmetry
We introduce electric open boundary conditions as in Fig. 2, with dynamical links extending from the bulk to the vacuum which allow electric flux to pass through the boundary. We denote the link variables on the boundary link touching site by , with the convention that the link is oriented “in” from the vacuum to site . We have minimal open Wilson strings going around the boundary plaquettes, which we can write as
(41) |
Using these, we define the boundary action for our theory directly analogous to the Abelian case as
(42) |
Electric flux can thus enter and leave the system, but matter charges cannot. Under gauge transformations the fields transform as
(43) |
Notice that each boundary link only receives a transformation from the inside end, where it terminates on a matter field. In addition to this gauge symmetry, the boundary action Eq. 42 has a physical global symmetry acting on the boundary,
(44) |
where every boundary link is translated from the outside end.
This boundary symmetry is similar to the Abelian case, but with a subtle distinction. The total color charge of the system, including boundary links, is
(45) |
which rotates all matter fields in the fundamental representation and all gauge fields, including boundary links, in the adjoint representation. But this is not the generator of global gauge transformations, Eq. 43, under which the boundary links only rotate from the inside. The generator of global gauge transformations is the “bulk charge”
(46) |
where the second sum contains only the bulk links, and in the last sum the boundary links are oriented inwards. By the Gauss law, this operator must be zero on the physical gauge-invariant Hilbert space. On the other hand, the generator of the global symmetry is the “boundary charge”
(47) |
again with inward orientation. Together these make up the total charge of the system,
(48) |
where we have assumed global gauge invariance to identify . Thus the boundary symmetry Eq. 44 may be viewed, by the Gauss law constraint, as being generated by the total color charge of the system. This is analogous to the Abelian case, where the total charge of the system is just the matter charge, Eq. 15, since the links do not carry any charge. Note, however, that in the non-Abelian case the charge of the boundary links is “fractionalized” into a piece that contributes to the bulk charge and a piece that contributes to the boundary charge.
III.2 Boundary Symmetry Breaking in the SU(2) Higgs Model
We focus first on the case where the gauge group is . In this case the fundamental and group-valued representations are equivalent. In the fundamental representation, each is a vector with unit length, and the configuration space is a 3-sphere. Notice that is also topologically the 3-sphere, meaning that every configuration of the Higgs vector can be written as a unique fundamental-representation matrix times a fixed vector, for example as
(49) |
where and are complex numbers. Denoting the matrix Eq. 49 as , note that the determinant of this matrix is the length of the rotor. The Wilson link for the vector Higgs model, Eq. 28, can then be written as
(50) |
which is exactly equivalent to the group-valued Higgs definition, Eq. 29. This makes the SU(2) Higgs rotors special, since they may be viewed as either vector-valued or group-valued (which is also the case for U(1) rotor).
III.2.1 Bulk Phase Diagram
To map the phase diagram, we perform classical Monte Carlo simulations for the 4D model defined by the Fradkin-Shenker action, Eq. 7, with periodic boundary conditions on an hypercubic lattice. The phase diagram can be mapped out by measuring the Wilson plaquette and link observables along with their variances (i.e. susceptibilities), which are shown in Fig. 5. There is a roughly horizontal first-order transition line extending from a critical endpoint at [36] (cyan box) towards . This line is clearly visible in the link susceptibility, Fig. 5(b).
The phase diagram exhibits only one clearly distinct bulk phase—the confined-Higgs phase—a deconfined phase being absent in non-Abelian theories. The phase diagram is, however, roughly separated into two regions indicated by the behavior of the Wilson link expectation value, shown in Fig. 5(a), with indicating the confining regime and indicating the Higgs regime. To the left of this critical endpoint is a supercritical region (a Widom line [38]) extending to smaller and larger , a rapid crossover from the Higgs to the confined regimes. The location of this supercritical region is most evident in the Wilson plaquette susceptibility, Fig. 5(d), which shows a pronounced intensity emanating from the critical endpoint. We expect this phase diagram to be qualitatively consistent with those for general non-Abelian Higgs models with either fundamental vector- or group-valued Higgs fields. In the group-valued case the continuity of the Higgs and confined regimes was proven by Fradkin and Shenker [12].
III.2.2 Boundary Symmetry Breaking
We now consider open boundaries, with boundary action Eq. 42 in a slab geometry. As in Section II.3, we start by considering the limiting behavior when . In this limit every bulk link satisfies the constraint . In this section we will resolve this constraint by fixing a gauge. Gauge-invariant formulations for the group-valued representation are presented in Section III.3, and for the fundamental representation in Section III.4.
From Eq. 50, if we fix to unitary gauge where all the Higgs rotors are aligned globally, or , the bulk constraint becomes , or , which can only be satisfied if on every bulk link. Thus in the limit of the SU(2) Higgs model, the bulk is completely frozen and has no remaining degrees of freedom, as in the Abelian case. The boundary links, however, have no constraint and remain fluctuating. In unitary gauge where the bulk links are set to the identity, the boundary action becomes
(51) |
This boundary action may be viewed as a lattice discretization of a nonlinear -model with target space the SU(2) group manifold, i.e. a principal chiral model [33].
The boundary model, Eq. 51, has an SU(2)SU(2)O(4) symmetry. To see this, we can re-express it as an O(4) Heisenberg model by representing the SU(2) group-valued link variables as unit quaternions,
(52) |
where the are real numbers, and , , and are Pauli matrices. The boundary action then becomes an O(4) Heisenberg model,
(53) |
Therefore, in the limit , the system exhibits a boundary phase transition in the 3D O(4) universality class, at a critical coupling [37], with order parameter,
(54) |
The gauge-invariant object which reduces to in the unitary gauge is
(55) |
where the former is an matrix which decomposes according to Eq. 52, and the latter is explicitly a 4-component unit-length vector. For large but finite , the bulk fluctuations are strongly gapped and the boundary should behave as quasi--dimensional, and we expect the boundary phase transition to persist as in the Abelian case.
To test this prediction, we perform Monte Carlo simulations with open boundaries and measure the order parameter, Eq. 54 (defined in terms of gauge invariant observables Eq. 55), at finite values of . In Fig. 6 we show the evolution of the boundary order parameter as a function of , for in (a) and in (b), for different system sizes. These reveal a transition from a disordered, symmetric boundary on the confined side (small ) to an ordered, symmetry-broken boundary on the Higgs side (large ). This value of is quite close to the bulk critical point (), demonstrating that the boundary phase transition persists far into the phase diagram where the bulk is quite strongly fluctuating. In Fig. 6(c) and (d), we show the Binder cumulant for the same cuts, showing crossing behavior at a -dependent critical coupling. In (e) and (f) we have rescaled using the 3D O(4) critical exponent [37], demonstrating a clean scaling collapse, thus verifying that the transition remains second order and in the same universality class even for relatively small values of . The transition line appears to terminate at the bulk critical point, which can be seen in the red line in Fig. 7, though verifying this numerically is difficult as the bulk correlation length grows larger than the finite width of the open boundaries as the system approaches bulk criticality.
III.2.3 Tuning the Boundary Coupling
We now consider varying the parameter , which shifts the location of the transition. Figure 7 shows the behavior of the boundary order parameter and its susceptibility for different values of , along different constant- cuts at fixed system size. The corresponding behavior of the bulk link susceptibility is shown in black in the background for reference. The bulk transition line moves as is varied, but appears to remain second-order throughout. Figure 8 summarizes the results by showing the approximate location of the boundary transition line for different values of , with results very similar to the Abelian case (Fig. 4). For the transition moves to smaller values of and appears to terminate at the bulk critical endpoint (cyan box). No boundary transition is detected for small values of below the bulk first-order line. For the location of the boundary transition line moves to larger values of , and appears to terminates on the line of bulk first-order transitions, at least for sufficiently small .
III.3 Boundary Symmetry Breaking in Non-Abelian Group-Valued Higgs Models
Having verified the existence of a boundary phase transition in the SU(2) Higgs model, which is both a fixed-length-rotor Higgs model and a group-valued Higgs model, we now consider how these results generalize to these two types of models separately for a general gauge group . From the point of view of the boundary action, the group-valued Higgs is the simpler case, so we consider it first. The action (up to an overall normalization convention for the trace) is
(56) |
The boundary action is given by Eq. 42. In addition to the global color charge symmetry acting on the boundary, Eq. 44, it also has a global symmetry given by right multiplication of the Higgs field
(57) |
The global symmetry group is therefore . We expect the bulk phase diagram to be qualitatively similar to the SU(2) case, Fig. 5, with a single thermodynamic phase, a first-order line terminating at a critical endpoint. These models were considered by Fradkin and Shenker [12], who showed that the Higgs and confining regimes are contiguous, as in the Abelian models.
We now consider taking the limit. Maximizing the trace in Eq. 56 yields the constraint , which implies that the bulk links can be expressed in terms of the matter field as
(58) |
The bulk of the system is completely frozen in this limit, which can be most easily seen in unitary gauge where . Substituting Eq. 58 into Eq. 41, the boundary action can then be expressed in terms of the gauge-invariant observables
(59) |
where was defined with the boundary link oriented “in”, which are short half-open Wilson strings coming from the vacuum and ending at site . The boundary action becomes
(60) |
which is a lattice chiral model. Compare this to the equivalent results for the Abelian case, Eq. 22 and Eq. 51, to which it reduces when . The symmetry acts on this chiral model by left and right multiplication of the , respectively. Such a model is known to exhibit symmetry breaking to the diagonal subgroup [39]. We discuss boundary phase transitions in these models further in Section III.5
Note that if is Abelian there is no distinction between left and right group multiplication, therefore and are not independent symmetries of the system. corresponds to global rotation of the Higgs field generated by the total matter charge, Eq. 15, while is generated by the net electric flux through the boundary links, Eq. 17. These two generators are the same operator on the physical gauge-invariant Hilbert space by the Gauss law. In contrast, in the group-valued non-Abelian Higgs models these are really distinct symmetries generated by different physical operators.
III.4 Boundary Symmetry Breaking in Non-Abelian Fundamental-Higgs Models
We now consider non-Abelian gauge groups SU() (and by trivial generalization SO()) with fundamental Higgs fields. As a convenient shorthand, in this section we will represent the Higgs vectors as kets, , which should not be confused with quantum states. Furthermore, we suppress the subscript “f” on the fundamental representation matrices . The action we write as
(61) |
The boundary action is given by Eq. 42. For the Higgs rotors can no longer be identified as group elements. Because of this, the constraint,
(62) |
does not completely trivialize the bulk. To see this, note that the constraint enforces that nearest-neighbor Higgs rotors are parallel relative to the gauge field. However, since there is still freedom to perform rotations about the colinear axis of the remaining components of the Higgs field, the gauge field can continue to fluctuate so long as it does not rotate the Higgs field away from the parallel axis. This is explicitly seen by rotation to unitary gauge where the Higgs field is parallel in a global frame, which then fixes one diagonal element of the gauge field to unity, i.e.
(63) |
This constraint forces the link matrices to take the form
(64) |
Thus in this limit the bulk theory becomes (gauge-equivalent to) an gauge theory. Note that is trivial, making the case special, as discussed in Section III.2.
III.4.1 Gauge-Invariant Resolution of the Infinite Kappa Limit
To resolve the constraint in a fully gauge-invariant fashion analogous to Eqs. 22 and 60, we first note that the constraint Eq. 62 actually implies that
(65) |
which follows from the fact that real part of a Hermitian inner product on is the Euclidean product when the vector space is viewed as . In other words, and have the same real and imaginary components, and so are the same complex vector. We can think of and as unit vectors spanning a two-dimensional complex vector space, in which case must act within this subspace as the unique SU(2) rotation rotating to .444 Within the two-dimensional subspace this rotation is given by , using the notation of Eq. 49. That this rotation is unique follows from the fact that the two vectors have unit norm within this subspace, thus they live on a 3-sphere, and SU(2) is isomorphic to the 3-sphere, so each point on the 3-sphere corresponds to a unique SU(2) rotation. Therefore we can resolve the constraint as
(66) |
where is an matrix which preserves , i.e. an rotation in the subspace orthogonal to . Note that the ’s are independent for every link, i.e. they are associated to the ends of the links and are independent on different links touching the same site. Furthermore, there is only one independent degree of freedom on each link, because is equivalent to . Thus the constraint reduces each link variable to an degree of freedom, and the whole theory reduces to an gauge theory.
While Eq. 66 demonstrates the reduction of the gauge group, it is not that useful for formulating the boundary theory. A more useful way is to decompose each by sandwiching it between two resolutions of the identity decomposed into the parallel and perpendicular subspace of and . Namely, for site
(67) |
where is the projector to the orthogonal complement of . Note that under a gauge transformation . Inserting this identity on either side of a link variable decomposes it into four pieces,
(68) |
In the limit , this simplifies significantly. Firstly, the two cross terms (the last line) are exactly zero by Eq. 66, i.e. because . Second, in the second term, Eq. 65 reduces it to . In summary, in the limit every link variable in the bulk can be expressed as
(69) |
It follows that a product of two consecutive link variables is
(70) |
where we used that and . Therefore gauge-invariant closed Wilson loops have a projector to the orthogonal subspace inserted between each consecutive link,
(71) |
which is another manifestation of the Higgsing down to an gauge theory.
Now consider the three-legged plaquettes appearing in the boundary action, Eq. 42. Inserting the identities at the two bulk sites, we obtain
where we have defined the gauge-invariant variables
(72) |
which are vectors corresponding to the short half-open Wilson strings at the boundary. This highlights an important distinction for fundamental Higgs compared to group-valued Higgs models—here the half-open Wilson string is a vector degree of freedom, not group-valued. The theory then can be expressed in the gauge-invariant form
(73) |
Thus the bulk Higgses down to an gauge theory while the boundary decomposes into a pure part plus a part which may be viewed as an ferromagnet. The global boundary symmetry acts as and simultaneously . We generically expect this symmetry to spontaneously break above a critical down to , and the short Wilson string rotors on the boundary to exhibit long range order and Goldstone modes. We discuss the nature of this phase transition further in Section III.5.
III.4.2 Formulation for General Gauge Groups
For a general gauge group with fundamental Higgs, the large- limit Higgses it down to a subgroup . We consider here cases where the residual gauge group is non-Abelian. In the limit the bulk fluctuates as a pure gauge theory with gauge group governed by the Wilson action, while the boundary links continue to explore the full gauge group . The generic picture (which can be obtained in unitary gauge) is that the action Higgses down to
(74) |
where is the link variable restricted to the subgroup, and is the corresponding Wilson plaquette for these -valued bulk links. Taking the traces in any faithful representation of the group should yield the same physics.
Assuming the bulk is gapped with a finite correlation length, as it must be if is non-Abelian, the boundary is quasi--dimensional with some finite correlation length extending into the bulk. The system retains a global symmetry rotating all of the boundary links, together with the bulk gauge symmetry. The appropriate boundary theory is therefore expected to be a gauged nonlinear -model with target space with subgroup gauged, or equivalently, a nonlinear -model with target space the quotient space . For example, Higgses down to and the quotient space is , which agrees with the finding in Eq. 73 of a boundary theory of rotors, whose configuration space is a sphere in dimensions. Similarly, Higgses to , with quotient space .
III.4.3 Hamiltonian Perspective
Similar considerations apply to the Hamiltonian formulation of the non-Abelian Higgs theory described in Section III.1.2. The analog of the large limit in the Hamiltonian formulation (Eq. 33) is first of all to drop the conjugate variables leaving only the variables in the matter sector. This essentially renders the Higgs fields classical and they may be gauge fixed without loss of generality along some fixed direction . There is now a Hamiltonian constraint
(75) |
that breaks gauge fluctuations from down to .
Now if we consider only the boundary plaquettes, this constraint acts only on the bulk links parallel to the boundary while those extending out of the boundary have no such constraint. Therefore the boundary theory of the four dimensional bulk in the large limit is a three dimensional chiral model that is partially gauged by an gauge group where the gauge theory permeates the bulk.
III.5 Boundary Phase Transition Order and Universality Class
So far we have demonstrated that Higgs models with group-valued or fundamental Higgs fields have large- limits with well-defined boundary degrees of freedom that may exhibit symmetry breaking. Having discussed the boundary actions for the family of fundamental Higgs SU() models we are in a position to say something about the phase transitions at the boundary. As the bulk is gapped for finite we should expect that the boundary theory has a finite correlation length into the bulk making it a quasi-()-dimensional boundary theory so that statements at hold also for finite .
In the light of Refs. [19, 20] and our analysis above, the quasi-()-dimensional boundary theory typically should have a symmetry breaking phase transition as the boundary coupling is tuned. But what is the nature of the phase transition? When there are only global symmetries, the symmetry group and the spacetime dimension determine the type of the transition. When some symmetries are gauged, the full set of global symmetries (including higher-form symmetries originating from the gauging [3]) are expected to determine the nature of the phase transition, while the gauge redundancy only serves to reduce to the quotient space. For the coupled bulk-boundary models considered here, the bulk gap ensures the integrity of the boundary model across the phase diagram well away from bulk critical points. Having identified the physical gauge-invariant boundary variables on the lattice that can go critical, in this part we write down the corresponding Landau boundary theories.
III.5.1 Group-Valued Higgs
The boundary degrees of freedom, Eq. 59, are gauge-invariant composites transforming under the global symmetry. Considering the case , the coarse-grained (generally complex) matrix-valued fields are denoted [40], where are color indices, which transform as under the symmetry. The Landau theory is
(76) |
The symmetry is susceptible to breaking down to the diagonal subgroup. This set of models has been studied in Ref. [40] to which we refer for more details. One finds, in the case , that the determinant contributes to the quadratic term. Numerically one finds a continuous transition consistent with a quartic term stabilizing the free energy. In the case , the determinant is cubic implying that the transition is first order. In the case , the determinant contributes a quartic term but with a negative sign that is expected to drive the transition first order. We therefore expect a continuous transition for and first order for and for .
Chiral models on a lattice, such as the one in Eq. 60, have been studied for many years especially in two dimensions at large , where they are integrable [41]. If one is interested in boundaries of four dimensional gauge theories, the three dimensional analogs of such models are of interest. One early work on in three dimensions, Ref. [40], contains Monte Carlo results for . For (Section III.2) the symmetry group is and the transition is continuous, consistent with O() criticality. For the numerical results reveal the transition to be first order in agreement with the mean field theory predictions.
III.5.2 Fundamental Higgs and SU(3) Numerical Results
The boundary model in the case of a fundamental Higgs has gauge-invariant degrees of freedom of the form where are the color indices. There is a single global symmetry that acts from the left . This will be broken spontaneously for sufficiently large (modulo Mermin-Wagner restrictions). The coarse-grained version of is denoted and the Landau theory for SU() is
(77) |
which is invariant under U() transformations. This may instead be viewed as a theory of real variables invariant under the enlarged O() symmetry group. This Landau theory gives the impression that, for the three dimensional boundary of a four dimensional SU() gauge theory, there is a phase transition in the O() universality class. In the case of SU() this is O() criticality as shown above at the level of the microscopic model both analytically and numerically. Note that the coupling of the rotors in Eq. 73 is invariant under O() rotations, even though the microscopic action manifestly only has a global symmetry.
It may seem surprising that the SU() invariant model exhibits O() criticality. One simple check is that SU() has generators and that SU() SU() symmetry breaking therefore has broken generators (corresponding to the number of Goldstone modes) which matches the count for O() O() symmetry breaking where O() has generators. But what about the terms that are SU() invariant but not O() invariant? These terms are certainly present. An example is given by taking an operator which transforms as and considering its determinant. This is only U() invariant, but it is also an irrelevant operator for in . More precisely, the couplings in the action originating from the determinant have mass dimension which, in dimensions higher than two, is negative for all but in three dimensions and in four dimensions. As we have seen, the case of SU() is special because it is identical to a problem with a group valued Higgs. So we have treated it separately. Therefore the remaining puzzle relates to SU() in three dimensions where the determinant coupling goes like and is marginal by power counting. As with the problem of scalar field theory in three dimensions [43], among other cases, we expect this sixth order term to be marginally irrelevant. Taking this together with our results for SU() we surmise that the SU() boundary theory phase transition is in the O() universality class since terms breaking O() down to SU() are irrelevant or marginally irrelevant.
The discussion above shows that we should set aside the three dimensional boundary theory of the SU() Higgs model in four dimensions for further examination. We anticipate a boundary phase transition in the O() universality class for this case. This model is also interesting because, unlike the SU() Higgs model, it has a dynamical bulk as the Higgs mechanism in this case leaves residual bulk degrees of freedom. We have therefore performed simulations of SU() gauge theory with a fundamental Higgs field. The Higgs field is a -component complex vector at each site, and we define the average of the Higgs field measured from the vacuum end of the boundary links,
(78) |
The numerical results are presented in Fig. 9. Figure 9(a) shows the behavior of this order parameter for various system sizes, while Fig. 9(b) shows the behavior of its variance (susceptibility). The behavior is qualitatively consistent with an order parameter at a second order phase transition, approaching zero at small and continuously increasing from , with a diverging susceptibility. Figure 9 shows the Binder cumulant scaled by the 3D O() critical exponent [42], which shows a clear collapse to a universal scaling function. We therefore conclude that the SU() Higgs model in indeed demonstrates a boundary phase transition in the D O() universality class, even though a single fundamental Higgs field does not remove all bulk degrees of freedom when , in agreement to the general picture of boundary O() criticality on the boundary of fundamental Higgs models.
III.6 Summary and Discussion: Non-Abelian Case
In this section we have discussed global boundary symmetries in non-Abelian Higgs models and their spontaneous symmetry breaking with both fundamental Higgs fields and group-valued Higgs fields. For general non-Abelian gauge group there is a global symmetry acting on the boundary of the system, which is equivalent by the Gauss law to the total color charge of the system, Eq. 48. For the group-valued case, there is an additional bulk matter symmetry given by right multiplication of the matter field. If is an Abelian group, these two symmetries exactly coincide, since there is no distinction between left and right multiplication. For , the two types of models are equivalent due to fact that is equivalent to the unit quaternions. For other gauge groups the two types of models are inequivalent. In all cases, the location of the boundary transition that we identify can be shifted by tuning the boundary coupling.
In the group-valued case, the infinite limit freezes the bulk of the system. The resulting boundary theory can be expressed in terms of gauge-invariant half-open Wilson strings at the boundary, Eq. 60. This may be viewed as a lattice discretization of a principal chiral model, with symmetry coming from the boundary color flux and bulk matter symmetries. This global symmetry is expected to be broken to the diagonal subgroup at large . Since the boundary is frozen in the limit and should remains strongly gapped at large , the boundary transition is expected to persist also for finite . Our SU(2) Monte Carlo simulations verify this, and the boundary transition line appears to terminate at the bulk critical endpoint.
For Higgs fields in the fundamental representation, there is only a single global symmetry, generated by the total color charge of the system and acting on the boundary. The general picture is that the bulk gauge symmetry is Higgsed down to a subgroup . If is trivial, then the bulk is completely frozen in the limit. If is Abelian, there may be a bulk phase transition even at infinite (though that does not preclude a boundary symmetry breaking). If is non-Abelian then the bulk reduces to a gapped gauge theory as . The boundary theory is expected to be a gauged nonlinear -model with target space , Eq. 74. Note that by placing the Higgs in other representations (e.g. adjoint) or adding addition Higgs fields, one may target different subgroups [44] and thus obtain different boundary theories. For example, by starting from an gauge theory with Higgs fields in the fundamental representation, the general expectation is that it reduces to an gauge theory in the bulk, and the boundary target space is the Stiefel manifold , which have recently attracted significant interest [45]
We considered in particular the case (and by direct extension SO()). In the limit, the bulk Higgses down to a gauge theory, while the boundary decomposes into a pure part and an rotor part. The global symmetry, under most circumstances, will be broken spontaneously at large . We argued that this SU() fundamental Higgs boundary phase transition lies in the O() universality class, since there are no relevant operators differentiating from . We tested this prediction numerically for the lattice gauge theory with fundamental Higgs, demonstrating a clean scaling collapse with O( critical exponents.
One interesting exceptional case is SO() gauge theory with fundamental Higgs. This is equivalent to the Georgi-Glashow electroweak theory, with gauge group and the Higgs field in the adjoint representation. In the limit we can fix to unitary gauge which Higgses the bulk gauge group down to . In and , the bulk has a confinement-deconfinement transition as a function of , while we predict an additional O() boundary critical point in the Higgs regime. This may have implications for the physics of domain walls or cosmic strings formed in the early universe.
IV Boundary Symmetry Breaking in Higher-Form Abelian Higgs Models
Thus far we have discussed boundary symmetry breaking in Higgs phases of both Abelian and non-Abelian gauge theories with 1-form gauge fields, drawing a general picture of a global charge symmetry realized on the boundary due to the Gauss law constraint of a gauge-invariant system. Here we consider a further extension, to higher-form gauge fields. A -form gauge field describes the parallel transport of charged -dimensional objects [46, 47, 48, 49]. Two-form gauge fields are often called Kalb-Ramond fields in string theory literature [50], they appear in the dual descriptions of superfluids and superconductors [51, 52] and may be realized in certain spin models on frustrated lattices [53]. The gauge group for is generically Abelian, because a unique path ordering only exists on 1-dimensional contours [46, 47].
IV.1 Higher-Form Abelian Higgs Models
Here we consider the case of -form U(1) gauge theory, though restriction to subgroups follows. Let be a -form field and a -form field, each taking values in in Euclidean spacetime dimensions, with .555 If there is no gauge field. If then the gauge field is not dynamical and can be completely integrated out using the Gauss law, yielding long-range interactions for the Higgs field. If then identically. By a -form we mean a function on oriented -dimensional cells of the lattice, such that , where denotes the cell with the opposite orientation. See Appendix A for a more detailed discussion of the discrete differential forms notation used throughout this section. Let denote the collection of -cells of the lattice, each with a fixed orientation. The fields are governed by the generalized Fradkin-Shenker action
(79) |
The exterior derivative of a -form is a -form whose value is defined by the discrete Stoke’s theorem,
(80) |
where the sum is over the -cells forming the oriented boundary of the -cell . This is a straightforward generalization of the Abelian Higgs model, Eq. 8, to a theory of extended -dimensional charged objects attached to -dimensional electric flux branes [46].
The action is invariant under higher-form gauge transformations,
(81) |
for an arbitrary -form . We will see that this enforces the Gauss law attaching electric branes to the charged objects. When , this gauge invariance includes “gauge-of-gauge” transformations
(82) |
for arbitrary -form . This enforces an additional Gauss law which states that the -dimensional electrically charged objects are closed.
When , Eq. 79 reduces to a pure -form U(1) gauge theory, which has a global electric -form symmetry given by
(83) |
corresponding to conservation of global electric flux. In spacetime dimensions, it also admits magnetic homotopy defects, whose cores trace out -dimensional worldsheets in spacetime [46, 47, 54], where
(84) |
In the limit it has a -form magnetic symmetry, discussed further in Section IV.2.
The phase diagram of this model is expected to be similar to that of the 1-form U(1) gauge theory, sketched in Fig. 1, as long as (). In the marginal case (), the magnetic defect is an instanton (zero dimensional in spacetime) and is expected to destabilize the deconfined phase [55, 56, 57]—a generalization of the Polyakov mechanism for 1-form U(1) gauge theory in [58], which itself may be viewed as a higher-form generalization of the Mermin-Wagner theorem for 0-form symmetries [3, 4]. In those marginal cases, the bulk phase diagram will generically be similar to that of the non-Abelian models, as in Fig. 5 [55, 59, 60].
IV.1.1 Gauss Law and Higher-Form Matter Symmetry
Before we introduce the Hamiltonian formulation and discuss matter symmetry, we point out that the Higgs charges for behave slightly differently than for . When , the zero-dimensional point charges must come in pairs at the two ends of oriented electric strings, and are either positive or negative depending on which end of the string they sit. When the charges are extended oriented objects (e.g. strings when ) living on the edges of electric -branes. Such a brane for can have a single edge, meaning that it is perfectly valid to have a single charged object in the system. As such the extended charged objects are net-charge-neutral if they are contractible [46, 47]. It is only if they wrap around periodic boundaries that they have non-trivial global charge.
The Hamiltonian formulation follows the exact arguments laid out in Section II.1, with conjugate operators on all -cells and on all -cells. Invariance of physical states under the gauge transformations, Eq. 81, enforces the Gauss law
(85) |
where is a -cell and the sum is over its coboundary, the set of -cells containing in their positively-oriented boundary. This simply says that the number of electric -branes emanating from is equal to the amount of electric charge on . For , the charges carry a sense of orientation. For example, when the charges are nothing but the electric strings of a 1-form gauge field. The “gauge-of-gauge” symmetry, Eq. 82, enforces the constraint
(86) |
which says that the charged objects are closed.666 This constraint obviously follows from the Gauss law Eq. 85 (following from ), just as Eq. 82 is already implied by Eq. 81, but it is worth spelling out for those unfamiliar with this point.
In the theory discussed in Section II with 0-form Higgs field, the net charge of the Higgs field generates a global 0-form symmetry, which is pure gauge with periodic boundaries but becomes physical with the choice of electric boundary conditions. The natural extension of this global matter symmetry for general is a -form symmetry,
(87) |
Letting be the dimension of space, the generators of these transformations are Gukov-Witten operators [3] supported on -dimensional closed surfaces in the dual lattice
(88) |
where is the -cell in the direct lattice corresponding to in the dual lattice, and is a -form Poincaré dual to .777 For our purposes, the defining property of the Poincaré dual of a -dimensional closed surface in the dual lattice is that it is a -form in the direct lattice acting as a generalized delta function, for example the unit -form supported on the direct lattice -cells piercing with the appropriate orientation. Note that is only defined up to an exact form because is closed, i.e. it is a cohomology class. These operators generate the symmetry, Eq. 87,
(89) |
where with a constant. Because is a closed surface its Poincaré dual is a closed form, .
These operators are topological, i.e. they only depend on the homology class of , owing to the matter Gauss law Eq. 86. Consider replacing it with another surface such that for some -volume . Then . Plugged into Eq. 88, we have
(90) |
where we used the matter Gauss law, Eq. 86. Therefore there is one -form charge generator for each homology class in . Deforming without changing its homology class corresponds to shifting by an exact form, which are just the gauge transformations of Eq. 82.
The operators charged under are the charge creation/annihilation operators, i.e. Wilson branes supported on open -dimensional surfaces which insert an electric membrane with charge on its boundary,
(91) |
where we have defined the intersection number between the -dimensional in the direct lattice and the -dimensional in the dual lattice,
(92) |
Thus the operator simply counts the intersection number of the closed charged objects with .
In the Maxwell case, , the matter charges are point particles, is a closed -dimensional volume in the dual lattice, and the only non-trivial choice is to take it to be all of space. The associated charge then simply counts the number of positive minus the number of negative charges. For a less trivial example, consider the case and , so that the matter charges are 1-dimensional strings and is a closed two-dimensional surface in the dual lattice, intersecting a collection of links in the direct lattice. A choice of its Poincaré dual is a 1-form which is zero everywhere except on the intersected links, on which it has unit value in the direction normal to the surface. The charge operator, Eq. 88, then counts (with signs, because ) the number of strings piercing the surface , i.e. the intersection number between closed strings and the surface. Because the charges are closed strings, this intersection number must be zero unless winds around a periodic boundary and intersects a charge which also winds around the periodic boundaries in the transverse direction.
In a closed system, gauge invariance guarantees that the matter charge operators are exactly zero because of the Gauss law, Eq. 85, which when inserted into Eq. 88 yields . Equivalently, the intersection numbers, Eq. 92, are exactly zero because has trivial homology class, i.e. since is closed. Consider for example the case , where the charge operators count the number of charged strings wrapped around a periodic boundary. We cannot, in a closed system, have a single non-contractible charged string, because it is attached to an electric membrane which must end somewhere inside the system. The only possibility is that it ends on another non-contractible charged string going the other direction, such that the two strings constitute the boundary of the electric membrane. This is the sense in which a closed system is charge-neutral when .
IV.1.2 Open Boundaries and Global Matter Symmetry
To make the charge operators non-trivial, we must impose boundary conditions such that either can terminate on the boundary, and is not a closed form, or such that the contour can terminate at the surface and thus not be closed. The latter case, however, means that the charged objects can pass through the boundary, meaning charge is not conserved. Therefore, extending our results from Section II, we take the former case. We consider open boundaries in the same form as Fig. 2. To be precise, the boundary consists of a layer of -dimensional hypercubes, with all cells in the boundary layer which do not touch a bulk cell removed (or, equivalently, and on those cells). In other words, we remove the cells on the vacuum side, creating links missing one end, plaquettes missing one edge, cubic cells missing one face, etc. The boundary action is given by
(93) |
where denotes the set of -cells in the boundary layer. As in Section II, this boundary action does not allow the matter field to tunnel through the boundary meaning that all the charged objects are closed and contained inside the bulk.
In the presence of these boundary conditions, the matter charge operators, Eq. 88, can generate a physical symmetry. For we must take care to consider how the surfaces , which live in the dual lattice, can terminate at the boundary. Because we chose a “rough” boundary, as in Fig. 2, with a layer of cells sticking out from the edge of the system, the dual lattice boundary is flat. This is illustrated in Fig. 10(a) which shows the direct lattice bulk and boundary as in Fig. 2, along with the dual lattice in red and purple. The bulk part of the dual lattice is colored purple, while the boundary layer of the dual lattice is colored red and can be seen to form a flat surface without any protruding cells. The charge operators can terminate on this flat dual boundary layer.
For concreteness, consider the case in , in which case the matter charges are strings and is a two-dimensional membrane in the dual lattice. Due to our boundary conditions the charge strings cannot terminate at the boundary, but the electric membranes can exit the system through the boundary. Referring to Fig. 10(b), consider a state with a single charged string wrapping around a periodic direction, as shown in Fig. 10(b) by the red line, attached to an electric membrane which exits through the boundary, illustrated by the yellow surface. This charged string is detected by taking the membrane to intersect it transversely, as illustrated by the cyan surface. The charge associated to the surface, Eq. 88, is then related to the electric flux through the boundary via the Gauss law,
(94) |
In the figure, is shown as a dark green line which pierces a collection of boundary plaquettes (green squares with black borders). According to Eq. 94, the amount of charge measured by is equal to the number of electric branes exiting through the boundary measured by the set of plaquettes pierced by . In summary, the physical matter symmetry is generated by charge operators supported on which terminate at the boundary and, by the Gauss law, acts on the gauge field on the boundary elements pierced by .
IV.1.3 Boundary Higher-Form Symmetry Breaking
Returning now to the higher-form gauge-Higgs action, Eq. 79, let us now see how the physical charge symmetry is spontaneously broken at the boundary. In the limit we have the constraint on every bulk -cell. This completely freezes the bulk degrees of freedom, as is seen by rotating to unitary gauge, , resulting in . A covariant way to see this is that the gauge field operators for -dimensional closed surfaces trivialize—if is contained entirely within the bulk, then since is closed. This constraint is not imposed on the field variables on the cells touching the vacuum, however. As a result, if exist through the boundary, then we can decompose it into two pieces, , where is the part of supported on bulk cells and is the part supported on boundary cells. Since , we have . As a result, the gauge field operators decompose as
(95) | ||||
(96) |
In other words, for closed surfaces exiting through the boundary, the constraint reduces to a half-open Wilson operator terminating on Higgs operators as soon as it touches the bulk. Similarly, longer half-open Wilson surfaces ending on the Higgs field in the bulk, e.g. Eq. 16, are reduced to short half-open Wilson surfaces ending on the Higgs field at the boundary.
This is depicted in Fig. 11 for the cases and , where is taken to be the smallest Wilson surfaces which are closed in the bulk, for . In the case , depicted in Fig. 11(a), is an oriented boundary plaquette (green), and consists of three links (yellow), one of which is in the bulk. We decompose into the two links in the boundary, denoted , and the one link in the bulk, denoted . The gauge field on the link in the bulk trivializes by the constraint into evaluated at its two endpoints, i.e. , depicted as red spheres in Fig. 11(b). These are the degrees of freedom appearing in the boundary action in Eq. 22. In the case , Fig. 11(c), is a three-dimensional cube and is a set of five plaquettes (yellow). In Fig. 11(d), it decomposes into four plaquettes in the boundary (yellow), , and one plaquette in the bulk, . The constraint turns the gauge field on this bulk plaquette into a Higgs string operator on its boundary (red). The extension to larger is obvious, but can’t be illustrated.
Using this, we can see that the boundary action Eq. 93 contains precisely these minimal operators. In the general case, the action reduces in this limit to
(97) |
We have already seen in Section II how in the case this can be recast as a 0-form XY model, Eq. 22, which we can now extend to the case . Note that every boundary -cell is associated to a unique bulk -cell (the one which trivializes under the constraint). We define composite degrees of freedom for each such pair,
(98) |
where is evaluated on the unique boundary cell corresponding to , which we treat as a -form rather than a -form. Note that are not gauge invariant when , because they would create open charged objects, but we can combine them to construct the gauge-invariant degrees of freedom appearing in Eq. 97. For example, in the case , each consists of on a boundary plaquette and on the bulk link it touches, four of which combine to form the gauge-invariant composite object in Fig. 11(d). Let us denote the layer of bulk cells which touch the boundary layer, consisting of cells up to dimension , as , illustrated in Fig. 10 by the dark layer between and . Each -cell is associated to a unique -cell . We treat as a -form gauge field in , so that its exterior derivative is given by
(99) |
where we identified with . We can then write the boundary action Eq. 97 as
(100) |
Equation 100 is precisely the action of a -form U(1) gauge theory defined at the boundary of the system in terms of composite half-open Wilson operators. In the case it reduces to the 0-form XY model identified in Section II.
We conclude that a -form Abelian-Higgs model in dimensions Higgses in the limit down to a -form gauge theory on the dimensional boundary. If (), this implies there will be a boundary phase transition. For this will be in the -dimensional XY universality class, with a symmetry-broken phase at large . For it will be a confinement-deconfinement transition (expected to be first-order [61]), with a confined phase at small and a deconfined phase at large . This reduces to the results reported in Section II for and . The symmetry that is spontaneously broken at large is the -form matter symmetry, corresponding to the electric symmetry of the boundary gauge theory when , under which (as a half-open Wilson operator) is charged. A summary of this result is given in the table Table 1 for and , which may be extended straightforwardly to all and . In the cases where (), a generalized Mermin-Wagner theorem [4, 6] (equivalently, the Polyakov mechanism or magnetic instanton proliferation [58]) will prevent the symmetry from breaking.
IV.2 Electric-Magnetic Dual Picture
We can gain further insight into these Abelian models and the boundary symmetry breaking by reformulating them in terms of dual magnetic variables. We will do so for gauge group U(1). The duality transformation is well-established in a variety of forms [27, 28, 29, 26]. We derive it here in the presence of our open boundary conditions, which provides a concrete picture of both the electric and magnetic sectors of the theory near the boundary.
IV.2.1 Open Boundary Duality Transformation
The partition function of the original action for the -form gauge theory is
(101) |
This can be turned into a theory of electric strings and particles by utilizing the identity
where are Bessel functions. We introduce integer -form coupled to and -form coupled to , to rewrite the partition function exactly as
(102) |
The restriction arises from the fact that we did not include half-open Wilson strings coupling the bulk to the vacuum. Utilizing the adjointness relation , this reduces to
(103) |
We interpret the -form as the worldlines swept out by the -dimensional Higgs excitations (carrying integer electric charge), and as the worldsheets swept out by the -form electric field. The first delta function is just the Gauss law—it tells us that the worldsheets of electric flux must terminate on the worldlines of the Higgs charges. The second delta function enforces that the the worldlines are “divergence-free,” i.e. they form closed -dimensional surfaces, corresponding to global conservation of charge. Note that, because the electric worldsheets can end, the electric 1-form symmetry , present when , is explicitly broken, though it can be restored as an emergent symmetry at energies below the charge gap.
With the choice of “electric” boundary conditions (Fig. 2), we have the constraint on all links extending from the bulk to the vacuum in Eq. 103, because there is no term for these links, while can be non-zero on the boundary plaquettes. This means that electric charge cannot enter or exit the system, but electric flux can. The fluxes in the bulk form closed surfaces unless they terminate on charges. In the presence of the boundary these fluxes may be cut off at the boundary without terminating on charges.
We resolve the two constraints by writing888Note that and can have a harmonic component corresponding to the electric winding sectors which we neglect. These give rise to ground state degeneracies in the topological Coulomb phase.
(104) |
for integer ()-form and ()-forms , respectively. Note that these are not gauge-invariant and can be shifted by co-exact forms. We utilize the large-argument expansion of the Bessel functions,
to approximate the partition function by999 We use the shorthand to denote .
(105) |
where the delta function comes from the constraint and we dropped the prefactor. Finally, we apply Poisson resummation to turn these into real-valued fields,
We promote the integer fields to real fields,
(106) |
coupled respectively to integer currents and via Poisson resummation. Lastly, we move to the dual lattice, replacing for each -form , where is a ()-form in the dual lattice, the discrete Hodge dual of . We thus obtain the dual partition function
(107) |
where and are forms, while and are forms, with the constraint that is zero within the dual boundary layer coming from the delta function in Eq. 105. Finally, let us rescale the fields by absorbing the into the definition of and , to obtain the dual action
(108) |
with dual couplings and .
The currents and correspond to the winding defects of the U(1) gauge field and the Higgs field , respectively. The former are the magnetic monopole worldlines, while the latter are the worldsheets of the vortices of the Higgs field. Under a gauge transformation of the -valued gauge fields , , the action is shifted by
(109) |
If we integrate over all gauges, i.e. over all the generators , we obtain delta functions that yield the constraint
(110) |
This says that the magnetic monopole worldlines form the boundaries of the Higgs vortex worldsheets. This is precisely the magnetic Gauss law, i.e. it says that magnetic monopoles are the sources of magnetic strings (compare to the electric Gauss law Eq. 103). This is familiar from superconductor phenomenology—vortex cores carry magnetic flux. It also follows from this that , i.e. that the magnetic worldlines are closed and magnetic charge is conserved, which follows from the “gauge-of-gauge” invariance for arbitrary .
This theory can be recast as a U(1) gauge theory as follows. Summing over undoes one of the Poisson resummations and forces , where the integer field is the Hodge dual of in Eq. 105. The residual gauge symmetry is and , where is an integer shift, meaning that the gauge-invariant configuration space for is actually and for is . The resulting theory is therefore a Villainized -form U(1) gauge theory [62, 20] (Eq. 84) coupled to magnetic currents, the dual of the original Abelian Higgs model model, Eq. 101, which was coupled to electric currents. The integer gauge field measures the winding numbers of the compact , and its fluxes are the homotopy defects of , which are the electric charges of the original theory. They act as sources for the fluxes , which are the electric strings in the direct lattice.
Let us briefly review how the dual bulk behaves in the various limits in the Maxwell case, and . First, consider the () limit: in the electric formulation the gauge field is turned off, , and the remaining Higgs sector is a gauged XY model. In the dual theory the magnetic charges are turned off (integrating sets ) and the magnetic 1-form symmetry is restored, resulting in a gas of closed membranes interacting via their coupling to the 2-form gauge field. Performing the Gaussian integration over we obtain the action , i.e. the gauge field generates Coulomb interactions among the membranes.
Next, consider the () limit. Electric charges are turned off, restoring the electric 1-form symmetry and reducing to a pure U(1) gauge theory. In the dual theory the gauge field is turned off, . The theory reduces to a Coulomb gas of magnetic monopole worldlines [26, 58], which has a transition separating the deconfined phase (low temperature condensate) and the confined phase (high temperature gas).
Lastly, consider the () limit, the strong coupling limit of the original theory. In the electric theory, we can fix unitary gauge to remove and obtain the action on every link independently, so the system is fully disordered. Equivalently, if we use Eq. 103, setting forces all of the Bessel functions to vanish except when , which reduces the partition function to . In the dual theory, , the weak coupling limit, and we have the constraint , which further implies . The resulting theory then just has Lagrange multipliers that force the charge loops and membranes to vanish, so the theory trivializes completely.
IV.2.2 Dual Boundary Symmetry Breaking
We now consider setting in Eq. 108 (). The variables in play are and , defined on all links of the dual lattice, and and on all plaquettes. In the bulk, the means that the kinetic term for drops out and so the fluxes of (corresponding to electric charges of the original action) are completely unconstrained. In other words, the bulk is in the strong coupling limit. We may integrate out , to obtain
(111) |
For a closed system, we can integrate out and express the bulk partition function in the form
(112) |
The sum is over all possible configurations of the (open or closed) worldsheets with a bare surface tension . An intuitive picture in the Maxwell case, and , is that in a time slice this corresponds to magnetic monopole pairs attached by a magnetic string with linearly rising potential, i.e. the magnetic charges are confined in this limit, as expected for a bulk electric condensate which collimates the magnetic field into flux tubes. The characteristic size (“Debye” screening length) of the neutral monopole pairs tends to zero as tends to . Alternatively, we may view this as monopole strings (worldlines) interacting electrostatically through the membranes of the field. For large the membranes are short, meaning that the strings are bound into charge-neutral pairs. This phase persists to all because the entropic gain of dipole strings outweighs their energetic cost at all effective temperatures.
Now consider an open boundary. Recall that the boundary of the dual lattice is “flat”, as shown in Fig. 10, i.e. it has no cells extending into the vacuum. This means that no magnetic charge or magnetic flux can exit the system. More concretely, the constraint , in Eq. 102, in the electric variables is reflected in the dual constraint on every dual boundary link, which enforces that is pure gauge in the dual boundary layer. This means that the boundary action has no dependence. We resolve the boundary constraint as , so that
(113) |
We then define a composite field , which is gauge invariant under the gauge transformations and . This allows us to rewrite the boundary action as
(114) |
To proceed from here we need to be careful about how the and can move between the bulk and boundary layers. We do so in the limit, by combining this with the bulk action Eq. 111. We can then integrate out to generate the magnetic Gauss law, which is enforced on every link. This leaves us with the total action
(115) |
The boundary portion describes monopoles moving in the boundary layer interacting via a non-compact gauge field.
One should be concerned here as to how the from the bulk (edges of ) couple to the boundary. The key to understand what happens here is that (i) the membranes can lie in the boundary layer where they cost zero action, and (ii) this action to be gauge-invariant under shifts of , we must have an additional boundary Gauss law,
(116) |
This implies that, in the limit , magnetic monopoles cannot move between the boundary and the bulk, i.e. there is a -form symmetry on the boundary (0-form in the Maxwell case) corresponding to conservation of boundary magnetic charge. The action Eq. 115 is precisely the dual of the boundary -form U(1) gauge theory Eq. 100, a 3D XY model in the Maxwell case.
The in the boundary layer must be coupled to membranes by the magnetic Gauss law, but these have zero tension if they lie entirely within the boundary layer. This means that the Higgs vortices (magnetic field lines) are effectively not present within the boundary layer. This was by construction, since there were no terms in the original action which would allow for vortices of the Higgs field within the boundary layer. As a result, in the partition function to leading order the boundary and bulk are effectively decoupled from each other. Given any configuration of monopole worldlines , which are closed in the boundary and closed in the bulk, the dominant contribution to the partition function will be for the boundary to be connected to tensionless membranes in the boundary, rather than to be connected to a bulk monopole by a tensionful one. As a result, the monopoles on the boundary can condense at small , which is the dual to the boundary symmetry breaking transition we found in the electric formulation. One can then consider turning on small and doing a strong coupling expansion. This will have the effect of renormalizing the bulk length scale enabling boundary monopoles to extend further into the bulk by extending a magnetic flux tube while preserving the quasi--dimensional nature of the boundary.
IV.3 Summary and Discussion: Higher-Form Case
In this section, we have generalized our results from Section II on -form Abelian-Higgs models to higher form Abelian-Higgs models in spacetime dimensions. In particular, we found that -form gauge field coupled to a -form Higgs field reduces at infinite to a -form gauge theory on the boundary whose dynamical degrees of freedom are half-open Wilson branes terminating on the Higgs field as soon as it enters the bulk. This boundary theory exhibits a confinement-deconfinement transition when (), which spontaneously breaks the -form global matter symmetry at large . Table 1 summarizes this pattern of boundary symmetry breaking. In the marginal cases, (), a generalized Mermin-Wagner theorem prevents the symmetry from breaking, except in the case and , where we predict a BKT boundary transition. There are also precisely the cases where the same mechanism destabilizes the deconfined phase in the bulk [56, 55]. As in the 1-form Abelian and non-Abelian cases studied numerically in this paper, we expect that this boundary phase transition extends into the phase diagram when bulk fluctuations are restored and will end at a bulk critical point, demarcating a boundary between the Higgs and confining regimes.
The general mechanism for this emergent boundary theory at large identified by considering higher-form Abelian-Higgs models revolves around the constraint enforced exactly at infinite . This constraint implies that the Wilson operators which create electric charges attached to electric membranes, for open surfaces , act as the identity. This naïvely indicates that the system is a condensate of electric charge in this limit, i.e. the ground state is a coherent state of the charge annihilation operators. As a consequence, the electric-brane insertion operators for closed surfaces trivialize in the bulk. In the presence of electric-flux-permeable open boundaries, however, operators inserting electric flux through the boundary are immediately screened as soon as they enter the bulk, terminating on the Higgs field, as shown in Fig. 11, and become charged under the matter symmetry. These operators are the dynamical degrees of freedom at play at the boundary which exhibit the matter symmetry breaking. Presumably when is reduced the electric flux can penetrate further into the system before being screened by the electric charge condensate, forming a quasi--dimensional boundary. It would be interesting to explore whether this mechanism can be extended to higher-form non-Abelian theories described by higher-categorical gauge groups [63, 64, 65, 66, 67].
In the second part of this section, we studied the dualized version of these theories, identifying the boundary degrees of freedom and accounting, in the -form case and , for the existence of the dual to the D XY model that we found in terms of direct variables. Magnetic charge moves in the boundary layer, and there is an extra Gauss law in the limit which originates from the constraint that electric charge cannot leave the system. This enforces that magnetic charge cannot leave the boundary into the bulk, and thus the charges on the boundary can condense, leading to the dual phase transition. It would be of interest to study the dual theory in more detail, since it gives a clearer picture of the boundary symmetry breaking in the large- (small ) limit. In particular, it would be worthwhile to pursue Monte Carlo simulations in the dual representation, for which efficient algorithms have been developed [68, 69, 70]. It is also worth re-emphasising the importance of our choice of boundary conditions, which prevented electric charge from leaving the system while allowing electric flux to leave. In the dual theory this led to a flat dual boundary, meaning magnetic charge and flux is always contained in the system and cannot leave. It follows that if one started with flat boundaries, which keep all electric charge and flux inside the system, the dual boundaries would be open, i.e. magnetic charges are kept in the system but magnetic flux can leave. This implies a physical magnetic charge symmetry which can spontaneously break in the confined regime rather than the Higgs regime, with a phase transition on the axis instead of the axis.
V Discussion and Conclusion
V.1 Summary and Outlook
In this work we have explored a variety of models of charged Higgs fields coupled to gauge fields in the fundamental representation, under the imposition of boundary conditions which allow electric flux, but not charge, to exit the system. In a closed system, the gauge field does not have a physical global charge symmetry which can spontaneously break, because a charge is always attached to electric flux, which must end inside the system on another charge. With the “electric-flux-permeable” boundary conditions we consider, the charge sectors and global symmetry become physical, as non-zero bulk charge can be compensated by non-zero boundary flux, and can in principle spontaneously break. This work, focusing on the boundary degrees of freedom, complements work exploring the interplay of global and gauge symmetries in the bulk. See, for example, recent work in Refs. [71, 72, 73, 74, 75].
We have considered the Abelian-Higgs model, two types of non-Abelian Higgs models (with fundamental representation and group-valued Higgs fields), and higher-form Abelian-Higgs models. In terms of the inverse gauge coupling and the matter coupling , all of these models share the common feature of a continuity between an electric-charge confining regime at small , and a Higgs regime at large , which are two ends of one continuous thermodynamic phase, as proven by Fradkin and Shenker [12]. We have demonstrated through a combination of analytical argument and numerical investigation that, under all but a few marginal cases, there is a boundary phase transition which indicates the spontaneous breaking of the matter symmetry at large , i.e. in the Higgs regime.
Ref. [19] predicted a boundary phase transition in the case of gauge group , and verified it using DMRG in dimensions. Ref. [20] predicted the boundary transition for a magnetic monopole-free Abelian U(1) Higgs model without numerics. Using lattice gauge theory Monte Carlo simulations, we have explored a wider range of models, both for Abelian and non-Abelian gauge groups. In Section II, we have numerically verified and studied the boundary transition for gauge group U(1) (with monopoles) in dimensions, which exhibits a boundary XY transition. In Section III, we extended this to non-Abelian gauge groups, performing numerical simulations for both and gauge theories coupled to fundamental Higgs fields in , and considered generalizations to , , and general gauge groups. Lastly, in Section IV we studied higher-form generalizations of Abelian-Higgs models, and demonstrated that the corresponding higher-form matter symmetry can spontaneously break at the boundary. In all of our numerical simulations, the nature of the boundary phase transitions deduced from the numerical data conforms to predictions obtained by studying the limit.
We expect the considerations in this work to extend further to all gauge groups and different Higgs representations. Notably we have not discussed higher-rank gauge theories coupled to scalar matter that are connected to fractonic excitations, but here too one may preserve the Gauss law for tensor fields on the boundary and expect a U(1) global symmetry that, for certain classes of such models, can be broken spontaneously. It would also be of interest to extend these results to discrete non-Abelian gauge groups, and non-Abelian higher-form gauge theories with higher-categorical gauge groups.
Already in the seminal work of Fradkin and Shenker [12] it was understood that in gauge theories with fundamental Higgs matter, under most circumstances thermodynamic quantities exhibit no singularities along paths between Higgs and confined regimes, meaning that they form the same bulk phase of matter. The boundary spontaneous symmetry breaking investigated in this paper does not contradict this result, because it does not define a precise bulk phase boundary between the Higgs and confined regimes. Indeed, as we observed, one may tune the boundary coupling while preserving all the bulk properties to shift the location of the boundary phase transition.
V.2 Higgs=SPT
It would be remiss of us to conclude this paper without a description of some of the work that motivated this study—namely the papers [18, 19, 20] discussing the relationship of Higgs phases to symmetry-protected topological (SPT) phases. We briefly summarize the main findings therein and comment on how our results bear on the generality of this relationship.
An SPT phase is a state of matter that cannot be adiabatically connected to a trivial phase under local symmetry-preserving perturbations without closing a gap. In general, a trivial phase can be reached without gap closure only if those symmetries are broken. A classic example of an SPT phase in an interacting lattice model is the spin one-half chain—the cluster model—whose nontrivial topology is protected by symmetry. A consequence of the topological nature of the phase is the presence of gapless states localized at the boundary of an open chain. In [18] it was realized that the cluster SPT order emerges in the Higgs phase of the -gauged Ising chain in the absence of the electric string tension term. Moreover, the fermionic SPT order of the gauged Kitaev chain in the Higgs regime was identified.
Many 2+1 dimensional SPT phases are known and one of the latest to be added to the inventory of phases lives in the phase diagram of the dimensional gauge theory coupled minimally to Ising matter. Using the notation of this paper, where the gauge coupling is denoted by and the matter coupling by , the famous deconfined toric code model lives at and . It was argued in Ref. [19] that deep in the Higgs limit and (where both gauge and matter degrees of freedom are frozen in the bulk and the theory has a magnetic 1-form Ising symmetry) the model maps to the two-dimensional cluster SPT protected by the form magnetic symmetry and a -form Ising matter symmetry 101010 For earlier work on SPTs protected by generalized symmetries see for example Refs. [67, 76] . The Ising matter symmetry is trivialized in the bulk by gauging, but survives as a global symmetry in presence of a boundary. This then has gapless surface states localized on symmetry-preserving boundaries. It was further shown that this phase persists to large but finite at , where both protecting symmetries remain intact. Going to finite breaks the magnetic -form symmetry yet at only the boundary degrees of freedom can fluctuate and are governed by an effective one-dimensional transverse field Ising model that is in the spontaneously broken phase for large . The gapless edge modes thus originate from the Ising degeneracy of this phase. It is expected [77] that the form symmetry survives in the infrared even though it is explicitly broken away from . In other words one can think about it as a low-energy emergent symmetry. Protected by the Ising matter -form symmetry and the emergent magnetic -form symmetry, the low-energy edge modes thus should survive within some window of parameters even away from limit. It was verified numerically that there is indeed a phase with a gapless boundary within the entire Higgs regime of the model. The phase diagram of this model bears a great deal of resemblance to our Fig. 1 with a phase transition line originating at a triple point at the innermost corner of the deconfined phase and ending on a critical endpoint. Within the resolution of the DMRG numerics of Ref. [19] the phase boundary to the boundary gapless phase extends from the critical endpoint.
A follow-up work extended the results to Higgs phases [20]. In 3+1D the relevant symmetries are the matter U() -form symmetry associated with charge conservation and an exact magnetic -form symmetry, whose conserved current . We may couple the currents to background gauge fields—one-form for the matter current and two-form for . Deep in the Higgs regime, the topological response was found to be
(117) |
This SPT response is a consequence of a mixed anomaly between the matter and magnetic symmetries, which is cured by adding boundary degrees of freedom. For the D Higgs phase the boundary theory takes the form
(118) |
for the conjugate pair: a compact scalar field and a compact gauge field . This describes a boundary superfluid phase. The continuum discussion was supplemented in Ref. [20] with a lattice Villain formulation, where monopoles are under control and the magnetic -form symmetry is preserved.
We have briefly reviewed the results of [18, 19, 20] on connections between SPT and Higgs phases. But what bearing, if any, do the results we have reported have on this scenario? After all, at no stage of this work did we invoke bulk topological arguments to understand the existence of the boundary phase transition. It was possible to understand the global symmetry breaking purely at the level of the boundary theory. Even so, given the arguments sketched above, one might seek a description in terms of an SPT order since it would shed light on how robust the boundary symmetry breaking actually is.
In the 4D compact Wilson model coupled to the fundamental Higgs discussed in Sec. II the 1-form magnetic symmetry is broken explicitly by monopoles. However, in the limit , where electromagnetic fields freeze, monopole creations/annihilations are suppressed and the magnetic 1-form symmetry emerges. Based on general arguments of [20], the interplay of the 0-form matter and 1-form magnetic symmetries should give rise to a gapless symmetry-preserving boundary. Sure enough we have established the presence of such a boundary phase numerically. Moreover, in agreement with arguments of [19, 20], the gapless nature of the boundary survives beyond the symmetric limit, as is clear in the limit where the boundary maps to a three dimensional XY model. This much follows from ideas put forward in [19, 20]. However, in contrast to the models studied in [19, 20], our bulk model in the limit becomes gapless for large since it can be thought of (in the unitary gauge) as the ordered phase of a four-dimensional XY model.111111 In contrast, approaching from the limit apparently preserves a finite gap. Obviously, it cannot support a gapped SPT phase. It may be interesting to explore this regime of the phase diagram more closely and determine whether it exhibits a gapless SPT order. It is also important and interesting to understand better the role of the coupling from the point of view of the bulk SPT.
The physics of the non-Abelian gauge theories with Higgs matter gives rise to another puzzle. In Sec. III we have argued that in the limit boundary is expected to be gapless rather generically. If the Higgs=SPT scenario is general enough to encompass the non-Abelian cases too, one should be able to identify the protecting (higher-form) bulk symmetries. Pure gauge theory has a center 1-form symmetry and the electric field (Wilson) lines in transform non-trivally under the center of the group. The center symmetry, however, is broken explicitly in the presence of the fundamental Higgs matter. This should be analogous to the breaking of electric one-form symmetry in the Abelian case discussed above. But, in the limit the Abelian theory also has a magnetic -form symmetry which, together with the matter symmetry, protects the SPT order. In contrast, ’t Hooft lines transform trivially under the center of the group. Currently it is unclear to us whether there is a generalized symmetry present at for non-Abelian gauge group that is an analog of the magnetic one-form symmetry emerging in the Abelian case. Whether SPT order is responsible for the gapless nature of the symmetry-broken boundary phase in the non-Abelian gauge theories coupled to Higgs matter is an open question.
Taking a broader perspective, it would be interesting to search for mutual anomalies between the bulk matter (equivalently, boundary flux) symmetry and any emergent higher-form symmetries which may be present in the limit arising from suppressing homotopy defects of the gauge field.
Acknowledgements.
PM and PR acknowledge useful discussions with Pedro Bicudo and Nuno Cardoso as a part of a related collaboration. KTKC acknowledges Chris Hooley for useful discussion. S.M. is supported by Vetenskapsrådet (grant number 2021-03685), Nordita and STINT. This work was in part supported by the Deutsche Forschungsgemeinschaft under the cluster of excellence ct.qmat (EXC-2147, project number 390858490).Appendix A Discrete Differential Calculus for Abelian Fields
We utilize notation which mimics continuum differential forms on a lattice, borrowed from algebraic topology [78] and used extensively in Section IV. Pedagogical treatments can be found in [29, 48]. Fields are described (locally) as differential forms (i.e. anti-symmetric tensors), whose primary property is that they can be integrated over surfaces. For our purposes we can think of them simply as functions of such surfaces, i.e. if is a differential -form and an oriented -dimensional surface, acts on by integration
(119) |
which results in some number. There are two main properties that we wish to preserve on the lattice: reversing the orientation of changes the sign of the integral, and if is divided into a collection of smaller parts the total integral is the sum of the integrals over the parts, i.e.
(120) |
where denotes the reversed orientation. Lastly, we can naturally define the derivative of a -form to be a -form whose value on a -dimensional surface is defined by Stoke’s theorem,
(121) |
where denotes the boundary of .
In this paper, our space(time) is a -dimensional cubical cell complex—a collection of -cells for , i.e. vertices, links, plaquettes, cubes, hypercubes, etc., where -cells are glued along their -dimensional boundary cells. We denote the collection of all cells by and the collection of all -cells by . Because we also consider open boundaries in the form of Fig. 2, we let denote just the bulk -cells and denote the -cells in the boundary layer which touch the vacuum.
The cells naturally provide the integration surfaces once equipped with an orientation. It is natural to define the possible integration surfaces therefore as integer weighted linear combinations of oriented -cells, call -chains. The signs of the integer coefficients determine the orientations and their magnitudes determine how many times to integrate over each cell. Since we can formally add such chains together, they form an Abelian group, . The structure of the cell complex is contained in the boundary relation,
(122) |
which distributes over the linear combinations on -cells, and sends each oriented -cell to the linear combination of its oriented boundary -cells.
The sensible lattice analog of a differential form, i.e. a discrete -form, also called a -cochain, is a function which (i) maps chains to numbers, Eq. 119, and (ii) disributes over linear combinations, Eq. 120. In full generality, a discrete -form is a linear map from chains to elements of any Abelian group ,
(123) |
In practice what this means is that a disrete -form is defined by its value on each oriented -cell , and satsfies , where is understood as the inverse operation in the Abelian group . The space of discrete -forms is denoted . Since we already have a natural notion of the boundary operation on chains, we can define a natural exterior derivative operation on cochains, , according to Stoke’s theorem, i.e.
(124) |
where is a -form, is a -form, and is a -chain.
Denoting the coefficients in a -chain by
(125) |
where each -cell is summed once with a fixed orientation, we can define a natural inner product on -chains as
(126) |
Using this, we can define an adjoint of the boundary operator, which we call the coboundary,121212 Note that this differs from the algebraic topology terminology, where the discrete exterior derivative is often called the coboundary.
(127) |
In particular, the coboundary of a single -cell is
(128) |
By choosing to orient all the in the sum such that or , we can read this to say that the coboundary of an oriented -cell is the sum of all oriented -cells containing in their positively oriented boundary. For example, the coboundary of a point is the set of oriented links which terminate at it, the coboundary of a link is the set of plaquettes touching it, oriented so that circulates in the same direction as is oriented, etc.
The couboundary defines a co-exterior derivative, via a “co-Stoke’s theorem”,
(129) |
which reduces the degree of forms. In the case that is , , or , we can also define an inner product for -forms,
(130) |
where denotes complex conjugation. It is easy to see that the codifferential is adjoint to the differential,
(131) |
Lastly, we review a bit of basic algebraic topology terminology used in Section IV. The boundary operation is nilpotent, , and thus defines an exact sequence of maps,
(132) |
This allows us to define the homology groups, The interpretation of these quotient groups is that they classify the non-contractible -dimensional surfaces. Here, is generated by the set of surfaces without boundaries (called cycles), while is generated by the set of surfaces which are boundaries of a -dimensional volume, and are therefore contractible. The elements of are equivalence classes of surfaces which differ by a boundary, i.e. homology classes.
The dual of the homology classes on the differential form side are the cohomology classes, which are defined by the discrete equivalent of the de Rham complex:
(133) |
as
(134) |
Here is the set of locally-constant -forms (called closed forms), and is the set of -forms which are gradients of -forms (called exact forms). Every exact form is closed, but there can exist closed forms which are not exact, which are classified by cohomology. Cohomology classes are equivalence classes of closed forms which differ by an exact form, just as homology classes are equivalence classes of closed surfaces which differ by a boundary. Note that exact forms integrate to zero by Stoke’s theorem on any closed surface, while closed forms need only integrate to zero on closed surfaces that are boundaries. Forms with non-trivial cohomology class integrate non-trivially on non-contractible surfaces, defining a pairing between cohomology and homology classes.
Appendix B U(1) Bulk Phase Diagram: Limits, Symmetries and Monte Carlo
Here we provide a brief review of the bulk phase diagram of the 4D Abelian-Higgs model, described by the action Eq. 8, the Hamiltonian Eq. 10, and dual action Eq. 108, by considering the various limits and undertaking numerical Monte Carlo simulation in the action formulation. The phase diagram is sketched in Fig. 1, which conveniently summarizes the discussion that follows. In particular, it is well-known that there are only two distinct phases [12], the Coulomb phase and the Higgs-confined phase. According to Eq. 102, this is a theory of electric strings terminating on electric point charges. The dual magnetic description is in terms of magnetic strings terminating on magnetic point charges (vortices and monopoles of the Higgs and gauge fields, respectively). The basic structure of the phase diagram can be deduced by consider each of the four limits.
Pure gauge limit (): In the limit the gap of the electric point charges diverges, and the theory reduces to 4D gauge theory. In this limit, the system has a global 1-form symmetry, with , called electric symmetry. This symmetry yields corresponds to electric strings forming closed loops. This theory has two phases, a confined phase at small (strong coupling), where the system is gapped and electric strings cost energy proportional to their length; and a deconfined phase at large (weak coupling), where the electric strings condense, the electric symmetry is spontaneously broken, and the system has a gapless photon excitation. Turning on a small explicitly breaks the 1-form symmetry by introducing gapped electric point charges at which open electric strings end. Because the charges are strongly gapped, qualitatively speaking we expect the 1-form symmetry to re-emerge at low energies below the charge gap, thus allowing the deconfined phase to extend to finite .
Frozen gauge limit (): In this limit the gauge fields are completely trivialized by the constraint , and equivalently the mass of the magnetic monopoles diverges. We can choose a gauge where and the action turns into that of a XY model, whose dual description is a gas of vortex strings with Coulomb interactions. This theory has a 1-form symmetry corresponding to the closure of these magnetic strings and the corresponding absence of magnetic monopoles, called magnetic symmetry. From the XY model we deduce a gapless superfluid phase at large and a gapped phase at small separated by a second-order phase transition, but this description is not gauge-invariant. The gauge-invariant statement is that at small the magnetic strings (disorder operators o the XY model) condense, spontaneously breaking the magnetic symmetry. For large but finite the magnetic monopoles explicitly break the magnetic symmetry, though one expects it to be effectively restored at low energies below the monopole gap. The result is that the magnetic symmetry is spontaneously broken in the gapless Coulomb phase, while the superfluid at large is gapped out by the Higgs mechanism.
Strong coupling limit (): In this limit the curvature of the gauge field is not penalized, and we can say that the magnetic monopoles are maximally proliferated. If we fix to unitary gauge (), the action becomes , reducing to a set of completely disconnected link variables. Thus the system is in a trivial phase for all , and there is no bulk phase transition on this line.
Infinite Higgs coupling limit (): In this limit we have the constraint . This implies that the electric charge creation and annihilation operators, Eq. 14, act as the identity, and the system can be described as an electric condensate. The bulk has no dynamics and is completely frozen, which can be seen by fixing to unitary gauge (), in which and the Higgs field is frozen. There is therefore no bulk phase transition along this line. For large but finite the bulk action can be roughly understood as the proca-type action, Eq. 23, describing a massive 1-form field.
This completes the general outline of the phase diagram in the vicinity of the edges in Fig. 1. The only question that remains is how the two transitions at small and large reach each other to separate the confined phase (in which the electric and magnetic symmetries are emergent and spontaneously broken) from the Higgs-confined phase. Figure 12 shows results from Monte Carlo simulations measuring the average Wilson plaquette and Wilson link, which fills in the remainder of the phase diagram schematically shown in Fig. 1. The two transitions extend as first-order transition lines and meet at a triple point in the vicinity of an . A third first order line extends from the this triple point towards smaller and larger which ends at a critical endpoint. We show evidence for this first-order line in Fig. 12(c).
References
- Wegner [1971] F. J. Wegner, Duality in Generalized Ising Models and Phase Transitions without Local Order Parameters, J. Math. Phys. 12, 2259 (1971).
- Wilson [1974] K. G. Wilson, Confinement of quarks, Phys. Rev. D 10, 2445 (1974).
- Gaiotto et al. [2015] D. Gaiotto, A. Kapustin, N. Seiberg, and B. Willett, Generalized global symmetries, J. High Energ. Phys. 2015 (2), 172.
- Lake [2018] E. Lake, Higher-form symmetries and spontaneous symmetry breaking (2018), arxiv:1802.07747 .
- Delacrétaz et al. [2020] L. V. Delacrétaz, D. M. Hofman, and G. Mathys, Superfluids as higher-form anomalies, SciPost Physics 8, 047 (2020).
- McGreevy [2023] J. McGreevy, Generalized Symmetries in Condensed Matter, Annual Review of Condensed Matter Physics 14, 57 (2023).
- Anderson [1963] P. W. Anderson, Plasmons, Gauge Invariance, and Mass, Phys. Rev. 130, 439 (1963).
- Englert and Brout [1964] F. Englert and R. Brout, Broken Symmetry and the Mass of Gauge Vector Mesons, Phys. Rev. Lett. 13, 321 (1964).
- Higgs [1964] P. W. Higgs, Broken Symmetries and the Masses of Gauge Bosons, Phys. Rev. Lett. 13, 508 (1964).
- Guralnik et al. [1964] G. S. Guralnik, C. R. Hagen, and T. W. B. Kibble, Global Conservation Laws and Massless Particles, Phys. Rev. Lett. 13, 585 (1964).
- Osterwalder and Seiler [1978] K. Osterwalder and E. Seiler, Gauge field theories on a lattice, Annals of Physics 110, 440 (1978).
- Fradkin and Shenker [1979] E. Fradkin and S. H. Shenker, Phase diagrams of lattice gauge theories with Higgs fields, Phys. Rev. D 19, 3682 (1979).
- Greensite [1983] J. Greensite, An Introduction to the Confinement Problem, 2nd ed. (Springer Cham, 1983).
- Caudy and Greensite [2008] W. Caudy and J. Greensite, Ambiguity of spontaneously broken gauge symmetry, Phys. Rev. D 78, 025018 (2008).
- Cherman et al. [2019] A. Cherman, S. Sen, and L. G. Yaffe, Anyonic particle-vortex statistics and the nature of dense quark matter, Phys. Rev. D 100, 034015 (2019).
- Cherman et al. [2020] A. Cherman, T. Jacobson, S. Sen, and L. G. Yaffe, Higgs-confinement phase transitions with fundamental representation matter, Phys. Rev. D 102, 105021 (2020).
- Cherman et al. [2024] A. Cherman, T. Jacobson, S. Sen, and L. G. Yaffe, Line operators, vortex statistics, and Higgs versus confinement dynamics (2024), arxiv:2401.17489 .
- Borla et al. [2021] U. Borla, R. Verresen, J. Shah, and S. Moroz, Gauging the Kitaev chain, SciPost Physics 10, 148 (2021).
- Verresen et al. [2022] R. Verresen, U. Borla, A. Vishwanath, S. Moroz, and R. Thorngren, Higgs Condensates are Symmetry-Protected Topological Phases: I. Discrete Symmetries (2022), arxiv:2211.01376 .
- Thorngren et al. [2023] R. Thorngren, T. Rakovszky, R. Verresen, and A. Vishwanath, Higgs Condensates are Symmetry-Protected Topological Phases: II. U(1) Gauge Theory and Superconductors (2023), arxiv:2303.08136 .
- Fredenhagen and Marcu [1986] K. Fredenhagen and M. Marcu, Confinement criterion for QCD with dynamical quarks, Phys. Rev. Lett. 56, 223 (1986).
- Gregor et al. [2011] K. Gregor, D. A. Huse, R. Moessner, and S. L. Sondhi, Diagnosing deconfinement and topological order, New J. Phys. 13, 025009 (2011).
- Kennedy and King [1986] T. Kennedy and C. King, Spontaneous symmetry breakdown in the abelian Higgs model, Commun.Math. Phys. 104, 327 (1986).
- Kogut [1979] J. B. Kogut, An introduction to lattice gauge theory and spin systems, Rev. Mod. Phys. 51, 659 (1979).
- Creutz [1984] M. Creutz, Quarks, Gluons and Lattices (Cambridge University Press, Cambridge, 1984).
- Banks et al. [1977] T. Banks, R. Myerson, and J. Kogut, Phase transitions in Abelian lattice gauge theories, Nuclear Physics B 129, 493 (1977).
- Savit [1980] R. Savit, Duality in field theory and statistical systems, Rev. Mod. Phys. 52, 453 (1980).
- Savit [1982] R. Savit, Duality transformations for general abelian systems, Nuclear Physics B 200, 233 (1982).
- Drühl and Wagner [1982] K. Drühl and H. Wagner, Algebraic formulation of duality transformations for abelian lattice models, Annals of Physics 141, 225 (1982).
- Ranft et al. [1983] J. Ranft, J. Kripfganz, and G. Ranft, Phase structure, magnetic monopoles, and vortices in the lattice Abelian Higgs model, Phys. Rev. D 28, 360 (1983).
- Gottlob and Hasenbusch [1993] A. P. Gottlob and M. Hasenbusch, Critical behaviour of the 3D XY-model: A Monte Carlo study, Physica A: Statistical Mechanics and its Applications 201, 593 (1993).
- Hansson et al. [2004] T. H. Hansson, V. Oganesyan, and S. L. Sondhi, Superconductors are topologically ordered, Annals of Physics 313, 497 (2004).
- Polyakov [2017] A. M. Polyakov, Gauge Fields and Strings (Routledge, London, 2017).
- Kogut and Susskind [1975] J. Kogut and L. Susskind, Hamiltonian formulation of Wilson’s lattice gauge theories, Phys. Rev. D 11, 395 (1975).
- Susskind and Kogut [1976] L. Susskind and J. Kogut, XIV. New ideas about confinement, Physics Reports 23, 348 (1976).
- Bonati et al. [2010] C. Bonati, G. Cossu, M. D’Elia, and A. Di Giacomo, Phase diagram of the lattice SU(2) Higgs model, Nuclear Physics B 828, 390 (2010).
- Kanaya and Kaya [1995] K. Kanaya and S. Kaya, Critical exponents of a three-dimensional O(4) spin model, Phys. Rev. D 51, 2404 (1995).
- Sordi and Tremblay [2023] G. Sordi and A.-M. S. Tremblay, Introducing the concept of Widom line in the QCD phase diagram (2023), arxiv:2312.12401 .
- Chandrasekharan and Wiese [2004] S. Chandrasekharan and U. J. Wiese, An introduction to chiral symmetry on the lattice, Progress in Particle and Nuclear Physics 53, 373 (2004).
- Kogut et al. [1982] J. B. Kogut, M. Snow, and M. Stone, Mean field and Monte Carlo studies of SU(N) chiral models in three dimensions, Nuclear Physics B 200, 211 (1982).
- Faddeev and Reshetikhin [1986] L. D. Faddeev and N. Y. Reshetikhin, Integrability of the principal chiral field model in 1 + 1 dimension, Annals of Physics 167, 227 (1986).
- Kos et al. [2014] F. Kos, D. Poland, and D. Simmons-Duffin, Bootstrapping the O(N) Vector Models, J. High Energ. Phys. 2014 (6), 91.
- Binney et al. [1992] J. J. Binney, N. J. Dowrick, A. J. Fisher, and M. E. J. Newman, The Theory of Critical Phenomena: An Introduction to the Renormalization Group (Oxford University Press, Oxford, New York, 1992).
- Li [1974] L.-F. Li, Group theory of the spontaneously broken gauge symmetries, Phys. Rev. D 9, 1723 (1974).
- Zou et al. [2021] L. Zou, Y.-C. He, and C. Wang, Stiefel Liquids: Possible Non-Lagrangian Quantum Criticality from Intertwined Orders, Phys. Rev. X 11, 031043 (2021).
- Henneaux and Teitelboim [1986] M. Henneaux and C. Teitelboim, P-Form electrodynamics, Foundations of Physics 16, 593 (1986).
- Teitelboim [1986] C. Teitelboim, Gauge invariance for extended objects, Physics Letters B 167, 63 (1986).
- Wise [2006] D. K. Wise, P-form electromagnetism on discrete spacetimes*, Class. Quantum Grav. 23, 5129 (2006).
- Omero et al. [1983] C. Omero, P. A. Marchetti, and A. Maritan, Gauge differential form theories on the lattice, J. Phys. A: Math. Gen. 16, 1465 (1983).
- Kalb and Ramond [1974] M. Kalb and P. Ramond, Classical direct interstring action, Phys. Rev. D 9, 2273 (1974).
- Davis and Shellard [1989] R. L. Davis and E. P. S. Shellard, Global strings and superfluid vortices, Phys. Rev. Lett. 63, 2021 (1989).
- Franz [2007] M. Franz, Vortex-boson duality in four space-time dimensions, EPL 77, 47005 (2007).
- Chung and Gingras [2023] K. T. K. Chung and M. J. P. Gingras, 2-Form U(1) Spin Liquids: Classical Model and Quantum Aspects (2023), arxiv:2310.17607 .
- Nepomechie [1985] R. I. Nepomechie, Magnetic monopoles from antisymmetric tensor gauge fields, Phys. Rev. D 31, 1921 (1985).
- Rey [1989] S.-J. Rey, Higgs mechanism for Kalb-Ramond gauge field, Phys. Rev. D 40, 3396 (1989).
- Orland [1982] P. Orland, Instantons and disorder in antisymmetric tensor gauge fields, Nuclear Physics B 205, 107 (1982).
- Pearson [1982] R. B. Pearson, Phase structure of antisymmetric tensor gauge fields, Phys. Rev. D 26, 2013 (1982).
- Polyakov [1975] A. M. Polyakov, Compact gauge fields and the infrared catastrophe, Physics Letters B 59, 82 (1975).
- Quevedo and Trugenberger [1997a] F. Quevedo and C. A. Trugenberger, Condensation of p-Branes and Generalized Higgs/Confinement Duality, Int. J. Mod. Phys. A 12, 1227 (1997a).
- Quevedo and Trugenberger [1997b] F. Quevedo and C. A. Trugenberger, Phases of antisymmetric tensor field theories, Nuclear Physics B 501, 143 (1997b).
- Iqbal and McGreevy [2022] N. Iqbal and J. McGreevy, Mean string field theory: Landau-Ginzburg theory for 1-form symmetries, SciPost Physics 13, 114 (2022).
- Sulejmanpasic and Gattringer [2019] T. Sulejmanpasic and C. Gattringer, Abelian gauge theories on the lattice: Theta-terms and compact gauge theory with(out) monopoles, Nuclear Physics B 943, 114616 (2019).
- Pfeiffer [2003] H. Pfeiffer, Higher gauge theory and a non-Abelian generalization of 2-form electrodynamics, Annals of Physics 308, 447 (2003).
- Baez and Huerta [2011] J. C. Baez and J. Huerta, An invitation to higher gauge theory, Gen Relativ Gravit 43, 2335 (2011).
- Rey and Sugino [2010] S.-J. Rey and F. Sugino, A Nonperturbative Proposal for Nonabelian Tensor Gauge Theory and Dynamical Quantum Yang-Baxter Maps (2010), arxiv:1002.4636 .
- Lipstein and Reid-Edwards [2014] A. E. Lipstein and R. A. Reid-Edwards, Lattice gerbe theory, J. High Energ. Phys. 2014 (9), 34.
- Kapustin and Thorngren [2017] A. Kapustin and R. Thorngren, Higher Symmetry and Gapped Phases of Gauge Theories, in Algebra, Geometry, and Physics in the 21st Century: Kontsevich Festschrift, Progress in Mathematics, edited by D. Auroux, L. Katzarkov, T. Pantev, Y. Soibelman, and Y. Tschinkel (Springer International Publishing, Cham, 2017) pp. 177–202.
- Mercado et al. [2013a] Y. D. Mercado, C. Gattringer, and A. Schmidt, Surface worm algorithm for abelian Gauge–Higgs systems on the lattice, Computer Physics Communications 184, 1535 (2013a).
- Mercado et al. [2013b] Y. D. Mercado, C. Gattringer, and A. Schmidt, Dual Lattice Simulation of the Abelian Gauge-Higgs Model at Finite Density: An Exploratory Proof of Concept Study, Phys. Rev. Lett. 111, 141601 (2013b).
- Schmidt et al. [2012] A. Schmidt, Y. D. Mercado, and C. Gattringer, Monte Carlo simulation of abelian gauge-Higgs lattice models using dual representation (2012), arxiv:1211.1573 .
- Pelissetto and Vicari [2019] A. Pelissetto and E. Vicari, Multicomponent compact Abelian-Higgs lattice models, Phys. Rev. E 100, 042134 (2019).
- Bonati et al. [2021a] C. Bonati, A. Pelissetto, and E. Vicari, Global symmetry breaking in gauge theories: The case of multiflavor scalar chromodynamics (2021a), arxiv:2110.05341 .
- Bonati et al. [2021b] C. Bonati, A. Pelissetto, and E. Vicari, Breaking of the gauge symmetry in lattice gauge theories, Phys. Rev. Lett. 127, 091601 (2021b).
- Bonati et al. [2022a] C. Bonati, A. Pelissetto, and E. Vicari, Critical behaviors of lattice U(1) gauge models and three-dimensional Abelian-Higgs gauge field theory, Phys. Rev. B 105, 085112 (2022a).
- Bonati et al. [2022b] C. Bonati, A. Pelissetto, and E. Vicari, Scalar gauge-Higgs models with discrete Abelian symmetry groups, Phys. Rev. E 105, 054132 (2022b).
- Yoshida [2016] B. Yoshida, Topological phases with generalized global symmetries, Phys. Rev. B 93, 155131 (2016).
- Pace and Wen [2023] S. D. Pace and X.-G. Wen, Exact emergent higher-form symmetries in bosonic lattice models (2023), arxiv:2301.05261 .
- Hatcher [2002] A. Hatcher, Algebraic Topology (Cambridge University Press, 2002).