arxiv: v3 [quant-ph] 23 Jun 2011

Transcription

1 Feasibility of self-correcting quantum memory and thermal stability of topological order Beni Yoshida Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA arxiv: v3 [quant-ph] 23 Jun 2011 Abstract Recently, it has become apparent that the thermal stability of topologically ordered systems at finite temperature, as discussed in condensed matter physics, can be studied by addressing the feasibility of self-correcting quantum memory, as discussed in quantum information science. Here, with this correspondence in mind, we propose a model of quantum codes that may cover a large class of physically realizable quantum memory. The model is supported by a certain class of gapped spin Hamiltonians, called stabilizer Hamiltonians, with translation symmetries and a small number of ground states that does not grow with the system size. We show that the model does not work as selfcorrecting quantum memory due to a certain topological constraint on geometric shapes of its logical operators. This quantum coding theoretical result implies that systems covered or approximated by the model cannot have thermally stable topological order, meaning that no system can be stable against both thermal fluctuations and local perturbations simultaneously in two and three spatial dimensions. Keywords: self-correcting quantum memory, thermal stability of topological order, quantum coding theory, topological phase, stabilizer formalism, topological quantum field theory Contents 1 Introduction 2 2 Stabilizer code with physical constraints Stabilizer code Stabilizer code with Translation and Scale symmetries Feasibility of self-correcting quantum memory Self-correction, energy barrier and storage time Dimensional duality of logical operators Thermal stability of topological order Classical equivalence Existence of topological order and quantum codes Quantum equivalence Emergence of topology in logical operators 25 6 Summary and open questions 31 address: rouge@mit.edu (Beni Yoshida) Preprint submitted to Elsevier June 24, 2011

2 Appendix A Topological order at finite temperature 33 Appendix A.1 Topological order parameters and logical operators Appendix A.2 Topological order at finite temperature Appendix A.3 Thermal stability of topological order and energy barrier Appendix B Topology and quantum codes 38 Appendix B.1 Dimensional duality as a corollary of continuous deformability Appendix B.2 D-dimensional STS model and D + 1-dimensional TQFT Appendix C Decomposition of logical operators 44 Appendix C.1 Sketch of proof of theorem Appendix C.2 Identity generating matrix Appendix C.3 Existence of an odd matrix for m = Appendix C.4 Characteristic vectors Appendix C.5 Proof of lemma Appendix C.6 The existence of an odd matrix for m > Appendix D Derivation of logical operators Introduction In recent years, ideas from quantum information science have become increasingly useful in condensed matter physics where quantum information theoretical viewpoints give an important insight on studies of many-body entanglement arising in ground states or quasi-particle excitations of correlated spin systems [1 11]. In particular, it has been realized that many interesting physical systems in condensed matter physics may be described in the language of quantum codes such as stabilizer codes and its generalization [12 16]. The emerging closeness between two fields has made it possible to study several problems concerning many-body correlated spin systems through quantum coding theoretical tools [17]. In the present paper, we shall explore such a capacity of quantum coding theory for studying problems in condensed matter physics further, by making use of the correspondence between the following two open problems: (a) Feasibility of self-correcting quantum memory. (b) Thermal stability of topological order. Feasibility of self-correcting quantum memory is an important open problem in quantum information science concerning reliable storage of qubits. Thermal stability of topological order at finite temperature is an open problem of fundamental importance in condensed matter physics which concerns whether topological order may survive at finite temperature or not. While these two problems may look very different from each other, they are fundamentally akin to each other, in a sense that by searching for self-correcting quantum memory, one can search for topological ordered spin systems which are stable at finite temperature [18 22]. This surprising correspondence enable us to address a condensed matter theoretical question, thermal stability of topological order, by analyzing a quantum information theoretical question, feasibility of self-correcting quantum memory. The main contributions of the present paper toward these two open problems are the followings; We propose a certain model of quantum codes which may cover a large class of physically realizable quantum codes, which is called a stabilizer code with translation and scale symmetries. We solve the D-dimensional model exactly by determining its coding properties completely and analyze the feasibility of self-correcting quantum memory for D = 1, 2, 3. We establish the connection between self-correcting quantum memory and stable topological order at finite temperature, and analyze thermal stability of topological order arising in the model. 2

3 Figure 1: The correspondence between the feasibility of self-correcting quantum memory and the thermal stability of topological order, and the plan of the paper. The plan of the paper is schematically summarized in Fig. 1. We begin the introduction of the paper by reminding why these two problems are important in each field. Then, we give brief summaries of our results in the rest of the introduction. Feasibility of self-correcting quantum memory: Quantum entanglement decays easily. This underlying difficulty in quantum information science gave birth to the beautiful art of protecting qubits from decoherence; quantum coding theory. The central idea of quantum error- correcting codes is to encode a qubit in many-body entangled states, and perform error-corrections so that encoded qubits are not lost. After discoveries of first examples of quantum codes [12, 23 26] which culminated in stabilizer codes [27], a large number of quantum codes have been found. Now, quantum coding theory constitutes one of the most important building blocks for realizing fault-tolerant quantum computation [28]. Yet, there still remain important gaps between theoretical constructions of quantum codes and their physical realizations as quantum memory devices. First of all, most quantum coding schemes encode qubits dynamically by applying a large number of logical gates, and encoding is discussed only in terms of the Hilbert space. However, since encoded qubits will be eventually lost in the presence of interactions with the external environment, it is desirable to store logical qubits statically in some physical subspaces which are naturally protected from decoherence. One plausible approach may be to store logical qubits in the gapped ground space of some quantum many-body system. In this light, stabilizer codes constructed with geometrically local generators are promising candidates for physical realizations of quantum codes since such local stabilizer codes can be realized as the ground space of gapped Hamiltonians by using their local generators as interaction terms. Indeed, from purely theoretical viewpoint, it is known that sufficiently accurate and frequent error-corrections can store logical qubits reliably [18]. However, one may still dream of having a quantum memory device which would work without active error-corrections, given the difficulties and inefficiencies of performing fast and accurate error-corrections in reality 1. Self-correcting quantum memory is an ideal memory device which corrects errors by itself [16, 18 21, 29]. Due to the large energy barrier separating degenerate ground states, natural thermal dissipation processes restore the system into the original encoded states by correcting errors automatically without any active error-correction. If such a memory device could exist, it will be a perfect quantum information storage device which may be used commercially in the future. Also, the reliable storage of qubits seems to be the starting point for building scalable quantum computers. 1 To the best of our knowledge, there have been no convincing experimental demonstration of error-corrections. All the demonstrations are limited only to phase-type errors, and do not correct bit-type errors. Also, it seems very difficult to reach the frequency threshold for reliable storage of logical qubits [18]. Finally, it is extremely inefficient to keep performing error-corrections during the time one is storing a qubit if one stores a qubit for a long duration such as months or years. 3

4 There has been significant progress toward construction of self-correcting quantum memory. It has been pointed out that the Toric code defined on a four-dimensional system (D = 4) serves as self-correcting quantum memory [18, 30]. Yet, such a four-dimensional code cannot be embedded in a three-dimensional system. While there have been several proposals for three-dimensional self-correcting memory [16, 31, 32], validities of none of these proposals have been verified yet 2. In addition, it has been shown that a self-correcting local stabilizer code cannot exist in two-dimensional systems [20, 33]. It seems that two-dimensional spin systems are not promising as resource for self-correcting quantum memory. Now, the question we want to address can be summarized as follows: Is it possible to have three-dimensional self-correcting quantum memory? Thermal stability of topological order: The feasibility of self-correcting quantum memory is closely related to another important open problem in condensed matter physics; the thermal stability of topological order. Studies of topologically ordered systems [7, 8, 12, 34 39] have been frontiers of researches in condensed matter physics community, as systems with topological order are beyond the description of the Landau s symmetry-breaking paradigm which was once considered almost as theory of everything for studies of many-body systems. Topologically ordered systems are also of practical importance in quantum information science since many-body entanglement arising in ground states and quasi-particle excitations of topologically ordered systems is a primary resource for realizing various quantum information processing tasks [8, 12]. The notion of topological order was originally introduced in order to characterize the stability of ground states of many-body quantum systems against local perturbations [35]. Loosely speaking, a system is said to have topological order when its ground state properties do not change significantly under any types of small, but finite local perturbations. This stability of ground states against local perturbations is also valuable for quantum information processing since topologically ordered spin systems can be used as good quantum codes with macroscopic code distances [40]. However, the situation changes completely when one considers the effect of thermal fluctuations on topologically ordered systems. In fact, it is known that topological order in a two-dimensional Toric code is not stable at any finite temperature which may be quantitatively seen from the fact that topological order parameters such as topological entanglement entropy vanish at any non-zero temperature at the thermodynamic limit [41]. A similar result is obtained in a recent numerical work on topological entanglement entropy in a spin liquid model at finite temperature [42]. It seems that topological order in a two-dimensional system is not stable at finite temperature according to general studies on the ground state properties of two-dimensional frustration-free Hamiltonians [20, 40]. Now, the question concerning the stability of topological order can be summarized as follows 3 : Is there any system with thermally stable topological order, which is stable against both thermal fluctuations and local perturbations? Main results: Interestingly, this condensed matter theoretical question on the stability of topological order can be addressed through quantum coding theory. In fact, it has been pointed out that when topological order in correlated spin systems is stable at finite temperature, such a system can be used as self-correcting quantum memory [20, 43, 44]. This correspondence between self-correcting quantum memory and the thermal stability of topological order may be better understood by identifying thermal fluctuations as random errors acting on a quantum memory. Therefore, one may take a slight liberty and say that a search for self-correcting quantum memory is equivalent to a search for a novel quantum phase with topological order which is stable at finite temperature. In this paper, with this correspondence in mind, we analyze coding and physical properties of a certain model of local stabilizer codes with physically reasonable constraints [17]. The model, which is called Stabilizer code with Translation and Scale symmetries (STS model), is constrained to the following physical conditions. Qubits are defined on a D-dimensional square (hypercubic) lattice with periodic boundary conditions. The Hamiltonian consists only of geometrically local interaction terms with translation symmetries. The number of logical qubits does not grow with the system size (scale symmetries). 2 We shall discuss these proposals in detail in section The definition of topological order at finite temperature will be clearly stated in section 4 4

5 We show that the model has a certain dimensional duality on geometric shapes of logical operators, as summarized in the following informal theorem: Theorem (Dimensional Duality). In a D-dimensional STS model (D 3), m-dimensional and (D m)-dimensional logical operators always form anti-commuting pairs where m is an integer. Based on this dimensional duality on logical operators, we give answers to open questions (a) and (b), as summarized below: (a) The three-dimensional STS model does not work as self-correcting quantum memory since the energy barrier is finite for the encoding with respect to a two-dimensional logical operator. (b) Systems covered or approximated by the model cannot have thermally stable topological order, meaning that systems cannot be stable against both thermal fluctuations and local perturbations simultaneously in two and three spatial dimensions. Code distance: Our results on STS models also provide a partial answer to a long-standing open problem in quantum coding theory concerning the upper bound on the code distance of local stabilizer codes. The code distance d is a measure of the robustness of quantum codes against errors, and one of the ultimate goals in quantum coding theory is to find a quantum code with a large code distance for a fixed system size N (the total number of qubits). While an upper bound on the code distance of stabilizer codes is roughly known [45], the upper bound for local stabilizer codes is currently not known yet. In other words, despite the fact that the stabilizer formalism is a canonical framework in quantum coding theory, we still do not know how robust the best local stabilizer code can be! For a long time, it had been believed that the code distance of local stabilizer codes with N qubits is upper bounded by O( N): d O( N) at N since all the examples of local stabilizer codes ever found satisfied this upper bound [12, 18]. Later, an example of a local stabilizer code whose code distance scales as O( N log N) was found [46]. While the code distance of this local stabilizer code exceeds the previously believed upper bound O( N) logarithmically, an example of a local stabilizer code whose code distance exceeds O( N) polynomially has not been found yet. Therefore, the code distance of local stabilizer codes seemed to be upper bounded by O(N 1 2 +ɛ ) at N where ɛ is an arbitrary small positive number, although no analytical result was known on the upper bound. It was recently proven that the code distance of local stabilizer codes is upper bounded by O(L D 1 ) where L is a linear length of the system and D is the spatial dimension (N O(L D )) [20]. While this work does not rule out the possibilities for the existence of a local stabilizer code whose code distance exceeds O(N 1 2 ) polynomially, this bound was proven to be tight only for D = 1, 2. Thus, whether the tight upper bound is d O(N 1 2 +ɛ ) or d O(L D 1 ) at N (1) for D > 2 seems to be one of the most important open questions concerning coding properties of local stabilizer codes. Our analysis on STS models provides the following tight upper bound on the code distance of STS models: A three-dimensional STS model has a code distance which is tightly upper bounded by O(L) where L is the linear length of the system. Thus, in a three-dimensional system, the upper bound on code distances turns out to be more strict than not only O(L 2 ), but also O( N). Organization: The paper is organized as follows. In section 2, we give a brief review of stabilizer codes and introduce STS models. In section 3, we discuss the feasibility of self-correcting quantum memory. In section 4, we discuss the thermal stability of topological order at finite temperature. In section 5, we present a certain universal property of logical operators shared among all the STS models. In section 6, we discuss possible future problems. The definition of topological order at finite temperature is discussed further in Appendix A. A possible relevance to topological quantum field theory is discussed in Appendix B. The proof of the dimensional duality of logical operators is presented in Appendix C and Appendix D. The main part of the paper is self-consistent and accessible to readers without previous knowledge on quantum coding theory and topological order. However, the proof part is rather technical, and relies heavily on theoretical tools developed previously in [17, 47] by the author. In particular, we owe a lot of arguments to [17] which introduced and solved the two-dimensional STS model originally. 5

6 2. Stabilizer code with physical constraints While local stabilizer codes are physically realizable as quantum memory devices in principle, realistic physical systems are often constrained to not only the locality of interaction terms, but also various physical symmetries. In this section, we give the definition of the STS model which are local stabilizer codes with translation and scale symmetries [17]. In section 2.1, we give a brief review of stabilizer codes. In section 2.2, we describe the definition of the STS model Stabilizer code Here, we give a brief review of stabilizer codes which are quantum codes possessing Hamiltonians to support logical qubits in the ground space with a finite energy gap [27]. Some notations which will be used throughout this paper are also fixed here. Note that we shall use the notations {} for a set and for a group. Stabilizer formalism: The main idea of stabilizer codes is to encode k logical qubits into N physical qubits (N > k) by using a subspace V S spanned by states ψ that are invariant under the action of the stabilizer group S: V S = { ψ (C 2 ) N : U ψ = ψ, U S }. (2) Here, the stabilizer group S is an arbitrary Abelian subgroup of the Pauli group S P = ii, X 1, Z 1,..., X N, Z N (3) such that I S. The elements in S are called stabilizers. The logical subspace V S can be realized as the ground space of the following Hamiltonian (see Fig. 2) H = S j, S = S 1, S 2, (4) j since the energy eigenvalue is minimized for states satisfying S j ψ = ψ for all j. There are k logical qubits encoded in V S where k N G(S). Here, G(S) represents the number of independent generators in S. The ground space is separated from excited states by a finite energy gap since eigenstates are simultaneously diagonalized with respect to eigenvalues ±1 of S j. Figure 2: Energy spectrum of a stabilizer Hamiltonian. Logical qubits are encoded in the degenerate ground space. Logical operators: In analyzing properties of logical qubits stored in the ground space, operators called logical operators play central roles. Logical operators are Pauli operators which commute with the Hamiltonian, but not inside the stabilizer group S. Logical operators can be found inside the centralizer group: C = { U P : [U, S j ] = 0, for all j } (5) 6

7 which is a group of Pauli operators commuting with all the stabilizers. Then, a set of logical operators is L = { U C : U 2 = I, U S }. (6) Logical operators may transform encoded qubits since they act non-trivially inside the ground space V S. Equivalence relation: One may introduce an equivalence relation between logical operators by seeing how they act inside the ground space. Two logical operators l and l are said to be equivalent if and only if l and l act in the same way inside the ground space: l l l ψ = l ψ, ψ V S (7) ll S. (8) Therefore, logical operators remain equivalent under multiplications of stabilizers. Canonical form: It is often convenient to represent a set of 2k independent logical operators in the following canonical form [47]: { } l1,, l k. (9) r 1,, r k Here, l p and r p are independent logical operators whose commutation relations are {l p, r p } = 0, [l p, r q ] = 0 for p q, [l p, l q ] = 0 and [r p, r q ] = 0. Thus, only the operators in the same column anti-commute with each other. Note that choices of logical operators are not unique. Code distance: The code distance is a measure of the robustness of a quantum code, which is quantified by the minimal weight of logical operators: d = min(w(u)) where U L. (10) Here, w(u) denotes the number of non-trivial Pauli operators constituting U. The code distance corresponds to a minimal number of single Pauli errors necessary to destroy an encoded qubit. Roughly speaking, a quantum code with a large code distance can securely protect logical qubits. Bi-partition: It is often convenient to split the entire system of qubits into two complementary subsets of qubits in studying coding properties of stabilizer codes [48, 49]. Let us recall a useful formula to study geometric shapes of logical operators in stabilizer codes through a bi-partition. For a stabilizer code in a bi-partition, the following theorem is known to hold [47] (Fig. 3). Theorem 1 (Bi-partition). For a stabilizer code with k logical qubits, let the number of independent logical operators supported by a subset of qubits R be g R. Then, for an arbitrary bi-partition into two complementary subsets of qubits R and R, the numbers of logical operators supported by R and R obey the following constraint: g R + g R = 2k. (11) Figure 3: A bi-partition of a stabilizer code. Each dot represents a qubit. 7

8 This bi-partition theorem is useful for analyzing geometric sizes and geometric shapes of logical operators. For example, if we find a region R where there is no logical operator: g R = 0, we immediately know that all the logical operators can be supported inside R since g R = 2k [20]. Thus, one can restrict geometric regions of qubits where logical operators are supported. For an extension of the theorem to subsystem codes, see [50, 51] Stabilizer code with Translation and Scale symmetries Here, we describe the definition of STS models, which are local stabilizer codes with translation and scale symmetries. (1) Locality of interaction: Physically realistic systems must have only geometrically local interaction terms. To introduce the notion of locality to stabilizer codes, we consider a system of qubits defined on a D-dimensional square lattice (hypercubic lattice) which consists of N = L 1 L D qubits where L m is the total number of qubits in the ˆm direction for m = 1,, D. Therefore, qubits are distributed in the physical space with a metric. Here, the entire system is separated into a collection of hypercubes which consists of v = v 1 v D qubits without overlaps by assuming that n m L m /v m are integer values (see Fig. 4). We consider a block of v = v 1 v D qubits as the single unit block which constitutes the entire system. In particular, we consider these unit blocks as single composite particles with a larger Hilbert space (C 2 ) v (Fig. 4). Thus, the entire system is viewed as a hypercubic lattice of n 1 n D composite particles. Now, we assume that interaction terms of the stabilizer Hamiltonian are defined locally: H = S j (12) where S j are supported inside some regions with 2 2 composite particles (Fig. 4). (Otherwise, we coarse-grain the system). In this paper, instead of qubits, we consider composite particles as the smallest building blocks of the system. j Figure 4: An illustration of the STS model. A two-dimensional example is shown where a unit block of 3 2 qubits is considered as a composite particle with a larger Hilbert space. Interaction terms S j are defined locally inside a region of 2 2 composite particles. The Hamiltonian is invariant under unit translations of composite particles. (2) Translation symmetries: Physically realistic systems often have not only local interactions, but also some physical symmetries. Here, we assume that the stabilizer Hamiltonian possesses translation symmetries: T m (H) = H (m = 1,, D) (13) where T m represent unit translations of composite particles in the ˆm direction (Fig. 4). 8

