Arxiv:1707.05490V2 [Quant-Ph] 1 Nov 2017

Universal Quantum Computation with Gapped Boundaries

Iris Cong,1, 2, 3 Meng Cheng,4, 3 and Zhenghan Wang3, 5 1Department of Computer Science, University of California, Los Angeles, CA 90095, U.S.A. 2Department of Physics, Harvard University, Cambridge, MA 02138, U.S.A. 3Microsoft Station Q, University of California, Santa Barbara, CA 93106-6105 U.S.A. 4Department of Physics, Yale University, New Haven, CT 06520-8120, U.S.A. 5Department of Mathematics, University of California, Santa Barbara, CA 93106-6105 U.S.A. This Letter discusses topological quantum computation with gapped boundaries of two- dimensional topological phases. Systematic methods are presented to encode quantum information topologically using gapped boundaries, and to perform topologically protected operations on this encoding. In particular, we introduce a new and general computational primitive of topological charge measurement and present a symmetry-protected implementation of this primitive. Throughout the Letter, a concrete physical example, the Z3 toric code (D(Z3)), is discussed. For this example, we have a qutrit encoding and an abstract universal gate set. Physically, gapped boundaries of D(Z3) can be realized in bilayer fractional quantum Hall 1/3 systems. If a practical implementation is found for the required topological charge measurement, these boundaries will give rise to a direct physical realization of a universal quantum computer based on a purely abelian topological phase.

Introduction. The quantum model of computation study and experimental realization of TQC. We intro- strikes a delicate balance between classical digital and duce a new computational primitive—topological charge analog computing models, as its stability lies closer to measurement (TCM), which extends topological charge digital models, while its computational power is closer to projection [4]. We propose a physical realization of analog ones. Still, a major obstacle to developing quan- symmetry-protected TCM in a fractional quantum spin tum computers lies in the susceptibility of qubits to de- Hall state, while leaving a fully topologically protected coherence. One elegant theoretical solution to this prob- one to the future because which measurement is possible lem is topological quantum computation (TQC) [14–16]. in gauge theory is an open fundamental question [5].

TQC is a paradigm that information is encoded in topo- Our universal gate set for D(Z3) is close to experi- logical degrees of freedom of certain quantum systems, mental technology in bilayer quantum Hall liquids. If a thereby protected from local decoherence. While the practical implementation is found for our TCM primitive, standard implementation uses (non-abelian) anyons in this gate set is a direct physical realization of a universal topological phases of matter, recent studies revealed that quantum computer. certain topological phases also support gapped boundaries. It is hence natural to study TQC with gapped Realization of Z3 toric code by bilayer ν = 1/3 boundaries [4,7,8, 13]. fractional quantum Hall liquids. The Z3 toric code D(Z3) can be realized in bilayer fractional quantum Hall Real samples of topological phases of matter such as systems: Ref. [2] considers an electron-hole bilayer FQH fractional quantum Hall liquids and topological insula- system, with a 1/3 Laughlin state of opposite chirality tors have boundaries, which are usually conducting (gap- in each layer. The topological order in this system is less) even though the bulk are insulating (gapped). How- 1 SU(3)1 SU(3)1 , which is equivalent to the Z3 toric ever, they can be modified to realize Dijkgraaf-Witten × code D(Z3). Hence, we will recycle many of the results (DW) gauge theories, which are also given by Kitaev’s of Ref. [2]. quantum double Hamiltonian [16]. These theories sup- We briefly summarize the basic data for D( ). Math- port gapped boundaries in the sense that the extensions Z3 arXiv:1707.05490v2 [quant-ph] 1 Nov 2017 ematically, a topological phase is described by a modu- of the Hamiltonians to spaces (surfaces) with boundaries lar tensor category (MTC) B [1]. The anyon types2 are are still gapped; the Hamiltonian and algebraic frame- eamb, a, b = 0, 1, 2, where e and m are unit gauge works are developed in Refs. [8,9]. These frameworks Z3 charge and flux respectively (so e2 = e, m2 = m). The show that a gapped boundary effectively behaves as a braiding statistics of the anyons is encoded in the mod- non-abelian anyon. However, while the existence of non- abelian anyons is still uncertain, gapped boundaries of abelian phases are much more routine and support topologically protected degeneracies even on the plane. 1 In this Letter, we apply our theory to a concrete phys- Together with physical electrons, SU(3)1 is topologically equiv- ical example—the toric code D( )—to obtain a uni- alent to a 1/3 Laughlin state Z3 Z3 2 There are many terms in the literature referring to the same versal gate set, which is a striking example of the ex- thing: a simple quasiparticle, an anyon, or a simple object of tra computational power from gapped boundaries. This D(G). An anyon type, topological charge, or superselection sec- new direction opens up new vistas in both the theoretical tor is an isomorphism class of the above. 2

