Electric-magnetic duality in the quantum double models of topological orders with gapped boundaries

We generalize the Electric-magnetic (EM) duality in the quantum double (QD) models to the case of topological orders with gapped boundaries. We also map the QD models with boundaries to the Levin-Wen (LW) models with boundaries. To this end, we Fourier transform and rewrite the extended QD model with a finite gauge group $G$ on a trivalent lattice with a boundary. Gapped boundary conditions of the model before the transformation are known to be characterized by the subgroups $K \subseteq G$. We find that after the transformation, the boundary conditions are then characterized by the Frobenius algebras $A_{G,K}$ in $\mathrm{Rep}_G$. An $A_{G,K}$ is the dual space of the quotient of the group algebra of $G$ over that of $K$, and $\mathrm{Rep}_G$ is the category of the representations of $G$. The EM duality on the boundary is revealed by mapping the $K$'s to $A_{G,K}$'s. We also show that our transformed extended QD model can be mapped to an extended LW model on the same lattice via enlarging the Hilbert space of the extended LW model. Moreover, our transformed extended QD model elucidates the phenomenon of anyon splitting in anyon condensation.


Introduction
Two-dimensional phases of matter with intrinsic topological orders, or topological phases for short, can be well studied by effective topological field theories, whose Hamiltonian extensions are exactly solvable lattice models. Two major families of such models are the quantum double (QD) models [1] and the string-net models or Levin-Wen (LW) models [2]. The QD models have been generalized to be the twisted QD models [3,4], and the LW models have also been generalized similarly [5]. In this paper, we shall not deal with such generalizations.
A QD model is a lattice gauge theory with a finite gauge group G as its input data. The Hamiltonian of the model converts the input data to an output data-the quantum double D(G) of good quantum numbers. The quasiparticle excitations of the topological phase described by the model, namely the anyons, carry representations ([c], µ) of D(G), where [c] labels the conjugacy classes of G, and µ labels the irreducible representations of the centralizer of [c] in G. Anyons of the type ([e], µ) in which e ∈ G is the identity element are the pure charges, which live at the vertices of Γ; the anyons of the type ([c], 0) with 0 being the trivial representation of the centralizer of [c] are the pure fluxes, which live in the plaquettes of Γ; the anyons of the mixed type are the dyons. Hence, intuitively, the QD models would exhibit an explicit EM duality. Indeed, to every QD model on a lattice Γ, there corresponds a QD model on the dual lattice Γ * , such that the charges (fluxes) on Γ are the fluxes (charges) on Γ * , and vice versa [6]. This EM duality is immediately understood in the cases of Abelian groups. In such a case, the irreducible representations of G are all 1-dimensional and form a group isomorphic to G itself. Denote the set of all irreducible representations of G by Rep G . While the QD model on Γ takes G as its input data, the dual model on Γ * takes Rep G ∼ = G as the input data. That is, the dual model is truly a QD model, which by definition has a finite group as its input data. If G is non-Abelian, however, Rep G cannot be a group and thus cannot serve as the input data of a QD model. In such cases, the QD models must be generalized to allow Hopf algebras as their input data to exhibit the EM duality [6].
On the other hand, a QD model with a finite group G can also be mapped via Fourier transform to an LW model with Rep G as its input data on the same lattice [7]. A subtlety is that when G is non-Abelian, a truncation of the Hilbert space of the QD model must be done to complete the mapping [7] unless one enlarges the Hilbert spaces of the LW model [6,8]. Via this mapping, the LW models with input data Rep G bear an EM duality as well.
The EM duality and the aforementioned mapping to the LW models of the QD models are restricted to topological phases on closed surfaces only. Nevertheless, topological phases on surfaces with boundaries are of more practical and theoretical importance because 1) materials with boundaries are much more available than the closed ones, 2)boundary modes are easier to measure experimentally, and 3) a dynamical theory of topological phases is incomplete if unable to encompass different boundary conditions. Recently the QD models and LW models have been extended to be defined on lattices with boundaries by adding appropriate boundary Hamiltonians and are called the extended QD models [9,10] and the extended LW models [11,12]. This motivates us to examine in this paper whether and how the extended QD models still possess an EM duality, in particular along the boundaries, and can be mapped to the extended LW models.
Because of topological invariance and for convenience, we consider a trivalent lattice Γ with a single boundary (see Fig. 1). The input data of an extended QD model on Γ is still a finite group G; however, the boundary Hamiltonian projects the boundary degrees of freedom into a subgroup K ⊆ G. Each K characterizes a gapped boundary condition. We first Fourier transform the Hilbert space on Γ. The Fourier-transformed Hilbert space basis begs us to rewrite the Fourier-transformed model on a slightly different latticeΓ, which modifies Γ by adding near to each vertex of Γ a dangling edge (see Fig. 2(d)). While the bulk degrees of freedom after the Fourier transform become Rep G , the boundary degrees of freedom would be projected by the Fourier-transformed boundary Hamiltonian into a Frobenius algebra A G,K = (C[G]/C[K]) * , the dual space of the quotient of the group algebra of G over that of a given K. When G is Abelian, the emergent Frobenius algebra A G,K happens to be an Abelian group too; hence, the boundary EM duality can be understood as one between the Fourier-transformed extended QD model onΓ with boundary condition specified by A G,K and the extended QD model with boundary condition specified by K. When G is non-Abelian, A G,K cease to being a group but truly an algebra.
We also show that our Fourier-transformed extended QD model onΓ can be mapped to an extended LW model on the same lattice. In doing so, instead of truncating the Hilbert space of the extended QD model, we enlarge the extended LW model. This enlargement is necessary because the Hilbert space of the original extended LW model is too small to contain the full spectrum of excited states with charges; it has already been done for the original LW model on a closed surface [8]. Since the extended QD model and the extended LW model are Hamiltonian extension of the extended Dijkgraaf-Witten and extended Turaev-Viro types of topological field theories, our results also offer a correspondence between the two types of topological field theories.
The EM duality on the boundary and the mapping to the extended LW model can be revealed by mapping a K ⊆ G to (C[G]/C[K]) * . Three cases of such mappings are listed in Table 1. We explain them in order. (c) For a nontrivial subgroup K, the corresponding Frobenius algebra is Another understanding of the boundary EM duality would require generalizing the entire extended model to one with input data being a Hopf algebra, similar to what is done only to the bulk as in Ref. [6]. We shall not discuss such generalizations in this paper but offer a justification of the boundary EM duality when G is non-Abelian.
According to the mechanism of anyon condensation [13], a gapped boundary condition of a topological phase can be accounted for by a condensate at the boundary formed by certain types of anyons from the bulk [14,15]. To condense at the boundary, certain types of bulk anyons would have to first split into a few pieces, several or all of which are allowed to condense at the boundary, depending on the structure of the condensate. Current understandings of this phenomenon are categorical and abstract. It would be interesting to understand the phenomenon of splitting and partial condensation based on concrete lattice models of topological phases. As we will show, our Fourier transform of the extended QD model explains this phenomenon in terms of solely the input data of the model.
To facilitate our studies in the paper, we introduce also a graphical tool of group representation theory. We provide concrete examples, one for Abelian groups G and one for the non-Abelian group S 3 , to illustrate our results.
Our paper is organized as follows. Section 2 reviews the extended QD model. Section 3 Fourier-transforms and rewrites the extended QD model. Section 4 verifies the emergent Frobenius algebra structure on the boundary and elucidates the phenomenon of anyon splitting in boundary anyon condensation. Section 5 illuminates the boundary EM duality. Section 6 maps the Fourier-transformed extended QD model to the extended LW model on the same lattice. Section 7 provides two concrete examples of our results. Finally, the appendices collect a review of the extended LW model and certain details to avoid clutter in the main text.