9 For simplicity of discussion and in order to accommodate translation symmetries, we set the periodic boundary conditions. Then, the entire system may be viewed as a D-dimensional torus: T D = S 1 S 1 where S 1 is a circle. Thus, the entire system has a topologically non-trivial geometric shape a priori. (3) Scale symmetries: In this paper, we are interested in coding properties at the limit where the system size goes to the infinity (in other words, at the thermodynamic limit). So far, we have considered the cases where the system size n (n 1,, n D ) is fixed. Here, we consider changes of the number of composite particles n m while keeping interaction terms S j the same. It is commonly believed that there is a tradeoff between the number of logical qubits k and the code distance [40, 51] where the code distance d decreases as the number of logical qubits k increases for a fixed N. Then, it may be legitimate to limit our considerations to the cases where the number of logical qubits k remains small when the system size increases. We assume that stabilizer codes have scale symmetries by requiring that the number of logical qubits k n is independent of the system size n: k n = k, n. (14) Here, we emphasize that, in a system with scale symmetries, the number of logical qubits k remains constant under not only global scale transformations: n c n where c is some positive integer, but also arbitrary changes of n m. One might think that scale symmetries are too strong as physical constraints. However, through appropriate coarse-graining, a fairly large class of local stabilizer codes with translation symmetries can be considered as the STS model. For example, let us consider the cases where the number of logical qubits k n is small: k n k 0, n (15) where k n does not grow with the system size n. Then, there always exists some finite coarse-graining such that a coarse-grained system has scale symmetries, as proven in [17]. As a result, by analyzing coding properties of stabilizer codes in the presence of scale symmetries, one can easily deduce coding properties of stabilizer codes which satisfy Eq. (15). Therefore, solutions of STS models are sufficient to discuss coding properties of translation symmetric stabilizer codes with a small number of logical qubits. Note that STS models cover the Toric code and its generalizations to D-dimensional systems. Translation equivalence of logical operators: There is a certain property of logical operators which emerges naturally as a result of translation and scale symmetries. For STS models, the following theorem holds [17] (Fig. 5). Theorem 2 (Translation equivalence). For each and every logical operator l in an STS model, a unit translation of l with respect to composite particles in any direction is always equivalent to the original logical operator l: T m (l) l, l L n (m = 1,, D) (16) where L n is a set of all the logical operators for an STS model defined with the system size n. Figure 5: The translation equivalence of logical operators. Each square represents a composite particle. Shaded regions represent translated logical operators which are equivalent to each other. Here, we only give an intuition on why this theorem holds. Let us consider the case where the system size n is large. Then, since the number of logical qubits k does not depend on the system size, k is relatively small compared 9

10 with the system size n. Now, due to the translation symmetries of the system Hamiltonian, translations of a given logical operator l are also logical operators. However, there are only 2k independent logical operators. Then, there must be a finite integer a m such that l T a m m (l) for all the logical operators l. (Otherwise, there would be so many independent logical operators). It turns out that a m = 1 for any l and m. While we have used only the condition that the number of logical operators k is small, due to scale symmetries (k is constant), one can prove the above theorem by showing a m = 1 for any l, m and n. 3. Feasibility of self-correcting quantum memory In this section, we present one of our main results in this paper, concerning coding the feasibility of self-correcting quantum memory. In section 3.1, we begin by giving brief discussion on the relation between the qubit storage time and the energy barrier. In section 3.2, we describe geometric shapes of logical operators in STS models, and discuss whether a three-dimensional STS model works as self-correcting quantum memory or not Self-correction, energy barrier and storage time In subsection, we review how self-correcting quantum memory works by establishing the connection between the energy barrier and the qubit storage time. Self-correcting classical memory: For simplicity of discussion, let us start by analyzing an example of selfcorrecting classical memory. Consider two-dimensional Ising model: H = Z i, j Z i+1, j Z i, j Z i, j+1 (17) i, j which consists of L L qubits with periodic boundary conditions. The model works as a classical code since one can encode a classical bit in the ground space by labeling 0 0 as 0 and 1 1 as 1. Now, let us see why this model works as self-correcting classical memory. Assume that the system is originally 0 0. Then, in order for errors to change a ground state 0 0 into another ground state 1 1, errors must flip all the spins from 0 to 1. However, during these spin flips, the excitation energy becomes at least O(L) because there is a domain wall separating the regions with 0 s and 1 s (Fig. 6). In other words, ground states 0 0 and 1 1 are separated by a large energy barrier. Then, before errors accumulate, natural thermal dissipation processes restore the system into the original encoded state 4. Therefore, the system corrects errors by itself. One may estimate the bit storage time of two-dimensional Ising model by using the so-called Arrhenius law: i, j τ exp( E/T) (18) up to some polynomial corrections where τ is the storage time, E is an energy barrier and T is the temperature. According to the Arrhenius law, the bit storage time can be estimated as τ EXP(L) since the energy barrier is E O(L). While it should be noted that this is an empirical law, it is commonly believed that the law correctly estimates the bit or qubit storage time of spin systems (for systems with well defined energy barrier) below a critical temperature. Indeed, the law is rigorously verified in various models of classical and quantum memory [19, 21, 22, 52, 53]. However, it should be noted that verification of the Arrhenius law for classical and quantum memory is usually a very difficult task, involving an evaluation of time evolution of the system in a presence of interactions with an external environment. Next, let us consider an example of a classical code which is not self-correcting; one-dimensional Ising model: H = Z j Z j+1. (19) j 4 Precisely speaking, the system does not return to the original state, but return to a state which is sufficiently close to the original state with a probability approaching to unity at the thermodynamic limit. Therefore, one can reliably read out the encoded bit from such a state. 10

11 One may see that this model does not work as self-correcting classical memory since one can change from 0 0 to 1 1 by costing only a finite energy. In other words, one can create a kink by costing only a finite energy: (20) which may propagate the lattice freely without costing any extra energy. Then, according to the Arrhenius law, the bit storage time is τ O(1) which is independent of the system size since the energy barrier is E O(1). Figure 6: How self-correcting classical memory works in two-dimensional Ising model. Energy barrier and logical operators: One can associate the self-correcting property of two-dimensional Ising model with geometric shapes of logical operators by viewing the model as a stabilizer code. Note that D-dimensional Ising model satisfies the definition of STS models since interaction terms are local and translation symmetric, and there is a single logical qubit (k = 1 and d = 1) regardless of the system size. Two-dimensional Ising model has the following pair of zero-dimensional and two-dimensional logical operators: l = Z 1,1, r = X i, j. (21) Then, a classical bit is encoded in eigenstates of a zero-dimensional logical operator l. In order to change the encoded bit, one needs to apply a two-dimensional logical operator r since 1 1 = r 0 0. Then, an intermediate state during the change from 0 0 to 1 1 may be represented as r 0 0 where r is some subpart of the original two-dimensional logical operator r (see Fig 7(a)). Since interaction terms anti-commute with Pauli operators at the boundary of r, the excitation energy associated with r 0 0 is proportional to the perimeter of r. Thus, during the change from 0 0 to 1 1, the excitation energy must become O(L) since the perimeters of subparts of a two-dimensional logical operator r are always one-dimensional. One can also understand why one-dimensional Ising model does not work as self-correcting classical memory through geometric shapes of logical operators. One-dimensional Ising model has the following pair of zerodimensional and one-dimensional logical operators: l = Z 1 and r = i X i. Since a subpart r of a one-dimensional logical operator r is always zero-dimensional, the energy barrier is E O(1) (see Fig 7(b)). Self-correcting quantum memory: Next, let us discuss how self-correcting quantum memory works. It is known that many of good local stabilizer codes with macroscopic code distances do not have self-correcting properties since they have string-like logical operators which lead to O(1) energy barrier. As an example, let us consider 11 i, j

12 Figure 7: (a) A subpart of a two-dimensional logical operator. The excitation energy is O(L). (b) A subpart of a one-dimensional logical operator. The excitation energy is O(1). two-dimensional Toric code: H Toric = A s B p (22) s p where qubits live on edges of L L square lattice. Interaction terms A s and B p are shown in Fig. 8. Two-dimensional Toric code can be viewed as an STS model since the number of logical qubits is constant: k = 2 regardless of the system size. Figure 8: (a) A subpart of a two-dimensional logical operator. The excitation energy is O(L). (b) A subpart of a one-dimensional logical operator. The excitation energy is O(1). One can understand why two-dimensional Toric code does not work as self-correcting quantum memory through geometric shapes of logical operators. Since the Toric code has pairs of one-dimensional logical operators as shown in Fig. 8, the energy barrier is E O(1). As a result, the qubit storage time is τ O(1) according to the Arrhenius law. A similar discussion holds for three-dimensional Toric code too. As two-dimensional and three-dimensional Toric codes have O(1) qubit storage time, they are not reliable quantum memories. Then, what kinds of quantum codes may serve as reliable quantum memory devices? It is worth noting that two-dimensional Toric code works as a reliable quantum memory device if one perform sufficiently frequent and accurate error-corrections. In particular, it has been shown [18] that the qubit storage time can be exponentially long: τ EXP(L) in the presence of active error-corrections; however, it seems very difficult to reach such accuracy and frequency which are necessary for reliably storing qubits. In the same paper [18], it has been also pointed out that four-dimensional Toric code may have exponentially long storage time τ EXP(L) below the critical temperature since the model has a large energy barrier which scales as O(L). This remarkable insight has been further investigated in [30], and later, rigorously verified in [22]. One can associate the self-correcting property of four-dimensional Toric code with geometric shapes of logical operators, as 12

13 the model has pairs of two-dimensional logical operators which lead to O(L) energy barrier. Yet, since one cannot embed four-dimensional Toric code in a three-dimensional space, one hopes to have three-dimensional self-correcting quantum memory whose qubit storage time grows as the system size increases. Previous proposals: So far, there have been several interesting proposals of three-dimensional self-correcting quantum memory. First, a certain model of three-dimensional subsystem codes [16], called three-dimensional quantum compass model, was proposed as a candidate of self-correcting quantum memory. Subsystem codes may be considered as generalizations of stabilizer codes which are supported by Hamiltonians whose interaction terms are Pauli operators, but may anti-commute with each other. While the model opened a new possibility of quantum memory devices supported by frustrated Hamiltonians and initiated studies of Hamiltonian realizations of subsystem codes, its validity as self-correcting quantum memory has not been verified since it is difficult to solve the frustrated Hamiltonian, and its properties are hard to determine. Also, it seems difficult for such frustrated systems to have the thermal stability since the model supports gapless excitations. Furthermore, its validity has been denied in several literature including [19]. Later, an interesting model of the mixture of two-dimensional Toric code and three-dimensional Bosonic gas, called the Toric-Boson model, was proposed [31] where Bosonic gas induces confining potential between anyonic excitations. While the model opened new capacities of quantum codes constructed in the so called mixed-dimensional configurations, which are of particular interest in ultracold atom physics community, the model itself has two serious drawbacks as a candidate of self-correcting quantum memory. First, its storage time is only polynomial in L: τ POLY(L) since an effective energy barrier is: E LOG(L). However, it is known that it takes at least O(d) gate operations, which cannot be implemented simultaneously, to read out the encoded qubit, where d is the code distance of the system. Then, since the code distance d scales polynomially in L in the model, polynomially long storage time is not sufficient for a model to work as efficient quantum memory device [21]. Second, the model needs a very strong coupling between the Toric code and Bosonic gas whose strength scales polynomially in L. Therefore, the model has a problem in the scalability. Recently, an interesting proposal of three-dimensional spin glass models [32] has been made whose relaxation dynamics is very slow due to the existence of a large number of energy local minima. The model seems to have a polynomially long storage time: τ POLY(L) with logarithmically large energy barrier: E LOG(L), though it needs verifications (so, it could be longer or shorter). The system undergoes a phase transition at T = 0, which mat imply a potential thermal instability of the system properties. Finally, it seems difficult o find an efficient decoding scheme for the model. Here, it should be noted that the model does not have scale symmetries since the number of logical qubits is highly sensitive to the system size, and there is no finite upper bound on it. See section 6 for discussion on stabilizer codes without scale symmetries too. At this moment, validities of none of these proposals of three-dimensional self-correcting quantum memory has been verified yet, and the feasibility of self-correcting quantum memory remains open. Below, we summarize the qubit storage time and the energy barrier in these proposals of self-correcting quantum memory. Code Energy Barrier Storage Time Comment 2-dim Toric Code O(1) O(1) 3-dim Toric Code O(1) O(1) 2-dim Toric Code O(1) EXP(L) with error-correction 4-dim Toric Code O(L) EXP(L) self-correcting 3-dim subsystem code?? gapless excitation? Toric-Boson model LOG(L) POLY(L) POLY(L) coupling Spin glass model LOG(L) POLY(L) T = 0 phase transition 3.2. Dimensional duality of logical operators In this subsection, we describe geometric shapes of logical operators in STS models, and discuss the feasibility of self-correcting quantum memory. One-dimension: We begin by reviewing geometric shapes of logical operators in one-dimensional STS models. In a one-dimensional chain of n 1 composite particles, we denote jth composite particle as P j. Then, the following theorem holds [17]. 13

14 Theorem 3 (Dimensional Duality). For a one-dimensional STS model, there exists a canonical set of logical operators: { } l1,, l k (23) r 1,, r k whose pair of anti-commuting operators l j and r j has one of the following property ( j = 1,, k). l j is a zero-dimensional logical operator defined inside P 1, while r j is a one-dimensional logical operator defined in a periodic way: T 1 (r j ) = r j. It is worth presenting geometric shapes of logical operators graphically (Fig. 9). The code distance of onedimensional models is upper bounded by v which is the inner dimension of a composite particle. Since there always exist zero-dimensional logical operators, such systems cannot have topological order. Figure 9: Logical operators in a one-dimensional system. Two-dimension: Next, we analyze geometric shapes of logical operators for two-dimensional STS models. Let us first introduce some regions of composite particles in order to define geometric shapes of logical operators concisely. A square region of x 1 x 2 composite particles is denoted as P(x 1, x 2 ) (Fig. 10(a)): P(x 1, x 2 ) { P r1,r 2 : 1 r 1 x 1, 1 r 2 x 2 } (24) where 1 x 1 n 1 and 1 x 2 n 2. Note that a composite particle at the position (r 1, r 2 ) is denoted as P r1,r 2. Then, the following theorem holds [17]. Theorem 4 (Dimensional Duality). For a two-dimensional STS model, there exists a canonical set of logical operators: { } l1,, l k (25) r 1,, r k whose pair of anti-commuting operators l j and r j has one of the following two properties ( j = 1,, k). l j is a zero-dimensional logical operator defined inside P(1, 2v), while r j is a two-dimensional logical operator defined in a periodic way: T 1 (r j ) = r j and T 2 (r j ) = r j. l j is a one-dimensional logical operator defined inside P(1, n 2 ) in a periodic way: T 2 (l j ) = l j, while r j is a one-dimensional logical operator defined inside P(n 1, 1) in a periodic way: T 1 (r j ) = r j. It is worth presenting geometric shapes of logical operators graphically (Fig. 10(b)). There is a dimensional duality on geometric shapes of logical operators as follows: { } 0 dim, 1 dim. (26) 2 dim, 1 dim 14

15 Logical operators are periodic in the directions in which they are stretched. The code distance is upper bounded by O(L) which is consistent with the upper bound O(L D 1 ) obtained in [20]. The existence of one-dimensional logical operators gives rise to topological order. Figure 10: Dimensional duality in a two-dimensional system. (a) A region of x 1 x 2 composite particles is denoted as P(x 1, x 2 ). (b) Geometric shapes of logical operators in a two-dimensional STS model. Three-dimension: Finally, let us proceed to coding properties of a three-dimensional STS model. A region with x 1 x 2 x 3 composite particles is denoted as P(x 1, x 2, x 3 ): P(x 1, x 2, x 3 ) { P r1,r 2,r 3 : 1 r m x m, m = 1, 2, 3 } (27) where P r1,r 2,r 3 represents a composite particle at (r 1, r 2, r 3 ). Then, for logical operators in a three-dimensional STS model, the following theorem holds. Theorem 5 (Dimensional Duality). For a three-dimensional STS model, there exists a canonical set of logical operators: { } l1,, l k (28) r 1,, r k whose pair of anti-commuting operators l j and r j has one of the following four properties. l j is a zero-dimensional logical operator defined inside P(1, 2v, (2v) 2 ), while r j is a three-dimensional logical operator defined in a periodic way: T 1 (r j ) = r j, T 2 (r j ) = r j and T 3 (r j ) = r j. l j is a one-dimensional logical operator defined inside P(n 1, 2v, 1) in a periodic way: T 1 (l j ) = l j, while r j is a two-dimensional logical operator defined inside P(1, n 2, n 3 ) in a periodic way: T 2 (r j ) = r j and T 3 (r j ) = r j. l j is a one-dimensional logical operator defined inside P(1, n 2, 2v) in a periodic way: T 2 (l j ) = l j, while r j is a two-dimensional logical operator defined inside P(n 1, 1, n 3 ) in a periodic way: T 1 (r j ) = r j and T 3 (r j ) = r j. l j is a one-dimensional logical operator defined inside P(2v, 1, n 3 ) in a periodic way: T 3 (l j ) = l j, while r j is a two-dimensional logical operator defined inside P(n 1, n 2, 1) in a periodic way: T 1 (r j ) = r j and T 2 (r j ) = r j. This is the main technical result of the present paper. We present the proof of the theorem in appendices Appendix C and Appendix D. It is worth representing geometric shapes of logical operators graphically (Fig. 11). Note that 15

16 logical operators are periodic in the directions in which they are stretched. There is a dimensional duality on geometric shapes of logical operators as follows: { } 0 dim, 1 dim (29) 3 dim, 2 dim Figure 11: Dimensional duality in a three-dimensional system. As a result of this dimensional duality, one may find the upper bound on the code distance. When n 1 = n 2 = n 3 = L, the code distance of a three-dimensional STS model is upper bounded as follows: d 2vL O(L). (30) Note that this bound is tight for the three-dimensional Toric code. The system may have topological order when it is free from zero-dimensional logical operators. Higher-dimensions: Though our primary interests are in coding properties of three-dimensional STS models, it is worth extending the analysis to higher dimensions. For a D-dimensional STS model (D 4), we make the following conjectures: In a D-dimensional system, m-dimensional and (D m)-dimensional logical operators form anti-commuting pairs where m is an integer. The code distance is tightly upper bounded by O(L D 2 ) when D is even and by O(L D 1 2 ) when D is odd. Note that generalizations of the Toric code to D-dimensional systems have the above dimensional duality for arbitrary integer m. See Appendix B.2 for constructions of D-dimensional Toric code. Such a dimensional duality arising in geometric shapes of logical operators may be viewed as a manifestation of Poincaré duality on coding properties of stabilizer codes supported on a D-torus, as very briefly discussed in Appendix B. However, we report that straightforward generalizations of analysis tools presented in this paper did not work for STS models with D 4, and the conjecture above needs to be verified. Feasibility for a three-dimensional STS model: Finally, let us show that three-dimensional STS models do not work as self-correcting quantum memory. First, in order for the system to work as a quantum code, there must be an anti-commuting pair of one-dimensional and two-dimensional logical operators. We denote them as l and r 16