FIG. 1. Ground state for n gapped boundaries Ai on a plane and total charge vacuum. All edges are directed to point downward. 2 FIG. 2. Braiding of two gapped boundaries (σ2 ). Solid lines ular S = [S ] and T = [T ] matrices [1]: indicate tunneling operators from the basis vectors (i.e. not ab ab motion of the holes), while dotted lines indicate how the holes a2b1 a1b2 a1a2 move in the braiding process. Sab = ω− − ,Tab = ω δab. (1) Here, ω = e2πi/3. another (A ), where a (resp. a) condenses on A (resp. Gapped boundaries, degeneracy, and topological 2 1 A ). Physically, this corresponds to applying the a string operations. Let us first review the physics of gapped 2 operator [16] along a path γ connecting the two gapped boundaries, as they will encode our topological qudits. boundaries. This operation, known as a Wilson line op- A convenient physical description for a gapped bound- erator, is denoted by W (γ). For the D( ) theory, it ary type is the consistent collection of (bosonic) anyons a Z3 can be represented as follows: that can condense to vacuum to the boundary at no energy cost. Mathematically, this is given by a Lagrangian W (γ) b = a b . (2) algebra A in the MTC B3, which can be represented as a | i | × i a direct sum of all condensed anyon types. For D( 3), Z Expressing Wa(γ) as a matrix that acts on the ground there are two gapped boundary types: e-boundary (resp. state subspace, we see that W (γ) implements the single- 2 2 e m-boundary) where e, e (resp. m, m ) condense. The qutrit Pauli-X gate σx. corresponding Lagrangian algebras are 1 + e + e2 and 3 Loop-a operations. Analogously, one can create a pair 1 + m + m2. In the bilayer 1/3 Laughlin state descrip- of anyons a, a in the bulk, loop one of them around a tion, the e/m-boundary types correspond to holes with gapped boundary, and annihilate the pair. When we loop the two layers connected via electron pairing (i.e. super- a counter-clockwise around the boundary, this is known conducting) or tunneling. as the Wilson loop operator W (α ) where α is the loop Multiple gapped boundaries support a degenerate a i i encircling boundary A . Appendix B shows that ground state manifold 4. Consider a closed system with i n gapped boundaries (Fig.1). Refs. [8,9] show that the Sab ground state of the system is given by the different ways Wa(α2) b = b . (3) | i db | i we can create n anyons out of vacuum, and condense all of them onto the boundaries as a fusion tree (Fig.1). Braiding gapped boundaries. Another topological oper- This fusion tree also specifies a choice of basis states for ation is to braid gapped boundaries around each other. the ground state manifold. For example, if we have two This gives multiple-qudit operations that can produce en- e-boundaries in a planar D(Z3) theory, the GSD is 3, la- tangling gates. Physically, braiding corresponds to mov- beled by a1 =c ¯, a2 = c, c = 1, e, e¯. We denote the basis ing gapped boundaries around each other, e.g. by tuning elements by c and encode our qutrit in this space. the Hamiltonian H of Refs. [8,9] adiabatically. | i G.B. We now discuss the topological operations on gapped We may arbitrarily braid n gapped boundaries with to- boundaries, which induce unitary transformations in the tal charge vacuum around each other to obtain a unitary degenerate subspace. We focus on the D(Z3) example transformation on the ground state, so long as we return and leave the general results to the appendix. each boundary to its original position. Mathematically, Tunnel-a operations. The first topological operation is this means that the braiding matrices form a representa- to tunnel an anyon a from one gapped boundary (A1) to tion of the (spherical) n-strand pure braid group Pn [12]. They can be computed using the diagrammatic rules of anyon models and the basis states of gapped boundaries. For most purposes of quantum computation, it is suffi- 3 See Refs. [6,8,9, 17, 18] and references therein for precise defi- cient to consider 2-qudit encodings, where n = 4 (Fig. nitions. 4 The degeneracy is exponentially protected in all length scales, 2). In general, one must compute all 6 generators of P4. including distance between boundaries as well as lengths of the As an example, we derive the formula for the generator boundaries. 2 σ2 in Appendix B. 3