Extended quantum double model
An extended QD model [9,10] is an extension of the QD model to the case with boundaries by adding boundary Hamiltonians to the QD Hamiltonian. The model is a Hamiltonian extension of the Dijkgraaf-Witten topological gauge theory with a finite gauge group. The model can be defined on an arbitrary lattice with one or multiple boundaries. Topological invariance allows the model to be defined on a fixed lattice for computational convenience. In this paper, we consider an oriented trivalent lattice Γ, part of which is shown in Fig. 1.
The input data of the model is a finite gauge group G, whose elements are assigned to the edges of Γ. The total Hilbert space is spanned by all possible configurations of the group elements of G on the edges of Γ and is the tensor product of the local Hilbert spaces respectively on the edges. Namely, where e is an edge in Γ. Note that reversing the orientation of an edge graced with a group element g changes the group element toḡ := g −1 ; however, we work with a fixed lattice with a fixed orientation. The Hamiltonian of the model is the sum of a bulk Hamiltonian and a boundary Hamiltonian: The bulk Hamiltonian consists of two sums of local operators: where the two sums are respectively over all vertices and all plaquettes in the bulk of Γ. A local vertex operator A QD v acts locally on the three edges incident at the vertex v as follows.
which is understood as a discrete gauge transformation averaged over G.
and it projects out any such local states that are not invariant under the gauge transformation. In gauge theory terminologies, A QD v imposes a local Gauss constraint. A local plaquette operator acts locally on the six edges outlining the plaquette p as follows.
which is also a projection. In gauge theory terminologies, B QD p imposes a local flatness condition in p. All plaquette operators and vertex operators commute. A gapped boundary condition is specified by a subgroup K ⊆ G. The boundary Hamiltonian comprises boundary local operators: where the two sums are respectively over all vertices and all virtual plaquettes (to be defined shortly) on the boundary of Γ. We define where the trace is taken over the total Hilbert space (2.1). Quasiparticle excitations of the model are charges on the vertices, fluxes through plaquettes, and dyons as bound-states of charges and fluxes. A charge at a vertex v arises when the local Gauss constraint is violated; a flux through a plaquette p occurs when the local flatness condition is violated; when both constraints are violated in p, a dyon shows up in p. Other properties of the model and topological phases that are classified by this model can be found in Ref. [10] and references therein.

Fourier transforming and rewriting the extended QD model
In this section, we Fourier-transform the basis of the Hilbert space of the extended QD model from the group space to the representation space, and as urged by this transformation, rewrite the extended QD model on a slightly different lattice.