17 respectively. Let us assume that the system is initially in the eigenstate of l = 1 denoted as ψ(l = 1). Then, ψ(l = 1) = r ψ(l = 1) and r is a two-dimensional logical operator. Since the boundary of a subpart of r is onedimensional, the excitation energies associated with intermediate states are O(L). Thus, the encoding with respect to l is self-correcting. Next, let us assume that the system is initially in the eigenstate of r = 1 denoted as ψ(r = 1). Then, ψ(r = 1) = l ψ(r = 1) and l is a two-dimensional logical operator. Since the boundary of a subpart of l is zero-dimensional, the excitation energies associated with intermediate states are O(1). Thus, the encoding with respect to r is not selfcorrecting, and the qubit storage time is τ O(1), independent of the system size. Therefore, while such a system is a good quantum code with the code distance O(L), it works only as self-correcting classical memory, with the bit storage time τ EXP(L). However, in order for the system to work as self-correcting quantum memory, there must be a pair of anti-commuting two-dimensional logical operators. Summary and discussion: We summarize coding properties of STS models based on dimensions of pairs of logical operators: Spatial dim Logical operators Code distance Bit storage time Qubit storage time Self-correction 1 dim 0 dim + 1 dim O(1) O(1) O(1) 2 dim 0 dim + 2 dim O(1) EXP(L) O(1) classical 2 dim 1 dim + 1 dim O(L) O(1) O(1) 3 dim 0 dim + 3 dim O(1) EXP(L) O(1) classical 3 dim 1 dim + 2 dim O(L) EXP(L) O(1) classical 4 dim 2 dim + 2 dim O(L 2 ) EXP(L) EXP(L) quantum where, for D = 4, we presented coding properties of four-dimensional Toric code. 4. Thermal stability of topological order In this section, we present another main result of the present paper, concerning the thermal stability of topological order at finite temperature by establishing the connection between the feasibility of self-correcting quantum memory and the stability of topological order. In section 4.1, we begin by establishing the connection between self-correcting classical memory and the thermal stability of ferromagnetic order at finite temperature. In section 4.2, we establish the connection between quantum codes and topological order at zero temperature. In section 4.3, we establish the connection between self-correcting quantum memory and topological order at finite temperature, and analyze the thermal stability of topological order arising in STS models. Note that discussion in this section is rather heuristic, and more rigorous treatment follows in Appendix A 4.1. Classical equivalence In this subsection, we establish the connection between self-correcting classical memory and thermal stability of ferromagnetic order at finite temperature: Classical code Ferromagnetic order at T = 0 Self-correcting classical memory Ferromagnetic order at T > 0 Two-dimensional ferromagnet: We have seen that two-dimensional Ising model works as self-correcting classical memory since the energy barrier separating two degenerate ground states is O(L). This self-correcting property of two-dimensional Ising model is closely related to the thermal stability of ferromagnetic order as seen from the following thermodynamic argument. Consider two-dimensional Ising model: H(ɛ) = Z i, j Z i+1, j Z i, j Z i, j+1 ɛ Z i, j (31) i, j i, j i, j 17

18 with an initial bias (symmetry-breaking field); ɛ i, j Z i, j for ɛ 0. As a result of an initial bias, the system is not degenerate anymore, and the ground state is 0 0. The thermal stability of ferromagnetic order can be analyzed through the expectation value of the total magnetization m: m = at finite temperature, which can be computed as follows: 1 Z i, j (32) N i, j m = lim β ɛ 0 1 log Z ɛ (33) where the partition function is: Z(β, ɛ) = Tre βh(ɛ). Here, we evaluate the expectation value of m at the limit where ɛ 0 after we take the limit where N goes to infinity. (Otherwise, there is no use of introducing an initial bias). Then, we have: 1 Z i, j = 1 (T = 0) (34) N i, j ɛ > Z i, j > 0 (T c > T > 0) (35) N 1 N i, j Z i, j i, j ɛ 0 ɛ 0 = 0 (T > T c ) (36) where T c is some finite transition temperature, as plotted in Fig. 12(a). One may notice that the total magnetization has some non-zero value as long as the temperature is below the transition temperature T c. This indicates that the system is stable against thermal fluctuations for T < T c. In particular, the system properties for T c > T > 0 are close to the ground state properties at T = 0. However, the system properties change drastically at T = T c, and for T > T c, the total magnetization vanishes. This indicates that the system undergoes a phase transition at T = T c, and is unstable against thermal fluctuations for T > T c. Figure 12: Total magnetizations. (a) Two-dimensional Ising model. (b) One-dimensional Ising model. Now, one can establish the connection between self-correcting properties and the thermal stability. In a coding theoretical language, the computation of the magnetization at finite temperature may be interpreted as follows. We first apply an initial bias and encode a bit 0 to a ground state 0. Then, we let the system interact with the external environment at finite temperature, and wait until the system reaches the equilibrium while slowly turning off the initial bias. Finally, we decode the initially encoded bit by measuring the total magnetization. Therefore, selfcorrecting property of two-dimensional Ising model with exponentially long bit storage time is a direct consequence of the thermal stability of ferromagnetic order at finite temperature. In general, the total magnetization we have used above is called an order parameter in condensed matter physics since it can be used to distinguish different phases of matters. On the other hand, quantum coding theoretically, one 18

19 may view the total magnetization as a symmetric summation of a zero-dimensional logical operator Z i, j with which a classical bit is encoded. Later, we shall see that expectation values of logical operators can be used as order parameters to characterize the stability of topological order. One-dimensional ferromagnet: Next, let us see that one-dimensional Ising model, which does not work as self-correcting classical memory, is not stable against thermal fluctuations. Consider one-dimensional Ising model: H(ɛ) = Z j Z j+1 ɛ Z j (37) with an initial bias, and compute the expectation values of the total magnetization. Then, we have 1 Z j = 1 (T = 0) (38) N j ɛ 0 1 Z j = 0 (T > 0) (39) N j ɛ 0 as plotted in Fig. 12(b). One may notice that the total magnetization becomes zero at any finite temperature. This indicates that the system is not stable against thermal fluctuations at any finite temperature, and the system undergoes a phase transition at T = 0 as seen from the sudden change of the magnetization at T = 0. The system properties for T = 0 and for T > 0 are significantly different, and the ground state properties are not stable against thermal fluctuations. This is the underlying reason why one-dimensional Ising model has O(1) bit storage time, and does not work as self-correcting classical memory. Summary of the equivalence: With observations above, one may notice that a large energy barrier, which is essential to self-correcting properties, is the key to the thermal stability of ferromagnetic order. In Fig. 13, we give a summary of the equivalence concerning classical memories. However, it is not clear whether a large energy barrier is sufficient for the thermal stability of the system or not. We address this issue briefly in Appendix A.3. Figure 13: A summary of the classical equivalence. Here, it is worth noting that ferromagnetic systems, such as one-dimensional and two-dimensional Ising model are not stable against perturbations, as described in Fig. 13. If one chooses V as the initial bias, one can easily break the ground state degeneracy, and thus, the ground state property is not stable against perturbations. This indicates that ferromagnetic systems are not topologically ordered (see section 4.2 too) Existence of topological order and quantum codes Now, we return to the main discussion in this section, concerning the stability of topological order at finite temperature. In this subsection and the next, we establish the connection between self-correcting quantum memory and 19

20 the thermal stability of topological order: Quantum code Topological order at T = 0 Self-correcting quantum memory Topological order at T > 0 In this subsection, we describe the definition of topological order in spin systems on a lattice, and argue that quantum codes have topological order at zero temperature. Note that our discussion closely follows pioneering works [44, 54]. Phenomenological definition: Let us first begin by describing the phenomenological definition of topological order which is commonly used in physics community. A system is said to have topological order when its ground state properties do not change significantly under any types of small perturbations [35]. Definition 1 (Stability against perturbations). Consider a degenerate spin system defined on some closed geometric manifold governed by a Hamiltonian H. The system is said to be topologically ordered at T = 0 if and only if the system satisfies the following conditions: There exists some finite positive number δ such that, for any perturbations V: H = H + V, V = V j (40) where V j are locally defined and V j δ, the ground state degeneracy is not broken at the thermodynamic limit (Fig. 14). The perturbed ground space G can be approximated by the original ground space G 0 through some local unitary transformation: which can be represented as j UG 0 U G (41) U = 1 0 exp( iĥt)dt (42) where ĥ is some hermitian operator which is a summation of local terms with finite amplitudes. Figure 14: The stability of the ground state properties against local perturbations. It is worth noting that a classical ferromagnet (Ising model) is not topologically ordered while it works as selfcorrecting classical memory. First of all, the ground state degeneracy of a classical ferromagnet is broken if a local magnetic field term Z j is added as a perturbation: V = ɛ Z j. (43) 20 j

21 While the original Hamiltonian has 0 0 and 1 1 as ground states, the perturbed Hamiltonian has 0 0 as a single ground state for positive ɛ. Second, the ground state property significantly changes as a result of perturbations. While a classical ferromagnet has 1 2 ( ) as a ground state, the perturbed ground state is 0 0, and there is no local unitary transformation which transforms 1 2 ( ) into 0 0 [1]. As a straightforward extension of this discussion, one may notice that STS models with zero-dimensional logical operators are not topologically ordered, and topological order may exist only in two or higher-dimensional systems due to the dimensional duality of logical operators. Coding theoretical definition: While the definition above captures the stability of ground state properties in topologically ordered systems, one may wonder what kinds of systems are actually topologically ordered. Also, one may wonder if there is a universal feature commonly shared among all the topologically ordered systems. Somewhat surprisingly, a quantum coding theoretical viewpoint gives a hint to answer these questions by providing an alternative definition of topological order which may capture universal properties of topologically ordered systems [1]. In fact, a system of spins on a lattice may be said to be topologically ordered when the ground space is separated from excited states by a finite energy gap at the thermodynamic limit and the ground space realizes a quantum code whose code distance d is comparable to the system size [54]. One may rewrite this characterization of topological order through coding properties more explicitly in the following way. 5 Definition 2 (Quantum code). Consider a degenerate spin system defined on some closed geometric manifold. The system is said to be topologically ordered when it satisfies the following condition: Let the degenerate ground state space be G, which is protected by a finite energy gap. Consider reduced density matrices of ground states for a region R which is contractable to a single point: ρ R ( ψ ) Tr R( ψ ψ ). (44) So, R is a topologically trivial zero-dimensional region (Fig. 15). Then, all the reduced density matrices of degenerate ground states are the same at the thermodynamic limit ρ R ( ψ ) = ρ R ( ψ ) for all ψ, ψ G at N. (45) Let us see why the condition above ensures that the system serves as a quantum code with a macroscopic code distance. When all the reduced density matrices are the same for all the topologically trivial regions, any errors acting inside such regions cannot change the encoded logical qubits. In order to change logical qubits, some global operator acting on a region with a topologically non-trivial geometric shape is necessary. Such a global operator has a weight equal to the code distance. Therefore, the ground space of such a system works as a good quantum code with a macroscopic code distance. 6 Equivalence of two definitions: Now, let us see that definition 2 implies definition 1. Under a small perturbation V which is a summation of local terms, we perform a perturbative analysis. Since the ground states are separated by excited states by a finite energy gap, only the coupling between degenerate ground states is dominant. Since there is no local operator to couple two ground states, one needs to consider O(L)th order perturbation to have a coupling between degenerate ground states. Since such a coupling will be exponentially suppressed by a factor exp( L/L 0 ) where L 0 is some finite length scale, and the energy splitting is exponentially suppressed and goes to zero as L goes to infinity. Thus, the ground state degeneracy is protected at the thermodynamic limit. When the ground state degeneracy is protected and the ground space is separated from excited states by a finite energy gap, one can approximate the ground space from the original unperturbed ground space according to the theory of adiabatic continuation [54, 55]. A complete mathematical proof was presented in [54]. However, it is not known if definition 1 implies definition 2 or not. 5 Strictly speaking, this definition is not complete. See [54] for a complete definition. 6 The definition used here is slightly different from the one used in [44] where the existence of the energy gap is not required. The difference comes apparent whether a two-dimensional Bacon s subsystem code [16] is classified as topologically ordered system or not. The ground state space of this subsystem code has the same reduced density matrix for any zero-dimensional regions. However, it is not likely to have a finite energy gap since it is highly frustrated, and thus, its coding properties are likely to be unstable against perturbations. 21

22 Figure 15: A coding theoretical definition of topological order. Here, we briefly mention the usage of the expression topological order since it is sometimes used in a fuzzy and elusive way. In the present paper, by topological order, we mean the ground state properties of topologically ordered systems. By its definition, topological order is stable against any types of small local perturbations, meaning that the ground state degeneracy is protected and perturbed ground states can be approximated through local unitary transformations. Therefore, by the ground state properties, we mean any global features of ground states which are free from length scales and are not affected by local unitary transformations. We will return to discussion on the meaning of global features later. Topological order in STS models: Let us conclude this subsection on the definition of topological order by seeing that two-dimensional and three-dimensional STS models can have topological order according to definition 2. Consider the case when two-dimensional STS model has k pairs of one-dimensional logical operators. Let us split the entire system into two complementary regions P(n 1 1, n 2 1) and P(n 1 1, n 2 1) where P(n 1 1, n 2 1) is a region with n 1 1 n 2 1 composite particles. Note that P(n 1 1, n 2 1) is topologically trivial zero-dimensional region since one can shrink it into a single point continuously. Since there are 2k logical operators defined inside P(n 1 1, n 2 1), we have g P(n1 1,n 2 1) = 2k. However, this means that g P(n 1 1,n 2 1) = 0 from theorem 1, and there is no logical operator defined inside P(n 1 1, n 2 1). Then, all the reduced density matrices are the same, and such a two-dimensional STS model is topologically ordered since it satisfies definition 2, and thus, definition 1. A similar discussion holds for a three-dimensional STS model with pairs of one-dimensional and two-dimensional logical operators Quantum equivalence In this subsection, we establish the connection between self-correcting quantum memory and topological order at finite temperature, and analyze the thermal stability of topological order arising in STS models. Two-dimensional Toric code: Two-dimensional Toric code is topologically ordered since it has a macroscopic code distance: d O(1). Yet, it does not work as self-correcting quantum memory since the energy barrier is E O(1). This is closely related to the thermal instability of topological order as we shall see below. Let us perform 22

23 a thermodynamic analysis on two-dimensional Toric code in a way similar to Ising model by adding an initial bias: H l (ɛ) = H Toric ɛ H r (ɛ) = H Toric ɛ L j=1 L i=1 T j 2 (l) (46) T1 i (r) (47) where l and r are one-dimensional logical operators extending in the ˆ1 and ˆ2 directions respectively. Here, we define the following normalized logical operators : m l m r = 1 L = 1 L L j=1 L i=1 T j 2 (l) (48) T1 i (r) (49) by taking symmetric summations of logical operators, as we did for Ising model. Then, we have m l ɛ 0 = m r ɛ 0 = 1 (T = 0) (50) m l ɛ 0 = m r ɛ 0 = 0 (T > 0) (51) as plotted in Fig. 16(a) where m l is evaluated for H l (ɛ), and m r is evaluated for H r (ɛ). This indicates that the system is not stable against thermal fluctuations at any finite temperature, meaning that topological order arising in two-dimensional Toric code is thermally unstable. This implies that one cannot read out initially encoded qubit by measuring m l and m r. Three-dimensional Toric code: Next, let us consider thermodynamic properties of three-dimensional Toric code. While three-dimensional Toric code does not work as self-correcting quantum memory, it works as self-correcting classical memory. This is because it has pairs of one-dimensional and two-dimensional logical operators, and as a result, the bit storage time is τ EXP(L) while the qubit storage time is τ O(1). These coding properties are closely related to thermodynamic properties of three-dimensional Toric code, as seen from expectation values of logical operators: m l = 1 L L z=1 m r = 1 L 2 L x,y=1 T z 3 (l) (52) T1 x T y 2 (r) (53) where l is a two-dimensional logical operator extending in the ˆ1 and ˆ2 directions, while r is a one-dimensional logical operator extending in the ˆ3 direction. Their expectation values behave as follows: and m l ɛ 0 = 1 (T = 0) (54) m l ɛ 0 = 0 (T > 0) (55) m r ɛ 0 = 1 (T = 0) (56) 1 > m r ɛ 0 > 0 (T c > T > 0) (57) m r ɛ 0 = 0 (T > T c ) (58) 23

24 where T c is some finite transition temperature, as plotted in Fig. 16(b). This implies that three-dimensional Toric code undergoes phase transitions both at T = 0 and T = T c, and the ground state properties are not completely stable against thermal fluctuations at any finite temperature. Yet, the ground state properties partially survive at finite temperature as a direct consequence of being self-correcting classical memory. Figure 16: Expectation values of logical operators. (a) Two-dimensional Toric code. (b) Three-dimensional Toric code. (c) Four-dimensional Toric code. Four-dimensional Toric code: Finally, let us see that four-dimensional Toric code, which works as self-correcting quantum memory, is stable against thermal fluctuations. Here, we define m l = 1 L 2 L z,w=1 m r = 1 L 2 L x,y=1 T z 3 T w 4 (l) (59) T1 x T y 2 (r) (60) where l is a two-dimensional logical operator extending in the ˆ1 and ˆ2 directions, and r is a two-dimensional logical operator extending in the ˆ3 and ˆ4 directions. Then, we have m l ɛ 0 = m r ɛ 0 = 1 (T = 0) (61) 1 > m l ɛ 0 = m r ɛ 0 > 0 (T c > T > 0) (62) m l ɛ 0 = m r ɛ 0 = 0 (T > T c ) (63) where T c is some finite transition temperature, as plotted in Fig. 16(c). This implies that the ground state properties are stable against thermal fluctuations and topological order arising in four-dimensional Toric code is stable at finite temperature. Summary of the equivalence: With these observations, one may notice that large energy barrier inside the ground space, which is essential to self-correcting properties, is the key to the thermal stability of topological order. In Fig. 17, we give a summary of the equivalence concerning quantum memory. With this connection between the feasibility of self-correcting quantum memory and the thermal stability of topological order, we conclude that topological order arising in STS models is not stable against thermal fluctuations. In 24