FIG. 3. Topological charge projection (n = 2).

Topological charge measurement. For a DW theory, the gapped boundary braidings only generate a fi- FIG. 4. Braid for the σz gate. ∧ 3 nite group [12]. Inspired by the results of Ref. [4], we introduce topological charge measurement based on the Wilson operators. Before we discuss the general case, erators to the effective Hamiltonian of the ground state recall that topological charge projection can detect the subspace: total charge of a collection of quasiparticles inside a cer- H0 = tW (β) + h.c. (5) tain region, by e.g. sending a probe particle along a − a path enclosing the region and performing interferomet- Here, t is the (complex) tunneling amplitude. This effec- ric measurement. As a generalization, we can use similar tive Hamiltonian then projects the system to the desired methods to perform measurement of topological charge state space. through any loop, not just contractible ones, possibly on a higher-genus surface [4]. Universal gate set with D(Z3) gapped boundaries. Let us now specialize to D(Z3), or the bilayer ν = 1/3 Recall that D(Z3) splits into two theories B = C C FQH. Ref. [2] proposed to use superconducting (1+e+e) with C = SU(3)1 which do not interact in the bulk, but are “stuck together” at the original boundaries of B. The boundaries to encode qutrits, so the read out can be done with electric charge measurement. We follow this scheme, planar region Y also splits into two mirror layers, S+(Y ) and S (Y ), which are completely disjoint in the bulk but and occasionally use the other (m-boundary) encoding as “stuck− together” at the boundaries of Y . This way, we an ancilla. can view the system as a single layer of C on a higher- By Ref. [10], one universal qutrit gate set is the meta- genus surface. Similarly, each loop α in Y becomes a loop plectic gate set: 1. The single-qutrit Hadamard gate H3. lα in S+(Y ) or S (Y ), while an arc γ connecting two − 2. The two-qutrit entangling gate SUM3. boundaries lifts to a loop lγ going around both layers. Let β be one of these loops. Fig.3 illustrates this for 3. The single-qutrit generalized phase gate Q3 = n = 2. diag(1, 1, ω). 4. Any nontrivial single-qutrit classical (i.e. Clifford) Define Ox(β) = Wx(αi) (tunneling operator in C) if β gate not equal to H2. is the lifting of the line αi, and Ox(β) = Wxx(γi) (loop 3 5. A projection M of a state in the qutrit space operator in B) if β is the lifting of the loop γi. By Ref. [4], C3 to Span 0 and its orthogonal complement the projection measuring topological charge a through β {| i} can be expressed as Span 1 , 2 , so that the resulting state is coherent if{| projectedi | i} into Span 1 , 2 . (a) {| i | i} Pβ = S0aSxa∗ Ox(β). (4) We now discuss how each of these gates can be imple- x C mented from the aforementioned topological operations. X∈ First, we discuss the implementation of 1-4: The sum runs over the anyon labels x of C, and Sab is 1. H3 is equal to the modular S matrix of the anyon the modular S-matrix of C. The Wilson operators Wx(αi) theory SU(3)1, so it is in the representation of map- and Wxx(γi) are computed using the formulas (2) and (3) ping class group of the torus and can be imple- C B with the data of and , respectively. As shown in [4], mented via a sequence of topological charge pro- topological charge projections generate all mapping class jections. group representations VC(Y ) of a closed surface Y from 2. For SUM , consider braiding one hole of a e- C 3 the anyon theory . boundary target qutrit with another hole of a m- For our purpose, we generalize these projections to 2 boundary control qutrit (i.e. apply σ2, as shown in topological charge measurements (TCM) which are the Fig.4). This gives complements of topological charge projections (the more general definition is in Appendix C). Specifically, given 2 2 2 z σ2 = diag(1, 1, 1, 1, ω, ω , 1, ω , ω) = σ3 . (6) an anyon label a and the lifting β of a Wilson line/loop ∧ (a) z as above, we consider the projection 1 Pβ . Physically, Because we implemented the Hadamard and σ3 = this can be implemented by adding such− non-local op- (I H )SUM (I H ), conjugating the∧ target 3 ⊗ 3 3 3 ⊗ 3 4