A graphical tool for group representation theory
To facilitate the derivations in this paper, we introduce the following graphical tool for group representation theory [16]. Let L G be the set of all (representatives of equivalence classes of) irreducible representations of a finite group G. For simplicity, we shall define Here, the line is oriented from the right index n µ to the left index m µ of the representation matrix. The g insertion on the line indicates the group action by g. In this graphical presentation, multiplying two representation matrices is done by simply concatenating two such lines: A line with group action g = e reads as an identity matrix, as in Eq.  where we invoke that two bases belong to different representation spaces are orthogonal to each other. Using the above inner product, we can define the complex conjugate of an irreducible representation diagrammatically as where the arrow is reversed after the complex conjugation. This way, the complex-conjugated group action should be read upward as D µ mµnµ (g) * = D µ nµmµ (ḡ), withḡ := g −1 , because of the unitarity of µ. As such, the great orthogonality theorem is presented by (3.6) Every irreducible representation µ ∈ L G has a dual µ * ∈ L G , such that µ * is equivalent to (but not necessarily identical to) the complex conjugate representation of µ. There is an invertible duality map where the Ω µ mµn µ * is a complex matrix satisfies normalization Ω † Ω = I and maps µ to µ * by similarity transformation Graphically, the duality map and its inverse has presentations And the similarity transformation in Eq. (3.8) is presented by (3.10) For the uniqueness of duality map, the matrix Ω µ is either symmetric or antisymmetric depends on whether µ is pseudo real or not. This is an intrinsic property specified by a number β µ called the Frobenius-Schur (FS) indicator. We have (Ω µ * ) T = β µ Ω µ with β µ = ±1 if µ is real or pseudo real. Frequently in later derivations, we will need 3j-symbols to deal with the coupling of three representations of G. A 3j-symbol is a tensor w abc mam b mc that is defined as an intertwiner: where (a, V a ), (b, V b ), and (c, V c ) are three irreducible representations of G, and m j labels the basis e m j in the vector space V j . Graphically, a 3j-symbol is presented by (3.12) Being an intertwiner, by definition a 3j-symbol is invariant under group actions, namely, b a c mc m b ma (3.13) In this paper, we will also use a lot of Clebsch-Gordan coefficients, which are values of intertwiners that can be defined using the 3j-symbols (3.12) and duality map (3.9) as (3.14) For later convenience, we list a few properties of 3j-symbols as follows. For a generic finite group G, we can always construct such 3j-symbols satisfying the following properties [16].  G µνλ ηκρ = mµ,n µ * ,mν ,n ν * ,...,mρ,n ρ * Ω µ mµn µ * Ω ν mν n ν * Ω λ m λ n λ * Ω η mηn η * Ω κ mκn κ * Ω ρ mρn ρ * × w κλ * η mκn λ * mη w η * ν * ρ n η * n ν * mρ w ρ * µ * κ * n ρ * n µ * n κ * w µνλ mµmν m λ , (3.19) which can be presented either by a planar graph or a tetrahedron (3.20)

Fourier transform on the Hilbert space
Let us first Fourier-transform the Hilbert space of an extended QD model with a finite gauge group G. The total Hilbert space of the model is defined in Eq. (2.1). Since the degrees of freedom of the extended QD model live on the edges of the lattice, the total Hilbert space H is the tensor product of all local Hilbert spaces H e on the edges. The local basis state on an edge is |g , which can be Fourier transformed to the basis states in terms of representations |µ, m µ , n ν , called a rep-basis, by performing the following Fourier transform (FT): |G| is the order of G, and D µ mµnµ is the representation matrix in V µ . The local rep-basis and the group-basis have the same dimension because µ d 2 µ = |G|. Note again that the total Hilbert space of the model is defined regardless of the Hamiltonian, which when imposed would separate the Hilbert space into ground and excited states. Hence, the Fourier transformation of the group-basis of total Hilbert space on the trivalent lattice can be done by transforming the local basis of the H e on each individual edge independently. Figure 2(a) shows a small part of the lattice Γ, including both bulk and boundary, and a corresponding group-basis state on the part of Γ. Fourier-transforming the group-basis state in Figure 2(a) results in a linear superposition of the rep-basis states, shown in Figure  2(b). The explicit transformation will be dealt with a little later. Here we focus on the logic and physics of the basis transformations.
While the basis state in Figure 2(a) is directly defined on the original trivalent lattice of the model, the basis states in Figure 2(b) lack the attachment to an actual trivalent lattice. Here is the reason. The Fourier transform turns the three group elements on the three lattice edges meeting at a vertex independently into three linear combinations of all irreducible representations of G. Taking the vertex v 2 in Figure 2(a) as an example, the three group elements g 2 , g 3 and g 5 are respectively transformed into three independent representations µ 2 , µ 3 , and µ 5 in Figure 2(b); however, there is no a priori a reason the three lines graced with µ 2 , µ 3 , and µ 5 should simply just meet at v 2 . If we did so and removed the matrix indices n 2 , m 3 , and n 5 , one would misinterpret the vertex v 2 as an intertwiner of the three representations and would certainly be wrong because the dimension of the basis would be lower than that of the group-basis. Hence, we leave the vertex v 2 and all other vertices open in Fig. 2(b), and the basis becomes detached from the original lattice. Now the question is: Can we rewrite the basis in Fig. 2(b) to one that has the same dimension and still attached to an actual trivalent lattice? The answer is "Yes". Let us Rewrite all vertices again stare at the vertex v 2 in Fig. 2(b). The strategy is to fuse (i.e., couple) the three representations in an order at v 2 . Since coupling representations is an associative linear transformation, we can choose our convention. Let us first fuse µ 2 and µ 5 by contracting their indices n 2 and n 5 , resulting in a linear combination of irreducible representations {η}, each member of which is graphically a line with a free end labeled by m η . Then, we can fuse an η with µ 3 by contracting m η and m 3 , resulting in a linear combination of representations {s}, each member of which is a line with a free end labeled by m s . For each η and each s, we obtain a new basis state, as in Fig. 2(c). A pivotal point is that the degrees of freedom η and s are both local with respect to the vertex v 2 because they arise from fusing the three representations µ 2 , µ 5 , and µ 3 at v 2 . Hence, the line graced with s can be considered an induced degree of freedom at v 2 , and we place the line very close to v 2 , as a dangling edge. In this procedure of rewriting the basis, there is no loss of dimension because the degrees of freedom associated with v 2 in Fig. 2 , which agrees with that contributed by g 2 , g 3 , and g 5 in Fig. 2(a) and that contributed by µ 2 , µ 3 , and µ 5 in Fig. 2(b). Here we used the identities (A.3) and a d 2 a = |G|. We can go through the procedure above on all the other open vertices in Fig. 2(c) and obtain the basis states in Fig. 2(d), which we shall call a rep-basis state. Now the basis states are again defined on an actual trivalent latticeΓ, which differs from the original lattice Γ by having a tail attached to each of the original vertices. The tails are necessary to maintain the correct number of local degrees of freedom. On the new lattice, each vertex can indeed be interpreted as where three representations fuse. The new latticeΓ is in fact the right lattice for defining a LW model, such that the Hilbert space contains both ground states and charge excitations [8]. We will come back to this point later when we map our Fourier-transformed model to the extended LW model in Section 6.
Since the vertex and plaquette operators comprising the Hamiltonian of the extend QD model are local operators, the action of such an operator only affects a few local degrees of freedom, without affecting others. Consider a boundary vertex operator (2.7) acting on the vertex v 2 in Fig. 2(a), it affects only the three group elements meeting at v 2 . When transformed into the basis in Fig. 2(d), the vertex operator acting on v 2 would possibly affect at most the degrees of freedom µ 2 , µ 3 , µ 5 , η, s, m s , which are local at v 2 . In other words, the other degrees of freedom other than these six are diagonal indices when the vertex operator is represented in the Hilbert space. Therefore, when studying a vertex operator action at a vertex, we can simply single out a local basis state consisting of only the degrees of freedom local at the vertex.
Let us study the Fourier transform and basis rewriting depicted in Fig. 2 by focusing on a local state at a boundary vertex in detail to retrieve the linear transformations of basis with precise coefficients. Consider a local basis state |g, h, l drawn on the right hand side of Eq. (3.22). In our graphical presentation, Fourier-transforming this basis to one in the representation space is as trivial as in Eq. (3.22), and the inverse transformation in Eq. (3.23). (3.23) In the basis state on the left hand side of Eq. 3.22, for any fixed µ, ν, and λ, the degrees of freedom are only at the ends of the three lines. We then rewrite the local basis state on the left hand side of Eq. (3.22) by first fusing the two representations µ and ν via contracting their indices m µ and m ν , resulting in a set of representations {γ}, Then, we can fuse a γ with λ by contracting m γ and m λ , resulting in a set of representations {s}. This procedure yields two 3j-symbols with a pair of indices contracted, resulting in the coefficients of the expansion on the right hand side of (3.24). This procedure produces the linear combination of the new local basis states on the right hand side of Eq. (3.24) (also recall Fig. 2(d)). We can denote the new local basis states at the vertex by |Ψ sms and write down the inverse transformation in Eq. (3.25). (3.25) Each black dot in Eqs. (3.24) and (3.25) and in derivations hereafter represents a 3jsymbol assigned to each vertex to ensure that the three representations can couple to trivial representation. Explanations of the notation |Ψ sms must follow. The linear transformation in Eq. (3.24) rewrites only a subspace of the local Hilbert space spanned by the basis on the left hand side of the equation. The subspace is the one comprised by the degrees of freedom m µ , m ν , and m λ , which are transformed into γ, s, and m s . Hence, the two bases before and after the rewriting both have the same labels µ, ν, λ, n µ , n ν , and n λ . Moreover, staring at the right hand side of Eq. (3.24), one can see that the two degrees of freedom γ and s cannot be both independent. If we choose s independent, then γ is determined by µ, ν, λ, and s. Therefore, we can denote the new basis after the rewriting by |Ψ sms for simplicity, while keeping all the other labels in the graph inexplicit. This simplification causes no confusion because in actual calculations, e.g., in computing an inner product of two such local basis states, i.e., Ψ s m s |Ψ sms , the prime in Ψ implies that the hidden labels in Ψ should all be the primed version of those in Ψ. In the next subsection, we will see another advantage of this simplified notation.
Equations (3.25) and (3.23) lead to a direct linear transformation between the local basis states in the rep-basis and those in the group-basis at a boundary vertex v: (3.26) We verify in Appendix C that the local basis states |Ψ sms indeed form a well-defined local basis by showing that it is orthonormal and complete.