25 Figure 17: The quantum equivalence. other words, for D 3, there is no system which is stable against both thermal fluctuations and local perturbations simultaneously. While discussion here is rather heuristic, we give a more rigorous treatment on the definition of the thermal stability of topological order in Appendix A. We summarize physical properties of STS models based on dimensions of pairs of logical operators: Spatial dim Logical operators Local perturbations Thermal fluctuations Memory property 1 dim 0 dim + 1 dim Classical code 2 dim 0 dim + 2 dim stable Classical self-correction 2 dim 1 dim + 1 dim stable Quantum code 3 dim 0 dim + 3 dim stable Classical self-correction 3 dim 1 dim + 2 dim stable Quantum code 4 dim 2 dim + 2 dim stable stable Quantum self-correction where, for D = 4, we presented coding properties of four-dimensional Toric code. 5. Emergence of topology in logical operators So far, we have addressed specific questions concerning coding and physical properties of gapped spin systems on a lattice. While these two questions are of particular importance in quantum information science and condensed matter physics, the ultimate goal is to find universal properties of arbitrary gapped spin systems on a lattice and develop a unified theoretical framework to discuss their coding and physical properties. In particular, it would be very beautiful if arbitrary gapped spin systems can be analyzed universally through some unified theoretical tool which is yet to be discovered. In this section, as a necessary first step toward this goal, we make an attempt to find universal properties of logical operators which are commonly shared among all the STS models. Logical operators in STS models have a certain interesting topological property which is a direct consequence of physical constraints we have imposed on stabilizer Hamiltonians. Here, we show that one can always deform geometric shapes of logical operators continuously while keeping them equivalent by applying appropriate stabilizers as long as one does not change geometric shapes in a topologically non-trivial way. We prove this continuous deformability of logical operators in STS models for D = 1, 2, 3 with various examples of continuous deformations of logical operators. Some of these results were also presented in our previous papers [17, 47], especially for D = 2. The program of finding a universal theory of gapped spin systems is continued in Appendix B where the role of topology in analyzing coding and physical properties of stabilizer codes is further examined, while we concentrate on demonstrating topological properties of logical operators in this section. One-dimension: We begin by illustrating a topological property of logical operators for one-dimensional STS models. Recall that P(x) represents a region of x composite particles, and g P(1) = k in one-dimensional systems regardless of the system size since there always exist k zero-dimensional logical operators. Then, as a result of the 25

26 bi-partition theorem (theorem 1), we have g P(1) = g P(1) = k. (64) Figure 18: A deformation of a zero-dimensional logical operator. Now, let us interpret the equation above from a geometric viewpoint. Consider some logical operator l defined inside P(1). Then, since g P(1) = g P(1) = k, there must exist some equivalent logical operator l l which is defined inside P(1). In other words, if there exists a logical operator defined inside P(1), one can shrink its geometric shape into P(1) by applying some stabilizer (Fig. 18): P(1) P(1). (65) Here, one may notice that P(1) and P(1) are topologically equivalent regions since one can deform P(1) into P(1) continuously. On the other hand, if a logical operator is a one-dimensional logical operator defined all over the lattice (i.e. defined inside P(n) where n is the linear length of the system), one may not be able to deform the logical operator into P(1). Two-dimensions: Next, we analyze a topological property of logical operators for two-dimensional STS models. For the convenience of presentation, we assume that zero-dimensional logical operators defined inside P(1, 2v) can be actually defined inside P(1, 1) in a two-dimensional STS model. This may be done through some appropriate coarse-graining. We first list regions which serve as references to classify geometric shapes of logical operators in a two-dimensional system. We define topological unit regions as follows (Fig. 19(a)): Q(0, 0) P(1, 1), Q(1, 0) P(n 1, 1), Q(0, 1) P(1, n 2 ), Q(1, 1) P(n 1, n 2 ). (66) 1 and 0 represent whether a region extends in the corresponding direction or not. For example, Q(1, 0) and Q(0, 1) are one-dimensional unit regions which extend in the directions of ˆ1 and ˆ2 respectively. Q(0, 0) is a zerodimensional unit region with a single composite particle. Q(1, 1) is a two-dimensional unit region which consists of all the composite particles in the system. These topological unit regions are shown graphically in Fig. 19(a). We also denote a union of all the m-dimensional topological unit regions as R m : R 0 Q(0, 0), R 1 Q(1, 0) Q(0, 1), R 2 Q(1, 1). (67) We call R m m-dimensional concatenated topological unit regions. All the concatenated unit regions are graphically shown in Fig. 19(b). In a two-dimensional system, there are five different unions of topological unit regions: R 0, Q(1, 0), Q(0, 1), R 1 and R 2. We call these regions, except R 2, reference regions, whose set is denoted as R re f : R re f = {R 0, Q(1, 0), Q(0, 1), R 1 }. (68) Then, one can introduce equivalence relations between these reference regions and their complements in terms of continuous deformations. For example, as shown in Fig. 20(a), R 0 can be continuously deformed into R 1 by enlarging the hole of R 0 gradually. Also, as shown in Fig. 20(b), R 1 can be continuously deformed into R 0 since both R 1 and R 0 are zero-dimensional regions without any winding around the torus. Finally, as shown in Fig. 20(c), Q(1, 0) can be deformed into Q(1, 0) since both regions have a winding in the ˆ1 direction. In summary, we have the following 26

27 Figure 19: Reference regions. (a) Topological unit regions. (b) Concatenated unit regions. equivalence relations between reference regions and their complements: R 0 R 1, R 1 R 0, Q(1, 0) Q(1, 0), Q(0, 1) Q(0, 1). (69) Figure 20: The topological deformations of logical operators. Now, we discuss how geometric shapes of logical operators can be determined. A useful observation regarding geometric shapes of logical operators can be obtained by considering the number of independent logical operators defined inside a region R. Let the number of independent logical operators inside R be g R. Here, we consider the case where we have two regions R and R where R is larger than R, meaning that R includes all the composite particles inside R (Fig. 21). Then, if g R = g R, R and R support the same logical operators since all the logical operators defined inside R have equivalent representations which are supported inside R. This means that, for a given logical operators l defined inside R, one can always find another equivalent logical operator l defined inside R. In other words, one can deform the geometric shape of l into l by applying some appropriate stabilizer (Fig. 21). Thus, by finding two connected regions R and R where R is larger than R and g R = g R, one can conclude that logical operators defined inside R can be deformed into R. With the above observation on geometric shapes of logical operators and deformations in mind, let us describe 27

28 Figure 21: A shrinkage from R to R when g R = g R. a topological property of geometric shapes of logical operators. The following lemma summarizes the topological property of logical operators in two-dimensional STS models. Lemma 1 (Topological shrinkage). In two-dimensional STS models, the following equations hold: g R0 = g R1, g R1 = g R0, g Q(1,0) = g Q(1,0) = k, g Q(0,1) = g Q(0,1) = k. (70) In other words, one can shrink geometric shapes of logical operators by applying some appropriate stabilizers in the following ways: R 0 R 1, R 1 R 0, Q(1, 0) Q(1, 0), Q(0, 1) Q(0, 1). (71) Note that these shrinkages preserve topological properties of geometric shapes of logical operators. Proof. Let us begin with the proof of g R0 = g R1. Let k 0 be the numbers of pairs of anti-commuting zero-dimensional and two-dimensional logical operators. Let k 1 be the number of pairs of anti-commuting one-dimensional logical operators. Then, we notice that g R1 = 2k 1 + k 0 since R 1 supports both zero-dimensional and one-dimensional logical operators. On the other hand, since g R0 = 2k g R0 from theorem 1 and g R0 = k 0, we have g R0 = g R1. If we use theorem 1 to g R0 = g R1, we readily obtain g R0 = g R1 too. Next, let us show g Q(1) = g Q(1). Note that g Q(1) = k. Then, we have g Q(1) = k, and have g Q(1) = g Q(1). Three-dimensions: Let us continue our analysis on STS models for higher-dimensional cases (D > 2). Below, we analyze the topological property of logical operators in three-dimensional STS models. We assume that zerodimensional logical operators and one-dimensional logical operators in a three-dimensional STS model can be defined inside P(1, 1, 1), P(n 1, 1, 1), P(1, n 2, 1) and P(1, 1, n 3 ). We begin by finding reference regions for D-dimensional systems. Reference regions for D > 2 can be defined from topological unit regions in a way similar to two-dimensional cases. Let d be an arbitrary binary D component vector d = (d 1,, d D ) with d m = 0, 1. Then, topological unit regions are: Q( d) P( x), where x m = n d m m. (72) For example, Q(1, 1, 0) = P(n 1, n 2, 1) for D = 3. We denote the weight of the binary vector d as w( d) D m=1 d m, which represents the dimension of Q( d). Then, concatenated unit regions are defined as follows: R m Q( d). (73) w( d)=m A set of reference regions can be obtained by considering all the possible unions of Q( d), which is denoted as R re f. One may notice that topological unit regions Q( d) are like independent generators for mth homology group for a D-torus: H m (T D ). 28

29 Figure 22: Topological unit regions in a three-dimensional system. Recall that we set the periodic boundary conditions. It is worth presenting some examples here. topological unit regions. In a three-dimensional system (D = 3), we have the following 0 dim: Q(0, 0, 0) 1 dim: Q(1, 0, 0), Q(0, 1, 0), Q(0, 0, 1) 2 dim: Q(1, 1, 0), Q(0, 1, 1), Q(1, 0, 1) 3 dim: Q(1, 1, 1). (74) Some examples are shown in Fig. 22. Also, concatenated topological unit regions are: R 0 = Q(0, 0, 0) R 1 = Q(1, 0, 0) Q(0, 1, 0) Q(0, 0, 1) R 2 = Q(1, 1, 0) Q(0, 1, 1) Q(1, 0, 1) R 3 = Q(1, 1, 1) (75) which are described in Fig. 23. One can introduce the equivalence relations in terms of reference regions. Equivalence relations among them are 29

30 Figure 23: Concatenated unit regions in a three-dimensional system. shown as follows: R 0 R 2 Q(1, 0, 0) Q(1, 1, 0) Q(1, 0, 1) Q(0, 1, 0) Q(1, 1, 0) Q(0, 1, 1) Q(0, 0, 1) Q(1, 0, 1) Q(0, 1, 1) Q(1, 0, 0) Q(0, 1, 0) Q(1, 1, 0) Q(0, 0, 1) Q(0, 1, 0) Q(0, 0, 1) Q(0, 1, 1) Q(1, 0, 0) Q(0, 0, 1) Q(1, 0, 0) Q(1, 0, 1) Q(0, 1, 0) (76) R 1 R 1 Q(1, 1, 0) Q(1, 1, 0) Q(0, 1, 1) Q(0, 1, 1) Q(1, 0, 1) Q(1, 0, 1). Then, we have the following theorem. Lemma 2 (Topological shrinkage). For D-dimensional STS models (D = 1, 2, 3), let R and R be reference regions: R, R R re f. When R R, one can deform geometric shapes of logical operators continuously from R to R: g R = g R for R R. (77) Proof. The proof of the lemma is straightforward from theorem 1 and theorem 5, so we present a proof only for R 1 R 1. Since all the zero-dimensional and one-dimensional logical operators can be supported inside R 1, we have g R1 k. Similarly, we have g R1 k. Then, since g R1 + g R1 = 2k, we have g R1 = g R1 = k. Topological deformations of logical operators: So far, our discussion on a topological property of logical operators has been limited only to regions generated by m-dimensional unit regions while some connected regions (such as examples presented in Fig. 24) cannot be generated by taking unions of unit regions. Then, a naturally arising 30

31 question is whether one can perform similar continuous deformations for arbitrary connected regions of composite particles or not. Here, we provide a complete description on topological properties of logical operators in STS models by extending the notion of continuous deformations to arbitrary connected regions. It turns out that one can continuously deform a geometric shapes of a logical operator defined on any connected region of composite particles continuously for STS models with D = 1, 2, 3. Here, we summarize topological properties of logical operators in STS models as follows 7 : Theorem 6 (Continuous deformation). Consider STS models for D = 1, 2, 3. Consider two connected regions of composite particles R and R which are topologically equivalent: R R. Then, we have g R = g R for R R. (78) Therefore, for any given logical operator l defined inside R, one can always find an equivalent logical operator defined inside R as long as R can be continuously deformed into R. Figure 24: Examples of continuous deformations. Note that systems have periodic boundary conditions. For clarity of presentation, we skip the proof of theorem 6. One may verify it by using scale symmetries, a bipartition theorem (theorem 1), and the translation equivalence of logical operators (theorem 2). We present some examples of continuous deformations in Fig Summary and open questions In this paper, we have established the connection between the feasibility of self-correcting quantum memory and the thermal stability of topological order, and provided partial answers to these two open problems by solving a model that may cover a large class of physically realizable quantum codes. Our discussion is limited to stabilizer codes with translation and scale symmetries, and these two problems still remain open for an even larger class of gapped spin systems on a lattice. Yet, we hope that our analysis will provide an important insight and a useful guidance on studies 7 Mathematically inclined readers may want to have more precise definitions of geometric shapes and continuous deformations. Here, we make some comments on this issue to make our discussion more rigorous. For a D-torus with n 1 n D composite particles, we split the entire system into n 1 n D hypercubic regions where each hypercube contains each composite particle, and cover the entire system completely, but without any overlap. For a given set of composite particles R, we define its geometric shape as a union of hypercubes which contain all the composite particles inside R. After assigning geometric shapes to sets of composite particles in this way, one can define a continuous map between two sets R and R, and introduce an equivalence relation between them in a straightforward way. 31

32 of coding and physical properties of gapped spin systems on a lattice. We also hope that our work will contribute to stimulating the use of quantum coding theoretical concepts in studying many-body quantum systems further. Below, we discuss possible future problems and make some comments on them. Frustration-free systems with scale symmetries: While our discussion is limited to stabilizer codes, any gapped spin systems with degenerate ground states can be used as quantum memory devices in principle. It seems that, for an arbitrary gapped spin system defined on a lattice, there exists some frustration-free Hamiltonian which approximates the original system and serve as its low energy effective Hamiltonian (except some subtle properties such as chirality). In fact, this claim has been rigorously proven for one-dimensional gapped spin systems [56] by showing that any ground state of one-dimensional gapped spin Hamiltonians can be efficiently simulated through the matrix product state formalism. Now, this claim seems to be widely believed among the condensed matter physics community, and it is at the heart of recent progress in classifications of quantum phases in gapped spin systems [17, 57, 58]. Therefore, the analyses on coding properties of arbitrary frustration-free Hamiltonians with translation and scale symmetries may provide a useful insight on questions concerning coding properties and quantum phases arising in arbitrary gapped spin systems on a lattice. To the best of our knowledge, all the examples of frustration-free Hamiltonians with translation and scale symmetries, such as the quantum double model [8] and the string-net model [59], have the dimensional duality on geometric shapes of certain operators which may be considered as generalizations of logical operators. Also, these operators seem to have the continuous deformability in a way similar to logical operators in STS models. With these observations, we feel that our results are universal for all the frustration-free Hamiltonians with translation and scale symmetries, and thus, effectively true for arbitrary gapped spin systems defined on a lattice with a small number of ground states. Also, these observations imply that any gapped spin systems with scale symmetries (or a small number of ground states) may be described through TQFT, as discussed in Appendix B. Yet, the connection between TQFT and coding properties of gapped spin systems on a lattice with scale symmetries must be further established. Beyond scale symmetries: While our treatments in the present paper are limited to quantum codes with scale symmetries, there are several interesting models of quantum codes which do not have scale symmetries [32, 60 62]. Of particular interests are models proposed in [61] and in [32] which are classical and quantum memories respectively with partially self-correcting properties. These models do not have scale symmetries since the number of logical qubits k n is highly sensitive to the system size n, and there is no constant upper bound on k n. These models are known to have logarithmically large energy barrier E LOG(L) with a large number of energy local minima. As a result, these models seem to have POLY(L) relaxation time with slow dynamics, which may result in POLY(L) bit or qubit storage time. Whether these models with broken scale symmetries are useful as classical and quantum memories is a complicated problem. First of all, for a system to be a efficient quantum memory device, it is desirable to have exponentially long qubit storage time: τ EXP(L), since it needs at least d gate operations to write or readout a logical qubit, and these encoding and decoding processes takes at least polynomially long time τ POLY(L). Also, it seems difficult to find an efficient decoding algorithm with an efficient error-correction for these models. Finally, these models undergo phase transitions at T = 0, which may imply their potential thermal instability. At this moment, there are many issues to be analyzed concerning implementability of quantum codes with broken scales symmetries from engineering viewpoints. Beyond TQFT: Stabilizer codes with broken scale symmetries are remarkable examples which may be beyond the description of the standard TQFT. Such models do not have continuously deformable logical operators, and are not expected to be characterized by topological properties of logical operators as in STS models. This observation is consistent with the fact that systems described by TQFT, in a sense of the axiomatic formulation developed by Atiyah, are allowed to have only a finite number of degenerate ground states. Thus, finding an effective field theoretical description for stabilizer codes with broken scale symmetries may be an interesting future problem which may lead to discoveries of novel quantum phases that are currently unknown and are governed by some deeper mathematical formalism than topology. To the best of our knowledge, all the examples of stabilizer codes with broken scale symmetries still possess certain discrete scale symmetries. STS models have continuous scale symmetries since k n = k for all the system sizes n, while stabilizer codes with broken scale symmetries seem to have discrete scale symmetries since the number of logical qubits k n has good scaling properties under global scale transformations: (n 1, n 2, ) (cn 1, cn 2, ) where c is some appropriate integer with behaviors such as k c n = ck n. The distinction between continuous and discrete 32