follows: we tune the system such that the quasiparticle tunneling along the desired loop is enhanced, so that the system has the projected charge state as the ground state. This can be achieved by e.g. using gate configura- (a) tions to diminish the energy gap. We consider 1 Pγ − FIG. 5. Short circuit (generalizing Ref. [13]) to use three as a concrete example. The desired term in the Hamil- tonian we would like to create is H0 = tW (e) + h.c. SUM gates between 1+m+m and 1+e+e qutrits to implement − γ a SUM gate between 1 + e + e qutrits. All entangling gates where t is the (complex) tunneling amplitude and Wγ (e) drawn are SUM3. is the Wilson tunneling operator. Wγ has eigenvalues 1, ω, ω¯. The coherent projection requires that the eigenvalues of H0 split into two sets, one of which has two qutrit by Hadamards gives the SUM gate between degenerate eigenvalues. This puts a stringent constraint a 1 + e + e qutrit and a 1 + m + m qutrit. We then on the complex phase of t. The simplest choice is that t have a short circuit (Fig.5) using these SUM gates is real. to implement a SUM gate between two 1 + e + e The requirement that t is real is beyond topological qutrits. After the circuit, one must interpret the protection. Physically, such condition can be met in a measurement outcome of the ancilla qutrit. If we fractional quantum spin Hall state [20, 21], an interacting measure mj , we must apply (σx)j to the control- | i 3 analog of quantum spin Hall insulator enriched by time- out (e.g. by applying W j (γ)). e reversal symmetry. Topologically, this phase is identical 3. By Ref. [4], topological charge projections can be to bilayer ν = 1/3 Laughlin state, if the layer index is used to implement diag(1, ω, ω), the Dehn twist of actually identified as the electron spin up and down. In the SU(3) theory. We follow this by a generalized 1 such a state, the time-reversal symmetry exchanges the Pauli-Z gate to obtain Q . 3 two layers. The e anyon in this physical realization is the 4. By Eq. (2), the tunneling operator W (γ) imple- e bound state of the spin up/down quasiholes. Therefore, ments the single-qutrit Pauli-X gate σx. 3 the tunneling amplitude of e has to be real since e is time- The implementation of the coherent projection M is reversal invariant, and the TCM is symmetry-protected. the most challenging part of the proposal. First, we re- Conclusions. Gapped boundaries provide the missing late M to a TCM. A planar D( ) with two 1 + e + e Z3 π -gate for a universal gate set from the doubled Ising the- boundaries can be viewed as double layers of SU(3) con- 8 1 ory [4]. In this Letter, we use our symmetry-protected nected via two handles; the curve γ connecting the two TCM to obtain a coherent projection, which augments boundaries lifts to a loop in this perspective. By Eq. (4), the topological operations from Ref. [2] for the toric projecting to vacuum within this loop gives Z3 code to a universal gate set for a qutrit computational 1 1 1 model. The Z3 toric code is realized by bilayer fractional 1 P (1) = 1 1 1 . (7) quantum Hall liquids [2], whereas it is not yet clear how γ 3 1 1 1 to physically realize the doubled Ising theory. The chal-   lenge for a realistic implementation of our universal gate (1) set now lies in a practical realization of the coherent pro- The eigenvalues and eigenspaces of P are: γ jection. 1 1 1 Acknowledgment. The authors thank Maissam λ = 0 : Span ω , ω λ = 1 : Span 1 (8) Barkeshli, Shawn Cui, and CésarGalindo for answering       many questions. I.C. would like to thank Michael Freed-  ω ω   1        man and Microsoft Station Q for hospitality in hosting One then obtains the coherent projectionM by con- the summer internship and visits during which this work (1) jugating the orthogonal projector 1 Pγ with the was done. Z.W. is partially supported by NSF grant (1) − (1) DMS-1411212. Hadamard, i.e. H3†(1 Pγ )H3. While Pγ is a topolog- − (1) ical charge projection as in Ref. [4], 1 Pγ is a general TCM. − We now have universal quantum computation using gapped boundaries of D(Z3). This is very significant, as we achieve universal quantum computation using only an [1] B. Bakalov, A. A. Kirillov. Lectures on tensor categories abelian TQFT (all anyon braidings in D(Z3) are projec- and modular functors. Vol. 21. American Mathematical tively trivial), without using state injection, as in Ref. Soc. (2001). [13]. [2] M. Barkeshli. Charge 2e/3 superconductivity and topological degeneracies without localized zero modes in bi- Symmetry-protected realization. In physical realizations layer fractional quantum Hall states. Phys. Rev. Lett. such as bilayer FQH, the TCM can be implemented as 117, 096803 (2016). 5