Fourier transform of the vertex operators
We are now ready to study how a boundary vertex operator A QD v acts on a local basis state |Ψ sms in the rep-basis. In other words, we need to find the Fourier-transformed version acts on a group-basis local state |g, h, l ∈ H v by modifying the group elements g, h, and l on the three edges incident at v via gauge transformation. When the group-basis local state is expanded in terms of the local states |Ψ sms , the action of A QD v are expected to spread over all the degrees of freedom µ, ν, λ, γ, s, and m s local to the vertex v. Nevertheless, as we will see, the indices µ, ν, λ, γ are all readily diagonal with respect to the A QD v represented in this local Hilbert space. This further renders |Ψ sms a good notation of the local basis states for A QD v . Second, recall that a bulk vertex operator A QD v (2.4) acts as a gauge transformation averaged over the entire group G; it is a projector that projects out non-intertwiner states of µ, ν, and λ but keep the intertwiner states as its +1 eigenstates. A boundary vertex operator A QD v (2.7) however performs a gauge transformation averaged over a subgroup may not project out all non-intertwiner states of µ, ν, and λ.
To see this point, let us consider two extreme cases. When K = G, A QD v acts exactly the same as a bulk vertex operator and will project out all non-intertwiner states of µ, ν, and λ, and leaves no degrees of freedom for s. That is, the s on the dangling edge (tail) attached to v must be the trivial representation 0 of G. Note that the γ, as a degree of freedom on an internal edge, is not independent; hence, if s is trivial, γ = ν. When K = {e}, the action of A QD v is trivial and thus does not project out any states in the fusion space of µ, ν, and λ. That is, the s can be any irreducible representation of G. Now for a generic nontrivial subgroup K ⊂ G, the action of A QD v may render s taking value in certain subset of the set L G of all irreducible representations of G. Moreover, certain states in V s may also be projected out. Therefore, the states that survive the action of A QD v may be labeled by pairs (s, α s ) and collected into a set L A = {(s, α s )}, which is to be defined shortly. The intuitive argument above can be explicitly verified as follows, using Eqs. (2.7) and (3.26), and the completeness of the local rep-basis. Note that for certain s ∈ L G , P s K may project out the entire H s v ; therefore, the set L A collects all the +1 eigenstates of A QD v or its representation P s K for all possible s that is not annihilated by P s K . That is, The states |Ψ sms with (s,m s ) / ∈ L A are zero eigenstates of A QD v ; they have higher energy than the +1 eigenstates according to the Hamiltonian (2.6) and thus are excited states. The excitations appear at the end of the dangling edge graced with representations s and are the point-like charge excitations at the boundary.