33 scale symmetries becomes particularly important when one considers effects of RG transformations. As is pointed out in [63], systems with continuous scale symmetries correspond to fixed points of RG transformations, while systems with discrete scale symmetries correspond to limit cycles of RG transformations. Physical realizations of systems with discrete scale symmetries include the Effimov effect which is recently of particular interest in the ultracold atom physics community. This observation implies that stabilizer codes with discrete scale symmetries correspond to limit cycles of RG transformations. Yet, the connection between limit cycles and discrete scale symmetries for lattice systems must be further established since discussion in [63] is given primary for continuum systems. Spin glass behaviors and translation symmetry breaking: Another interesting property of stabilizer codes with broken scale symmetries is a glassy behavior with slow relaxation dynamics. This glassy behavior may be understood as a direct consequence of broken scale symmetries, and associated broken translation symmetries in the ground space. In general, even if a Hamiltonian possesses translation symmetries, its ground states may not possess translation symmetries. Let us assume that a stabilizer Hamiltonian is translation invariant under unit translations of composite particles: T m (H) = H for all m. Then, ground states of this stabilizer Hamiltonian may break translation symmetries: T m ( ψ ) ψ. Some examples of systems with broken translation symmetries are presented in [17, 38, 64], and translation symmetry breaking of ground states arising in topologically ordered systems are studied in [38, 64]. An interesting connection between scale symmetries and translation symmetries is that when a stabilizer code has scale symmetries: k n = k, all the ground states are invariant under unit translations of composite particles: T m ( ψ ) = ψ regardless of the system size, as proven in [17]. In other words, if a system is coarse-grained such that scale symmetries are satisfied, translation symmetries are protected inside the ground space. On the other hand, if a system does not have scale symmetries, translation symmetries may be broken inside the ground space. In particular, when the system has a large number of logical qubits, translation symmetries are strongly broken such that there is no finite translation which keep ground states equivalent. This strong breaking of translation symmetries seems to be the key to the glassy behavior observed in stabilizer codes with broken scale symmetries. In conventional spin glass models, a Hamiltonian consists of mutually commuting terms whose signs are random, and translations symmetries of the Hamiltonian are initially broken due to the randomness. This randomness in the Hamiltonian gives rise to spin configurations of ground states which are not uniform over real space, and leads to a complicated and slow relaxation dynamics in spin glasses. On the other hand, in stabilizer codes with broken scale symmetries, stabilizer Hamiltonians possess translation symmetries initially. However, spin configurations are not uniform since translation symmetries are strongly broken inside the ground space due to broken scale symmetries. Thus, systems with broken scale symmetries are expected to exhibit spin glass behaviors with slow relaxation dynamics and a large number of local minima. Yet, the connection between the glassy behavior and broken scale symmetries must be further established since the observation given here is very heuristic. Appendix A. Topological order at finite temperature In the main part of the paper, we have discussed the thermal stability of topological order by using expectation values of logical operators. In this appendix, we verify the use of expectation values of logical operators as topological order parameters, and define the thermal stability of topological order more rigorously. In Appendix A.1, we begin by verifying that expectation values of logical operators can be used for topological order parameter at T = 0. In Appendix A.2, we define topological order at finite temperature by extending the definition of topological order at zero temperature to finite temperature. In Appendix A.3, we discuss whether the existence of a large energy barrier is sufficient for the thermal stability of the system or not. Appendix A.1. Topological order parameters and logical operators While the ground state properties change only slightly under small perturbations in topologically ordered systems, the original ground state properties will be eventually lost under large perturbations. One may see whether topological order is lost by checking whether a new ground state can be approximated by the original ground state through a local unitary transformation. However, this naive approach will require a substantial amount of computations and cannot capture changes of coding and physical properties of the ground state under perturbations. Fortunately, the loss of the ground state properties under perturbations can be quantitatively characterized by some physical quantities, called topological order parameters. Here, we characterize the stability of topological order against perturbations through topological order parameters. 33

34 Topological order parameter: Topological order parameters are physical quantities with which one can judge if topological order is lost under perturbations or not. In this light, one may deduce the necessary properties of topological order parameters. By its definition, a topological order parameter must be a global function of the system which does not change under local unitary transformations. Also, a topological order parameter must change only slightly when topological order is protected, while it undergoes some non-analytic change when topological order is lost as a result of perturbations. There have been several proposals for such topological order parameters. The number of ground states serves as a topological order parameter since the ground state degeneracy is protected for topologically ordered systems at the thermodynamic limit [35]. Another interesting proposal is to use a certain entanglement measure, called topological entanglement entropy [6, 7], which corresponds to a constant correction to the entanglement area law. Since topological entanglement entropy is a non-local quantity which involves a large number of spins, it does not change under local unitary transformations. Below, we analyze the stability of topological order thorough expectation values of logical operators in order to further build the connection between topological order and quantum codes. Expectation values of logical operators: Let us first consider two-dimensional Toric code and denote a pair of anti-commuting logical operators as l and r where l extends in the ˆ1 direction and r extends in the ˆ2 direction. Then, one may use the expectation values of the summations of translations of logical operators as order parameters: U(l) = 1 T y L 2 (l), y U(r) = 1 T1 x L (r) where L is the linear length of the system. Here, we took the summation of logical operators in a symmetric way. Now, let us discuss if topological order is stable under a perturbation V by considering the following perturbed Hamiltonian with an initial bias: x (A.1) H l (ɛ) = H S TS ɛu(l) + V. (A.2) Note that, when V = 0, the small bias breaks the ground state degeneracy and U(l) ɛ 0 = 1. (A.3) Recall that, in evaluating the expectation value of U(l) in a presence of V, one must take the limit where ɛ goes to zero after taking the limit where L goes to infinity. Since the ground space is separated from excited states by a finite energy gap, one may neglect the effect of excited states. So, let us examine the effect of ɛu(l) and V on the ground space. As a result of ɛu(l), the ground state space is split into two sectors with l = 1 and l = 1 where the energy splitting between two sectors is ɛ. On the other hand, the energy splitting induced by a perturbation V is exponentially suppressed by a factor exp( L/L 0 ). Then, in considering the effect of V on a sector with l = 1, one may neglect the effect of a sector with l = 1 when ɛ is O(1). Now, we take the limit of L. Then, the ground space corresponds to a sector with l = 1, and we have U(l) ɛ 0 > 0 (A.4) even when ɛ goes to zero. However, if V is large enough, the ground state properties will be lost and U(l) ɛ 0 will be close to the value for the ground state of V. In a similar way, one may consider the following initial bias: H r (ɛ) = H S TS ɛu(r) + V (A.5) and have U(r) ɛ 0 > 0. (A.6) Therefore, the stability of topological order can be characterized by expectation values of U(l) and U(r). This argument implies that the loss of the ground state properties may be captured through the stability of logical 34

35 qubits encoded in the ground space. The initial bias of ɛu(l) represents the initial encoding of a logical qubit. The expectation value of U(l) under a perturbation V represents how much information is protected in a presence of perturbations. Stability of topological order: Based on these observations, let us formulate the stability of topological order in terms of logical operators. Let us first generalize the definition of logical operators. We call operators U which satisfy the following conditions generalized logical operators: [H 0, U] = 0 (A.7) U e iθ I and U max = 1 in the ground space (A.8) where the action of U is defined inside the ground space of the original unperturbed Hamiltonian H 0. In other words, generalized logical operators are any operators which does not change the energy of the system, but acts non-trivially inside the ground space. Generalized logical operators are not c-numbers inside the ground space. For example, a scaled summation of translations of logical operators is a generalized logical operator, while projectors onto the ground space are not generalized logical operators. Based on generalized logical operators, we define the stability of topological order in the following way. Definition 3 (Stability of topological order). The system has topological order when there exists a pair of generalized logical operators U(l) and U(r) which satisfy the following conditions: U(l) and U(r) do not commute with each other: [U(l), U(r)] 0 inside the ground space. This implies the existence of the ground state degeneracy, and U(l) and U(r) may characterize a logical qubit or qudit. For perturbed Hamiltonians: with any types of local perturbations V, H l = H 0 ɛu(l) + V, H r = H 0 ɛu(r) + V (A.9) U(l) ɛ 0, V 0 = 1 for H l (A.10) U(r) ɛ 0, V 0 = 1 for H r. (A.11) In evaluating U(l) ɛ 0,V 0, we first take the thermodynamic limit, and then, take the limit of V 0, and finally, take the limit of ɛ 0. It should be emphasized that we consider any types of perturbations V in analyzing the stability of topological order 8. Let us look at some examples here. For a one-dimensional classical ferromagnet, let l and r be a zero-dimensional and a one-dimensional logical operators: l = Z 1, r = X j. (A.12) j Then, for V = t j Z j, we have H l = Z j Z j+1 ɛ ( 1 ) Z j + t L j j = Z j Z j+1 + ( t ɛ ) Z j. L j Z j (A.13) At the thermodynamic limit, we have U(l) = 1 for H l (A.14) 8 The definition used here is similar to the one used in [44] except for taking the limit of V 0. This is because the definition in [44] is applicable to systems with gapless energy spectrum, while our definition is applicable only to systems with a finite energy gap. 35

36 regardless of t > 0 and ɛ > 0. Then, taking the limit of t 0 and ɛ 0, we have U(l) ɛ 0, V 0 = 1 for H l. (A.15) Therefore, the system is not topologically ordered. Topological order parameters and local unitary: Strictly speaking, expectation values of logical operators cannot be used as topological order parameters. In fact, while logical operators are non-locally defined for topologically ordered systems, their expectation values may change under local unitary transformations. However, topological order must be characterized by some physical quantities or objects which are not affected by local unitary transformations. In fact, both the ground state degeneracy and topological entanglement entropy do not change under local unitary transformations. In this light, one might think that the above definition of the stability of topological order is not legitimate. Here, we make some comments on this issue briefly. Despite the fact that expectation values of logical operators are not topological order parameters, we can formulate the stability of topological order through them. This is because we considered not only a specific type of local perturbations but all the types of local perturbations. Here, we demonstrate that, if we consider only one type of perturbations, expectation values of logical operators may fail to capture quantum phase transitions. In two-dimensional Toric code, let us consider a one-dimensional logical operator l which consists only of Pauli X operators. Then, consider the following perturbation: V = t X r, (A.16) where t is some positive parameter and r represents the position of each qubit. Then, when we increase t, there must be a quantum phase transition and the ground state properties will be close to the ones for V. However, such a perturbation may not change the value of U(l) since V commutes with the logical operator l. In general, if one carefully choose the types of perturbations, expectation values of logical operators remain unchanged. This implies that if one tries to study a quantum phase transition only through U(l), one may fail to detect the transition. 9 Appendix A.2. Topological order at finite temperature So far, we have focused on the ground state properties of topologically ordered systems at zero temperature and seen that topological order characterizes the ground state properties which are stable against any types of small perturbations. Here, we address the thermal stability of topological order. The stability of topological order against perturbations can be discussed by seeing if the perturbed ground state can be reached by local unitary transformations or not. However, the stability of topological order against thermal fluctuations cannot be formulated in a similar way since there is no unitary transformation which connects a pure state and a statistical ensemble. Here, we define the stability of topological order against thermal fluctuations by the changes of topological order parameters. Stability against thermal fluctuations: Below, we give the definition of the stability of topological order through expectation values of logical operators. Definition 4 (Topological order at finite temperature). The system is said to have topological order which is stable at finite temperature if and only if there exists a pair of non-commuting generalized logical operators U(l) and U(r), and some finite transition temperature T c which satisfy the following conditions: Consider perturbed Hamiltonians with initial biases r H l = H 0 ɛu(l) + V, H r = H 0 ɛu(r) + V. (A.17) 9 While expectation values of logical operators cannot be used as topological order parameters, geometric shapes of dressed logical operators [17], which are operators characterizing the transformations inside the ground space of the perturbed Hamiltonians [55], may be used. Since continuous deformations between topologically distinct logical operators are not allowed, geometric shapes of dressed logical operators must undergo some discontinuous change during the transition. Thus, some non-analyticities will be induced in ground states during quantum phase transitions between Hamiltonians with topologically distinct logical operators. 36

37 For T = 0, U(l) ɛ 0, V 0 = 1 for H l (A.18) U(r) ɛ 0, V 0 = 1 for H r, (A.19) and for T c > T > 0, U(l) ɛ 0, V 0 > 0 for H l (A.20) U(r) ɛ 0, V 0 > 0 for H r (A.21) for any types of V. One may see the connection between self-correcting quantum memory and the stability of topological order. In particular, if topological order in a spin system is stable at finite temperature, such a system works as self-correcting quantum memory. A logical qubit encoded with respect to U(l) and U(r) can be read out by measuring U(l) and U(r) even at finite temperature. This implies that encoded logical qubits are not lost in a presence of the interaction with the external environment, and there must be some self-correcting thermal dissipation processes. It should be noted that the expectation values of logical operators l and r usually vanishes at any finite temperature as one may see from a direct calculation. Only the symmetric summations of logical operators may give rise to nonvanishing expectation values of generalized logical operators. Comment on the definition of stability: Our definition of the stability of topological order relies on topological order parameters associated with a pair of non-commuting logical operators l and r. The definition of the stability of topological order in the present paper is motivated purely from quantum information theoretical viewpoints. However, the definition we have used in the present paper may be too strict. For example, if one is interested only in the ground state properties associated to the expectation value of l, it makes perfect sense to claim that a three-dimensional STS model has topological order which is stable at finite temperature. In fact, once the system is held at finite temperature, the system properties are stable against any types of small perturbations. Similarly, the stability of the ground state properties depends crucially on types of perturbations. For example, while there is no topological order in a one-dimensional system in a strict sense, the ground state properties of a onedimensional classical ferromagnet are stale against small perturbations if one declares to consider only perturbations which are products of X operators. Similarly, the ground state properties of the one-dimensional AKLT model is stable against small perturbations which do not break the time-reversal symmetry. Therefore, in discussing the stability of topological order, one needs to specify types of symmetries which are of interest. Appendix A.3. Thermal stability of topological order and energy barrier We have seen that a large energy barrier in classical or quantum memory is the key to the thermal stability of ferromagnetic order or topological order at finite temperature. Yet, it is not clear if the existence of a large energy barrier is sufficient for the thermal stability or not. While STS models always have energy barrier O(L a ) where a is an integer (a 0) due to the dimensional duality of logical operators, there are several examples of stabilizer codes with broken scale symmetries which have logarithmic energy barrier: E LOG(L). Let us look at an example of classical memory with broken scale symmetries [61]. The model is constructed on a square lattice with L L qubits: H = Z i, j Z i+1, j Z i, j+1. (A.22) i, j The model does not have scale symmetries, and there is no upper bound on the number of logical qubits (bits) k since k is highly sensitive to L. The model has logarithmic energy barrier E LOG(L), and as a result, the bit memory time is expected to be τ POLY(L) if one trusts the Arrhenius law. Despite a large energy barrier which scales with the system size, the model is known to be thermally unstable. This may be easily verified by computing the partition function: Z(β) = Tre βh. Since N k stabilizers are independent, 37

38 we have (e β + e β ) N k (e β ) k Z(β) (e β + e β ) N k (e β ) k (A.23) and, we have 1 lim N N log Z(β) = log(e β + e β ). (A.24) Therefore, the thermodynamic property of this model is equivalent to a single qubit in a magnetic field, and the model does not have the thermal stability. A similar discussion holds for stabilizer codes with N qubits, supported by N interaction terms, which include one-dimensional Ising model, two-dimensional Toric code, models proposed in [32, 60]. Note that M > N in two-dimensional Ising model and four-dimensional Toric code where M is the number of interaction terms. It is not clear whether polynomial energy barrier is necessary for thermal stability of ferromagnetic order or not. Another important characteristic of models with logarithmic energy barrier is the existence of a logarithmically large number of local minima. As a result, at finite temperature, the entropic term easily dominates the free energy function. So, it seems that the thermal instability results from both the logarithmic energy barrier and the logarithmically large number of local minima. It is worth mentioning that there is an interacting spin model with a logarithmic energy barrier, but is not defined on a lattice. The model is called two-dimensional XY model, which is known to undergo a thermal phase transition at finite temperature. The model is of particular interest since it was not expected to possess any thermal phase transitions as a result of the Mermin-Wagner theorem which states that continuous symmetries cannot be spontaneously broken at finite temperature for D 2. The reason why two-dimensional XY model undergoes a thermal phase transition at finite temperature is because the transition, known as Kosterlitz-Thouless transition, is induced by topological defects which are not directly related to continuous symmetries of the model. In two-dimensional XY model, it takes logarithmic energy to create a vortex: E log r where r is the size of a vortex. As a result, at low temperature, a configuration without vortices is favored, while at high temperature, a configuration with a large number of vortices is favored. Between these configurations, a phase transition occurs. From the observations above, it seems that the logarithmic energy barrier may lead to the thermal stability, while the existence of local minima may break the thermal stability. Yet, at this moment, the connection between the energy barrier and the thermal stability has not been completely established. Appendix B. Topology and quantum codes In section 5, we have shown that one can deform geometric shapes of logical operators continuously while keeping them equivalent. This continuous deformability of logical operators implies that one can introduce the notion of topology in analyzing and classifying coding and physical properties of STS models. In this appendix, we further discuss the role of topology in analyzing stabilizer codes. Many interesting properties of STS models, such as the dimensional duality of logical operators, naturally appear as corollaries of the continuous deformability of logical operators. In Appendix B.1, we begin by showing that the dimensional duality of logical operators can be derived only by assuming the continuous deformability of logical operators. In particular, we show that, for stabilizer codes with continuously deformable logical operators, m-dimensional and D m-dimensional logical operators always form anti-commuting pairs where m is an arbitrary positive integer (m D). As an example of D-dimensional stabilizer codes with continuously deformable logical operators, we present generalizations of the Toric code to D-dimensional systems with anti-commuting pairs of m-dimensional and D m-dimensional logical operators for arbitrary positive integers D and m. The topological property of logical operators naturally leads us to consider a possible relevance to another wellcelebrated theoretical framework equipped with the notion of topology; called topological quantum field theory (TQFT) [65 67]. In Appendix B.2, we establish the connection between quantum codes with the continuous deformability and TQFT further by demonstrating that the braiding of anyonic excitations in a D-dimensional stabilizer code is characterized by a topological invariant in a D+1-dimensional system. This implies that such a D-dimensional stabilizer code can be effectively described by D + 1-dimensional TQFT. 38

39 Appendix B.1. Dimensional duality as a corollary of continuous deformability The continuous deformability of logical operators implies that topology is the essential notion in analyzing coding and physical properties of STS models. Then, a naturally arising question concerns the role of topology in determining coding and physical properties of stabilizer codes. Here, we show that the dimensional duality of logical operators is a universal property for all the stabilizer codes with continuously deformable logical operators. We also give concrete examples of D-dimensional stabilizer codes which have continuous deformability of logical operators and the dimensional duality by generalizing the Toric code to D-dimensional systems. Dimensional duality from topological deformation: Let us begin by counting the number of independent logical operators defined inside m-dimensional regions. Recall that m-dimensional concatenated unit regions are obtained by taking unions of all the m-dimensional topological unit regions. Then, one may call logical operators which can be defined inside R m, but cannot be defined inside R m 1, m-dimensional logical operators. Definition 5. m-dimensional logical operators are logical operators which have representations supported inside R m, but do not have representations supported inside R m 1. Now, let us denote the number of independent m-dimensional logical operators as g m g Rm g Rm 1 where g 0 g R0 by setting g R 1 0. Then, there exists an interesting relation among the numbers of m-dimensional logical operators. In particular, the following lemma holds. Lemma 3. There are the same number of m-dimensional and D m-dimensional logical operators: g m = g D m for m = 0,, D. (B.1) The proof of this lemma can be obtained through a simple algebra by combining theorem 1 and the deformability of logical operators. Proof. Consider a bi-partition of the entire system into R m and R m. From the topological deformation of logical operators, we have since R m R D m 1. Thus, R m and R m support the following logical operators: g Rm = g RD m 1 (B.2) R m : R m : 0-dim, 1-dim,, m-dim 0-dim, 1-dim,, D m 1-dim. Therefore, we have g Rm = m g j, g Rm = j=0 D m 1 j=0 g j. (B.3) Recall that g R + g R = 2k as presented in theorem 1. Using this formula for R = R m, we have Then, we have m j=0 g Rm + g Rm = 2k. (B.4) D m 1 g j + g j = 2k, for m = 0,, D. (B.5) j=0 Since the total number of independent logical operators is D j=0 g j = 2k, we have g m = g D m for all m. 39