[3] M. Barkeshli, P. Bonderson, M. Cheng, Z. Wang. Sym- Group Representations from Braiding Gapped Bound- metry, defects, and gauging of topological phases. arXiv aries of Dijkgraaf-Witten Theories. arXiv preprint preprint arXiv:1410.4540 (2014). arXiv:1707.03884 (2017). [4] M. Barkeshli, M. Freedman. Modular transformations [13] A.G. Fowler, M. Mariantoni, J. M. Martinis, A. N. Cle- through sequences of topological charge projections. arXiv land. Surface codes: Towards practical large-scale quan- preprint arXiv:1602.01093 (2016). tum computation. Phys. Rev. A 86(3), 032324 (2012). [5] D. Beckman, D. Gottesman, A. Kitaev, J. Preskill. Mea- [14] M.H. Freedman. P/NP, and the quantum field computer. surability of Wilson loop operators. Phys. Rev. D 65, Proceedings of the National Academy of Sciences 95(1), 065022 (2002). 98-101 (1998). [6] S. Beigi, P. W. Shor, D. Whalen. The quantum dou- [15] M. Freedman, A. Kitaev, M. Larsen, Z. Wang. Topo- ble model with boundary: condensations and symmetries. logical quantum computation. Bulletin of the American Communications in Mathematical Physics 306(3), 663- Mathematical Society, 40(1), 31-38 (2003). 694 (2011). [16] A. Y. Kitaev. Fault-tolerant quantum computation by [7] S. B. Bravyi, A. Y. Kitaev. Quantum codes on a lattice anyons. Ann. Phys. 303(2) (2003). with boundary. arXiv:quant-ph/9811052 (1998). [17] A. Kitaev, L. Kong, Models for Gapped Boundaries and [8] I. Cong, M. Cheng, Z. Wang. Topological quantum Domain Walls. Commun. Math. Phys. 313, 351-373 computation with gapped boundaries. arXiv preprint (2012). doi: 10.1007/s00220-012-1500-5. arXiv:1609.02037 (2016). [18] T. Lan, J. C. Wang, X.-G. Wen. Gapped Domain Walls, [9] I. Cong, M. Cheng, Z. Wang. Hamiltonian and Al- Gapped Boundaries, and Topological Degeneracy. Phys. gebraic Theories of Gapped Boundaries in Topologi- Rev. Lett. 114(7), 076402 (2015). cal Phases of Matter. Commun. Math. Phys. (2017). [19] M. A. Levin, X.-G. Wen. String-net condensation: A doi:10.1007/s00220-017-2960-4 physical mechanism for topological phases. Phys. Rev. B [10] S. X. Cui, Z. Wang. Universal quantum computation with 71, 045110 (2005). metaplectic anyons. J. Math. Phys. 56, 032202 (2015). [20] M. A. Levin, A. Stern. Fractional Topological Insulators. [11] S. X. Cui, S.-M. Hong, Z. Wang. Universal quantum com- Phys. Rev. Lett. 103, 196803 (2009) putation with weakly integral anyons. Quantum Inf. Pro- [21] M. A. Levin, A. Stern. Classification and analysis of cess. 14:26872727 (2015). two dimensional abelian fractional topological insulators. [12] N. Escobar-Velásquez, C. Galindo, Z. Wang. Braid Phys. Rev. B 86, 115131 (2012) 1