Fourier transform of the plaquette operators
We can now proceed to check how a boundary plaquette operator B is represented in the Hilbert space. Hence, following the same logic as that of constructing the local basis states for the boundary vertex operators, we can denote the local basis states to be acted on by a boundary plaquette operator by |Ψ ηλ rmr;sms , defined as follows. (3.31) Here,ṽ j = d j withd j = β j d j , and the coefficients G abc def are the symmetric 6j-symbols of the irreducible representations of the group G; their properties can be found in Appendix A or in Ref. [16]. The full derivation of the equation above is found in Appendix C. The local basis states | Ψ ηλ rmr;sms may be +1 or zero eigenstate of the boundary vertex operators acting on the relevant vertices. It is useful to study the matrix elements of B where in the last equality we define We have successfully Fourier transformed and rewritten the extended QD model on a trivalent latticeΓ. In the following sections, we shall study the physical consequences.

Emergence of Frobenius algebras and anyon condensation
Interestingly, the Fourier transform and rewriting of the extended QD model leads to an emergent Frobenius algebra structure on the boundary ofΓ. Namely, the set L A defined in Eq. (3.29) together with the symbols f defined in Eq. (3.33) form a Frobenius algebra, as an object in the UFC Rep G -the category of linear representations of G, which is the tensor category generated by L G -the set of irreducible representations of G.
Before we show the emergent Frobenius algebra, let us review what Frobenius algebras are. A Frobenius algebra A is a pair (L A , f ), where L A is the set of elements of A, and f is the multiplication. An element of L A is a pair (s, α s ) (or just sα s for short), where s is a simple object of a UFC F. We denote the number of different pairs (s, α s ) with the same s by |s| and call it the multiplicity of s in A. The multiplication f is a map : L A × L A × L A → C that satisfies the following associativity and non-degeneracy.
where 0 is the unit element of A and has multiplicity 1, i.e. 0 = (0, 1). Here,ṽ c = d c , which is defined in Appendix A.
That the Frobenius algebra A defined above is an object of the corresponding UFC F is understood by writing A as A = s| (s,αs)∈L A s ⊕|s| , which is in general a non-simple object in F. For the sake of computation, one may also write A = (s,αs)∈L A s αs , explicitly taking different appearances of s as distinct elements of A.
We can show that the symbols f aαabα b cαc defined in Eq. 3.33 indeed satisfy the defining conditions (4.1) of a Frobenius algebra. We first prove the associativity (4.1a) as follows. In the first equality above, we simply rewrite the left hand side using the definition (3.33). The result is then graphically presented via the second equality. The third and fourth equalities are due to Eq. (3.15), while the fifth equality due to Eq. (3.16). The last equality is again rewriting using definition (3.33). To prove the non-degeneracy condition (4.1b), note that since ñ c * (w 0bc * 0α bñc * ) * (Ω c ) −1 n c * αc ∝ δ b,c δ α b ,αc due to Schur's lemma, det( ñ c * (w 0bc * 0α bñc * ) * (Ω c ) −1 n c * αc ) = 0 . Therefore, for any finite group G and a subgroup K ⊆ G, the set L A defined in Eq. Recall that according to Eq. (3.29), each pair (s, α s ) labels a local +1 eigenstate |Ψ sαs at vertex v of A QD v , which is a projector. All such eigenstates sharing the same hidden labels span a subspace in the d s -dimensional representation space V s . The dimension |s| of this subspace is the number of pairs (s, α s ) with the same s. A salient point of our recognition of an emergent Frobenius algebra out of the Fourier transform of the extended QD model with gauge group G is that it identifies this dimension |s| as the multiplicity |s| of s appearing in A (see above Eq. (4.1a)). This identification is intricately related to the mechanism of anyon condensation in topological phases. Here we briefly describe this relation and shall report the detailed studies elsewhere.
Anyon condensation has been extensively studied recently (see Refs. [13,[17][18][19] and references therein). In a topological phase C, certain types of anyons may condense and cause a phase transition that breaks the topological phase to a simpler child topological phase U. In an extreme case, U may be merely a vacuum, and the original topological phase is said to be completely broken. An alternative and equivalent perspective is that there is a gapped domain wall separating C and U. In particular in the case where U is a vacuum, we say certain types of anyons of C can move to and condense at the gapped boundary between C and the vacuum. An interesting phenomenon is that certain types of anyons do not condense at the boundary straightforwardly; rather, such an anyon may split into a number of pieces at the boundary, and not all of these pieces can and necessarily condense. If an anyon splits into two pieces, it is said to have multiplicity one (two) in the condensate if only one (both) of the two pieces condense. Anyon splitting only occurs to anyons with quantum dimension greater than or equal to 2. So far, the understanding of anyon splitting is categorical; hence, it would be important to understand such splitting in a concrete lattice model of topological phases with gapped boundaries in terms of the input data of the model.
For the extended QD model with input gauge group G, the gapped boundary conditions are specified by the subgroups of K ⊆ G. It is known that K = {e} corresponds to condensing all the G-charges at the boundary, K = G corresponds to condensing all the G-fluxes, and any K in between corresponds to condensing certain types of dyons.
Nevertheless, by the subgroup K alone one cannot immediately tell whether a type of condensed anyons should have a multiplicity greater than one in the condensate. In the Fourier-transformed picture, as we now show, anyon splitting and multiplicity becomes lucid. According to our discussion earlier in this subsection, a subgroup K on the boundary of Γ gives rise to a Frobenius algebra A G,K = (L A , f ) G,K at the boundary ofΓ. Consider the case where K = {e}, i.e., charge condensation at the boundary. Then by Eqs. where V αs s denotes the direct sum V s that appears in A G,K with the multiplicity index α s . The multiplication is also the 3j-symbol (up to some normalization factor) as defined in Eq. (3.33). By the definition of L A in Eq. (3.28), α s is given by a basis vector such that ρ s (k)|s, α s = |s, α s for all k ∈ K. The space A G,K spanned by |s, α s is thus identified with (C[G]/C[K]) * defined above.

EM duality in the bulk
The QD model exhibits an EM duality. Consider the QD models with finite Abelian groups G. In such a case, all the irreducible representations of G are 1-dimensional and form a group, with the group multiplication defined as the tensor product of representations. For example, the irreducible representations of G = Z n form the group Rep G = Z n . Let g ∈ {0, 1, . . . , n − 1} and g · h = g + h mod n. Then the irreducible representations are j ∈ {0, 1, . . . , n − 1} and j · k = j + k mod n in the sense that ρ j (g)ρ k (g) = ρ j·k (g) for all g ∈ G.  As reviewed in Section 2, the QD model on a closed trivalent graph Γ consists of vertex and plaquette operators. For example, in Fig. 3(a), we have an operator at vertex v = 1: where g 1 , g 2 , g 3 ∈ G are the group elements on the three edges meeting at v. At plaquette p = 2, we have B p |g 1 , g 2 , g 3 , g 4 , g 5 , g 6 = δ g 1 ·g 2 ·g 3 ·g 4 ·g 5 ·g 6 ,0 |g 1 , g 2 , g 3 , g 4 , g 5 , g 6 , where the g s are the group elements on the edges outlining the plaquette p.
We also Fourier transform the basis of Hilbert space H Γ on Γ from the group space to the representation space. Hence, the vertex and plaquette operators act on the new basis states as where j 1 , j 2 , j 3 ∈ RepG are on the three edges meeting at the vertex v = 1 in Γ, and where the j 's are graced on the six edges outlining the plaquette p = 2 on Γ.
On the other hand, we can draw the dual lattice Γ * (Fig. 3(b)) of Γ and place an element of Γ * on each edge of Γ * . Since Rep G ∼ = G, H Γ * ∼ = H Γ . The dual QD model with Rep G is naturally defined on Γ * . The basis states |j 1 , j 2 , j 3 v=1 and |j 1 , j 2 , j 3 , j 4 , j 5 , j 6 p=2 on Γ are the same as the basis states |j 1 , j 2 , j 3 p * =1 and |j 1 , j 2 , j 3 , j 4 , j 5 , j 6 v * =2 on Γ * . Hence, we have where B * p * and A * v * are the plaquette and vertex operators of the dual QD model. Comparing Eqs. (5.3) and (5.4) to (5.5) and (5.6), and since the above discussion applies to any vertex and plaquette on Γ, we conclude that the QD model with a group G on Γ is mapped to the dual QD model with a group Rep G on Γ * , with H Γ ∼ = H Γ * , and Since A v measures the gauge charges (which have an electric nature) and B p measures the gauge fluxes (which have a magnetic nature), the above duality is an electric-magnetic duality: an electric charge/magnetic flux on Γ is mapped to a magnetic flux/electric charge on Γ * .
For non-Abelian cases, Rep G is no longer a group. In a generalized context using quantum groups (Hopf algebra), we can still define the charges and fluxes of Rep G , and thus are able to study the EM duality [6] in the QD model on a closed surface. In Ref. [8], the authors has also checked the EM duality in the enlarged Hilbert space of the LW model on a closed surface.

EM duality on the boundary
The EM duality in the bulk can be extended to one that is on the boundary. We again consider the Abelian cases first. As reviewed in Section 2, on the boundary, given a subgroup K, there are two types of operators,Ā v acting on the boundary vertices, andB p on the boundary edges, as illustrated in Fig. 4(a).
TheĀ v measures the K-charges, whileB p measures the G/K-fluxes.
To reveal the EM duality on the boundary, we need to extend the lattice Γ in Fig  4(a) to the latticeΓ in 4(b), following the procedure illustrated in Fig. 2. As explained in Section 3.2, this extension preserves the Hilbert space of the extended QD model, i.e., H Γ = HΓ, and is merely a basis transformation. That is, the basis of the Hilbert space is transformed from the group space to the representation space.
Then according to Eqs. (2.7) and (2.8) and the basis transformation, a boundary vertex operatorĀ v and plaquette operatorB p acts respectively on the highlighted dangling edge and open plaquette onΓ in Fig. 4(b) as where j ∈ Rep G is graced on the red dangling edge in Fig. 4(b), and where j 1 , j 2 , j 3 , j 4 ∈ Rep G are respectively placed on the four blue edges outlining the open plaquette in Fig. 4(b), and L A is a subgroup of Rep G defined by We then draw in Fig. 4(c) the dual latticeΓ * of the extended latticeΓ. This dual latticeΓ * is a triangular lattice with a boundary and naturally defines a dual extended QD model with Rep G being the gauge group. The dual boundary operators B * p * and A * v * act on the dual boundary edges and vertices as (5.14) Comparing Eqs. (5.10) and (5.11) to (5.13) and (5.14), we conclude that the extended QD model with a group G on Γ is mapped to the dual extended QD model with the group Rep G onΓ * , with H Γ ∼ = HΓ * , and TheĀ v measuring K-charges is mapped to B * p * measuring (Rep G /L A )-fluxes, andB p measuring G/K-fluxes to A * v * measuring L A -charges. Although we consider the Abelian examples here, the duality (5.15) can be examined for an arbitrary finite group, in the language of Frobenius algebras. After the Fourier transformation, Eqs. (5.10) and (5.11) become Eqs. (3.28) and (3.32), expressed in terms of the Frobenius algebra A. They can be viewed as the operators that measure the (Rep G /L A )-fluxes and L A -charges in a more generalized context.

Mapping the Fourier-transformed extended QD model to the extended LW model
The Fourier transformed extended QD model with a finite gauge group G can actually be regarded as an extended LW model with a input unitary fusion category whose simple objects are the irreducible representations of G, commonly denoted by Rep G too. According  . Note that as defined in Eq. (3.28), here we choose the indicesm s i andm q i that diagonalize the vertex operators A v and A v . The original LW model [2] is known to have a Hilbert space smaller than that of the QD model [7] because the Hilbert space of the original LW model does not contain excited states with charges or dyons [8]. In Ref. [8], the authors generalized the original LW model by adding a tail to each vertex in the trivalent lattice of the model. The generalized LW model has a larger Hilbert space encompassing a full dyon spectrum due to the tails, which play the same role as the dangling edges in Fig. 5. Therefore, to map the Fourier-transformed extended QD model to the extended LW model with an enlarged Hilbert space, there is no need of reducing the Hilbert space of the extended QD model, in contrary to what was done in Ref. [7].
What acts on the basis states in Fig. 5 is the Fourier transformed version of the Hamiltonian of the extended QD model. Namely, where use is made of P s G = 1 |G| g D s msm s (g) = δ s,0 , with 0 being the trivial representation of G. Clearly, when s = 0, m s is not a degree of freedom and can be omitted. As to a bulk B QD p operator, its action on a bulk plaquette would be straightforwardly obtained and would be similar to that of the boundary B QD p but with seven more 6j-symbols because a bulk plaquette has seven more vertices on its perimeter. Besides, there is no restriction to L A in a bulk B QD p operator. By comparing to the boundary vertex and plaquette operators in the extended LW model reviewed in Appendix B or in Ref. [11,12] and the bulk operators in the enlarged LW model in Ref. [8], we can actually identify the model H QD,Γ L G ,A G,K defined onΓ as the extended LW model H LW,Γ Rep G ,A with A = A G,K defined on the same latticeΓ. That is, The two models have the same Hilbert space and Hamiltonians term by term.

Examples
In this section we offer two explicit examples, one for G being Abelian and one for G = S 3 , to aid the understanding of the results.

Abelian groups
Since an Abelian group is always isomorphic to the product of cyclic groups, let us consider Since the irreducible representations are 1-dimensional, the representation matrices are merely complex numbers; hence, all the matrix indices can be removed. The Fourier transform of the local group-basis states on an edge takes the simple form: where The rep-basis local states (3.25) of the vertex operators are now denoted simply by |Ψ s . The fusion of three irreducible representations µ, ν, and λ is determined by the delta function δ µνλ just introduced. Note that in Eq. (7.2), the group elements and irreducible representations are on equal footing, the fusion of irreducible representation is the same as the multiplication of group elements. Indeed, for an Abelian group G, Rep G has a group structure and is isomorphic to G itself. This fact results in the self duality under Fourier transforming the extended QD model with an Abelian gauge group G. This self duality is actually the EM-duality in Eq. (5.15).
To elucidate the self duality in this example, let us study how an A This result can be understood in the following way. Since G ∼ = Rep G in this case, the constraint δ t * s * s = 1 indicates the equality s = s·t * , where the dot is the group multiplication.
Then, a B QD p acts on its local basis states as  Table 2  For notation simplicity, we rename the three irreducible representations of S 3 by 0, 1, and 2, respectively corresponding to the 1, sign, and π in the table above. One set of irreducible representations matrices is listed as follows.
D 0 (e) = D 0 (r) = D 0 (r 2 ) = D 0 (sr) = D 0 (s) = D 0 (sr 2 ) = 1, The nonvanishing 3j-symbols of S 3 are w 000 111 = w 011 111 = 1, Plugging the 3j-symbols (7.8) into Eq. (3.33), we can obtain four emergent Frobenius f abc = f cab , f aa * 0 = 1 algebras at the boundary of the Fourier-transformed extended QD model, corresponding to the four distinct subgroups of S 3 , as recorded in Table 3. The four emergent Frobenius algebras can indeed be identified with the four Frobenius algebra objects in Table 4 obtained by solving the defining conditions (4.1) of a Frobenius algebra object in a UFC F, which in this case is Rep S 3 .
Here we can try to understand explicitly the relation between the emergent Frobenius algebras at the boundary and anyon condensation, in particular anyon splitting and multiplicity. Table 5. Correspondence between (column 1) anyon condensates of the D(S 3 ) model, emergent Frobenius algebras (column 2), and subgroups of S 3 (column 3).

Condensation type Boundary condensate Boundary
In the current example, there are four different gapped boundary conditions respectively characterized by the four subgroups 1 of S 3 , imposed by the boundary plaquette operators B QD p of the model. These four gapped boundary conditions respectively correspond to four different boundary condensates composed of different types of anyons condensed at the boundary. These correspondences are shown in Table 5. One can see in the table that the type C anyons, which are pure charges ([e], 2), appear with multiplicity 1 in the condensate A 3 but multiplicity 2 in the condensate A 4 .
Condensate A 3 corresponds to the emergent Frobenius algebra A 3 = (0, 1)+(2, 1). The element s = 2 in the pair (2, 1) in A 3 actually is the 2-dimensional irreducible representation of S 3 and thus can be directly identified with the pure charge C in the condensate A 3 . Note that for each K ⊆ S 3 in Table 3, the pairs (s, α s ) in the corresponding row each labels a +1 local eigenstate |Ψ sαs of the A QD v operator acting on the relevant vertex. In the case with K = Z 2 and s = 2, seen in Eq. (3.27), an operator A QD v is represented in the 2-dimensional local space V s by the projector where we didn't differentiate a matrix and its matrix elements labeled by m s m s to emphasize that the operator is not diagonalized yet. We can diagonalize the projector in the basis states em s labeled bym s and find the +1 eigenstate e 1 = ( √ 3, 1) T and 0 eigenstate e 2 = (− 1 √ 3 , 1) T . Hence, only α s=2 :=m s = 1 is allowed. That is, the pair (2, 1) ∈ L A labels the only +1 local eigenstate |Ψ 21 of the A QD v . Hence, only half of the 2-dimensional space V s or half of the charge-2 condenses at the boundary. As such, we can say that the charge-2 splits: 2 = (2, 1) ⊕ (2, 2) but only the half (2, 1) condenses. That is why only (2, 1) appears in the A 3 in Table 3, and type C anyons appear with unit multiplicity in A 3 in Table 5. If we denote (2, 1) by 2 1 and (2, 2) by 2 2 , we found that the Frobenius algebra object A 3 in Table 4 should be rigorously written as A 3 = 0 ⊕ 2 2 ; however, one often just write A 3 = 0 ⊕ 2 customarily. in V s are +1 eigenstates, which can be labeled by (2, 1) and (2,2). Note that these two eigenstates are not those of the operator in Eq. (7.9). As such, the charge 2 splits as 2 = (2, 1) ⊕ (2, 2), but both pieces condense at the boundary. Therefore, in the boundary condensate, the type C anyons count twice, rendering A 4 = A ⊕ B ⊕ 2C. If we write (2, 1) and (2, 2) as 2 1 and 2 2 , we can identify the emergent Frobenius algebra A 4 with the Frobenius algebra object A 4 = 0 ⊕ 1 ⊕ 2 1 ⊕ 2 2 in Table 4.
As a further corroboration of the mapping between the extended QD model and the extended LW model, we can compute the GSDs of both models on a cylinder, whose two boundaries may not necessarily possess the same boundary conditions. Using Eqs. (2.9) and (B.10), we obtain respectively Table 6 and Table 7. Matching the GSDs in these two tables futher confirms the correspondence between the gapped boundary conditions specified by the subgroups of S 3 in the extended QD model and those specified by the Morita inequivalent Frobenius algebra objects in the extended LW model.
For completeness of the data in this example, we collect in Table 8 the Modules of the Frobenius algebras in Rep S 3 . Table 6. The GSDs of the extended QD model with G = S 3 on a cylinder, whose two boundaries are specified by two subgroups of S 3 .

A Some properties of fusion categories
Simple objects in a UFC F obey a set of fusion rules: where N k ij are the fusion coefficients satisfying Here, j * is the dual of j defined by the second row in the equation above. Each object j also has a nonzero characteristic numberd j , called quantum dimension of j, such that In particular,d 0 = 1. Let β j = sgn(d j ), then β i β j β k = 1 is satisfied if N k ij > 0. When F = Rep G for a finite group G, we haved j = β j d j , and β j is identified with the FS indicator of the representation j of G. Now we can define tetrahedral-symmetric unitary 6j-symbols as the coefficients of the map G : L 6 → C that satisfies the following conditions G ijm kln = G mij nk * l * = G klm * ijn * = β m β n (G j * i * m * l * k * n ) * , nd n G mlq kp * n G jip mns * G js * n lkr * = G jip q * kr * G riq * mls * , nd n G mlq kp * n G l * m * i * pk * n = δ iq d i δ mlq δ k * ip .

B Extended LW model
In this section, we review the basic ingredients of extended LW model. The extended LW model is defined on an oriented trivalent latticeΓ , part of which is depicted in Fig. 6. The Hilbert space is spanned by all configurations of assigning of the simple objects of a UFC F on the oriented edges ofΓ . Conventionally, we call these simple objects string types and label them by integers in L = {0, 1, ..., N }. Each string types j has a dual j * , which is also an element of L and satisfies j * * = j. If we reverse the orientation of the edge and conjugate the string type j → j * associated with the edge at the same time, the state remains the same. There is always a trivial (unit) element 0 ∈ L satisfying 0 * = 0.
The set {N k ij ,d j , G ijk lmn } is the input data of the extended LW model, which can be derived from representation theory of a finite group or quantum group. In this work, we take the string types µ to be the irreducible representation of a finite group G and label them with Greek letters. The trivial representation is µ = 0. The set L is identified with L G -the set of irreducible representations of G. The fusion coefficients N λ µν are the numbers of trivial representations appearing in the decomposition of µ ⊗ ν ⊗ λ. In such cases, as mentioned earlier,d µ = β µ d µ with β µ the FS indicator. The dual representation µ * of µ satisfiesd µ * =d µ . The 6j-symbols G µνλ κηγ are the same as the symmetrized Racah 6j-symbols of G. Similar to the extended QD models, there are also two types of local operators in the bulk Hamiltonian of the extended LW model, namely where A LW v and B LW p are respectively the vertex and plaquette operators. The action of A LW v on a local trivalent state as where δ ijk is delta function that is equal to 1 if N k ij > 0 and otherwise 0. The action of B LW p is more involved and reads Note that on a cyliner, ∂Γ contains two components, on which the operators A LW n may project the boundary degrees of freedom into respectively two (not necessarily) different Frobenius algebras.