40 This lemma implies the existence of a dimensional duality in geometric shapes of logical operators. One may notice that this lemma is just a manifestation of Poincaré duality in a D-torus where the mth and D mth Betti numbers are equal. To completely establish relations between each logical operator with different dimensions, let us analyze their commutation relations. We have the following theorem. Theorem 7 (Dimensional duality). One can choose a set of 2k independent logical operators of D-dimensional systems with the deformability of logical operators in the following way: { } l1,, l k. (B.6) r 1,, r k where l p are m p -dimensional logical operators and r p are D m p -dimensional logical operators for some integer m p (0 m p D) for any p = 1,, k. In other words, one can choose logical operators such that the summation of dimensions of pairs of anti-commuting logical operators is always D. Theorem 7 follows immediately from the following lemma. Lemma 4. m-dimensional and m -dimensional logical operators commute with each other if m + m < D. Proof. Consider a m-dimensional logical operator l and a m -dimensional logical operator l which is defined inside R m and R m respectively. For m + m < D, there exists a translation of R m such that R m and R m have no overlap. Then, due to the translation equivalence of logical operators, some translation of l do not have an overlap with l, which leads to [l, l ] = 0. With this lemma, the proof of theorem 7 is immediate by using lemma 3. For example, from the lemma, zerodimensional logical operators may anti-commute only with D-dimensional logical operators. Since there are the same number of zero-dimensional and D-dimensional logical operators, there exists a canonical set of logical operators where D-dimensional logical operators can anti-commutes only with zero-dimensional logical operators. Similarly, one can show that there exists a set of 2k independent logical operators such that m-dimensional logical operators anti-commute only with D m-dimensional logical operators for all m. Generalization of the Toric code: Now, we give concrete examples of D-dimensional stabilizer codes which have continuously deformable logical operators. The model we present here is a straightforward generalization of two-dimensional Toric code to D-dimensional settings for arbitrary positive integer D. In particular, we illustrate the construction of D-dimensional Toric code with anti-commuting pairs of m-dimensional and D m dimensional logical operators for arbitrary positive integers D and m (m D). Figure B.25: An example of the construction for D = 3 and m = 2. We first consider a D-dimensional hypercubic lattice which consists of N = L L of D-dimensional unit hypercubes with periodic boundary conditions. We denote a set of p-dimensional unit hypercubes in this lattice as h p. Note that h p = D C p N since one needs to specify p directions from D directions in defining p-dimensional unit hypercubes. We put qubits at the centers of m-dimensional hypercubes. Then, the total number of qubits is D C m N, 40

41 and qubits are labeled by m-dimensional unit hypercubes in h m. Fig. B.25 shows a construction of the model for D = 3 and m = 2 where qubits live at centers of two-dimensional unit squares. The Hamiltonian consists of plaquette terms and star terms as in the conventional two-dimensional Toric code: H = H plaquette + H star (B.7) Here, we call (m + 1)-dimensional unit hypercubes in h m+1 plaquettes. Then, plaquette terms consist of Pauli Z operators which act on qubits included in plaquettes: H plaquette = B p, B p = Z r (B.8) p h m+1 r p where r represent m-dimensional unit hypercubes included inside a m + 1-dimensional unit hypercube p. Note that there are D C m+1 N plaquette terms. Next, we call (m 1)-dimensional unit hypercubes in h m 1 stars. Then, star terms consist of Pauli X operators which act on qubits neighboring to a star: H star = A s, A s = X r (B.9) s h m 1 s r where s r means that a star s is included inside m-dimensional unit hypercube r. Note that there are D C m 1 N star terms. Noting that [A s, B p ] = 0 since A s and B p share either zero or two qubits in common, the model is a stabilizer code. Fig. B.25 shows constructions of star terms and plaquette terms for D = 3 and m = 2. One can verify that the model has k = D C m logical qubits with k anti-commuting pairs of m-dimensional and (D m)-dimensional logical operators where a m-dimensional logical operator consists of Pauli Z operators supported on a m-dimensional hyperplane, while a (D m)-dimensional logical operator consists of Pauli X operators supported on a (D m)-dimensional hyperplane. Since they may share either zero or one qubit, they may commute or anti-commute with each other. One may easily see that the construction above reproduces two-dimensional, three-dimensional and four-dimensional Toric code for choices of (D, m) = (2, 1), (3, 1), (4, 2). For m = 0 and m = D, the model is reduced to the D-dimensional Ising model. One can verify that logical operators arising in this model can be deformed continuously by using the bi-partition theorem (theorem 1). The construction above can be easily generalized to arbitrary D-dimensional graph embedded in a D-dimensional geometric manifold in a way similar to two-dimensional Toric code. Appendix B.2. D-dimensional STS model and D + 1-dimensional TQFT The emergence of the notion of topology in geometric shapes of logical operators leads us to consider a possible relevance to topological quantum field theory (TQFT), which also deals with systems whose physical properties depend heavily on topological characteristics of the systems. Here, we make an attempt to establish the connection between stabilizer codes with continuously deformable logical operators and TQFT further. Roughly speaking, TQFT is a field theory which is invariant under diffeomorphism (continuous deformation), and particularly suited for describing topologically ordered systems. The most important characteristic in systems described by TQFT is the invariance of all the correlation functions under diffeomorphism. Consider an arbitrary diffeomorphism through a continuous change of space-time coordinates x x. Then, the correlation function of a scalar operator φ(x) satisfies 0 i φ(x 1 )φ(x 2 ) φ(x n ) 0 j = 0 i φ(x 1 )φ(x 2 ) φ(x n) 0 j (B.10) where 0 i represent degenerate ground states. Therefore, the vacuum expectation value of any products of scalar operators is invariant under differmorphism, and only the topological properties of products may characterize their expectation values. Below, in order to establish the connection between stabilizer codes and TQFT, we show that the braiding of anyonic excitations in D-dimensional stabilizer codes can be characterized by some topological invariant in a D + 1- dimensional system. In particular, by characterizing propagations of anyonic excitations in a D + 1-dimensional 41

42 system, we demonstrate that the braiding of anyonic excitations corresponds to a configuration of m-dimensional and D m-dimensional closed objects with a non-zero linking number. We begin by analyzing propagations of anyonic excitations in two-dimensional Toric code. Anyonic excitations can be created by applying a segment of a one-dimensional logical operator to a ground state of the Toric code Hamiltonian since endpoints of a segment of a logical operator l seg may not commute with interaction terms and create localized excitations (Fig. B.26(a)): l seg ψ gs, and one may make anyons propagate along a geometric shape of a one-dimensional logical operator. Since one can deform a geometric shape of a one-dimensional logical operator continuously in the Toric code, anyons can propagate freely on the lattice by applying a segment of a deformed onedimensional logical operator (Fig B.26(b)). Therefore, the continuous deformability of logical operators is the key to propagations of anyonic excitations. Figure B.26: Anyonic excitations created from segments of logical operators. (a) A one-dimensional logical operator. (b) A deformed onedimensional logical operator. Propagations of anyonic excitations can be characterized by a one-dimensional closed loop drawn in a threedimensional space. To illustrate this point, let us consider a process of creation, propagation, and annihilation of anyons, as described in Fig. B.27(a). By drawing propagations of anyonic excitations in a three-dimensional system by adding the time axis, the entire process can be represented as a one-dimensional closed loop as shown in Fig. B.27(b). In general, arbitrary one-dimensional closed loop in a three-dimensional system may characterize some propagations of anyonic excitations. See an example in Fig. B.28 which involves creations of two pairs of anyonic excitations. Figure B.27: The correspondence between anyonic excitations in two-dimensional Toric code and a closed loop in a three-dimensional system. (a) Creation, propagation and annihilation of anyonic excitations. (b) Anyonic excitations described in a three-dimensional system. Next, let us consider the braiding between anyonic excitations in the Toric code. One interesting property of topologically ordered spin systems is the non-trivial braiding property between anyonic excitations. There are two 42

43 Figure B.28: Propagations of anyonic excitations and an associated closed loop in a three-dimensional system. distinct type of anyons in two-dimensional Toric code which are created by a pair of anti-commuting logical operators respectively, and the braiding of different types of anyonic excitations may give rise to non-trivial change inside the ground space. In two-dimensional Toric code, the braiding of different anyons give rise to an additional phase 1 to the original ground state ψ gs ψ gs, which is a direct consequence of the anti-commutation between logical operators. One may understand this non-trivial braiding arising in two-dimensional Toric code as a topological invariant in a three-dimensional system. Let us consider the braiding of anyons described in Fig. B.29(a) where two types of pairs of anyonic excitations are created. One can characterize this braiding process as two closed loops in a threedimensional system where loops are linked as described in Fig. B.29(b). In a more technical language, the braiding of anyons occurs only if the linking number between two loops has non-zero value. In particular, the final state is ψ gs ( 1) N link ψ gs where N link is the linking number of a given configuration of one-dimensional closed loops. This observation indicates that the braiding of anyonic excitations can be characterized by topological invariants, such as the linking number, in a (2 + 1)-dimensional system. Figure B.29: The braiding as a topological invariant in a three-dimensional system. (a) A braiding of anyonic excitations. (b) Loops with non-zero linking number. A similar observation holds for D-dimensional Toric code with a pair of m-dimensional and D m-dimensional logical operators. The first type of anyonic excitations can be created by a segment of m-dimensional logical operator, while the second type of anyonic excitations can be created by a segment of D m-dimensional logical operator. 43

44 The propagation of anyonic excitations of the first type can be characterized by a m-dimensional object in a D + 1- dimensional systems, while the propagation of anyonic excitations of the second type can be characterized by a D m-dimensional object. These anyonic excitations are braided when the linking number between m-dimensional and D m-dimensional objects is non-zero. While we have discussed the braiding for D-dimensional Toric code, similar discussion holds for any stabilizer codes with continuously deformable logical operators. Therefore, one may expect that D-dimensional stabilizer codes with continuous deformability can be effectively described by D + 1- dimensional TQFT. Appendix C. Decomposition of logical operators We present the proof of theorem 5 in this and the next appendices. The goal of this appendix is to prove the following theorem which will be the key to the proof of theorem 5. Theorem 8 (Decomposition). Consider a three-dimensional STS model with the system size n 1 = 2 2 2n 2v!, n 2 = 2 m and arbitrary n 3 where m is an arbitrary positive integer. For a given logical operator l supported inside P(n 1, n 2, 1), one can decompose l as a product of the following centralizer operators where l l a l b, l a, l b C P(n1,n 2,1) (C.1) T β 1 (l b) = l b, β 2 2n 2v (C.2) and l a is defined inside P(2v, n 2, 1). Here, C R represents the restriction of the centralizer group C onto a region of composite particles R, meaning that C R is a subgroup of centralizer operators defined inside R. Therefore, l a, l b C P(n1,n 2,1) means that l a and l b are centralizer operators defined inside P(n 1, n 2, 1). We show the claim of the theorem graphically in Fig. C.30. The theorem claims that a two-dimensional logical operator defined inside P(n 1, n 2, 1) can be decomposed as a product of a one-dimensional centralizer operator l a and a two-dimensional centralizer operator l b which is periodic in the ˆ1 direction. As we shall see later, 2v comes from the number of independent generators for the Pauli group acting on a single composite particle. Before starting the proof of theorem 8, let us describe the entire sketch of the proof of theorem 5. As a simple extension of theorem 8, one can show that a one-dimensional logical operator defined inside P(2v, n 2, 1) can be further decomposed as a product of a one-dimensional and a zero-dimensional centralizer operators. After these decompositions, one can classify geometric shapes of logical operators according to their dimensions and can find commutation relations between them. Although theorem 8 is limited to some specially chosen system sizes: n 1 = 2 2 2n 2v! and n 2 = 2 m, one can construct logical operators for arbitrary system sizes from theorem 8. For example, due to scale symmetries, one can show that one-dimensional logical operators found in theorem 8 are also logical operators for the systems with arbitrary n 1. In fact, one can find logical operators in the forms described in theorem 5. These arguments will be presented in Appendix D. Appendix C.1. Sketch of proof of theorem 8 First, we note that theorem 8 was proven for n 2 = 1 (m = 0) in [17] since such a system with n 2 = 1 can be considered as a two-dimensional system which extends only in the ˆ1 and ˆ3 directions. For a two-dimensional STS model, we have the following lemma. Lemma 5. Consider a two-dimensional STS model where n 1 = 2 2 2v! and arbitrary n 2. For a given logical operator l supported inside P(n 1, 1), one can decompose l as a product of the following centralizer operators l l a l b, l a, l b C P(n1,1) (C.3) 44

45 Figure C.30: The claim of theorem 8. One can decompose a two-dimensional logical operator as a product of a one-dimensional centralizer operator l a and a two-dimensional centralizer operator l b. where and l a is defined inside P(2v, 1). T β 1 (l b) = l b, β 2 2v (C.4) We present the claim of the lemma graphically in Fig. C.31. The lemma claims that a one-dimensional logical operator defined inside P(n 1, 1) can be decomposed as a product of a zero-dimensional centralizer operator l a and a one-dimensional centralizer operator l b which is periodic. Figure C.31: The claim of lemma 5. One can decompose a one-dimensional logical operator as a product of a zero-dimensional centralizer operator l a and a one-dimensional centralizer operator l b. A three-dimensional STS model may be viewed as a two-dimensional system if one considers 1 n 2 1 composite particles as a single composite particle which consists of vn 2 qubits (see Fig. C.32). In other words, we view the entire system as a two-dimensional lattice of one-dimensional tubes. Then, as a direct consequence of the lemma above, we notice the following corollary. Corollary 1. A logical operator l considered in theorem 8 can be decomposed as a product of the following centralizer 45

46 operators l l a l b, l a, l b C P(n1,n 2,1) (C.5) where T β 1 (l b) = l b, where β 2 2n 2v (C.6) and l a is defined inside P(2vn 2, n 2, 1). We present the claim of the corollary graphically in Fig. C.32. The corollary claims that a two-dimensional logical operator defined inside P(n 1, n 2, 1) can be decomposed as a product of a one-dimensional logical operator l a and a two-dimensional logical operator l b which is periodic in the ˆ1 direction. However, a one-dimensional logical operator l a described in corollary 2 is not one-dimensional in a strict sense since it is defined inside P(2 m 2v, 2 m, 1), and its width 2 m 2v grows as m increases. On the other hand, l a described in theorem 8 is truly one-dimensional since its width is at most 2v. Therefore, we need to show that a logical operator defined inside P(2 m 2v, 2 m, 1) have an equivalent logical operator defined inside P(2v, 2 m, 1). Figure C.32: The claim of corollary 2. The width of a one-dimensional logical operator l a increases as m increases. The rest of this appendix is dedicated to the proof of the following lemma. Lemma 6 (Shrinkage). For system sizes considered in theorem 8, a logical operator operator l defined inside P(x, 2 m, 1) always has an equivalent logical operator l which is defined inside P(x 1, 2 m, 1) when 2v < x < n 1. By using this lemma, one can shorten the width of l a from 2 m 2v to 2v. Appendix C.2. Identity generating matrix Here, we discuss how to shorten the width of l a. In particular, we introduce a certain binary matrix which is essential in reducing the width of l a. Shrinkage in two dimensions: To give an intuition on how to shorten the width of l a, let us first consider the case where m = 0. (So, this is a two-dimensional system with n 2 = 1, and instead of the width, we use the length ). As an example, consider the case when m = 0, x = 4 and l is given by l = [ A, B, C, AB ] (C.7) 46

47 where l is defined inside P(4, 1, 1), and A, B and C are some Pauli operators. Then, consider the following logical operator: l lt1 2 (l)t 3 1 (l) (C.8) = [ A, B, AC, I, BC, ABC, AB ]. (C.9) Note that l l since l is a product of three logical operators which are equivalent to each other due to the translation equivalence of logical operators. Then, we notice that following operators are centralizer operators: l 1 = [ A, B, AC ], l 2 = [ BC, ABC, AB ] (C.10) where l = l 1 T 4 1 (l 2) since stabilizers in STS models are defined inside 2 2 composite particles and cannot overlap with l 1 and T 4 1 (l 2) simultaneously. Now, due to the translation equivalence of logical operators, we have l l l 1 l 2 = [ A, B, AC ] [ BC, ABC, AB ] = [ ABC, AC, BC ]. (C.11) Thus, a logical operator l, with the length 4, is shrunk into an equivalent logical operator l, with the length 3. An important observation is that one can form an identity operator I by taking a product of Pauli operators in l. In general, if the length x of l is larger than 2v, one can always form an identity operator I by taking a product of some Pauli operators in l since there are 2v independent generators for single Pauli operators acting on a single composite particle. Now, let us consider l represented as l = [ U 1, U 2,, U x ]. (C.12) which is defined inside P(x, 1, 1) where x > 2v. Then, there always exists a binary vector B = (B 1,, B x ) (0,, 0) which satisfies the following condition: Now, let us take the following product of translations of l: x j=1 x j=1 U B j j = I. (C.13) T x j 1 (l B j ). (C.14) Then, one may readily know that the xth entry of the above operator is I. From this operator, one can find two centralizer operators. By using them, one can readily shrink the length of l from x to x 1. This trick is the key to the proof of lemma 5. Although the argument above works only when B has an odd number of 1 entries, a slight modification makes the shrinkage of l possible when B has an even number of 1 entries. Identity generating matrix: Now, we consider more general cases with m > 0. Let us represent a logical operator l defined inside a region P(x, 2 m, 1) as a x 2 m matrix whose entries are single Pauli operators: l = U 1,1, U x, U 1,2 m, U x,2 m where each Pauli operator U i, j acts on each composite particle. We also represent each column of l as follows: U j = U j,1. U j,2 m 47 (C.15) ( j = 1,, x). (C.16)

48 Here, we denote a group of Pauli operators supported by a single column P(1, 2 m, 1) as P m col and call it the column operator group and its elements column operators. Note that U j P m col, and Pm col has 2m 2v independent generators: G(P m col ) = 2m 2v. When m = 0 and x > 2v, we found a binary vector B with x components which characterizes how to form an identity operator from l. When m > 0 and x > 2v, we can find an x 2 m binary matrix B which characterizes how to form an identity operator from l. Here, we introduce the identity generating matrices as follows. Definition 6. Consider a logical operator l defined inside P(x, 2 m, 1). For a x 2 m binary matrix B B = B 1,1, B x,1....., B i, j = 0, 1, (C.17) B 1,2 m, B x,2 m we define the following operations: l(b) l(b) x i, j x i=1 T x i 1 T j 1 2 (l B i, j ) (C.18) 2 m j=1 T j 1 2 (U B i, j i ) P m col (C.19) where l(b) is a product of translations of l taken according to B while l(b) x is the xth column of l(b). We call a binary matrix B identity generating matrix if and only if 0, 0 l(b) x = I and B , 0 (C.20) We assign parities to each column of a binary matrix B as follows: Par(B) i B i, j (mod 2) (C.21) j where i = 1,, x. We call a binary matrix B odd if and only if i s.t Par(B) i = 1. (C.22) Therefore, when we form an identity operator, we considered translations of l both in the ˆ1 and ˆ2 directions. The identity generating matrix is said to be odd when there exists a column with an odd parity. Shrinkage through odd matrices: Note that there always exists an identity generating matrix when x > 2v. Then, with the existence of identity generating matrices, one might hope that the width of l a can be reduced until it becomes 2v in a way similar to the cases where m = 0. However, there is a caveat. In fact, only identity generating matrices with some special properties can be used for shrinking. In particular, we have the following lemma. Lemma 7 (Shrinkage through odd matrices). If there exists an odd identity generating matrix for l defined inside P(x, 2 m, 1), l has an equivalent logical operator l defined inside P(x 1, 2 m, 1). Therefore, if there exists an odd identity generating matrix for any x with n 1 > x > 2v, one can complete the proof of lemma 6. Below, we present the proof of lemma 7. The existence of an odd identity generating matrix will be proven later. 48