Appendix A: Notations Appendix B: Deriving formulas for the topologically protected operations In this Appendix, we list all of the notations that are used throughout the paper. Tunnel-a operations The Drinfeld center of a category is denoted Z(C). When C is the representation category of a ﬁnite group G, Let us consider an arbitrary basis element Wb(γ) 0 of | i we have D(G) = Rep(D(G)) = Z(Rep(G)) = Z(VecG). the qudit. Diagrammatically, after applying the Wa(γ) We will adopt the following conventions for labeling operator, we have arrived in the following state: anyons, gapped boundaries, and their excitations:

1. Bulk excitations (a.k.a. anyons or topological charges), which are the simple objects within the modular tensor category B = Z(Rep(G)) or B = Z(C) will be labeled by a, b, c.... Their dual excita- (4) tions are labeled by a, b, c, ..., respectively.

2. The gapped boundary will be given as a Lagrangian algebra A which is an object in B. Here, and for the rest of the section, solid black lines are 3. Excitations on the boundary will be labeled as used to indicate a basis element of the hom-space that α, β, γ, .... When necessary, the local degrees of describes the ground state, while solid red lines are used freedom during condensation will be labeled as to denote Wilson operators. µ, ν, λ, .... To express this in terms of our original basis, we must convert the two anyon-tunneling ribbon operators into Furthermore, when using any F symbols and R sym- one. To do this, we can ﬁrst apply the M-3j opera- bols for a fusion category or a modular tensor category, tor and its Hermitian conjugate to the bottom and top we will adopt the following conventions for indices: boundaries of (4), respectively, to get1:

(1) (5)

ab Here, Mc (Ai) indicates that the M-3j symbol is for the gapped boundary given by the Lagrangian algebra Ai. . (2) We are now left with a bubble in the bulk. This can be eliminated using θ symbols of the bulk modular tensor category, by the following relation: arXiv:1707.05490v2 [quant-ph] 1 Nov 2017

Given a gapped boundary A and anyons a, b, c that condense to vacuum on the boundary, the associated M 3j symbols are deﬁned as in Refs. [8,9]: −

(6)

(3) Hence, we have the following equation:

These M symbols encode the associativity of bulk anyon 1 fusion and condensation to vacuum on the boundary. In this analysis, we will drop the multiplicity indices µ, ν, λ from the M symbols for concision; the generalization is obvious. 2

W (γ)W (γ) 0 = W (γ) a b | i a (9)

ab ab dadb The right hand side of Eq. (9) may be simpliﬁed using = M (A )[M ]†(A ) (7) c 1 c 2 d the deﬁnition of the S matrix as follows: Suppose c r c X

In general, Wa(γ) is unitary if and only if a cannot . (10) condense to any (non-vacuum) excitation on the boundary. Then, taking traces on both sides, we get

Loop-a operations

Another topological operation we can consider is to . (11) create a pair of anyons a, a in the bulk, move a counter- clockwise around a gapped boundary, and come back and annihilate the pair to vacuum. Physically, this corre- By deﬁnition of the S matrix of the modular tensor sponds to applying the a ribbon operator to a counter- category B [1], we have λ = Sab . Hence, we have the db clockwise closed ribbon encircling the gapped boundary. formula We will denote this operation by Wa(αi), where αi is the closed ribbon encircling boundary i. This operator is often known as the Wilson loop operator.

Suppose we have two gapped boundaries given by La- Sab grangian algebras A1, A2, which encode a qudit with or- Wa(α2) = (12) db thonormal basis as in Fig. ??. As before, we would like to compute the result of applying each Wa(αi) on each basis element of the ground state Hom(1, A A ), and 1 2 As before, W (α ) gives a d d matrix that acts on the express the result in terms of the original basis.⊗ a i ground state Hom(1, A A ×). Suppose we start as an arbitrary basis element 1 2 As in the case of the tunneling⊗ operator, the loop op- W (γ) 0 . Diagrammatically, the operator W (α ) trans- b a 2 erator W (α ) also need not be Hermitian or unitary. forms| thisi basis element into the following state: a i In general, it is just a Wilson loop operator, which is a holonomy.

Braiding gapped boundaries

2 (8) We now present the derivation for the formula of σ2. Consider an arbitrary basis element of Hom(1B, A1 2 ⊗ A2 A3 A4) as our start state. After applying σ2, we have⊗2: ⊗

Since we are working in a model where the total charge is 2 For the rest of this derivation, we assume for simplicity of il- vacuum, we may consider this picture as two holes on a lustration and computation that all anyons are self-dual. The sphere. We can hence push the anyon loop back through generalization is obvious, but one just needs to be more careful inﬁnity, to get in drawing orientations for each edge and using F and R symbols. 3

2 σ2 = (13)

FIG. 1. Topological charge projection (n = 2). As with the earlier cases, we would like to express the right hand side of Eq. (13) in terms of the original basis. 1. For each i = 1, 2, ...n, let αi be a simple loop en- We can apply an F -move to get: circling hole i (oriented counter-clockwise). Then α Γ(Y ). i ∈ 2. For each pair 1 i < j < n, let γij be a simple arc connecting hole≤ i and hole j (oriented to point a a b from i to j). Then γ Γ(Y ). = F 2 1 1 (14) ij ∈ b2;c0c Examples of loops and arcs in Γ(Y ) are shown in Fig. c X0 ??, when n = 2. Each knot diagram in Y can be resolved using graphical calculus to a linear combination of loops Next, applying two R-moves gives: and arcs in Γ(Y ). Next, we construct a basis LW(B, Ai ) for the Wilson operator algebra W(B, A ): { } { i} 1. For each anyon a B, and for each loop α Γ(Y ), ∈ i ∈ the Wilson loop operator Wa(αi) is a basis element. a1b1 b1a1 = R R (15) 2. For each pair 1 i < j < n, let Ai,j be the set of all c0 c0 anyon types a ≤B such that a (resp. a) condenses on the j-th (resp.∈ i-th) boundary. Then, for each a A , the Wilson line operator W (γ ) is a basis ∈ i,j a ij Finally, we apply one more F -move, which gives: element. We posit that any Hermitian operator O ∈ W(B, Ai ) = Span(LW(B, Ai )) can be measured. Such operators{ } O are called topological{ } charge measure-

a2a1b1 1 ment operators. The corresponding projective measure- = (Fb )c− c (16) 2 00 0 ments PO are called topological charge measurements c X00 (TCM). To be physical, we consider only operators which are monomials of basis operators. Hence, we see that action of the pure braid group gen- One special case of TCM, namely topological charge 2 erator σ on an arbitrary basis vector of Hom(1B, A projection [4], has been studied for doubled theories B. 2 1 ⊗ A2 A3 A4) is given by the following formula: In this case, B splits into two theories B = C C which ⊗ ⊗ do not interact in the bulk, but are “stuck together” at the original boundaries of B. The planar region Y also

2 a2a1b1 b1a1 a1b1 a2a1b1 1 splits into two mirror layers, S+(Y ) and S (Y ), which σ2 = F R Rc (F )− − b2;c0c c0 b2 c00c0 are completely disjoint in the bulk but “stuck together” Xc,c0 at the boundaries of Y . Fig.1 illustrates this for n = 2. (17) Notice that S+ and S together becomes a surface with genus 2. − Consider the Wilson operator algebra W(B, A ) in { i} this context. Each loop αi Γ(Y ) becomes a loop lαi Appendix C: General definition of topological charge ∈ in S+(Y ) or S (Y ), while each arc γi Γ(Y ) lifts to a measurement − ∈ loop lγi going around both layers. Let β be one of these loops. Define Ox(β) = Wx(αi) (tunneling operator in C) B If is a DW theory, gapped boundary braiding gen- if β is the lifting of the line αi, and Ox(β) = Wxx(γi) erates only a finite group [12]. Inspired by the results (loop operator in B) if β is the lifting of the loop γi. The of Ref. [4], we introduce topological charge measurement Wilson operator measuring charge a through β gives the (TCM) based on the Wilson operator algebra W(B, Ai ) TCM [4] for the symmetries of the theory at low energy. { } First, we construct a set Γ(Y ) of simple loops and arcs (a) Pβ = S0aSxa∗ Ox(β). (18) in Y : x C X∈ 4