49 Proof. Assume that B is an odd identity generating matrix for l. Assume that for some i (1 i x), we have Par(B) i = 1 and Par(B) i = 0 for i < i. (C.23) So, i is the smallest integer such that i th column has an odd parity. Here, we define the following binary matrix B (see Fig. C.33): B i, j B i, j ( i i ) (C.24) 0 ( i > i ). (C.25) Note that B consists of ith columns of B with i i. Based on B, we consider the following logical operator l : l l(b ) l. (C.26) Note that l is equivalent to l since l is a product of an odd number of translations of l. (Note that l(b) may not be equivalent to l since the number of 1 entries in B may be even). See Fig C.33 for graphical representations of B, B, l(b) and l(b ). Note that l(b) has an identity operator at xth column. Note that l(b ) has identity operators in the first x i columns since B i, j = 0 for i > i. Figure C.33: Constructions of B, B, l(b) and l(b ). Since l(b) has an identity operator at xth column, we can decompose it as a product of two centralizer operators whose lengths are at most x 1. Let us denote these centralizer operators as U a and U b : l(b) = U a U b (C.27) where U a is the centralizer on the left hand side and U b is the centralizer on the right hand side (See Fig. C.33). U a is defined from 1st column to x 1th column, and U b is defined from x + 1th column to 2x 1th column. Since B and B have the same entries from 1st column to i th column, we notice that l(b ) and U b have the same Pauli operators from 2x 2 i st column to 2x 1th column as shown in Fig. C.33. Then, by applying U b to l(b ), one can shrink the size of l(b ). In particular, U b l(b ) has the length at most x 1. Although U b l(b ) may not be equivalent to l(b ) as U b may not be a stabilizer, one may consider the following logical operator: U b l(b ) T i 1 (U b) l(b ) l (C.28) which is equivalent to l(b ) due to the translation equivalence of logical operators, and is defined inside a region with x 1 2 m 1 composite particles. Then, due to the translation equivalence of logical operators, there exists a logical operator l l which is defined inside P(x 1, 2 m, 1). This completes the proof. 49

50 Appendix C.3. Existence of an odd matrix for m = 1 Next, we present a proof of the existence of an odd identity generating matrix for x > 2v, in order to complete the proof of lemma 6 and theorem 8. In particular, we shall prove the following lemma. Lemma 8. When x > 2v, there always exist an odd identity generating matrix. We start by discussing cases with m = 0 and m = 1 before presenting general discussion. First of all, when m = 0, identity generating matrices are always odd since all the binary matrices are odd except B = (0,, 0). Therefore, we consider the cases where m = 1 below. Characteristic value: Recall that we represented l as a x 2 binary matrix: [ ] U1,1, U l = x,1 (C.29) U 1,2, U x,2 and each column of l as follows: U j = [ ] U j,1 U j,2 ( j = 1,, x). (C.30) Now, we suppose that there is no odd identity generating matrix for l, in order to use the contradiction for the proof of lemma 8. First, it is convenient to classify column operators in P 1 col into two types as follows. For a column operator U represented as [ ] U1,1 U =, (C.31) U 1,2 we assign a characteristic value b and characteristic operator V as follows: If U 1,1 = U 1,2, we assign a characteristic value b = 1 and a characteristic operator V = U 1,1. If U 1,1 U 1,2, we assign a characteristic value b = 0 and a characteristic operator V = U 1,1 U 1,2. One may easily understand this classification by representing U explicitly. A column operator U with b = 1 is [ V U = (C.32) V] and a column operator U with b = 0 is U = [ ] VU1,2. (C.33) U 1,2 So, a column operator with b = 1 is symmetric while a column operator with b = 0 is not. Note that a characteristic operator is not an identity operator I except when U = I. Next, we introduce the following x 2 binary matrices: [ ] E(i E1,1,, E ; 0) x,1 (C.34) E 1,2,, E x,2 such that E i,1 = 1 (i = i ) (C.35) E i, j = 0 otherwise (C.36) 50

51 and E(i ; 1) [ ] E1,1,, E x,1 E 1,2,, E x,2 (C.37) such that E i,1 = E i,2 = 1 (i = i ) (C.38) E i, j = 0 otherwise (C.39) For example, E(2; 0) = [ 0, 1, 0,, ] 0 0, 0, 0,, 0 (C.40) and E(2; 1) = [ ] 0, 1, 0,, 0. (C.41) 0, 1, 0,, 0 Now, let us represent characteristic values and characteristic operators for each column of l as follows: U 1 b 1, V 1 U 2 b 2, V 2... U x b x, V x. (C.42) Then one can establish a connection between binary matrices E(i; 0) and E(i; 1), and characteristic values b i and operators V j as follows: [ Vi ] l(e(i; 0)) x = V i [ ] Vi l(e(i; 1)) x = V i when b i = 1 (C.43) when b i = 0. (C.44) Proof of lemma 8: Now, let us proceed to the proof of lemma 8 for m = 1. Without loss of generality, we can assume that b 1 = = b x0 = 0, b x0 +1 = = b x = 1 (C.45) for x 0 x since permutations of columns do not affect the parities of identity generating matrices. (We will justify this later). We define the following sets of integers: b(0) {1,, x 0 }, b(1) {x 0 + 1,, x}. (C.46) We denote groups of Pauli operators generated by V j with b j = 0 and V j with b j = 1 as V 0 and V 1 : V 0 { V i : i b(0) } (C.47) V 1 { V i : i b(1) }. (C.48) Let us show that a set of characteristic operators {V i } for i b(1) is independent : G(V 1 ) = x x 0. For this 51

52 purpose, we suppose that there exists some set of integers A b(1) such that V i = I. (C.49) i A Then, the following binary matrix is an identity generating matrix: B = E(i; 0) i A (C.50) since l(b) x = i A [ ] Vi V i = [ I I]. (C.51) However, since B is odd with Par(B) i = 1 for i A, this leads to a contradiction. Next, let us analyze V 0. For simplicity of discussion, we first assume that {V i } for i b(0) are independent. We consider more general cases where {V i } for i b(0) are over complete later. Then, we have G(V 0 ) = x 0. Now, we define the following operators for i b(0): U i U i T 2 (U i ) [ ] Vi =. V i (C.52) Note that U i has a characteristic value b i = 1 and a characteristic operator V i. Notice that U i = l(e(i; 1)) x. (C.53) Since x > 2v, there exists a set of integer A such that V i = I, A b(1) and A b(0). (C.54) i A Note that A includes integers both from b(0) and b(1). Then, one notices that the following matrix B is an identity generating matrix: B = E(i; 1) + E(i; 0) (C.55) i A b(0) i A b(1) since l(b) x = = U i U i {i A b(0)} {i A b(1)} i A [ ] Vi V i = [ I I]. (C.56) However, for an integer i in A b(0), Par(B) i = 1, and B is an odd matrix. This leads to a contradiction. Note that discussion above is valid under the permutations of columns. Proof of lemma 8, continued: Next, let us consider the case where {V i } for i b(0) are not independent. Let us denote the number of generators for V 0 as x 0 = G(V 0 ) ( x 0 < x 0 ). Without loss of generality, we may assume that V 1,, V x0 are independent since permutations do not affect parities of identity generating matrices. Here, we change our notations slightly: b(0) {1,, x 0 }, b(0) { x 0 + 1,, x 0 }, b(1) {x 0 + 1,, x}. (C.57) 52

53 and V 0 { V i : i b(0) } G(V 0 ) = x 0 (C.58) V 1 { V i : i b(1) } G(V 1 ) = x x 0. (C.59) Since a set {V i } for b i = 0 is over complete, there are x 0 x 0 sets of integers A i (i = x 0 + 1,, x 0 ) such that V i = I where i i, i A i and i A i (C.60) i A i where the largest integer in A i is i. For i = x 0 + 1,, x 0, we form the following operator: [ ] V i U i j A i U j V i (C.61) which has a characteristic value b i = 1 and a characteristic operator V i. Here, we denote a group of {V i } for i b(0) as V 0 = { V i : i b(0) }. (C.62) Now, we show that {V i } for i b(0) are independent: G(V 0 ) = x 0 x 0. If there exists A b(0) such that V i = I, (C.63) i A we have the following identity generating matrix: B = E( j; 0) (mod 2) (C.64) i A j A i since l(b) x = [ ] V i V i A i = [ I I]. (C.65) Recall that the largest integer in A i is i. Let the largest integer in A be i max. Then, we have Par(B) imax = 1, and thus, B is odd. This leads to a contradiction. So far, we have shown that G(V 0 ) = x 0, G(V 0 ) = x 0 x 0, G(V 1 ) = x x 1. (C.66) Since x > 2v, there exists a set of integers A such that V i V i {i A b(1)} {i A b(1) } {i A b(0)} The following matrix is the identity generating matrix: B = E(i; 1) + E(i ; 0) + i A b(0) i A b(0) i A i V i = I. (C.67) E(i; 0) (mod 2) (C.68) i A b(1) 53

54 since l(b) x = {i A:i b(0)} [ ] Vi V i {i A:i b(0) } [ ] V i V i {i A:i b(1)} [ ] Vi V i = [ I I]. (C.69) Since A b(0), A has some element in b(0) b(1). Let the largest integer in A be i max. Then, we have Par(B) imax = 1, and thus, B is odd. This leads to a contradiction. Again, permutations of columns do not affect this discussion. This completes the proof of lemma 8 for m = 1. Appendix C.4. Characteristic vectors Let us proceed to the proof for the cases where m > 1. When m = 1, we assigned characteristic values 0 and 1 to each column operator according to its symmetry. For m > 1, we will assign a binary vector with m components to each column operator, which we will call a characteristic vector. We encode symmetries of a column operator on these characteristic vectors. Characteristic vector: We first define two maps f 0 and f 1 from P m column operator U P m col which is represented as we define f 0 (U), f 1 (U) P m 1 col as follows: U = U 1. U 2 m, col to 1 Pm col for m > 0 as follows. For a given (C.70) f 0 (U) j U j U j+2 m 1, f 1 (U) j U j ( j = 1,, 2 m 1 ). (C.71) One may represent f 0 (U) and f 1 (U) more explicitly: U 1 U 2 m 1 +1 f 0 (U) =., f 1(U) = U 2 m 1U 2 m U 1. U 2 m 1 (C.72) Note that f 0 and f 1 decrease the length of the column by half. Let us denote a set of all the m component binary vectors as B m vec. Now, for a column operator U P m col, we assign an m component binary vector b B m vec through the following rule. If f 0 (U) = I, take b m = 1 and define U (1) f 1 (U). If f 0 (U) I, take b m = 0 and define U (1) f 0 (U). and, iterate this procedure: If f 0 (U ( j) ) = I, b m j = 1 and define U ( j+1) f 1 (U ( j) ). If f 0 (U ( j) ) I, b m j = 0 and define U ( j+1) f 0 (U ( j) ). for 1 j m 1. We define a characteristic operator V of U as follows: V U (m) = f b1 f b2 f bm (U) P 0 col. (C.73) Note that V I when U I. It is worth presenting examples of characteristic vectors here (m = 2): V V (1, 1), V V V I (0, 1), V I V V (1, 0), I I V I (0, 0). I I (C.74) 54

55 Thus, symmetries of column operators are encoded in characteristic vectors. Next, we introduce an order between binary vectors in B m vec. We define g( b) m b j 2 j 1 j=1 (C.75) where b is like a binary representation of an integer g( b). For a given pair of m component binary vectors b and b, we denote b < b (C.76) if and only if g( b) < g( b ). (C.77) For example, for m = 3, we have the following relations between binary vectors: (0, 0, 0) < (1, 0, 0) < (0, 1, 0) < (1, 1, 0) < (0, 0, 1) < (1, 0, 1) < (0, 1, 1) < (1, 1, 1). (C.78) Below, we shall see that a column operator with a larger characteristic vector is more symmetric than a column vector with a smaller characteristic vector. Property of characteristic vectors: Let us briefly recall the proof of lemma 8 for m = 1. In the proof, we constructed a column operator with b = 1 from a column operator with b = 0. In particular, if U is a column operator with a characteristic value b = 0 and a characteristic operator V: [ ] VV U = V (C.79) where V is some Pauli operator, we have UT 2 (U) = [ V V] (C.80) which is a column operator with a characteristic value b = 1 and a characteristic operator V. Thus, we can create a column operator with a larger characteristic value from a column operator with a smaller characteristic value. In a way similar to this, one can construct a column operator with b from a column operator with b as long as b > b. Let us represent a binary column B as follows: B = B 1,1. B 1,2 m and denote a set of all the binary columns as B m col. Here, we define a parity of B as, (C.81) Par(B) Par(B) 1 (C.82) by viewing a binary column B as a binary matrix. For a column operator U P m col, we define U(B) as follows: U(B) 2 m j=1 T j 1 2 (U B 1, j ) P m col (C.83) just like l(b). Note that U(B) is a product of translations of U taken according to a binary column B. Then, the following lemma holds. 55

56 Lemma 9. Consider an arbitrary pair of m component binary vectors b, b B m vec such that b < b. For a given column operator U which has a characteristic vector b and a characteristic operator V, there always exists some binary column B B m col with an even parity Par(B) = 0 such that U(B) has a characteristic vector b and a characteristic operator V. In other words, from a column operator U with a characteristic vector b, one can always create a column operator U with larger characteristic vector b by taking a product of translations of U. On the other hand, it is impossible to create a column operator with a smaller characteristic vector from an operator with a larger characteristic vector. Therefore, one can create a column operator with higher symmetries (a larger characteristic vector), but cannot create a column operator with lower symmetries (a smaller characteristic vector). Appendix C.5. Proof of lemma 9 Now, we prove lemma 9 by explicitly finding a binary column B B m col for creating U(B) for every pair of b and b. In order to derive such binary matrices, we introduce a certain binary column B( b), called a characteristic column, which can be used to change the characteristic vector of a column operator. Characteristic column: Given an integer p Z 2 m, one may have its binary representation by considering the inverse of g denoted as g 1 : p (p 1,, p m ) g 1 (p) (C.84) where p = g( p) = m 1 j=1 p j2 j 1. Now, we define the following sets: J b = { a Bm vec : a j b j for all j }. (C.85) For example, J (1,0) = {(0, 0), (1, 0)} and J (0,0,1) = {(0, 0, 0), (0, 0, 1)}. Based on J b, we define the characteristic binary column B( b) B m col as follows: Here, we give some examples: B( b) 1,p+1 = 1 p J b (C.86) B( b) 1,p+1 = 0 otherwise. (C.87) J (0,0) = {(0, 0)}, J (1,0) = {(0, 0), (1, 0)}, J (0,1) = {(0, 0), (0, 1)} J (1,1) = {(0, 0), (1, 0), (0, 1), (1, 1)} (C.88) and B(0, 0) = 1 0, B(1, 0) = , B(0, 1) = , B(1, 1) = (C.89) 1 1 One may see the relation between a characteristic column and a characteristic vector: V I (0, 0), I I V V (1, 0), I I V I (0, 1), V I V V (1, 1). V V (C.90) Thus, if we replace 1 entries in B( b) with V and create a column operator, it has a characteristic vector b. In order to discuss changes of characteristic vectors, let us introduce the summation rule between binary vectors. 56

57 We denote a summation of binary vectors a, b B m vec as a + b B m vec and define it as follows: g( a) + g( b) = g( a + b) (C.91) when g( a + b) 2 m 1. Therefore, a + b is just like a summation of two binary numbers a and b. The characteristic columns defined above can be used to change a characteristic vector of a column operator, as summarized in the following lemma. Lemma 10. Let b < b. When U has a characteristic value b and a characteristic operator V, has a characteristic value b and a characteristic operator V. U = U(B( b)) where b + b = b (C.92) Below, we present a proof of this lemma by finding some property of characteristic columns and a certain rule on multiplications of column operators. Property of characteristic columns: There is a useful relation between a summation of vectors and characteristic columns, as summarized in the following lemma. Lemma 11. Let B( a) B( b) 2 m j=1 T 2 (B( a) B( b) 1, j ) j 1 (mod 2). (C.93) Then, B( a + b) = B( a) B( b). (C.94) The proof involves some exercises on elementary math. Proof. We begin by defining a summation of sets of binary vectors. For B 1, B 2,, B α positive integer, and for b B m vec, consider decompositions: B m vec where α is some b = b b α, b i B i for all i (C.95) and denote the number of different decompositions of b as N( b; B 1, B 2,, B α ). Then, we define the following summation: B 1 + B B α { b B m vec : N( b; B 1, B 2,, B α ) = odd }. (C.96) With the summation defined above, the claim of the lemma can be written as follows. By setting B = B( a) B( b), we may notice that from a direct calculation. Therefore, we need to show that The proof relies on the following sublemma. B 1,p+1 = 1 p J b + J b (C.97) B 1,p+1 = 0 otherwise (C.98) J a + J b = J a+ b. (C.99) 57

58 Sublemma 1. Consider αc β α! β!(α β)! (C.100) for (2 m > α β 1). Let the binary representations of α and β be where α = m i=1 α i 2 i 1 and β = m i=1 β i 2 i 1. Then, α C β is odd if and only if α = (α 1,, α m ), β = (β 1,, β m ) (C.101) β i α i for all i. (C.102) We suspect that the sublemma above has been proven somewhere else as it seems elementary. Yet, we could not find a reference, and thus, we present a proof here. Proof. For a given integer p Z 2 m, let h(p) be the largest integer such that speculations, one may notice that where p = (p 1,, p m ). Here, h(p) = p! 2 h(p) is an integer. Then, with some m (2 i 1 1)p i (C.103) i=1 αc β = α! β!(α β)!, (C.104) and h( α C β ) = h(α) h(β) h(α β) 0. (C.105) Then, α C β is odd if and only if h(β) + h(α β) = h(α). (C.106) Let α β γ. Then, we have m m (2 i 1 1)β i + (2 i 1 1)γ i = i=1 i=1 m (2 i 1 1)α i (C.107) i=1 and β i + γ i = α i for all i. This is true if and only if This completes the proof of the sublemma. β i α i for all i. (C.108) Now, let us return to the proof of the lemma. Below, we prove J a + J b show that = J a+ b. Let e 1 (1, 0,, 0). First, we J e1 + + J e1 } {{ } α = J α (C.109) 58

59 where α e e } {{ } 1. (C.110) α Here, notice that β J e1 + + J e1 } {{ } α (C.111) if and only if α C β is odd since J e1 has two elements: (1, 0,, 0) and (0,, 0). Thus, β J α from the sublemma, and we have Now, let us show that J e1 + + J e1 } {{ } α Note that c J a+ b if and only if a+bc c is odd. Here, we have = J α. (C.112) J a + J b = J a+ b. (C.113) a+bc c = c ac i bc c i i=0 (C.114) where x C y = 0 when y > x. Let us assume that a C i bc c i is odd for i = α 1,, α p. Notice that a C i bc c i is odd if and only if both a C i and b C c i are odd. Then, a+b C c is odd if and only if p is odd. Notice that p is the number of decompositions of c such that c = c 1 + c 2 ( c 1 J a and c 2 J b ). (C.115) Therefore, p is odd if and only if c J a + J b. This completes the proof of the lemma. Multiplication of column operators: Next, let us consider how characteristic vectors change under multiplications of column operators. Lemma 12. Consider the following column operators: where Then, If b > b, b = b and V = V. If b = b and V V, b = b and V = VV. If b = b and V = V, b > b. Below, we present the proof. U = UU I (C.116) U V, b U V, b. (C.117) U V, b 59

60 Proof. We start with the first claim. When b m = 1, b m = 1 since b > b. Then, f 0 (U) = f 0 (U ) = I, and f 0 (U ) = I. Thus, b m = 1. When b m = 0 and b m = 0, we have f 0 (U) I and f 0 (U ) I. Suppose f 0 (U ) = I. Then, f 0 (U) = f 0 (U ), and U (1) = U (1). This means b = b which contradicts with b > b. Thus, f 0 (U ) I, and b m = 0. When b m = 0 and b m = 1, we have f 0 (U) = I and f 0 (U ) I, and f 0 (U ) I, and b m = 0. In summary, b m = b m. One can repeat the same discussion and show b = b and V = V. The second claim is easy to prove, so we shall skip the proof. Let us move to the third claim. Since f b1 f bm (U) = f b1 f bm (U ) = V, (C.118) we have f b1 f bm (U ) = I. Let the largest integer i such that f bi f bm (U ) = I be i max. Then, f bimax+1 f bm (U ) I and b i max = 1. If b imax = 1, f bimax+1 f bm (U ) = I which leads to a contradiction. Thus, b imax = 0. Since f bimax+1 f bm (U ) I, we have Thus, b > b. b i = b i (i max + 1 i) (C.119) b imax = 0, b i max = 1. (C.120) Proof of lemma 10: Finally, let us finish the proof of lemma 10. Consider the following column operator V( b) which replaces 1 entries in B( b) with V and 0 entries with I: V( b) j V B( b) 1, j. (C.121) Then, one may easily notice that V( b) has a characteristic vector b with a characteristic operator V. Therefore, from lemma 11, we have V( b)(b( a)) = V( a + b). (C.122) For U with b and V, we decompose U with V( b) as follows: U = V( b)u. (C.123) When U V( b), U has a characteristic vector b with b > b from lemma 12. One can repeat the same decomposition and obtain: U = V( b)v ( b )V ( b ) (C.124) where b < b < b. Then, from lemma 11, we have U(B( a)) = V( a + b)v ( a + b )V ( a + b ). (C.125) Note that V( c)(b( a)) = I if a + c > 2 m 1. From lemma 12, U(B( a)) is a column operator with a characteristic vector a + b and a characteristic operator V. Note that B( a) is an even column when a (0,, 0). This completes the proof of lemma 10. Appendix C.6. The existence of an odd matrix for m > 1 Finally, let us proceed to the proof of lemma 8, to complete the proof of theorem 8. Procedure: Consider a logical operator l defined inside P(x, 2 m, 1) with x > 2v. Let us represent characteristic 60

61 values and characteristic operators for each column of l as follows: U 1 b 1, V 1 U 2 b 2, V 2... U x b x, V x. (C.126) Without loss of generality, we may assume that b i b i+1 for all i since permutations of columns do not affect parities of binary matrices. We define a set of integers such that b i = b as b( b): We denote a group generated by {V i } for i b( b) as V b : b( b) { i : b i = b }. (C.127) V b = { V i : i b( b) }. (C.128) Since x > 2v, there always exists a set of integers A such that V i = I. (C.129) i A We denote the largest vector in { b i } i A as b α and the largest integer i with b i = b α as i = α. Here, we define the following binary matrix E(i ) such that E(i ) E 1,1,, E x,1... E 1,2 m,, E x,2 m (C.130) E i,1 = 1 (i = i ) (C.131) E i, j = 0 otherwise. (C.132) For example, Notice that E(1) = 1, 0,, 0 0, 0,, 0, E(2) =.... 0, 0,, 0 0, 1, 0,, 0 0, 0, 0,, 0. (C.133) , 0, 0,, 0 Par(E(i)) j = δ i, j. (C.134) Now, we define the following operation between a binary column A B m col and a binary matrix B: B A 2 m j=1 T 2 (B A 1, j ) j 1 (mod 2). (C.135) 61

62 Then, consider the following binary matrix: E = E(i) B( b i ) i A (C.136) where b i + b i = b α. Note that E(i) B( b i ) is odd if and only if b i = (0,, 0). This matrix can generate the following column operator: l(e) x = l E(i) B( b i ) i A x ( = U i B( b i ) ) (C.137). i A Note that U i ( B( b i ) ) (C.138) has a characteristic vector b α and a characteristic operator V i from lemma 11. Here, we notice that Par(E) α = 1 (C.139) and E is an odd matrix since b α = (0,, 0). Then, l(e) x I since there is no odd identity generating matrix. Now, we notice that l(e) x is a column operator with a characteristic vector b α > b α from lemma 12. We summarize the discussion so far as follows. From A such that i A V i = I, one can form a column operator l(e) x which has a characteristic vector b α > b α and a characteristic operator V α where E is an odd matrix which satisfies Par(E) α = 1 and Par(E) j = 0 ( j > α). (C.140) Update: Next, we update U α, b α, V α and E(α) to l(e) x, b α, V α and E: U α l(e) x b α b α V α V α E(α) E (C.141) In other words, we replace U α, b α, V α and E(α) with E(α), l(e) x, b α, V α and E, and rename them as U α, b α, V α and E(α). Note that these updated E(i) satisfy Par(E(i)) i = 1 and Par(E(i)) j = 0 ( j > i). (C.142) Then, one may repeat the discussion above. Since x > 2v, there is a set A such that i A V i = I. (C.143) We denote the largest vector in { b i } i A as b β and the largest integer i with b i = b β as i = β. Then, consider the following binary matrix: E = E(i) B( b i ) (mod 2) (C.144) i A 62

63 where b i + b i = b β. This matrix E can generate the following column operator: l(e ) x = l E(i) B( b i ) i A x ( = B( b i ) ). i A U i (C.145) Here, we notice that Par(E ) β = 1 (C.146) since b β = (0,, 0) and E is an odd matrix. Then, l(e ) x I. Note that U i ( B( b i ) ) (C.147) has a characteristic vector b β and a characteristic operator V i from lemma 11. Then, we notice that l(e ) x is a column vector with a characteristic vector b β > b β. Note that Then, we obtain the following observation. Par(E ) β = 0 ( j > β). (C.148) From A such that i A V i = I for updated V i, one can form a column operator l(e ) x which has a characteristic vector b β > b β and a characteristic operator V β. E is an odd matrix which satisfies Par(E ) β = 1 and Par(E ) j = 0 ( j > β). (C.149) Here, we again update U β, b β, V β and E(β) to l(e ) x, b β, V β and E. Then, since updated E(i) always satisfy Par(E(i)) i = 1 and Par(E(i)) j = 0 ( j > i) (C.150) one can repeat the same discussion again. In each update, characteristic vectors b i increase, and at the end, one ends up with the following column operators U 1 1, V 1, E(1) U 2 1, V 2, E(2).... U x 1, V x, E(x) (C.151) where 1 (1,, 1) and Par(E(i)) i = 1 and Par(E(i)) j = 0 ( j > i) (C.152) Then, there exists a set A such that V i = I (C.153) i A 63

64 and, the following matrix is an identity generating matrix E = E(i) (mod 2). (C.154) i A Let the largest integer in A be i max. Then, E is odd since Par(E) imax = 1. However, this contradicts with our original assumption that there is no odd identity generating matrix. This completes the proof of lemma 8, lemma 6 and theorem 8. Appendix D. Derivation of logical operators Having showed that two-dimensional logical operators can be decomposed as a product of two-dimensional and one-dimensional centralizer operators, let us proceed to the proof of theorem 5. The proof owes a lot to arguments presented in [17]. For simplicity of presentation and in order to avoid making the paper unnecessarily long, we shall skip some parts of the derivation. However, we believe that interested readers can easily construct rigorous proofs. Preliminaries: We begin by providing some corollaries and lemma which are useful in the derivations of logical operators. Let us first generalize theorem 8 for any n 2. Corollary 2. Consider a three-dimensional STS model with the system size n 1 = 2 2 2n 2v!, arbitrary n 2 and n 3 > 1. For a given logical operator l supported inside P(n 1, n 2, 1), one can decompose l as a product of the following centralizer operators where l l a l b, l a, l b C P(n1,n 2,1) (D.1) T β 1 (l b) = l b, where β 2 2n 2v (D.2) and l a is defined inside P(2v, n 2, 1). The proof relies on the fact that one can make a logical operator l a quasi-periodic. Proof. Let us represent n 2 as n 2 = 2 m n 2 where n 2 is some odd integer. Then, for a logical operator l defined inside P(n 1, n 2, 1), one can see that the following logical operator l = n 2 j=1 T ( j 1)2m 2 (l) l (D.3) is equivalent to l since l is a product of an odd number of translations of l. Notice that l is periodic in the ˆ2 direction: T 2m 2 (l ) = l (D.4) with the periodicity 2 m. Because of this periodicity, one can form identity generating matrices in a way similar to the cases when n 2 = 2 m. Therefore, one can see that theorem 8 holds for any n 2. We have seen that a two-dimensional logical operator can be decomposed as a product of a one-dimensional centralizer operator and a two-dimensional centralizer operator, as summarized in theorem 8. One can further decompose a one-dimensional logical operator as a product of a one-dimensional centralizer operator and a zero-dimensional centralizer operator, as summarized in the following lemma. Lemma 13. Consider a three-dimensional STS model with the system size n 1 = 2 2 2n 2v!, n 2 = 2 2 (2v)2! and n 3 > 1 where m is an arbitrary positive integer. For a given logical operator l supported inside P(2v, n 2, 1), one can decompose l as a product of the following centralizer operators l l a l b, l a, l b C P(2v,n2,1) (D.5) 64

65 where T β 1 (l b) = l b, where β 2 (2v)2 (D.6) and l a is defined inside P(2v, (2v) 2, 1). We show the claim of lemma 13 graphically in Fig. D.34. One can prove the lemma through discussion similar to the one used in the proof of lemma 5. So, we shall skip the proof. Figure D.34: The claim of lemma 13. Finally, let us extend the claim of theorem 4 slightly. Corollary 3. Consider a two-dimensional STS model with even n 1 and n 2. Let l be a logical operator which is periodic in the ˆ1 and ˆ2 directions: T 1 (l) = T 2 (l) = l. (D.7) Then, l is a two-dimensional logical operator, and there exists a zero-dimensional logical operator r which is defined inside P(1, 2v) and anti-commutes with l. Let l be a logical operator which is defined inside P(1, n 2 ) and periodic in the ˆ2 direction: T 2 (l) = l. (D.8) Then, l is a one-dimensional logical operator, and there exists another one-dimensional logical operator r which is defined inside P(n 1, 1), anti-commutes with l: {l, r} = 0 and periodic in the ˆ1 direction: T 1 (r) = r. (D.9) In other words, if we find a logical operator which is periodic in both ˆ1 and ˆ2 directions, we readily know that it is a two-dimensional logical operator. Similarly, if we find a logical operator which is defined inside P(1, n 2 ) and periodic in the ˆ2 direction, we readily know that it is a one-dimensional logical operator. Note that the corollary holds only for the cases where both n 1 and n 2 are even. 65

66 Proof. Since l is periodic and system sizes are even, l commutes with all the one-dimensional logical operators and all the two-dimensional logical operators. Therefore, l can anti-commute only with zero-dimensional logical operators. This means that l is a two-dimensional logical operator. The second claim can be proven in a similar way. Now, we derive all the logical operators in a three-dimensional STS models, as described in theorem 5. Consider the system size analyzed in lemma 13: (n 1, n 2, n 3 ) = (2 2 2n 2v!, 2 2 (2v)2!, n 3 ) where n 3 > 1 is some fixed integer. Here, we view the entire system as a two-dimensional system by considering P(1, n 2, 1) as a single composite particle (Fig. D.35). (So, the entire system is viewed as a two-dimensional lattice of one-dimensional tubes). For a two-dimensional STS model, we already know the geometric shapes of all the logical operators, as summarized in theorem 4. When viewed as a two-dimensional system, zero-dimensional logical operators are defined inside P(2vn 2, n 2, 1). From theorem 8, we notice that these zero-dimensional logical operators can be actually defined inside P(2v, n 2, 1). Then, we have the following logical operators (Fig. D.35). Figure D.35: Viewed as a two-dimensional system. An anti-commuting pair of logical operators in Fig. D.35(a) where l is defined inside P(2v, n 2, 1) and r is periodic: T 1 (r) = T 3 (r) = r. (D.10) An anti-commuting pair of logical operators in Fig. D.35(b) where l is defined inside P(1, n 2, n 3 ) and periodic: T 3 (l) = l (D.11) and r is defined inside P(n 1, n 2, 1) and periodic: T 1 (r) = r. (D.12) Below, we analyze anti-commuting pairs of logical operators in Fig. D.35(a) and Fig. D.35(b), and derive logical operators. Pairs in (a): Below, we analyze properties of logical operators described above. We start with anti-commuting pairs described in Fig. D.35(a). We stop viewing the system as a two-dimensional system for the moment. Logical operators defined inside P(2v, n 2, 1) consists of periodic one-dimensional logical operators and zero-dimensional logical operators defined inside P(2v, (2v) 2, 1) due to lemma 13. Let us first analyze a zero-dimensional logical operator l defined inside P(2v, (2v) 2, 1). 66

67 From theorem 4, l is also a logical operators for arbitrary n 1 and n 3. Consider the case when n 3 = 1 (Fig. D.36). Then, by viewing the system as a two-dimensional system which extends only in the ˆ1 and ˆ2 directions, one notices that there exists a two-dimensional logical operator r which is periodic in both ˆ1 and ˆ2 directions: T 1 (r) = T 2 (r) = r and anti-commutes with l: {l, r} = 0. Next, let us consider the case when n 3 construction of r as follows (Fig. D.36): (D.13) > 1. Then, one may extend the r r = n 3 x=1 T x 3 (r). (D.14) In other words, we put r in a periodic way in the ˆ3 direction to form r. We shall call such an extension the periodic extension. The three-dimensional logical operator r obtained after the periodic extension of r is periodic in all the directions: T 1 (r ) = T 2 (r ) = T 3 (r ) = r (D.15) and anti-commutes with l : {r, l} = 0. Then, one may notice that r and l form a pair of anti-commuting logical operators for any system size n. From this discussion, we obtain the following observation (Fig. D.36). For zero-dimensional logical operators defined inside P(2v, (2v) 2, 1), there always exists a three-dimensional logical operator r which is periodic: T 1 (r) = T 2 (r) = T 3 (r) = r (D.16) and anti-commutes with l: {l, r} = 0. l and r are logical operators for any system size n. Figure D.36: Constructions of zero-dimensional and three-dimensional logical operators. Next, let us consider a one-dimensional logical operator l defined inside P(2v, n 2, 1) which is periodic in the ˆ2 direction: T β 2 (l) = l (D.17) 67

68 for β 2 (2v)2. Recall that r is periodic: T 1 (r) = T 3 (r) = r (D.18) and anti-commutes with l : {r, l} = 0. Then, one can decompose r as a product of two centralizer operators from lemma 5: r r a r b (D.19) where and r a is defined inside P(n 1, 2v, n 3 ), and for β 2 2v. Then, we notice that T 1 (r a ) = T 3 (r a ) = r a, T 1 (r b ) = T 3 (r b ) = r b (D.20) T β 2 (r b) = r b (D.21) [r b, l] = 0 (D.22) since T β β 2 (l) = l and T2 (r b) = r b, and n 2 /ββ is an even integer for n 2 = 2 2 (2v)2!. Thus, we have {r a, l} = 0. (D.23) Since r a is periodic in the ˆ1 and ˆ3 directions, one can periodically extend its construction for arbitrary n 1 and n 3. Now, let us consider the system size such that n 1 and n 3 are odd. Here, note that r a has some equivalent logical operator r a defined inside P(n 1, 1, n 3 ) [17]. Then, the following operator r a = i, j T i 1 T j 3 (r a) r a (D.24) is equivalent to r a (Fig. D.37). Note that r a is defined inside P(n 1, 1, n 3 ) and periodic in the ˆ1 and ˆ3 directions: T 1 (r a ) = T 3 (r a ) = r a. (D.25) Finally, we show that l can be periodic in the ˆ2 direction. Consider the case when n 1 and n 3 are even, and n 2 = 1. Then, r a is also a logical operators. Now, there always exists some logical operator l defined inside P(1, 2v, 1) which anti-commutes with r a from corollary 3. One can periodically extend its construction to arbitrary n 2. We denote it l. Then, l and r a are logical operators for any system size. From this discussion, we obtain the following observation (Fig. D.37). For one-dimensional logical operators defined inside P(2v, n 2, 1), there exists a two-dimensional logical operator r defined inside P(n 1, 1, n 3 ) which is periodic: T 1 (r) = T 3 (r) = r (D.26) and anti-commutes with l: {l, r} = 0. l can be also periodic: T 2 (l) = l, (D.27) and, l and r are logical operators regardless of the system size n. Pairs in (b): Let us proceed to the analysis on pairs of logical operators in Fig. D.35(b). We consider the following anti-commuting logical operators l and r. 68

69 Figure D.37: Constructions of one-dimensional and two-dimensional logical operators. l is defined inside P(1, n 2, n 3 ) and periodic: T 3 (l) = l. r is defined inside P(n 1, n 2, 1) and periodic: T 1 (r) = r. Since l is periodic in the ˆ3 direction, one can decompose l as follows: l l a l b, T 3 (l a ) = l a and T 3 (l b ) = l b (D.28) where l a is defined inside P(1, 2v, n 3 ), and l b is defined inside P(1, n 2, n 3 ) and periodic: T β 2 (l b) = l b (D.29) where β 2 2v. Thus, logical operators defined inside P(1, n 2, n 3 ) consist of two-dimensional logical operators and one-dimensional logical operators defined inside P(1, 2v, n 3 ). Let us analyze a one-dimensional logical operator l defined inside P(1, 2v, n 3 ) first. We decompose r defined inside P(n 1, n 2, 1) as follows: r r a r b (D.30) where r a is defined inside P(n 1, 2v, 1), and r b is defined inside P(n 1, n 2, 1) and periodic: where β 2 2v. Then, we notice that T β 2 (r b) = r b (D.31) [l, r a ] = 0 (D.32) since there exists a translation of r a which does not overlap with l. Thus, we have {l, r b } = 0. (D.33) Let us consider the case where n 1 = 1, and n 3 is even. Note that l and r b are both logical operators. Then, from corollary 3, there exists a one-dimensional logical operator r which is defined inside P(1, n 2, 1), anti-commutes with 69