The sum runs over the anyon labels x of C, and Sab is cal Phases of Matter. Commun. Math. Phys. (2017). the modular S-matrix of C. The Wilson operators Wx(αi) doi:10.1007/s00220-017-2960-4 [10] S. X. Cui, Z. Wang. Universal quantum computation with and Wxx(γi) are computed using the formulas (??) and (??) with the data of C and B, respectively. metaplectic anyons. J. Math. Phys. 56, 032202 (2015). [11] S. X. Cui, S.-M. Hong, Z. Wang. Universal quantum computation with weakly integral anyons. Quantum Inf. Pro- cess. 14:26872727 (2015). [12] N. Escobar-Velásquez, C. Galindo, Z. Wang. Braid Group Representations from Braiding Gapped Bound- aries of Dijkgraaf-Witten Theories. arXiv preprint [1] B. Bakalov, A. A. Kirillov. Lectures on tensor categories arXiv:1707.03884 (2017). and modular functors. Vol. 21. American Mathematical [13] A.G. Fowler, M. Mariantoni, J. M. Martinis, A. N. Cle- Soc. (2001). land. Surface codes: Towards practical large-scale quan- [2] M. Barkeshli. Charge 2e/3 superconductivity and topo- tum computation. Phys. Rev. A 86(3), 032324 (2012). logical degeneracies without localized zero modes in bi- [14] M.H. Freedman. P/NP, and the quantum field computer. layer fractional quantum Hall states. Phys. Rev. Lett. Proceedings of the National Academy of Sciences 95(1), 117, 096803 (2016). 98-101 (1998). [3] M. Barkeshli, P. Bonderson, M. Cheng, Z. Wang. Sym- [15] M. Freedman, A. Kitaev, M. Larsen, Z. Wang. Topo- metry, defects, and gauging of topological phases. arXiv logical quantum computation. Bulletin of the American preprint arXiv:1410.4540 (2014). Mathematical Society, 40(1), 31-38 (2003). [4] M. Barkeshli, M. Freedman. Modular transformations [16] A. Y. Kitaev. Fault-tolerant quantum computation by through sequences of topological charge projections. arXiv anyons. Ann. Phys. 303(2) (2003). preprint arXiv:1602.01093 (2016). [17] A. Kitaev, L. Kong, Models for Gapped Boundaries and [5] D. Beckman, D. Gottesman, A. Kitaev, J. Preskill. Mea- Domain Walls. Commun. Math. Phys. 313, 351-373 surability of Wilson loop operators. Phys. Rev. D 65, (2012). doi: 10.1007/s00220-012-1500-5. 065022 (2002). [18] T. Lan, J. C. Wang, X.-G. Wen. Gapped Domain Walls, [6] S. Beigi, P. W. Shor, D. Whalen. The quantum dou- Gapped Boundaries, and Topological Degeneracy. Phys. ble model with boundary: condensations and symmetries. Rev. Lett. 114(7), 076402 (2015). Communications in Mathematical Physics 306(3), 663- [19] M. A. Levin, X.-G. Wen. String-net condensation: A 694 (2011). physical mechanism for topological phases. Phys. Rev. B [7] S. B. Bravyi, A. Y. Kitaev. Quantum codes on a lattice 71, 045110 (2005). with boundary. arXiv:quant-ph/9811052 (1998). [20] M. A. Levin, A. Stern. Fractional Topological Insulators. [8] I. Cong, M. Cheng, Z. Wang. Topological quantum Phys. Rev. Lett. 103, 196803 (2009) computation with gapped boundaries. arXiv preprint [21] M. A. Levin, A. Stern. Classification and analysis of arXiv:1609.02037 (2016). two dimensional abelian fractional topological insulators. [9] I. Cong, M. Cheng, Z. Wang. Hamiltonian and Al- Phys. Rev. B 86, 115131 (2012) gebraic Theories of Gapped Boundaries in Topologi-